cdnihezhishupinghuafa.m

signal and system answer

clear
load que1
scd1=zeros(1,40);scd2=zeros(1,40);scd3=zeros(1,40);
scd1(1)=(cdnihemean(1)+cdnihemean(2)+cdnihemean(3))/3;
scd2(1)=(cdnihemean(1)+cdnihemean(2)+cdnihemean(3))/3;
scd3(1)=(cdnihemean(1)+cdnihemean(2)+cdnihemean(3))/3;
acd=zeros(1,40);bcd=zeros(1,40);ccd=zeros(1,40);Xcd=zeros(1,20);
Acd=0.15;
for i=2:44
    scd1(i)=Acd*cdnihemean(i-1)+(1-Acd)*scd1(i-1);
end
for i=2:44
    scd2(i)=Acd*scd1(i)+(1-Acd)*scd2(i-1);
end
for i=2:44
    scd3(i)=Acd*scd2(i)+(1-Acd)*scd3(i-1);
end
for i=1:44
    acd(i)=3*scd1(i)-3*scd2(i)+scd3(i);
    bcd(i)=Acd*0.5/(1-Acd)^2*((6-5*Acd)*scd1(i)-2*(5-4*Acd)*scd2(i)+(4-3*Acd)*scd3(i));
    ccd(i)=Acd^2/(1-Acd)^2*(scd1(i)-2*scd2(i)+scd3(i));
end
for i=40:50 
    Xcd(i)=acd(39)+bcd(39)*(i-39)+0.5*ccd(39)*(i-39)^2;
end
n=1:50;y1=[cdnihemean(1:40),Xcd(41:50)];plot(n,y1)