系统建模与仿真结课作业

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

第三章

1、

clear;

ST=900;Dt=2;NP=ST/Dt;R=1;

for i=1:1:6

if i==2

x10=0;x11=0;x20=0;x21=0;x30=0;x31=0;x40=0;x41=0;x50=0;x51=0;x60=0;x61=0 ;x70=0;x71=0;U10=0;U11=0;

dat=5.5;Ti=51;Td=7.6;

a1=Dt/(dat*Ti);a2=exp(-Dt/(0.1*Td));a3=exp(-Dt/10);

for j=1:NP

y(j)=x71;

U11=R-x70;

x11=x10+a1*U11;

x21=a2*x20+10*(a2-1)*U11/dat;

U21=1/dat*U11+x11+x21+10/dat*U11;

x31=a3*x30+(1-a3)*U21;

x41=a3*x40+(1-a3)*x31;

x51=a3*x50+(1-a3)*x41;

x61=a3*x60+(1-a3)*x51;

x71=a3*x70+(1-a3)*x61;

x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;x60=x61;x70=x71;

end

plot((1:NP)*Dt,y,'r')

hold on;

end

if i==3

x10=0;x11=0;x20=0;x21=0;x30=0;x31=0;x40=0;x41=0;x50=0;x51=0;x60=0;x61=0 ;x70=0;x71=0;U10=0;U11=0;

dat=5.5;Ti=43.2;Td=10.8;

a1=Dt/(dat*Ti);a2=exp(-Dt/(0.1*Td));a3=exp(-Dt/10);

for j=1:NP

y(j)=x71;

U11=R-x70;

x11=x10+a1*U11;

x21=a2*x20+10*(a2-1)*U11/dat;

U21=1/dat*U11+x11+x21+10/dat*U11;

x31=a3*x30+(1-a3)*U21;

x41=a3*x40+(1-a3)*x31;

x51=a3*x50+(1-a3)*x41;

x61=a3*x60+(1-a3)*x51;

x71=a3*x70+(1-a3)*x61;

x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;x60=x61;x70=x71;

end

plot((1:NP)*Dt,y,'r')

hold on;

end

if i==4

x10=0;x11=0;x20=0;x21=0;x30=0;x31=0;x40=0;x41=0;x50=0;x51=0;x60=0;x61=0 ;x70=0;x71=0;U10=0;U11=0;

dat=5.84;Ti=30;Td=10;

a1=Dt/(dat*Ti);a2=exp(-Dt/(0.1*Td));a3=exp(-Dt/10);

for j=1:NP

y(j)=x71;

U11=R-x70;

x11=x10+a1*U11;

x21=a2*x20+10*(a2-1)*U11/dat;

U21=1/dat*U11+x11+x21+10/dat*U11;

x31=a3*x30+(1-a3)*U21;

x41=a3*x40+(1-a3)*x31;

x51=a3*x50+(1-a3)*x41;

x61=a3*x60+(1-a3)*x51;

x71=a3*x70+(1-a3)*x61;

x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;x60=x61;x70=x71;

end

plot((1:NP)*Dt,y,'g')

hold on;

end

if i==5

x10=0;x11=0;x20=0;x21=0;x30=0;x31=0;x40=0;x41=0;x50=0;x51=0;x60=0;x61=0 ;x70=0;x71=0;U10=0;U11=0;

dat=5.2;Ti=45;Td=11;

a1=Dt/(dat*Ti);a2=exp(-Dt/(0.1*Td));a3=exp(-Dt/10);

for j=1:NP

y(j)=x71;

U11=R-x70;

x11=x10+a1*U11;

x21=a2*x20+10*(a2-1)*U11/dat;

U21=1/dat*U11+x11+x21+10/dat*U11;

x31=a3*x30+(1-a3)*U21;

x41=a3*x40+(1-a3)*x31;

x51=a3*x50+(1-a3)*x41;

x61=a3*x60+(1-a3)*x51;

x71=a3*x70+(1-a3)*x61;

x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;x60=x61;x70=x71;

end

plot((1:NP)*Dt,y,'g')

hold on;

end

0100200300400500600700800900

00.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

2、 clear;

ST=500;Dt=0.4;NP=ST/Dt;dat=0.43;dat2=0.05;Ig=1;Ti=98;R=1

x10=0;x11=0;x20=0;x21=0;x30=0;x31=0;x40=0;x41=0;x50=0;x51=0;U10=0;U11=0;U20=0;U21=0;

a1=Dt/(dat*Ti);a2=exp(-1/15*Dt);a4=exp(-1/30*Dt);

for j=1:NP

y1(j)=x51;

y2(j)=x31;

U11=Ig-0.1*x51;

x11=x10+a1*U11;

U21=(x11+U11/dat-0.1*x30)/dat2+R;

x21=a2*x20+(10-10*a2)*U21;

x31=a2*x30+(1-a2)*x21;

x41=a4*x40+(1-a4)*x31;

x51=a4*x50+(1-a4)*x41;

x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;

end

plot((1:NP)*Dt,y1)

hold on ;

plot((1:NP)*Dt,y2,'g')