连续系统离散化.ppt

连续系统的数字仿真连续系统的离散化连续系统的离散化连续系统的离散化连续系统的离散化信号重构器的特性及传递函数信号重构器的特性及传递函数信号重构器的特性及传递函数信号重构器的特性及传递函数信号重构器的特性及传递函数常用环节的离散相似模型常用环节的离散相似模型常用环节的离散相似模型常用环节的离散相似模型常用环节的离散相似模型常用环节的离散相似模型连续系统按结构图的离散相似法仿真连续系统按结构图的离散相似法仿真连续系统按结构图的离散相似法仿真连续系统按结构图的离散相似法仿真连续系统按结构图的离散相似法仿真k0=10;t1=5;t2=10;kc=2;ti=10;
ts=0.1;N=1000;r=1;
sp=r*ones(1,N);
e0=0;x0=0;u0=0;y0=0;t=0;
fori=1:N
e=r-y0;
u=kc*e+(ts/ti-1)*kc*e0+u0;
x=x0*exp(-ts/t1)+k0*u*(1-exp(-ts/t1));
y=y0*exp(-ts/t2)+x*(1-exp(-ts/t2));
t=t+ts;
tt(i)=t;
yy(i)=y;
u0=u;e0=e;x0=x;y0=y;
end
plot(tt,yy,tt,sp,'--')连续系统按结构图的离散相似法仿真时域解法求取离散系统差分方程时域解法求取离散系统差分方程时域解法求取离散系统差分方程时域解法求取离散系统差分方程时域解法求取离散系统差分方程时域解法求取离散系统差分方程时域解法求取离散系统差分方程例：某火力发电站汽包水位控制系统采用单级三冲量调节方式，系统框图如下页所示。对此系统进行仿真，绘制蒸汽量D扰动时，水位H的变化曲线。+0.15×0.083 状态方程0阶 状态方程0阶0阶0阶0阶 状态方程离散相似系统的差分方程clear;
DT=1;N=300
x10=0;x20=0;x30=0;x40=0;x50=0;H=0;D=2;
fori=1:N
e1=D*0.15*0.083-0.033*H;
x11=exp(-DT/12.02)*x10+43.89*(1-exp(-DT/12.02))*e1;
e2=40.08*e1+x11;
x21=x20+0.037*DT*e2;
x31=exp(-DT/30)*x30+(1-exp(-DT/30))*x21;
x41=exp(-DT/15)*x40+3.6*(1-exp(-DT/15))*D;
x51=x50+0.037*DT*D;
H=x31+x41-x51;
tt(i)=i*DT;
yy(i)=H;
x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;
end
plot(tt,yy);Matlab仿真结果Matlab仿真结果例：某火力发电站主汽温调节系统采用导前微分调节方式，系统框图如下页所示。对此系统进行仿真，并输出导前区汽温2和主汽温1的仿真结果。ig=01010 状态方程 状态方程 状态方程 状态方程 状态方程 状态方程 状态方程 状态方程 状态方程离散相似系统的差分方程离散相似系统的差分方程clear;
DT=1;N=300;
R=1;Delta=1.7;Ti=15;KD=1.7;TD=40;
sp=0*ones(1,N);
x10=0;x20=0;x30=0;x40=0;x50=0;x60=0;x70=0;E30=0;
fori=1:N
E10=(R+E30/Delta+x10)*10;
x21=exp(-DT/15)*x20+8*(1-exp(-DT/15))*E10;
x31=exp(-DT/15)*x30+(1-exp(-DT/15))*x21;
x41=exp(-DT/20)*x40+1.125*(1-exp(-DT/20))*x31;
x51=exp(-DT/20)*x50+(1-exp(-DT/20))*x41;
x61=exp(-DT/20)*x60+(1-exp(-DT/20))*x51;
E20=0.1*x31;
x71=exp(-DT/TD)*x70+KD*(1-exp(-DT/TD))*E20;
E30=-0.1*x61+x71-0.1*KD*x31;
x11=x10+DT/(Delta*Ti)*E30;
tt(i)=i*DT;
yy1(i)=x61;
yy2(i)=x31;
x10=x11;x20=x21;x30=x31;x40=x41;x50=x51;x60=x61;x70=x71;
end
plot(tt,yy1,tt,yy2,tt,sp);Matlab仿真结果例：某发电机励磁系统框图如下图所示，试求URFF为单位阶跃函数时系统的响应。URFF 状态方程 状态方程-30 状态方程 状态方程 状态方程 状态方程离散相似系统的差分方程离散