计算机建模与仿真心得,计算机建模与仿真.docx

实验项目名称：控制系统模型的建立与仿真 实验学时:4

同组学生姓名： 实验地点： —实验日期：

2018.4实验成绩：

批改教师：批改时间：

一、实验目的和要求

熟悉MATLA控制系统工具箱中线性控制系统传递函数模型的相关函数；

熟悉SIMULINK模块库，能够使用SIMULINK进行控制系统模型的建立及仿真。

二、实验仪器和设备

1、PC机1台并安装MATLAB7.0以上版本。

三、实验过程

1.熟悉线性控制系统传递函数模型的相关函数。

tf ()函数可用来输入系统的传递函数

该函数的调用格式为 G = tf ( num, den );

其中num, den分别为系统传递函数的分子和分母多项式系数向量。返回的 G 为系统的传递函数形式。

但如果分子或分母多项式给出的不是完全的展开的形式， 而是若干个因式的 乘积，则事先需要将其变换为完全展开的形式， 两个多项式的乘积在MATLA下借 用卷积求取函数conv()得出，其调用格式为：p=conv(p1,p2)

MATLA还支持一种特殊的传递函数的输入格式，在这样的输入方式下，应该

先用s=tf( ' s')定义传递函数算子，然后用数学表达式直接输入系统的传递函 数。

请自己通过下面两个例子来演示和掌握tf ()和s=tf( ' s')算子这两种输

入方式。

Gs3 5s2 3s 2

G 432

例1设系统传递函数s 2s 4s 3s 1

输入方式一：num = [1,5, 3, 2]; den = [1,2, 4, 3, 1]; %分子多项

式和分母多项式

G = tf ( num, den ) %这样就获得系统的数学模型G

输入方式二：s=tf( ' s' ); G=( sA3 + 5* sA2 + 3* s + 2)/( sA4 + 2*sA3

+ 4* sA2 + 3* s + 1)

任务一：将下列传递函数分别采用上面两种输入方式进行输入，并截图记录。

s 5

G ~3

① s4 2s3 3s 4

>>

>> s=TfC ) ;€=(s+6)/(s 4+2*s 3+3*s-i-4)

s + 5

s\ 十 2 $3 + 3 s -i- 4

Continuous-time transfet function.

G5(s 24)

G 2 2 2

② (s 1)(s 3s 4)(s1)提示：借助conv()函数获得分子和分母多项

式系数

>> s=tf f 才)：G= (5* (s+2. 4)) / ( (&+1 r2* (s * 2+3*s+4) * (s"2+L))

G =

5 s + 12

e*6 4- 5 s*5 + 12 s*4 + 16 s"3 ■+ L5 s2 + 1L s + 4

Cdntinuous-time transfer fuiLction.

zpk()函数可得到控制系统的零极点形式的传递函数

该函数的调用格式为G = zpk (乙P,K )

其中K为系统传函零极点形式的增益，乙P分别为系统传递函数的零点和极 点列向量。返回的G为传递函数形式。

― (s 1.539)(s 2.7305 2.8538j)(s 27305 2.8538j)

G

例2设系统传递函数(s 4)(s 3)(s 2)(s 1)

任务二：在MATLA命令窗口中验证下面的命令，并截图记录。

Q

Q

>> z=[-l, 539 ：-2. 7305+2-85381 ：-2- 7305-2. S538L]:

P=卜丄;-2;-3;-41;6-显忆*,1)

G z

(s+L. 639) (s^2 卜 0.451s 十 15-6)

{s+l) (s+2) (s+3) (s^t)

Continuous-tine zflro/pole/gain nodel.

已知传递数函的分子、分母多项式系数，可用tf2zp()函数求出传函的零

点向量、极点向量和增益。该函数的调用格式为[Z, P,K] = tf2zp ( num, den )

已知传递数函的零点、极点和增益，可用zp2tf ()函数求出传函的分子、

分母多项式系数。该函数的调用格式为[num, den]= zp2tf (乙P,K)

32

s 5s 3s 2

G ~432,.,

任务三：以S 2s 4s 3s 1为例，练习这两个函数的使用，并截图记

录。

? nuti=C J. 5, 3. 21 : den= Z- 4, 3. 11 : CZ. P, Kl = tf 2ac (ntiin, deni

ZntiB. denl = sp2ff (7, Pu K)

-4. 4241 亠 OOOOa.

-0.29€0 4 0.6O76i

-0.2GSO - 0.60761

F -

-0.500