实验二FFT频谱分析及应用.doc

实验二：FFT频谱分析及应用一、实验内容 （1）模拟信号x=2*sin(4*pi*t)+5*cos(8*pi*t);以t=0.01*n进行采样，求出N=40点FFT的幅度频谱，提高到N=128，再求该信号的幅度频谱；（2）一个连续信号含三个频谱分量，经采样得到以下序列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，分别求N=64，N=128时其频谱特性。二、实验代码代码一：模拟信号x=2*sin(4*pi*t)+5*cos(8*pi*t);以t=0.01*n进行采样，求出N=40点FFT的幅度频谱clc;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(single x(n);subplot(2,1,2);plot(w/pi,magX);title(FFT N=40);xlabel(f (unit:pi);ylabel(|X|);grid代码二：模拟信号x=2*sin(4*pi*t)+5*cos(8*pi*t);以t=0.01*n进行采样，N=128点FFT的幅度频谱：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(single x(n);subplot(2,1,2);plot(w/pi,magX);title(FFT N=128);xlabel(f (unit:pi);ylabel(|X|);grid代码三：一个连续信号含三个频谱分量，经采样得到以下序列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，N=64时其频谱特性clc;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(single x(n);subplot(2,1,2);plot(w/pi,magX);title(FFT N=64);xlabel(f (unit:pi);ylabel(|X|);gridclc;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(single x(n);subplot(2,1,2);plot(w/pi,magX);title(FFT N=64);xlabel(f (unit:pi);ylabel(|X|);grid代码四：一个连续信号含三个频谱分量，经采样得到以下序列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，N=128时其频谱特性。clc;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(single x(n);subplot(2,1,2);p