电力系统时域研究MATLAB.docVIP

PAGE / NUMPAGES 电力系统时域分析 MATLABA提供了一系列对于控制系统进行的时域分析的函数，本人基于学习到的MATLAB的相关知识，结合自己的专业学习，在电路零输入响应和阶跃响应的仿真实验中应用此软件进行分析，得到了极好的效果。 一、零输入响应 （一）一阶零输入响应 根据传递函数求解单一系统的零输入响应。 H(s)=(s+1)/(s*s*s+9s*s+26s+24) 初始状态下的状态变量为X0=[1,1,0] 创建M文件，在M文件中输入如下程序： 构造传递函数模型。 num=[0,0,1,1,] den=[1,9,26,24] 将传递函数模型转换为状态方程模型。 [a,b,c,d]=tf2ss(num,den) sys=ss(a,b,c,d); 输入初始状态下的状态变量X0。 x0=[1,1,0]; 输入采样时间矢量。 t=0:0.000001:4; 求解系统的零输入响应。 initial(sys,x0,t) 激活run命令，得到相应的状态方程模型矩阵： num=[0,0,1,1,]; den=[1,9,26,24]; [a,b,c,d]=tf2ss(num,den) sys=ss(a,b,c,d); x0=[1,1,0]; t=0:0.000001:4; initial(sys,x0,t) 运行得到： a = -9 -26 -24 1 0 0 0 1 0 b = 1 0 0 c = 0 1 1 d = 0 系统的零输入响应曲线图 系统的零输入响应曲线图 （二）二阶零输入响应 如图所示的二阶电路，如L=0.5H，C=0.02F。初始值uc（0-）=1V，iL(0-)=0，试研究R分别为1Ω，2Ω，3Ω，…，10Ω时，uc（t）和iL(t)的零输入响应，并画出波形。 二阶电路图.. 二阶电路图 1、 R=1时 clear,format compact L=0.5;R=1;C=0.02; uc0=1;iL0=0; alpha=R/2/L;wn=sqrt(1/(L*C)); p1=-alpha+sqrt(alpha^2-wn^2); p2=-alpha-sqrt(alpha^2-wn^2); dt=0.01;t=0:dt:1; num=[uc0,R/L*uc0+iL0/C]; den=[1,R/L,1/L/C]; [r,p,k]=residue(num,den); ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t); iLn=C*diff(ucn)/dt; figure(1),subplot(2,1,1), plot(t,ucn),grid subplot(2,1,2) plot(t(1:end-1),iLn),grid 2、 R=2时 3、 R=3时 4、 R=5时 5、 R=7时 6、 R=10时 通过MATLAB对电阻由小变大过程中电路中的电流、电容上的电压随时间变化的仿真，我更加直观而深刻认识了二阶动态电路的基本工作状态。 二、单位阶跃响应 已知系统框图如图所示： U(s) Y(s) 系统 系统框图系统框图 s 其中， 输入以下MATLAB命令，可得系统的单位阶跃响应曲线如图所示 num=[7 7]; den=[conv(conv([1 0],[1 3]),[1 4 5])]; g=tf(num,den); gg=feedback(g,1,-1); [y,t,x]=step(gg) plot(t,y) 运行结果为： y = 0 0.0111 0.0633 0.1535 0.2639 0.3768 0.4803 0.5679 0.6381 0.6922 0.7330 0.7638 0.7876 0.8069 0.8234 0.8383 0.8522 0.8653 0.8777 0.8895 0.9004 0.9105 0.9196 0.9279 0.9354 0.9420 0.9480 0.9533 0.9580 0.9622 0.9660 0.9694 0.9725 0.9752 0.9777 0.9800 0.9820 0.9838 0.9854 0.9