sizhongliaofanihe.m

signal and system answer

clear
load data1
clear A
A=zeros(1251,5)
a=zeros(1239,5);b=zeros(1251,5);c=zeros(1254,5);d=zeros(1292,5);
j=1;k=1;l=1;m=1;
for i=1:5036
    if B(i,2)==2
        A(j,:)=B(i,:);
        j=j+1;
    end
end
A(:,5)=exp(A(:,5))-1
aaa=[1239 1251 1254 1292];
cdnihe=ones(356,45)*1000;
j=1;
for i=1:1251
    if i~=1
    l=fix(A(i-1,1));
    end
     if i~=1&l~=A(i,1)
         j=j+1;
     end
    k=fix(A(i,4)+1);
    if isnan(j)==0&isnan(k)==0&k>0&k<=45
       cdnihe(j,k)=A(i,5);
    end
   
end
pcd=zeros(356,3);xiabiao=zeros(356,45);
for i=1:356
    for j=1:45
        if cdnihe(i,j)~=1000
            xiabiao(i,j)=j;
        end
    end
end
nihex=zeros(356,10);
nihey=ones(356,10)*1000;
for i=1:356
    k=1;
    for j=1:45
        if cdnihe(i,j)~=1000&xiabiao(i,j)~=0
            nihex(i,k)=xiabiao(i,j);
            nihey(i,k)=cdnihe(i,j);
            k=k+1;
        end
    end
end
cdnihe=cdnihe;nihey=nihey;
l=ones(1,356);
for i=1:356
    k=1;
    while nihex(i,k)~=0
       k=k+1;
    end
    k=k-1;
    x=nihex(i,1:k);
  
  
    while nihey(i,l(i))~=1000
       l(i)=l(i)+1;
    end
    l(i)=l(i)-1;
    y=nihey(i,1:l(i));
   
    pcd(i,:)=polyfit(x,y,2);
end
for i=1:356
    for j=1:45
        if cdnihe(i,j)==1000
            cdnihe(i,j)=polyval(pcd(i,:),j);
        else
            l(i)=l(i)-1;
        end
          if cdnihe(i,j)<=0|l(i)==0
            cdnihe(i,j)=0;
        end
      
    end
end


bilvcd=zeros(356,45);
for j=1:356
    k=1;
    while k<=45&cdnihe(j,k)~=0
        k=k+1;
    end
    k=k-1;
    for i=1:k
        if i==1
            bilvcd(j,1)=(-3*cdnihe(j,1)+4*cdnihe(j,2)-cdnihe(j,3))/2;
        elseif i~=k
            bilvcd(j,i)=(cdnihe(j,i+1)-cdnihe(j,i-1))/2;
        elseif k>2
            bilvcd(j,k)=(cdnihe(j,k-2)-4*cdnihe(j,k-1)+3*cdnihe(j,k))/2;
        end
    end
end
x=zeros(1,45);y=zeros(1,45);
for i=1:356
    for j=1:45
        if bilvcd(i,j)~=0
            x(j)=x(j)+1;
        end
        if cdnihe(i,j)~=0
            y(j)=y(j)+1;
        end
    end
end

bilvcdmean=sum(bilvcd)./x;
cdnihemean=sum(cdnihe)./y;
cdnihemean(41:45)=[];
bilvcdmean(41:45)=[];
month=1:40;
plot(month,bilvcdmean,'r.')
figure
plot(month,cdnihemean,'r.')