利用MATLAB进行信号分析.docVIP

利用MATLAB进行信号分析 数字滤波器的频率响应： 例：对一数字滤波器H(z)=，求其幅频特性和相频特性： a=[1 0.4 1]; b=[0.2 0.3 1]; freqz(b,a,128) 离散系统零极点图： b =[0.2 0.1 0.3 0.1 0.2]; a=[1.0 -1.1 1.5 -0.7 0.3]; zplane(b,a) 数字滤波器的冲激响应： b=[0.2 0.1 0.3 0.1 0.2]; a=[1.0 -1.1 1.5 -0.7 0.3]; impz(b,a,50) 计算离散卷积： x=[1 1 1 1 0 0]; y=[2 2 3 4]; z=conv(x,y); stem(0:length(x)+length(y)-2,z) 系统函数转换： H(z)= y=[0.1 0.5 1]; x=[0.3 2 5]; k=1; b=k*poly(y); a=poly(x); --------------------------------------------------- b b = 1.0000 -1.6000 0.6500 -0.0500 a a = 1.0000 -7.3000 12.1000 -3.0000 则H（z）= ------------------------------------------------------ a=[1 –7.3 12.1 –3]; b=[ –1 6 0.65 –0.05]; y=roots(b); x=roots(a); d=b(1)/a(1); --------------------------------------------- y y = 1.0000 0.5000 0.1000 x x = 5.0000 2.0000 3000 ----------------------------------------------- b=[8 -4 11 -2]; a=[1 -1.25 0.75 -0.125]; y=roots(b); x=roots(a) ---------------------------------------------- y = 0.1550 + 1.1367i 0.1550 - 1.1367i 0.1900 x = 0.5000 + 0.5000i 0.5000 - 0.5000i 0.2500 计算离散信号的能量和功率： 序列的能量：E== 序列的功率：P== 例：x(n)=0.8nu(n)，计算前10点的能量占总能量的百分比。 前十点能量： N=10; n=0:N-1; x=(0.8).^n; e=sum(abs(x).^2) E=1/(1-0.8.^2) e/E 用FFT计算周期序列的频谱：X(ejω)= X(k) ， N为序列x(n)的周期。 x(n)=cos(n+)，x(n)=[1 2 3 4],y(n)=[5 6 7 8],求z(n)=x(n)*y(n)。 x=[1 2 3 4];y=[5 6 7 8]; z=conv(x,y) x=[1 2 3 4];y=[5 6 7 8]; X=fft(x)；Y=fft(y); Z=X.*Y；z=ifft(Z) x=[1 2 3 4];y=[5 6 7 8];L=7; X=fft(x,L);Y=fft(y,L); Z=X.*Y;z=ifft(Z,L) x=[1 2 3 4];y=[5 6 7 8];L=8; X=fft(x,L);Y=fft(y,L); Z=X.*Y;z=ifft(Z,L) 求逆Z变换： 求X(z)= ﹥0.9的逆Z变换。 b=1; a=poly([0.9 0.9 -0.9]); [R P C]=residuez(b,a) 执行结果为： R = 0.2500 0.5000 0.2500 P = 0.9000 0.9000 -0.9000 C = [] 则：X(z)= + 相应的逆Z变换为：x(n)=0.25(0.9)nu(n)+ (n+1)(0.9)n+1u(n+1)+ 0.25(-0.9)nu(n) =0.75(0.9)nu(n)+ 0.5n(0.9)n+1u(n)+ 0.25(-0.9)nu(n) FFT的应用 试用FFT计算