hivzhishupinghuafa.m

signal and system answer

clear
hivnihemean=[5.095820783341429,4.452483209581492,4.189913281945378,3.876701375187596,3.431933896729329,3.477797999814873,3.381813771892126,3.200390410226402,3.064234861011267,2.952680941706319,2.863186994331568,2.755253288942034,2.662849431311355,2.578974239966391,2.510102974111348,2.454728765291937,2.414402683154488,2.386519203783996,2.379641636977733,2.379072502648362,2.392709208639219,2.417585292565604,2.453317547941533,2.470899075643746,2.551853804278938,2.508943938900803,2.386126781280675,2.401407098847278,2.401921349556324,2.386963801780501,2.416437391239251,2.453469509869026,2.458700421496229,2.4879552391846,2.543260621171543,2.605332646386495,2.674171314829457,2.736242320747744,2.793008756258855,2.894207179694991,3.166003561365927,3.342090379785534,3.40695839646896,2.687299538344168];
shiv1=zeros(1,40);shiv2=zeros(1,40);shiv3=zeros(1,40);
shiv1(1)=(hivnihemean(1)+hivnihemean(2)+hivnihemean(3))/3;
shiv2(1)=(hivnihemean(1)+hivnihemean(2)+hivnihemean(3))/3;
shiv3(1)=(hivnihemean(1)+hivnihemean(2)+hivnihemean(3))/3;
ahiv=zeros(1,40);bhiv=zeros(1,40);chiv=zeros(1,40);Xhiv=zeros(1,20);
Ahiv=0.18;
for i=2:40
    shiv1(i)=Ahiv*hivnihemean(i-1)+(1-Ahiv)*shiv1(i-1);
end
for i=2:40
    shiv2(i)=Ahiv*shiv1(i)+(1-Ahiv)*shiv2(i-1);
end
for i=2:40
    shiv3(i)=Ahiv*shiv2(i)+(1-Ahiv)*shiv3(i-1);
end
for i=1:40
    ahiv(i)=3*shiv1(i)-3*shiv2(i)+shiv3(i);
    bhiv(i)=Ahiv*0.5/(1-Ahiv)^2*((6-5*Ahiv)*shiv1(i)-2*(5-4*Ahiv)*shiv2(i)+(4-3*Ahiv)*shiv3(i));
    chiv(i)=Ahiv^2/(1-Ahiv)^2*(shiv1(i)-2*shiv2(i)+shiv3(i));
end
for i=41:50
    Xhiv(i)=ahiv(39)+bhiv(39)*(i-39)+0.5*chiv(39)*(i-39)^2;
end
n=1:50;y1=[hivnihemean(1:40),Xhiv(41:50)];plot(n,y1)