新·Matlab与信号与系统的应用2.pptVIP

* 注： x可以是向量、矩阵和多维数组；也可以带参数N即X=fft（x，N）。如果X长度小于N，则会自动补0，大于N则自动截断。 * P26 辅助现代工程 * 2011-2012-1 * 模拟系统: 函数[R,P,K]=residue(B,A) [B,A]=residue(R,P,K) 数字系统 函数[R,P,K]=residuez(B,A) [B,A]=residuez(R,P,K) 二(3)、系统描述与转换 * 2011-2012-1 * Num=[2 0 9 1]; Dec = [1 1 4 4]; [R P K] = residue(Num,Dec); 二(3)、系统描述与转换 * 2011-2012-1 * 结果: p= 0.0000+2.0000i 0.0000-2.0000i -1.0000 k= 2 r= 0.0000-0.2500i 0.0000+0.2500i -2.0000 二(3)、系统描述与转换 * 2011-2012-1 * 注意点 以上模型结构都是针对连续时间系统的，对于离散时间系统而言，描述方式类似，不同的是离散时间系统的传递函数中z的幂次都是负值，若要使用连续时间系统的函数描述，必须注意z的阶次（多项式的表达方式）。 * 2011-2012-1 * 注意点 计算上述系统的零、极点时，系统的描述不是如下所示 num = [1 -1] den = [1 0 -1] 而是首先将其转换为正数幂次形式，再列出传递函数模型。 这样传递函数模型为： num = [1 -1 0]; den = [1 0 -1]; [z p k] = tf2zp(num,dec) × √ * 2011-2012-1 * 系统(连续、离散)的单位脉冲响应、阶跃响应 系统对任意输入信号的响应 二(4)、时域分析 * 2011-2012-1 * 求解系统的单位阶跃响应和脉冲响应 单位阶跃响应：step()、dstep()-离散信号 冲激(脉冲)响应：impulse()、dimpulse()、impz() 这里的调用参数都是系统的描述方式 二(4)、时域分析 * 2011-2012-1 * 求系统的阶跃响应曲线 %传递函数描述 num=[20]; den=[1 8 36 0 40 ]; %绘制系统的阶跃响应曲线 t=0:0.1:10; step(num, den, t); 二(4)、时域分析 * 2011-2012-1 * 2输入2输出系统： 求系统的冲激响应 二(4)、时域分析 * 2011-2012-1 * %系统状态空间描述 a=[-2.5 -1.22 0 0;1.22 0 0 0;1 -1.14 -3.2 -2.56; 0 0 2.56 0]; b=[4 1;2 0;2 0;0 0]; c=[0 1 0 3;0 0 0 1]; d=[0 -2;-2 0]; %绘制闭环系统的冲激响应 impulse(a,b,c,d) title(‘impulse response’) xlabel(time-sec) ylabel(amplitude) 二(4)、时域分析 * 2011-2012-1 * %求解系统的脉冲响应 num = 1; den = [1 -0.6 -0.16]; n = 0:1:20; y = dimpulse(num,den,n); figure(1) stem(n,y) Title( The impulse response) 二(4)、时域分析 * 2011-2012-1 * 对于impz()与dimpulse()不同：该函数是针对数字信号系统的，描述系统结构时采用 impz(b,a,n)的形式 其中b、a都是按照z的负值幂次形式，并且从z的零次幂开始，因此对于以上离散系统H(z), b、a分别为： b = [ 0 0 1]; a = [1 -0.6 -0.16]; 二(4)、时域分析 * 2011-2012-1 * 注意：这里的b不能写成b = 1，否则就成为如下的系统： b = [ 0 0 1]; a = [1 -0.6 -0.16]; n =0:1:20; y = impz(b,a,n); figure(2) stem(n,y) Title( The impulse response) 二(4)、时域分析 * 2011-2012-1 * 对任何输入信号，系统响应的求解。 已知系统脉冲响应的，利用输入信号与冲激信号的卷积得到系统的输出，基于MATLAB的conv()命令 已知系统的传递函数，利用滤波原