fig342.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理（第3版）》图3.4.2程序: fig342.m
% 电子工业出版社出版  高西全  丁玉美 合著 2006年8月
xn=ones(1,5);                    %设置x(n)=R10(n)
hn=ones(1,4);
L=length(hn);nh=0:L-1;%设置h(n)=R6(n)
subplot(5,2,1),stem(nh,hn,'.');axis([0,10,0,1.5])
ylabel('h(n)');xlabel('n');title('(a) h(n)波形')
L=length(xn);nx=0:L-1;
subplot(5,2,5),stem(nx,xn,'.');axis([0,10,0,1.5])
ylabel('x(n)');xlabel('n');title('(b) x(n)波形')
%调用conv计算线性卷积
yn=conv(hn,xn);
L=length(yn);ny=0:L-1;
subplot(5,2,9),stem(ny,yn,'.');axis([0,10,0,5])
ylabel('y(n)');xlabel('n');title('(c) y(n)=h(n)*x(n)')
%循环卷积等于线性卷积的条件:循环卷积区间长度L>=M+N-1
%L=15,按图3.2.3用DFT计算序列线性卷积
L1=6;
Xk=fft(xn,L1);                    %L1点FFT[x(n)]
Hk=fft(hn,L1);                    %L1点FFT[h(n)]
Yck=Xk.*Hk;                     %频域相乘得Yc(k)
ycln=ifft(Yck,L1);                 %L1点IFFT得到15点循环卷积结果yc(n)
nycl=0:L1-1;
subplot(5,2,2),stem(nycl,ycln,'.');axis([0,10,0,5])
ylabel('yc(n)');xlabel('n');title('(d) 6点循环卷积')
%L=8,按图3.4.1用DFT计算序列循环卷积
L2=8;
Xk=fft(xn,L2);                      % L2点FFT[x(n)]
Hk=fft(hn,L2);                      % L2点FFT[h(n)]
Yck=Xk.*Hk;                      %频域相乘得Yc(k)
yc2n=ifft(Yck,L2);                  % L2点IFFT得到13点循环卷积结果yc(n)
nyc2=0:L2-1;
subplot(5,2,6),stem(nyc2,yc2n,'.');axis([0,10,0,5])
ylabel('yc(n)');xlabel('n');title('(e) 8点循环卷积')
%L=10,按图3.4.1用DFT计算序列循环卷积
L3=10;
Xk=fft(xn,L3);                      % L3点FFT[x(n)]
Hk=fft(hn,L3);                      % L3点FFT[h(n)]
Yck=Xk.*Hk;                      %频域相乘得Yc(k)
yc3n=ifft(Yck,L3);                  % L2点IFFT得到13点循环卷积结果yc(n)
nyc3=0:L3-1;
subplot(5,2,10),stem(nyc3,yc3n,'.');axis([0,10,0,5])
ylabel('yc(n)');xlabel('n');title('(f) 10点循环卷积')