cdbadzhishu.m

signal and system answer

clear
cdbad=[4.108103919247526,7.339397661278048,8.562264442955064,9.813659098928667,12.120090934015199,12.00209019726287,12.98003606423615,14.195253704579523,15.518602250589257,16.04701221649483,16.215722783747026,16.93384749380301,17.738664978179674,18.346228033678397,18.95800154837301,19.54376566753861,20.170241713011233,20.501551167613588,20.91881073241269,21.287441809737896,21.60744439958924,21.791922254533315,22.316655832291975,22.18302078016475,22.15861659039198,22.601107912011283,23.6105351149787,23.409250054812894,23.652976294932213,23.71704609786545,23.54346391747302,23.85220213133485,23.847472594376935,23.986598207266006,24.207107916733353,24.17398349757664,24.08753393091195,24.213130584832108,24.45545074825179,23.613105799491667,21.34231387892763,17.577709553058178,18.94129985151815,24.924910293133735];
scd1=zeros(1,40);scd2=zeros(1,40);scd3=zeros(1,40);
scd1(1)=(cdbad(1)+cdbad(2)+cdbad(3))/3;
scd2(1)=(cdbad(1)+cdbad(2)+cdbad(3))/3;
scd3(1)=(cdbad(1)+cdbad(2)+cdbad(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*cdbad(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=[cdbad(1:40),Xcd(41:50)];plot(n,y1)