clear;clc;
% 因为数据很少,可以列举出年份
year_1=[1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,...
2002,2003,2004,2005,2006,2007,2008,2009,2010,2011,2012,2013,2014];
% 来判断
[~,col]=size(year_1);
DOY=[];
for i=1:col
if ((rem(year_1(i),100)==0 && rem(year_1(i),400)==0) || (rem(year_1(i),100)~=0 && rem(year_1(i),4)==0))
days=366;
else
days=365;
end
DOY=[DOY days];
end
% 将结果转为列的形式
DOY_Re=DOY';
year=year_1';
Date(:,1)=year;
Date(:,2)=DOY_Re;