数字信号处理 实验二 用DFT、FFT进行频谱分析

【实验内容及基本要求】：    

1.任务1：一个连续信号包含三个频谱分量，经采样得到以下序列x(n)：

1）N=64，df分别为1/16，1/64，观察其频谱；

2）N=64及128时，df=1/64，做128点FFT，观察其频谱，并与前述结果做比较。

2.任务2：对模拟信号x(t)以t=0.01n(n=0:N-1)进行采样，分析其在下面两种情况下的幅度频谱：

1）N=40点FFT的幅度频谱，从图中能否观察出信号的两个频谱分量？

2）提高采样点数，比如N=128，再求该信号的幅度频谱，此时幅度频谱发生了什么变化？原因是什么？

【实验过程及实验步骤】：（要求实验前作好预习，拟出实验步骤。）

  1.任务1：一个连续信号包含三个频谱分量，经采样得到以下序列x(n)：

  1）N=64，df分别为1/16，观察其频谱； 

%df为1/16代码如下：
N=64;n=0:N-1;
df=1/16;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
k=0:N/2;w=2*pi/N*k;
X=fft(x,N);
magX=abs(X(1:N/2+1));
subplot(2,1,1);stem(n,x,'.');
title('signal x(n)');
subplot(2,1,2);
plot(w/pi,magX);
title('FFT N=64');
xlabel('f(unit :pi)');
ylabel("TX1");
grid on;

结果如图：

 2）N=64，df分别为1/64，观察其频谱；

%df为1/64代码如下：
N=64;n=0:N-1;
df=1/64;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
k=0:N/2;w=2*pi/N*k;
X=fft(x,N);
magX=abs(X(1:N/2+1));
subplot(2,1,1);stem(n,x,'.');
title('signal x(n)');
subplot(2,1,2);
plot(w/pi,magX);
title('FFT N=64');
xlabel('f(unit :pi)');
ylabel("TX1");
grid on;

结果如图：

 3）N=128时，df=1/64，做128点FFT，观察其频谱。

  代码如下：

N=128;n=0:N-1;
df=1/64;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
k=0:N/2;w=2*pi/N*k;
X=fft(x,N);
magX=abs(X(1:N/2+1));
subplot(2,1,1);stem(n,x,'.');
title('signal x(n)');
subplot(2,1,2);
plot(w/pi,magX);
title('FFT N=128');
xlabel('f(unit :pi)');
ylabel("TX1");
grid on;

结果如图： 

N=64，df=1/64和N=128，df=1/64点FFT相比，其频谱由于采样点数的增加，频谱的细节也会增加，但同时也会降低频率分辨率。

 2.任务2：对模拟信号x(t)以t=0.01n(n=0:N-1)进行采样，分析其在下面两种情况下的幅度频谱：

 1）N=40点FFT的幅度频谱，代码如下：

N=40;n=0:N-1;t=0.01*n;
x=2*sin(4*pi*t)+5*cos(8*pi*t);k=0:N/2;w=2*pi/N*k;
X=fft(x,N);
magX=abs(X(1:N/2+1));
subplot(2,1,1);stem(n,x,".");
title('signal x(n)');
subplot(2,1,2);plot(w/pi,magX);
title("FFTN=40');
xlabel('f(unit:pi)');
ylable("|XI”);grid

 结果如图：

 2）提高采样点数，比如N=128，再求该信号的幅度频谱，此时幅度频谱发生了什么变化？原因是什么？

N=128;n=0:N-1;t=0.01*n;
x=2*sin(4*pi*t)+5*cos(8*pi*t);k=0:N/2;w=2*pi/N*k;
X=fft(x,N);
magX=abs(X(1 :N/2+1));
subplot(2,1,1);stem(n,x,".");
title('signal x(n)');
subplot(2,1,2);
plot(w/pi,magX);
title('FFT N=128');
xlabel('f(unit:pi)');
ylable('|XI');grid

结果如图：

提高采样点数会使幅度频谱更加精确和详细，但同时也会降低频率分辨率。由奈圭斯特采样定理所决定的，采样频率需要至少是信号中最高频率的两倍才能无失真地重建该信号