cdgoodzhishu.m

signal and system answer

clear
cdgood=[34.27310489207324,37.61588811710571,39.380411522702396,41.67148020343229,44.20607799366955,43.403540209057525,45.24560551825326,46.98563558905411,47.278706602935415,49.962572069198515,50.172668743778125,50.89165234591165,51.67593815504748,52.3633867449288,52.956003665171295,53.716019909658876,54.17056701819695,55.099268675708345,55.443738738754504,55.702170116945034,55.87456281027987,55.54707879055555,55.82722086540281,54.3237351738607,55.009893381652894,54.65863813538845,55.88154884082104,56.82492138113307,56.95883138752763,57.44737097235684,57.526119648197636,57.56296609194122,57.557910303587526,57.51095228313657,57.42209203058836,57.29132954594288,57.11866482920018,57.04674494193328,56.95654230486423,53.79311330301158,43.39480473389877,29.805912099668618,51.65955935535524];
scd1=zeros(1,40);scd2=zeros(1,40);scd3=zeros(1,40);
scd1(1)=(cdgood(1)+cdgood(2)+cdgood(3))/3;
scd2(1)=(cdgood(1)+cdgood(2)+cdgood(3))/3;
scd3(1)=(cdgood(1)+cdgood(2)+cdgood(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*cdgood(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=[cdgood(1:40),Xcd(41:50)];plot(n,y1)
scd1=zeros(1,40);scd2=zeros(1,40);scd3=zeros(1,40);
scd1(1)=(cdgood(1)+cdgood(2)+cdgood(3))/3;
scd2(1)=(cdgood(1)+cdgood(2)+cdgood(3))/3;
scd3(1)=(cdgood(1)+cdgood(2)+cdgood(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*cdgood(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=[cdgood(1:40),Xcd(41:50)];plot(n,y1)