信号与系统-matlab使用教程.ppt

一、微分方程(wēi fēn fānɡ chénɡ)到状态方程的转换 num,den 分别表示系统函数H(s)的分子(fēnzǐ)和分母多项式； A,B,C,D 分别为状态方程的矩阵。 [A,B,C,D]= tf2ss(num,den) 第八十九页，共100页。 二、状态方程系统函数矩阵(jǔ zhèn)H (s)的计算 A, B, C, D 分别表示状态方程的矩阵。 K 表示函数ss2tf计算的与第k个输入(shūrù)相关的系统函数，即H(s)的第k列。 num 表示H(s)第k列的m个元素的分子多项式 den 表示H(s)公共的分母多项式。 [num,den]= ss2tf (A,B,C,D,k) 第九十页，共100页。 三、MATLAB求解(qiú jiě)连续系统的状态方程 sys 由函数(hánshù)ss构造的状态方程模型 t 需计算的输出样本点, t=0:dt:Tfinal x(:,k) 系统第k个输入在t上的抽样值 q0 系统的初始状态(可缺省) y(:,k) 系统的第k个输出 to 实际计算时所用的样本点； q 系统的状态 获得连续系统(xìtǒng)状态方程的计算机表示模型 sys = ss(A,B,C,D) 求解状态方程 [y,to,q]=lsim(sys,x,t,q0) 第九十一页，共100页。 四、MATLAB求解(qiú jiě)离散系统的状态方程 sys 由函数ss构造的状态方程模型 x(:,k) 系统(xìtǒng)第k个输入序列 q0 系统(xìtǒng)的初始状态(可缺省) y(:,k) 系统(xìtǒng)的第k个输出 n 序列的下标; q 系统(xìtǒng)的状态 获得离散状态方程的计算机表示(biǎoshì)模型 sys = ss(A,B,C,D,[ ]) 求解状态方程 [y,n,q]=lsim(sys,x,[ ],q0) 或直接利用 [y, q]=dlsim(A,B,C,D,x,q0) 第九十二页，共100页。 例1 写出系统(xìtǒng) 的状态方程。 由 [A,B,C,D]=tf2ss([1],[1 5 10]) 可得 所以(suǒyǐ)系统的状态方程为 C=[0 1] D=0 第九十三页，共100页。 例2 已知某连续系统(xìtǒng)的状态方程和输出方程为 其初始状态和输入(shūrù)分别为 求该系统的系统函数(hánshù)矩阵H(s) 和输出。 第九十四页，共100页。 二、用数值积分分析连续非周期(zhōuqī)信号频谱 数值函数积分(jīfēn)quad8可用来计算非周期信号频谱 F 是一个字符串，它表示(biǎoshì)被积函数的文件名。 a,b 分别表示定积分的下限和上限 quad8的返回是用自适应Simpson算法得出的积分值 y = quad8(F,a,b) 第五十七页，共100页。 例3 试用(shìyòng)数值方法近似计算三角波信号的频谱 X(jw)= Sa2(w / 2) 解: 图示三角(sānjiǎo)波可表示为 三角(sānjiǎo)波信号频谱的理论值为 第五十八页，共100页。 例3 试用数值(shùzí)方法近似计算三角波信号的频谱 function y=sf1(t,w);y=(t=-1 t=1).*(1-abs(t)).*exp(-j*w*t); w=linspace(-6*pi,6*pi,512); N=length(w);X=zeros(1,N); for k=1:N X(k)=quad8(sf1,-1,1,[],[],w(k)); end figure(1); plot(w,real(X));title() xlabel(\omega);ylabel(‘X(j\omega)); figure(2); plot(w,real(X)-sinc(w/2/pi).^2); xlabel(\omega);title(计算误差); 第五十九页，共100页。 例3 试用数值方法近似计算三角(sānjiǎo)波信号的频谱 运行(yùnxíng)结果 第六十页，共100页。 三、离散(lísàn)周期信号频谱的MATLAB的实现 函数fft可用来计算离散周期(zhōuqī)信号频谱 x 是离散周期信号0~N-1 一个(yī ɡè)周期的序列值 Xm = fft(x) Xm 是离散周