华北电力大学热工系统建模实验作业.doc

科 技 学 院 名 称： 热工系统建模实验作业 院 系： 动力工程系 班 级： 自动化08K3班 学 号： 081912010116 学生姓名： 李方义 日期： 2011 年 6 月 11 日1、单容水箱水位的DT仿真——显示欧拉公式实验。 实验程序：clear; L=100; F=900;K1=1800; K2=90; RO=1; u1=0.6; ST = 2000; DT = 10;NP = ST / DT; for i=1:NP W1=K1*u1; W2=K2*sqrt(L); Y(i)=L; L=L+DT*(W1-W2)/(RO*F); end plot((1:NP)*DT,Y); title(单容水箱仿真); 实验曲线： 更改初值后的程序：clear; L=100; F=900;K1=1800; K2=90; RO=1; u1=0.6; ST = 10000; DT = 500;NP = ST / DT; for i=1:NP W1=K1*u1; W2=K2*sqrt(L); Y(i)=L; L=L+DT*(W1-W2)/(RO*F); end plot((1:NP)*DT,Y); title(单容水箱仿真); 试验曲线： 实验程序：clear; L=100; F=900;K1=1800; K2=90; RO=1; u1=0.6; ST = 20000; DT = 300;NP = ST / DT; for i=1:NP W1=K1*u1; W2=K2*sqrt(L); Y(i)=L; L=L+DT*(W1-W2)/(RO*F); end plot((1:NP)*DT,Y); title(单容水箱仿真); 试验曲线： 隐式欧拉公式仿真：clear; L=100; F=900;K1=1800; K2=90; RO=1; u1=0.6; ST = 2000; DT = 500;NP = ST / DT; for i=1:NP B2=K2/sqrt(L); L=( RO*F/DT*L+K1*u1)/(RO*F/DT+B2); Y(i)=L; end plot((1:NP)*DT,Y); title; 2.COS函数的仿真实验： 实验程序：clear; x=0:0.01*pi:4*pi; y=cos(x); plot(x,y); grid on; 实验曲线： 3. 单容水箱水位的PID仿真实验。 试验程序：clear all L=0;F=900;K1=1800;K2=90;RO=1;U0=0.5;Lsp=10;Kp=0.05;Err0=0;Ti=30; ST=300;DT=0.5;NP=ST/DT; for i=1:NP Err=Lsp-L; U1=U0+Kp*(Err-Err0)+DT*Err/Ti; % L=L+DT*(K1*U1-K2*(L))/(RO*F); L=(RO*F/DT*L+K1*U1)/(RO*F/DT+K2); Err0=Err; U0=U1; Y1(i)=U1;Y2(i)=L;Y3(i)=Err; end figure(1) %plot((1:NP)*DT,Y1);title(?????÷ê?3?); %figure(2) plot((1:NP)*DT,Y2);title(????); %figure(3) %plot((1:NP)*DT,Y3);title(??2?); grid on; 仿真曲线： 4.弹性阻尼系统的仿真实验。 试验程序：clear; M=10; K=200; B=30; U=75; ST = 16; DT = 0.1;NP = ST / DT-1; for i=1:NP X1(i)=0; X2(i)=0; X(i)=i; Y(i)=i; end for i=1:NP X(i+1)=i; F1(i)=X2(i); F2(i)=(-1.0)*K/M*X1(i)-B/M*X2(i)+U/M; X1(i+1)=X1(i)+DT*F1(i); X2(i+1)=X2(i)+DT*F2(i); end F1(NP+1)=F1(NP); F2(NP+1)=F2(NP); % plot(X,X1); % plot(X,X1,X,X2,X,F2); title(μˉD?×è?á?μí3); plot(X,X1,X,X2); 仿真曲线： 五、PID实验。 试验