3.1.4用MATLAB计算序列的DFT ■xn=[11111111];%输入时域序列向量xn=R8() ■Xk32=fft(xn,32; %计算xn的32点DFT ■ Xk64=t(xn,64): %计算xn的64点DFT ■ %以下为绘图部分 ■k=0:31;wk=2*k/32; %产生32点DFT对应的采样点频率（关于π归一化值） ■subplot(2,2,1);stem(wk,abs(Xk32),'.)；%绘制32点DFT的幅频特性图 tite('(a)32点DFT的幅频特性图)；xlabel(('wApi);ylabel(幅度) subplot(2,2,3);stem(wk,angle(Xk32),');%绘制32点DFT的相频特性图 title((b)32点DFT的相频特性图)； xlabel('wpi')ylabel('相位)；axis(0,2,-3.5,3.5]) k=0:63;wk=2*k/64;%产生64点DFT对应的采样点频率（关于π归一化值） subplot(2,2,2);stem(wk,abs(Xk64),');%绘制64点DFT的幅频特性图 title(c)64点DFT的幅频特性图)；xlabel('wpi门ylabel(幅度) subplot(2,2,4);stem(wk,angle(Xk64),'.);%绘制64点DFT的相频特性图 title('(d)64点DFT的相频特性图) xlabel('w/pi)ylabel('相位')；axis(0,2,-3.5,3.5]) 24 3.1.4用MATLAB计算序列的DFT  xn=[1 1 1 1 1 1 1 1]; %输入时域序列向量xn=R8(n)  Xk32=fft( ) xn,32 ; %计算xn的32点DFT  Xk64=fft(xn,64); %计算xn的64点DFT  %以下为绘图部分  k=0:31;wk=2*k/32; % k=0:31;wk=2*k/32; %产生32点DFT对应的采样点频率 对应的采样点频率(关于π归一化值)  subplot(2,2,1);stem(wk,abs(Xk32),'.'); %绘制32点DFT的幅频特性图  title('(a)32点DFT的幅频特性图');xlabel('ω/\pi');ylabel('幅度')  subplot(2,2,3);stem(wk,angle(Xk32),'.'); %绘制32点DFT的相频特性图  title('(b)32点DFT的相频特性图');  xlabel( xlabel('ω/\pi );ylabel( ');ylabel('相位');axis([0,2, );axis([0,2,-3.5,3.5]) 3.5,3.5])  k=0:63;wk=2*k/64; %产生64点DFT对应的采样点频率(关于π归一化值)  subplot(2,2,2);stem(wk,abs(Xk64),'.'); %绘制64点DFT的幅频特性图  titl ('( )64 title('(c)64点DFT的幅频特性图') l b l(' ');xlabel('ω/\pi') l b l(' i');ylabel('幅度')  subplot(2,2,4);stem(wk,angle(Xk64),'.'); %绘制64点DFT的相频特性图  title('(d)64点DFT的相频特性图') 24 (( )  xlabel('ω/\pi');ylabel('相位');axis([0,2,-3.5,3.5])