实验 1:
x1n=[ones(1,4)]; %产生序列向量 x1(n)=R4(n)
M=8;xa=1:(M/2); xb=(M/2):-1:1; x2n=[xa,xb]; %产生长度为 8 的三角波序列 x2(n)
x3n=[xb,xa];
X1k8=fft(x1n,8); %计算 x1n 的 8 点 DFT
X1k16=fft(x1n,16); %计算 x1n 的 16 点 DFT
X2k8=fft(x2n,8); %计算 x1n 的 8 点 DFT
X2k16=fft(x2n,16); %计算 x1n 的 16 点 DFT
X3k8=fft(x3n,8); %计算 x1n 的 8 点 DFT
X3k16=fft(x3n,16); %计算 x1n 的 16 点 DFT
%以下绘制幅频特性曲线
subplot(2,2,1);n=linspace(0,2-2/8,8);stem(n,abs(X1k8),'.'); %绘制 8 点 DFT 的幅频特性图
title('(1a) 8 点 DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X1k8))])
subplot(2,2,3);n=linspace(0,2-2/16,16);stem(n,abs(X1k16),'.'); %绘制 16 点 DFT 的幅频特性图
title('(1b)16 点 DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X1k16))])
figure(2)
subplot(2,2,1);n=linspace(0,2-2/8,8);stem(n,abs(X2k8),'.');%绘制 8 点 DFT 的幅频特性图
title('(2a) 8 点 DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X2k8))])
subplot(2,2,2);n=linspace(0,2-2/16,16);stem(n,abs(X2k16),'.'); %绘制 16 点 DFT 的幅频特性图
title('(2b)16 点 DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X2k16))])
subplot(2,2,3);n=linspace(0,2-2/8,8);stem(n,abs(X3k8),'.'); %绘制 8 点 DFT 的幅频特性图
title('(3a) 8 点 DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X3k8))])
subplot(2,2,4);n=linspace(0,2-2/16,16);stem(n,abs(X3k16),'.'); %绘制 16 点 DFT 的幅频特性图
title('(3b)16 点 DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,1.2*max(abs(X3k16))])
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;F=1/Tp; %频率分辨率 F
k=-N/2:N/2-1;fk=k*F; %产生 16 点 DFT 对应的采样点频率(以零频率为中心)
subplot(3,1,1);stem(fk,abs(X6k16),'.');box on %绘制 8 点 DFT 的幅频特性图
title('(6a) 16 点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度');
axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k16))])
N=32;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)32 点采样
评论0