利用MATLAB在信号与系统中的应用案例.ppt

1、部分分式展开的MATLAB实现 [r,p,k]=residuez(num,den) num,den分别为F(z)分子多项式和分母多项式的系数向量。 r为部分分式的系数，p为极点，k为多项式的系数。 若为真分式，则k为零。 2、H(z)的零极点与系统特性的MATLAB计算 利用tf2zp函数计算H(z)的零极点，调用形式为 [z,p,k]=tf2zp(b,a) H(z)零极点分布图可用zplane函数画出，调用形式为 zplane(b,a) b和a分别为H(z)分子多项式和分母多项式的系数向量。 返回值z为零点、p为极点、 k为增益常数。 [例1] 将F(z)用部分分式展开。 %program8_1 num = [18]; den = [18 3 -4 -1]; [r,p,k] = residuez(num,den) 运行结果为 r =0.3600 , 0.2400 , 0.4000 p =0.5000 , -0.3333 , -0.3333 k =[] 故F(z)可展开为 [例2]试画出系统 的零极点分布图，求其单位冲激响应h[k]和频率响应H(ejW)。 b =[1 2 1];a =[1 -0.5 -0.005 0.3]; figure(1);zplane(b,a); num=[0 1 2 1]; den=[1 -0.5 -0.005 0.3]; h=impz(num,den); figure(2);stem(h) xlabel(k) title(Impulse Respone) [H,w]=freqz(num,den); figure(3);plot(w/pi,abs(H)) xlabel(Frequency \omega) title(Magnitude Respone) 运行结果 八、利用MATLAB进行系统状态变量分析 微分方程到状态方程的转换 状态方程系统函数矩阵H(s)的计算 MATLAB求解连续系统的状态方程 MATLAB求解离散系统的状态方程 微分方程到状态方程的转换 num,den 分别表示系统函数H(s)的分子和分母多项式 A,B,C,D 分别为状态方程的矩阵。 [A,B,C,D]=tf2ss(num,den) 状态方程系统函数矩阵H(s)的计算 A,B,C,D 分别表示状态方程的矩阵。 K 表示函数ss2tf计算的与第k个输入相关的系统函数，即H(s)的第k列。 num 表示H(s)第k列的m个元素的分子多项式 den 表示H(s)公共的分母多项式。 [num,den]= ss2tf (A,B,C,D,k) MATLAB求解连续系统的状态方程 sys 由函数ss构造的状态方程模型 t 需计算的输出样本点, t=0:dt:Tfinal f(:,k) 系统第k个输入在t上的抽样值 x0 系统的初始状态(可缺省) y(:,k) 系统的第k个输出 to 实际计算时所用的样本点； x 系统的状态 获得连续系统状态方程的计算机表示模型 sys = ss(A,B,C,D) 求解状态方程 [y,to,x]=lsim(sys,f,t,x0); MATLAB求解离散系统的状态方程 sys 由函数ss构造的状态方程模型 f(:,k) 系统第k个输入序列 x0 系统的初始状态(可缺省) y(:,k) 系统的第k个输出 n 序列的下标; x 系统的状态 获得离散状态方程的计算机表示模型 sys = ss(A,B,C,D,[]) 求解状态方程 [y,n,x]=lsim(sys,f,[],x0) [例1]写出系统 的状态方程 由 [A,B,C,D]=tf2ss([1],[1 5 10]) 可得 所以系统的状态方程为 C=0 D=0 [例2]已知某连续系统的状态方程和输出方程为 其初始状态和输入分别为 求该系统的系统函数矩阵H(s) 和输出。 计算系统函数矩阵H(s) A=[2 3;0 -1];B=[0 1; 1 0]; C=[1 1; 0 -1];D=[1 0; 1 0]; [num1,den1]=ss2tf(A,B,C,D,1) [num2,den2]=ss2tf(A,B,C,D,2) 运行结果 num1 =1 0 -1 1 -2 0 den1 =1 -1 -2 num2 =0 1 1 0