MATLAB系统模型建立和动态特性实验教程分析.docVIP

实验 MATLAB系统模型建立实验 实验目的 1．掌握如何使用MALAB进行系统模型的建立； 2．实验仪器 计算机 实验内容 1.多项式的向量表示 MATLAB中多项式用向量表示，行向量元素依次为降幂排列的多项式各项系数。 例如：多项式，表示为： P=[1,3,0,2,5] 注意：尽管s2项系数为0，但输入P时不可缺省0 2.多项式乘法 MATLAB中多项式乘法处理函数调用格式为： C=conv（A，B） 例如：给定两个多项式A(s)=s+3和,求C（s）=A（s）B（s）， 则应先构造多项式A(s)和B(s)，然后再调用conv( )函数来求C(s)。 A=[1,3];B=[10,20,3]; C=conv(A,B) 输出结果： C = 10 50 63 9 或直接写成 C=conv([1,3],[10,20,3]) 例如： 3.传递函数 在MATLAB中直接用矢量组表示传递函数的分子、分母多项式系数，即： num = [bm bm-1…b0]； 表示传递函数的分子多项式系数 den = [an an-1…a0]； 表示传递函数的分母多项式系数 sys = tf(num,den) tf命令将sys变量表示成传递函数模型。例用MATLAB表示传递函数为??的系统。解：在MATLAB环境下输入 ng=[1 1]; dg=[1 3 2]; sys=tf(ng,dg)则执行后得到如下结果：sys = s + 1 ------------- s^2 + 3 s + 2 Continuous-time transfer function. 3.2零极点增益模型 在MATLAB中用z、p、k矢量组分别表示系统的零点、极点和增益，即： z = [ z1 z2…… zm ]； p = [ p1 p2…… pn ]； k =[ k ]； sys = zpk(z,p,k) zpk命令将sys变量表示成零极点增益模型。例2 用MATLAB表示传递函数为??的系统。解：在MATLAB环境下输入 z=-1; p=[0 -1 -2]; k=1.5; Gz=zpk(z,p,k) Zero/pole/gain: 1.5 (s+1) ------------- s (s+1) (s+2) 3.3、状态空间模型 x = ax + bu y = cx + du 在MATLAB中用（a、b、c、d）矩阵组表示，然后 sys = ss(a,b,c,d) ss命令将sys变量表示成状态空间模型。 程序： num=[12 24 12 20]; den=[2,4,6,2,2]; G=tf(num,den); GG=zpk(G) 运行结果： Zero/pole/gain: 6 (s+1.929) (s^2 + 0.07058s + 0.8638) ------------------------------------------------- (s^2 + 0.08663s + 0.413) (s^2 + 1.913s + 2.421) 2.零、极点形式转换为多项式形式 例： z=[-1,-2]; p=[0,-5,-10]; k=10; Gz=zpk(z,p,k); GG=tf(Gz) 运行结果为： Transfer function: 10 s^2 + 30 s + 20 ------------------- s^3 + 15 s^2 + 50 s 3.5传递函数的部分分式展开 当： ?时 在MATLAB中直接用分子/分母的系数表示时有????num=[b0，b1，…，bm]；????den = [a0，a1，…，an]；则命令????[r,p,k] = residue(num,den)将求出两个多项式Y(s)和X(s)之比的部分分式展开的留数、极点和直接项。Y(s)/X(s)的部分分式展开由下式给出： 例 考虑下列传递函数： 命令 [r,p,k] = residue(num,den)将给出下列结果： [r,p,k]=residue(num,den)r=-6.000-4.0003.000p=-3.000-2.000-1.000k=2 留数为列向量r，极点位置为列向量p，直接项是行向量k。以下是Y(s)/X(s)的部分分式展开的MATLAB表达形式： 命令[num,den] = residue(r,p,k)执行后得到