仿真技术与应用实验报告2.docVIP

仿真技术与应用 实验报告 电气工程及自动化 班 级： 学 号： 姓 名： 哈尔滨工业大学 2016年3月 实验二 利用根匹配法对RLC串联电路的仿真 一．实验目的 熟悉MATLAB的工作环境； 掌握在MATLAB命令窗口调试运行程序； 掌握M文件编写规则及在MATLAB命令窗口运行程序； 掌握利用根匹配法构造离散模型的方法。 二．实验内容 系统电路如图2.1所示。电路元件参数：直流电压源，电阻，电感，电容。电路元件初始值：电感电流，电容电压。系统输出量为电容电压。 连续系统输出响应的解析解为： （2-1） 其中，，。 三．实验过程 （1)建立根匹配法离散数据模型 对于题中的RLC电路，对它的数学模型进行拉布拉斯变换，可得到： 近一步求得它的传递函数为： 传递函数 两个极点分别为 将映射到z平面，得 写出离散传递函数，，根据终值相等的原则，求出，系统的差分方程 经计算，可获得的差分模型为： （2）仿真程序 %---------ＲＬＣ串联电路零输入阶跃响应----------------% % R,L,C,Us: ＲＬＣ串联电路的元件和电源参数 % t,T: 动态过程时间和离散时间间隔 % 实验目的：１）比较用根匹配法和双线性替换法得到的与连续系统相似的离散时间模型准确性 % ２）改变离散时间间隔Ｔ，观察对离散模型准确性的影响 function []=RLC(R,L,C,Us,t,h) %---------- 系统及仿真参数 -------------% R=10; L=0.01; C=1.0e-5; Us=1; t=0.02; T = 1.0e-3; NNN = fix(t/T);% 计算离散时间序列长度 %---------------------利用状态方程参数求G(z)------------- % AA=[-R/L,-1/L;1/C,0];BB=[1/L;0];CC=[0,1];DD=[0]; % [Az,Bz,Cz,Dz]=c2dm(AA,BB,CC,DD,T,matched); % [Zzs,Pzs,Kzs]=ss2zp(AA,BB,CC,DD); %[Zzx,Pzx,Kzx]=ss2zp(Az,Bz,Cz,Dz); %Zzs,Pzs,Kzs %------- 解析解 ------------------- tad = R/(2*L); W=sqrt(1/(L*C)-(R/(2*L))^2); for k=1:1:NNN y(k) = Us*(1-exp(-tad*(k-1)*T) * ( cos(W*(k-1)*T) + sin(W*(k-1)*T)*tad/W)); end %----------传递函数 G(s) = 1/(L*C*S^2+R*C*S+1)---------% a=[1]; b=[L*C R*C 1]; Ps = roots(b); % G(s)的极点；多相式求根 Ps1 = Ps(1); Ps2 = Ps(2); Ks =1/(L*C); % 放大系数 G(s) = Ks/(s-Ps1)(s-Ps2) % --------- 离散传递函数 G(z) = (Kz*(z-Qz1)(z-Qz2))/((z-Pz1)(z-Pz2))---------% Pz1 =exp(T*Ps1); Pz2 =exp(T*Ps2); Qz1 = 0; Qz2 = 0; % 终值 Ys = Ks/(Ps1*Ps2); Kz = Ys*(1-Pz1)*(1-Pz2)/((1-Qz1)*(1-Qz2)); %====================================== %ZZ=[0;0];PP=[Pz1;Pz2];KK=Kz; %[num,den]=zp2tf(ZZ,PP,KK); %num,den, % 差分方程 a2y(n+2)+a1y(n+1)+a0(y(n)=b2us(n+2)+b1us(n+1)+