MATLAB实验报告3-控制系统仿真.docVIP

MATLAB实验报告3控制系统仿真

1、一个传递函数模型：

将该传递函数模型输入到MATLAB工作空间。

num=6*[1，5];

den=conv(conv([1，3，1]，[1，3，1])，[1，6]);

tf(num,den)

2、若反馈系统为更复杂的结构如图所示。其中

，，

则闭环系统的传递函数可以由下面的MATLAB命令得出：

G1=tf([1,7,24,24],[1,10,35,50,24]);

G2=tf([10,5],[1,0]);

H=tf([1],[0.01,1]);

G_a=feedback(G1*G2,H)

得到结果：

Transferfunction:

0.1s^5+10.75s^4+77.75s^3+278.6s^2+361.2s+120

--------------------------------------------------------------------

0.01s^6+1.1s^5+20.35s^4+110.5s^3+325.2s^2+384s+120

图复杂反馈系统

3、设传递函数为：

试求该传递函数的部分分式展开

num=[2,5,3,6];

den=[1,6,11,6];

[r,p,k]=residue(num,den)

4、给定单位负反馈系统的开环传递函数为：

试画出伯德图。

利用以下MATLAB程序，可以直接在屏幕上绘出伯德图如图20。

num=10*[1,1];

den=[1,7,0];

bode(num,den)

5、已知三阶系统开环传递函数为：

画出系统的奈氏图，求出相应的幅值裕量和相位裕量，并求出闭环单位阶跃响应曲线。

G=tf(3.5,[1,2,3,2]);

subplot(1,2,1);

%第一个图为奈氏图

nyquist(G);

grid

xlabel(RealAxis)

ylabel(ImagAxis)

%第二个图为时域响应图

[Gm,Pm,Wcg,Wcp]=margin(G)

G_c=feedback(G,1);

subplot(1,2,2);

step(G_c)

grid

xlabel(′Time(secs)′)

ylabel(′Amplitude′)

6．绘制下列各系统根轨迹图。