matlab控制系统传递函数模型

合集下载
相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

MATLAB及

控制系统仿真

实验

班级:智能0702

姓名:刘保卫

学号: 06074053(18)

实验四控制系统数学模型转换及MATLAB实现

一、实验目的

熟悉MATLAB 的实验环境。

掌握MATLAB 建立系统数学模型的方法。

二、实验内容

(注:实验报告只提交第2 题)

1、复习并验证相关示例。

(1)系统数学模型的建立

包括多项式模型(Transfer Function,TF),零极点增益模型(Zero-Pole,ZP),状态空间模型

(State-space,SS);

(2)模型间的相互转换

系统多项式模型到零极点模型(tf2zp),零极点增益模型到多项式模型(zp2tf),状态空间模

型与多项式模型和零极点模型之间的转换(tf2ss,ss2tf,zp2ss…);

(3)模型的连接

模型串联(series),模型并联(parallel),反馈连接(feedback)

2、用MATLAB 做如下练习。

(1)用2 种方法建立系统的多项式模型。

程序如下:

%建立系统的多项式模型(传递函数)

%方法一,直接写表达式

s=tf('s')

Gs1=(s+2)/(s^2+5*s+10)

%方法二,由分子分母构造

num=[1 2];

den=[1 5 10];

Gs2=tf(num,den)

figure

pzmap(Gs1)

figure

pzmap(Gs1)

grid on

运行结果:

易知两种方法结果一样

Transfer function:

s

Transfer function:

s + 2

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

s^2 + 5 s + 10

Transfer function:

s + 2

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

s^2 + 5 s + 10

(2)用2 种方法建立系统的零极点模型和多项式模型。

程序如下:

%方法一

s=tf('s')

Gs1=10*(s+1)/((s+1)*(s+5)*(s+10))

% zpk模型

ZPK=zpk(Gs1)

%方法二

% tf模型

num=[10 10];

den=conv([1 1],conv([1 5],[1 10]));

Gs2=tf(num,den)

% zpk模型

ZPK=zpk(Gs2)

figure

pzmap(Gs1)

figure

pzmap(Gs1)

grid on

运行结果:

易知两种方法结果一样

Transfer function:

s

Transfer function:

10 s + 10

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

s^3 + 16 s^2 + 65 s + 50

Zero/pole/gain:

10 (s+1)

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

(s+10) (s+5) (s+1)

Transfer function:

10 s + 10

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

s^3 + 16 s^2 + 65 s + 50

Zero/pole/gain:

10 (s+1)

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

(s+10) (s+5) (s+1)

(3)如图,已知G(s)和H(s)两方框对应的微分方程是:

且初始条件为零。试求传递函数C(s)/R(s)及E(s)/R(s)。

程序如下:

%求微分方程的传递函数C(s)/R(s),E(s)/R(s)

%求Gs=Cs/Rs

n1=[20];

d1=[6 10];

Gs=tf(n1,d1)

%求Hs=Bs/Cs

n2=[10];

d2=[20 5];

Hs=tf(n2,d2)

% C(s)/R(s)

sys=feedback(Gs,Hs)

% E(s)/R(s)=(Es/Cs)*(Cs/Rs)

ER=sys/Gs

运行结果:

Transfer function:

20

-------- % Gs=Cs/Rs

6 s + 10

Transfer function:

10

-------- % Hs=Bs/Cs

20 s + 5

Transfer function:

400 s + 100

--------------------- % C(s)/R(s)

120 s^2 + 230 s + 250

Transfer function:

2400 s^2 + 4600 s + 1000

------------------------ % E(s)/R(s)=(Es/Cs)*(Cs/Rs) 2400 s^2 + 4600 s + 5000

相关文档
最新文档