que1cdniheright.m

signal and system answer

clear
load quecdnihe
cdnihe=ones(356,45)*1000;
for i=1:1665
    j=fix(A(i,1));k=fix(A(i,2)+1);
    if isnan(j)==0&isnan(k)==0&k>0&k<=45
       cdnihe(j,k)=A(i,3);
    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,8);
nihey=ones(356,8)*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=0.2*cdnihe;nihey=0.2*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*0.2
       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*0.2;
            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;
month=1:45;
plot(month,bilvcdmean,'r.')
figure
plot(month,cdnihemean,'r.')