数字信号处理Matlab实现实例.pdfVIP

实用标准文案 数字信号处理 Matlab 实现实例 第１章 离散时间信号与系统 例 1-1 用 MATLAB计算序列 {-2 0 1 –1 3} 和序列 {1 2 0 -1} 的离散卷积。 解 MATLAB程序如下： a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n); ylabel( 幅度 ); 图 1.1 给出了卷积结果的图形，求得的结果存放在数组 c 中为： {-2 -4 1 3 1 5 1 -3} 。 例 1-2 用 MATLAB计算差分方程 当输入序列为 时的输出结果 。 解 MATLAB程序如下： N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; 精彩文档 实用标准文案 k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel(n);ylabel( 幅度 ) 图 1.2 给出了该差分方程的前 41 个样点的输出，即该系统的单位 脉冲响应。 例 1-3 用 MATLAB计算 例 1-2 差分方程 所对应的系统函数的 DTFT。 解 例 1-2 差分方程所对应的系统函数为： 1 2 3 0.8 0.44z 0.36 z 0.02z H ( z) 1 2 3 1 0.7 z 0.45z 0.6 z 其 DTFT为 j j 2 j 3 j 0.8 0.44e 0.36e 0.02 e H (e ) j j 2 j 3 1 0.7e 0.45e 0.6e 用 MATLAB计算的程序如下： k=256; num=[0.8 -0.44 0.36 0.02]; den=[1 0.7 -0.45 -0.6]; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1); 精彩文档 实用标准文案 plot(w/pi,real(h));grid title( 实部 ) xlabel(\omega/\pi);ylabel( 幅度 ) subplot(2,2,2); plot(w/pi,imag(h));grid title( 虚部 ) xlabel(\omega/\pi);ylabel(Amplitude)