系统建模实验报告.doc

文档介绍：数学模型的MATLAB描述实验报告
一实验目的：
学****使用MATLAB软件对数学模型进行描述
二实验设备：
个人电脑
MATLAB 软件
三实验步骤与方法：
一 MATLAB中数学模型的表示
1请将下面的传递函数摸型输入到MATLAE
s? + 4s + 2
(a) ⑻-s3(s2 + 2)[(s2 + 1)3 + 2s + 5]' (b)———
s? + 2s +1
(a)解：
输入：s=tf(,s');G=(sA3+4*s+2)/(sA3*(sA2+2)*((sA2+l)A3+2*s+5))
输出：
Transfer function:
sA3 + 4 s + 2
sAll +5 sA9 + 9 sA7 + 2 sA6 + 12 sA5 + 4 sA4 + 12 sA3 (b)解：
输入:num=[l,l]；
dei匸[1,2,1]；
G=tf(num,den)
输出：
Transfer function:
s+ 1 sA2 + 2 s + 1
2请将下而的零极点模型输入到MATLAB
(a) G(s)=
8(s + 1 - j)($ + 1 + j)
s2(s + 5)(s + 6)(s2 + 1)''
解：输入:G=zpk([-l+li;-l-li],[0;0;-5;-6;li;li],8)
输出： Warning: Not all complex roots come in conjugate pairs (transfer function has complex coefficients).
In at 200
Zero/pole/gain:
8 (sA2 + 2s + 2) sA2 (s+5) (s+6) (s-l)A2
3、状态空间刚L」
P74~.2~.194-1
假设输出方程为：y = [l 0
输入:
ss([0,l;-l,-2],[0,l]',[0,l],0)
输出：
a =
xl
x2
xl
0
1
x2
-1
-2
b =
ul
xl
0
x2
1
xl x2
yl 0 1
1 =
ul
yl 0
Continuous-time model.
4、复杂传谨诙求取——conv (,) 磁用于多项式磁」
输入：num=5*[l,l,l];
» den=conv(conv(conv([l,3,l],[l,3,l]),[l,6,5,3]),[l,2]);
» G=tf(num,den)
输出：
Transfer function:
5 sA2 + 5 s + 5
sA8 + 14 sA7 + 76 sA6 + 209 sA5 + 320 sA4 + 289 sA3
+ 161 sA2 + 49 s + 6
二模型之间的转换
传递函数到零极点
输入： num=[l ,2] ;den=[ 1,2,1 ] ;[z,p,k]=tf2zp(num,den); » G=tf(num,den)
输出:Transfer function:
s + 2 sA2 + 2 s + 1
零极点到传递函数
» [num,den]=z