第三章 离散时间信号的变换域分析.pptVIP

? 3.5.1 两个实序列 DFT 的计算 3.5 实序列 DFT 的计算 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? N N N N k N X k X k X k X j k H k X k X k G n jh n g n x n h n g ? ? ? ? ? ? ? ? ? ? ? * * * * 2 1 2 1 其中 ，则 为实序列，令 、 设 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 2 0 1 2 2 1 0 1 2 2 2 2 1 0 2 1 0 1 0 2 1 0 1 0 1 2 2 1 0 2 2 1 2 0 2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? N k k H W k G W n h W W n g W W n h W n g W n v W n v W n v n v N n n v n h n v n g N n v N k N N N n nk N k N N n nk N N n k N nk N N n nk N N n k n N N n nk N N n nk N ， ， ， 点实序列， 为 设 ? 3.5.2 2N 点实序列 DFT 的计算 3.6 使用 DFT 计算线性卷积 ? 3.6.1 两个有限长序列的线性卷积 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n h L n g n y n h n g n y L n M M n n h n h L n N N n n g n g N M L M N n h n g e e C L e e ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 则 ，定义 ，令 和 长为 、 设 1 0 1 0 1 0 1 0 1 例 3.20 x=input( 输入第一个序列 =); h=input( 输入第二个序列 =); L=length(x)+length(h)-1; XE=fft(x,L); HE=fft(h,L); y1=ifft(XE.*HE); k=0:1:L-1; subplot(2,1,1); stem(k,y1); title( 基于 DFT 的线性卷积结果 ); xlabel( 时间序号 n);ylabel( 振幅 ); y2=conv(x,h); error=y1-y2; subplot(2,1,2); stem(k,error); title( 误差序列 ); xlabel( 时间序号 n);ylabel( 振幅 ); ? 3.6.2 有限长序列与无限长序列的线性卷积 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 点） 均为 和 点重叠（ 处有 在 与 注： 其中 ，则 其中 点序列 为因果序列，切分为 设 ） 、重叠相加（ 为无限长序列 点有限长序列， 为 设 1 1 2 0 1 0 1 1 1 1 0 1 0 0 1 0 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? M N n y n y M M rN n rN n x n h n y n x n h n y l n x n h n y mN n y l n x l h n y otherwise N n mN n x n x mN n x n x n x N n x method add overlap n x n h l n x l h n y n x M l h r r r r r r m m m M l m m m m M l 例 3.21 R=64; d=rand(R,1)-0.5; for m=1:1:R; s(m)=2*(m-1)*((0.9)^(m-1)); x(m)=s(m)+d(m); end k=0:1:R-1; M=input( 滑动平均滤波器长度 =); h=ones(1,M)/M; y=fftfilt(h,x,4); plot(k,s,r-,k,y,b*); legend(r-,s[n],b*,y[n]); xlabel( 时间