xn=[1 2 3 3 2 1];
M=1024;
N=32;
Xjw=fft(xn,M); %用1024个点近似xn的傅里叶变换
Xk32=fft(xn,N); %计算xn的32点的DFT
k=0:M-1;wk=2*k/M;
k32=0:N-1;wk32=2*k32/N;
%3.4题第(1)(2)问,绘制xn的傅里叶变换和DFT的幅频特性和相频特性图像
figure(12);
subplot(2,2,1);plot(wk,abs(Xjw)); %绘制xn的傅里叶变换的幅频特性
title('(a)FT[x(n)]的幅频特性');xlabel('w/\pi');ylabel('幅度');
subplot(2,2,2);plot(wk,angle(Xjw)); %绘制xn的傅里叶变换的相频特性
line([0,2],[0,0]);
title('(b)FT[x(n)]的相频特性');xlabel('w/\pi');ylabel('相位');
subplot(2,2,3);stem(wk32,abs(Xk32),'.'); %绘制Xk32的DFT的幅频特性
title('(c)32点DFT[x(n)]的幅频特性');xlabel('w/\pi');ylabel('幅度');
subplot(2,2,4);stem(wk32,angle(Xk32),'.'); %绘制Xk32的DFT的相频特性
line([0,2],[0,0]);
title('(d)32点DFT[x(n)]的相频特性');xlabel('w/\pi');ylabel('相位');
%%3.4题第(3)(4)问,将xn的傅里叶变换和DFT的幅频特性和相频特性图像分别绘制在同一幅图中
figure(34);
%绘制xn的傅里叶变换和Xk32的DFT的幅频特性
subplot(2,2,1);plot(wk,abs(Xjw)); %绘制xn的傅里叶变换的幅频特性
hold on;
stem(wk32,abs(Xk32),'.'); %绘制Xk32的DFT的幅频特性
title('(e)FT[x(n)]和32点DFT[x(n)]的幅频特性');xlabel('w/\pi');ylabel('幅度');
%绘制xn的傅里叶变换和Xk32的DFT的相频特性
subplot(2,2,2);plot(wk,angle(Xjw)); %绘制xn的傅里叶变换的相频特性
hold on;
stem(wk32,angle(Xk32),'.'); %绘制Xk32的DFT的相频特性
title('(f)FT[x(n)]和32点DFT[x(n)]的相频特性');xlabel('w/\pi');ylabel('相位');
subplot(2,2,3);stem(xn32,'.'); %绘制X(k)的32点IDFT
title('(g)32点IDFT[x(n)]波形');xlabel('n');ylabel('x(n)');