实验二利用根匹配法对RLC串联电路的仿真..docVIP

实验二利用根匹配法对RLC串联电路的仿真 一．实验目的 熟悉MATLAB的工作环境； 掌握在MATLAB命令窗口调试运行程序； 掌握M文件编写规则及在MATLAB命令窗口运行程序； 掌握利用替换法构造离散模型的方法 二．实验内容 电路如图１所示电路进行仿真试验。元件参数：，，，。初始值：，。输出量电容电压。 实验步骤 求解传递函数 根据系统的传递函数是在零初始条件下输出量的拉普拉斯变换与输入量的拉普拉斯变换之比，可得该系统的传递函数： ，其中输入，， 系统的离散化处理 G（s）有两个极点，分别是： 将p映射到z平面上，得 写出离散传递函数 求连续系统的单位阶跃响应的终值 确定离散系统单位阶跃响应的终值 根据终值相等的原则，得 = 取附加零点为，，则 确定差分模型 进行z逆变换得 准确解的求解 根据电路理论知识可得 由初始条件可得 A==arctan U=1 ==500 ==10000 = 带入数据可得该系统的精确解，此解可用于对仿真结果进行分析，比较和调试。 用MATLAB进行仿真 1）根据离散方程和精确方程编写MATLAB程序,程序段如下： y(1)=0; y(2)=0; a=(-1e-5+sqrt(1e-10-4e-8))/2e-8; b=(-1e-5-sqrt(1e-10-4e-8))/2e-8; T=1e-6; %步长设定 q1=0; q2=0; %附加零点取值 p1=exp(T*a) p2=exp(T*b) %初始常量计算 for n=1:13000 u(n)=1; end for n=1:13000 t=T*(n-1) z(n)=-[1e4/sqrt(1e8-250000)]*exp(-500*t)*sin(sqrt(1e8-250000)*t+87pi/180)+u(n) %精确解计算 end for n=1:12998 y(n+2)=(p1+p2)*y(n+1)-p1*p2*y(n)+(1-p1)*(1-p2)*(u(n+2)+(q1+q2)*u(n+1)+q1*q2*u(n)) %根匹配计算 end plot(u); hold on; plot(y); %用蓝线画出根匹配曲线 hold on; plot(z,r) %用红线画出精确解 clear %清除变量 取合适的T和n的值，进行仿真。 下面分别取四组代表数据和整体及局部截图为例进行分析。 当 时 q1=q2=0时 跟匹配法 双线性 当 q1=q2=0时 跟匹配 双线性 当 q1=q2=0时 跟匹配 双线性 当 时,改变q1，q2使图形相位相同，经实验仿真当 q1= q2= 时较好，但极值相差加大 实验结论 1）难易性：替换法是对于给定的函数G(s),找到S域到Z域的特殊映射关系，通过一些基本变换以及简化处理将G（s）变换成G（z）。而根匹配法是构造一个相应于系统传递函数的离散传递函数，使两者的零点，极点相匹配，并且使二者具有相同的动态响应值。其过程中涉及到较复杂的运算与过程，明显较之于简单替换法复杂。所以替换法比根匹配法简单。 2）模型的稳定性：在离散时间间隔T取不同值时，简单替换法的曲线变化很大，且在步长为时即为发散说明其稳定性较差，而双线性替换法和根匹配法的曲线变化很小，稳定性好。 3）模型的精度：根匹配法的精度比替换法好。 4）离散时间间隔： 根匹配法的离散时间间隔比前两种替换法要求的离散时间间隔要小。