信号分析处理第二次实验.doc

第二次实验 例6：已知序列x(n)=cos(0.82pi*n)+2sin(0.43pi*n) ,0≦n≦100，绘制想x(n)及其离散傅里叶X（k）的幅度，相位图。 原程序如下：clf; clc; N=100; n=[0:N-1]; xn=cos(0.82*pi*n)+2*sin(0.43*pi*n); Xk=dft(xn,N); magXk=abs(Xk); angleXk=angle(Xk); figure(1); plot(xn); xlabel(n); ylabel(x(n)); title(x(n)N=100); figure(2); k=0:length(magXk)-1; plot(k,magXk); xlabel(k); ylabel(|X(k)|); title(|X(k)|N=100); figure(3); plot(k,angleXk); xlabel(k); ylabel(angle(X(k))); title(angle(X(k))N=100); 源程序结果如下： 例15 用impz函数求离散时间系统y[k]+3y[k-1]+2y[k-2]=f[k]的单位脉冲响应好h[k],并与理论值h[k]=-(-1)^k+2(-2)^k,k≧0相比较。 源程序：k=0:10; a=[1 3 2]; b=[1]; h=impz(b,a,k); subplot(1,2,1) stem(k,h) hk=-(-1).^k+2*(-2).^k; subplot(1,2,2) stem(k,hk) 实验结果如下： 例17：分别用ztrans 和iztrans 函数求（1）x(k)=sin(ak)u(k)的z变换 （2）X(Z)=1/((1+Z)*(1+Z))的z反变换 原程序如下：x=sym(sin(a*k)); X=ztrans(x) X=sym(1/(1+z)^2); x=iztrans(X) 结果如下： 例18：使用matlab计算X(Z)=1/(1+3Z^(-1)+Z^(-2))的部分分式展开式。 原程序如下：num=[1]; den=[1,3,1]; [r,p,k]=residuez(num,den) 结果如下： 例19：使用matlab，求系统传函H（z）=（1+2z）/(1+2z+3z^2+z^3)的 零极点并画出零极点图 系统的单位脉冲响应好h（n）和频率特性函数H(Ω). 原程序如下：b=[2,1]; a=[1,3,2,1]; [r,p,k]=tf2zp(b,a) Zplane(b,a) 结果如下： 原程序如下：num=[0,0,2,1]; den=[1,3,2,1]; h=impz(num,den); figure(1); stem(h) xlabel(‘k’); title(‘Impulse Respone’) [H,w]=freqz(num,den); figure(2) plot(w/pi,abs(H)) xlabel(‘Frequency\omega’) title(‘Magnitude Response’) 结果如下：