cdbad.m

signal and system answer

clear
load cdbad
nihex1=zeros(140,8);
nihey1=ones(140,8)*1000;
for i=1:140
    k=1;
    for j=1:45
        if cdfen1nihe(i,j)~=1000&xiabiao1(i,j)~=0
            nihex1(i,k)=xiabiao1(i,j);
            nihey1(i,k)=cdfen1nihe(i,j);
            k=k+1;
        end
    end
end
cdfen1nihe=0.2*cdfen1nihe;
nihey1=0.2*nihey1;
l=ones(1,140);
for i=1:140
    k=1;
    while nihex1(i,k)~=0
       k=k+1;
    end
    k=k-1;
    x1=nihex1(i,1:k);
  
  
    while nihey1(i,l(i))~=1000*0.2
       l(i)=l(i)+1;
    end
    l(i)=l(i)-1;
    y1=nihey1(i,1:l(i));
   
    pcd1(i,:)=polyfit(x1,y1,2);
end
for i=1:140
    for j=1:45
        if cdfen1nihe(i,j)==1000*0.2;
            cdfen1nihe(i,j)=polyval(pcd1(i,:),j);
        else
            l(i)=l(i)-1;
        end
          if cdfen1nihe(i,j)<=0|l(i)==0
            cdfen1nihe(i,j)=0;
        end
      
    end
end


bilv1cd=zeros(140,45);
for j=1:140
    k=1;
    while k<=45&cdfen1nihe(j,k)~=0
        k=k+1;
    end
    k=k-1;
    for i=1:k
        if i==1
            bilv1cd(j,1)=(-3*cdfen1nihe(j,1)+4*cdfen1nihe(j,2)-cdfen1nihe(j,3))/2;
        elseif i~=k
            bilv1cd(j,i)=(cdfen1nihe(j,i+1)-cdfen1nihe(j,i-1))/2;
        elseif k>2
            bilv1cd(j,k)=(cdfen1nihe(j,k-2)-4*cdfen1nihe(j,k-1)+3*cdfen1nihe(j,k))/2;
        end
    end
end
numx1=zeros(1,45);numy1=zeros(1,45);
for i=1:140
    for j=1:45
        if bilv1cd(i,j)~=0
            numx1(j)=numx1(j)+1;
        end
        if cdfen1nihe(i,j)~=0
            numy1(j)=numy1(j)+1;
        end
    end
end
sc=sum(cdfen1nihe);
cdfen1nihemean=sc(1:44)./numy1(1:44);
month=1:44;
figure
plot(month,cdfen1nihemean,'r.')