用FFT对信号作频谱分析 (2).doc

河北联合大学轻工学院
信息科学与技术部
数字信号处理
实验名称：用FFT对信号作频谱分析
班 级： 通信三班
学 号： 201124440317
姓 名： 刘 芳
,'.'); %绘制8点DFT的幅频特性图
title('(3a) 8点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
subplot(2,2,4);stem(abs(X3k16),'.'); %绘制16点DFT的幅频特性图
title('(3b)16点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
%实验内容(2) 周期序列谱分析==================================
N=8;n=0:N-1; %FFT的变换区间N=8
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k8=fft(x4n); %计算x4n的8点DFT
X5k8=fft(x5n); %计算x5n的8点DFT
N=16;n=0:N-1; %FFT的变换区间N=16
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k16=fft(x4n); %计算x4n的16点DFT
X5k16=fft(x5n); %计算x5n的16点DFT
figure(3)
subplot(2,2,1);stem(abs(X4k8),'.'); %绘制8点DFT的幅频特性图
title('(4a) 8点DFT[x_4(n)]');xlabel('ω/π');ylabel('幅度');
subplot(2,2,3);stem(abs(X4k16),'.'); %绘制16点DFT的幅频特性图
title('(4b)16点DFT[x_4(n)]');xlabel('ω/π');ylabel('幅度');
subplot(2,2,2);stem(abs(X5k8),'.'); %绘制8点DFT的幅频特性图
title('(5a) 8点DFT[x_5(n)]');xlabel('ω/π');ylabel('幅度');
subplot(2,2,4);stem(abs(X5k16),'.'); %绘制16点DFT的幅频特性图
title('(5b)16点DFT[x_5(n)]');xlabel('ω/π');ylabel('幅度');
%实验内容(3) 模拟周期信号谱分析===============================
figure(4)
Fs=64;T=1/Fs;
N=16;n=0:N-1; %FFT的变换区间N=16
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)16点采样
X6k16=fft(x6nT); %计算x6nT的16点DFT
X6k16=fftshift(X6k16); %将零频率移到频谱中心
Tp=N*T;