现代控制理论实验01.docxVIP

实验1：求传递函数 自动实1301 杨晓丹 实验内容及目的 求出指定电路的传递函数并做仿真。熟悉求出传递函数的方法，学会用MATLAB仿真。 实验方案内容 求出传递函数 用阻抗法求传递函数 R1’和R2’分别表示C1和C2的阻抗 Uo=R2’/(R2+R2’)*U1 R=R1’*(R2+R2’)/(R1’+R2+R2’) U1=R/(R1+R)*Ui 其中R1’=1/C1/S,R2’=1/C2/S 带入解得Uo=1/(R1*R2*C1*C2*S^2+(R1*C1+R1*C2+R2*C2)*S+1)*Ui 用方框图求传递函数 根据电路图列出方程 (Ui-U1)/R1=I →(Ui-U1)/R1=I (U1-Uo)/R2=I2 →(U1-Uo)/R2=I2 I1=I-I2 →I1=I-I2 Uo=1/C2*∫I2dt → Uo=1/C2/S U1=1/C1*∫I1dt → U1=1/C1/S 带入后化简得 U1=(Ui-R1*C2*S*Uo)/(1+R1*C1*S) Uo=1/(R1*R2*C1*C2*S^2+(R1*C1+R1*C2+R2*C2)*S+1)*Ui 用方框图求传递函数 由方程画出方框图 通过变换后得到 Uo=1/(R1*R2*C1*C2*S^2+(R1*C1+R1*C2+R2*C2)*S+1)*Ui 用信号流图求传递函数 将方框图转化为信号流图 P1=1/(R1*C1*S*R2*C2*S) DT1=1 DT=1+1/R1/C1/S+1/C1/S/R2+1/R2/C2/S+1/R/C1/S/R2/C2/S G(S)=1/(R1*R2*C1*C2*S^2+(R2*C2+R1*C2+R1*C1)*S+1) 用MATLAB进行仿真 代码如下 clc clear all; close all; c1=1; c2=1; r1=1; r2=1; num1=[1]; den1=[r1*c1*r2*c2 r1*c1+r2*c2+r1*c2 1]; sys1=tf(num1,den1); [y]=step(sys1); plot(y); Simulink连接图如下 实验结果及分析 用求出的传递函数G(S)=1/(R1*R2*C1*C2*S^2+(R2*C2+R1*C2+R1*C1)*S+1) 得到了两条相似的曲线。 实验中存在的问题 两条曲线并不是完全相同的，似乎每次simulink仿真的曲线开头一段都比较平。 仿真产生的曲线是Uo(t)的曲线还是Uo(s)的曲线？ 若不用step，该如何仿真？编程时输入的信号是否也需要进行拉普拉斯变换？