基于matlab的控制系统与仿真-3

举报
AXYZdong 发表于 2022/01/21 10:37:45 2022/01/21
【摘要】 基于matlab的控制系统仿真及应用(第二版),张聚,习题3。

Author:AXYZdong 自动化专业 工科男
有一点思考,有一点想法,有一点理性!
定个小小目标,努力成为习惯!在最美的年华遇见更好的自己!


习题3


3.3 绘制以下传递函数模型的单位阶跃响应曲线

G ( s ) = 5 s + 8 s 4 + 4 s 3 + 6 s 2 + 3 s + 3 G(s)=\frac{5s+8}{s^4+4s^3+6s^2+3s+3}

>>sys=tf([5 8],[1 4 6 3 3])

sys =
 
             5 s + 8
  -----------------------------
  s^4 + 4 s^3 + 6 s^2 + 3 s + 3
 
Continuous-time transfer function.

>> step(sys)

运行结果:

在这里插入图片描述

3.7 绘制以下系统的根轨迹曲线
(1)

K s ( s 2 + 2 s + 2 ) ( s 2 + 6 s + 13 ) \frac{K}{s(s^2+2s+2)(s^2+6s+13)}

(2)

K ( s + 12 ) ( s + 1 ) ( s 2 + 12 s + 100 ) ( s + 10 ) \frac{K(s+12)}{(s+1)(s^2+12s+100)(s+10)}

G1= tf([1],conv([1 0],conv([1 2 2],[1 6 13])))

G1 =
 
                   1
  ------------------------------------
  s^5 + 8 s^4 + 27 s^3 + 38 s^2 + 26 s
 
Continuous-time transfer function.

>> rlocus(G1)
>> G2= tf([1 12],conv([1 1],conv([1 12 100],[1 10])))

G2 =
 
                  s + 12
  --------------------------------------
  s^4 + 23 s^3 + 242 s^2 + 1220 s + 1000
 
Continuous-time transfer function.

>> rlocus(G2)

在这里插入图片描述

3.8 已知反馈系统的开环传递函数为:

G ( s ) = K ( s 2 + 2 s + 4 ) s ( s + 4 ) ( s + 6 ) ( s 2 + 1.4 s + 1 ) G(s)=\frac{K(s^2+2s+4)}{s(s+4)(s+6)(s^2+1.4s+1)}

试画出系统的根轨迹和根轨迹渐近线。

>> G3= tf([1 2 4],conv([1 0],conv([1 4],conv([1 6],[1 1.4 1]))))

G3 =
 
                s^2 + 2 s + 4
  -----------------------------------------
  s^5 + 11.4 s^4 + 39 s^3 + 43.6 s^2 + 24 s
 
Continuous-time transfer function.

>> rlocus(G3)
>> sgrid

在这里插入图片描述

在这里插入图片描述


3.9 已知单位反馈控制系统的开环传递函数为:

G ( s ) = K ( s + 2 ) s ( s + 4 ) ( s + 8 ) ( s 2 + 2 s + 5 ) G(s)=\frac{K(s+2)}{s(s+4)(s+8)(s^2+2s+5)}

试画出下列两种情形的根轨迹图;
(1)负反馈控制系统的根轨迹图;
(2)正反馈控制系统的根轨迹图;

>> G4= tf([1 2],conv([1 0],conv([1 4],conv([1 8],[1 2 5]))))

G4 =
 
                   s + 2
  ---------------------------------------
  s^5 + 14 s^4 + 61 s^3 + 124 s^2 + 160 s
 
Continuous-time transfer function.

>> rlocus(G4)
>> rlocus(-G4)

在这里插入图片描述

▲ 负反馈

在这里插入图片描述

▲ 正反馈


3.10 已知某控制系统的开环传递函数

G ( s ) = K s ( s + 1 ) ( s + 2 ) G(s)=\frac{K}{s(s+1)(s+2)}

K = 1.5 K=1.5 时,试绘制系统的开环频率特性曲线,并求出系统的幅值裕量与相位裕量。

G5=zpk([],[0 -1 -2],1.5)

G5 =
 
       1.5
  -------------
  s (s+1) (s+2)
 
Continuous-time zero/pole/gain model.

>> margin(G5)
>> [Gm,Pm,Wcg,Wcp] = margin(G5)

Gm =

    4.0000


Pm =

   41.5340


Wcg =

    1.4142


Wcp =

    0.6118

在这里插入图片描述


3.11 已知一个典型的二阶环节传递函数为

G ( s ) = ω n 2 s 2 + 2 ζ ω + ω n 2 G(s)=\frac{\omega_n^2}{s^2+2\zeta\omega+\omega_n^2}

其中 ω = 0.7 \omega=0.7 ,试分别绘制 ζ = 0.1 0.4 1.0 1.6 2.0 \zeta=0.1、0.4、1.0、1.6、2.0 时的 Bode 图。

wn=0.7;
s=tf('s');
n=[0.1,0.4,1.0,1.6,2.0];
for i=n
figure
G=wn^2/(s^2+2*i*wn*s+wn^2);
bode(G);
end

在这里插入图片描述

在这里插入图片描述

▲ 从上到下$\zeta=0.1、0.4、1.0、1.6、2.0$


3.13 已知系统开环传递函数为

G k ( s ) = 3 ( s + 1 ) ( s + 0.8 + 1.6 j ) ( s + 0.8 1.6 j ) G_k(s)=\frac{3(s+1)}{(s+0.8+1.6j)(s+0.8-1.6j)}

试求出 Nyquist 曲线。

>> G=zpk([-1],[-0.8-1.6*j,-0.8+1.6*j],3)

G =
 
       3 (s+1)
  ------------------
  (s^2 + 1.6s + 3.2)
 
Continuous-time zero/pole/gain model.

>> nyquist(G)

在这里插入图片描述

3.14 已知二阶系统传递函数为:

G ( s ) = 1 s 2 + 2 ζ s + 1 G(s)=\frac{1}{s^2+2\zeta s+1}

试绘制阻尼系数 ζ \zeta 分别为 0.4 0.7 1.0 1.3 0.4、0.7、1.0、1.3 时的系统的 Nyquist 曲线。

n=[0.4,0.7,1.0,1.3];
for i=n
figure
G=tf([0 1],[1 2*i 1]);
nyquist(G);
end

在这里插入图片描述

在这里插入图片描述

▲ 从上到下$\zeta=0.4、0.7、1.0、1.3$


  本次的分享就到这里


11

好书不厌百回读,熟读自知其中意。将学习成为习惯,用知识改变命运,用博客见证成长,用行动证明努力。
如果我的博客对你有帮助、如果你喜欢我的博客内容,请 “点赞” “评论” “收藏” 一键三连哦!
听说 👉 点赞 👈 的人运气不会太差,每一天都会元气满满呦!^ _ ^
**码字不易,大家的支持就是我坚持下去的动力。点赞后不要忘了👉关注👈我哦!


如果以上内容有任何错误或者不准确的地方,欢迎在下面👇留个言。或者你有更好的想法,欢迎一起交流学习~~~

【版权声明】本文为华为云社区用户原创内容,转载时必须标注文章的来源(华为云社区)、文章链接、文章作者等基本信息, 否则作者和本社区有权追究责任。如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。