实验二用fft对信号作频谱分析

《实验二用fft对信号作频谱分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验二用FFT对信号作频谱分析一、实验目的学习用FFT对连续倌号和时域离散倌号进行谱分析的方法，了解可能出现的分析误差及K原因，以便正确应用FFT。二、实验步骤及内容1.对以下序列进行谱分析：Xi(n)=R4(n)n+10<77<3x2(n)=J8-z?4

2、,8);xl6k=fft(xl,16);yl=abs(x8k);y2=abs(xl6k);y21=angle(x8k);y22=angle(xl6k);subplot(2,2,l);stem(yl,'.’);程序:%产生矩形信号R4(n)%求系统8位FFT变换%求系统8位FFT变换%8点幅频特性%16点幅频特性%8点相频特性%16点相频特性subplot(2,2,2);stem(y21/.');title('8点相频特性*);subplot(2,2,3)；stem(y2/.');title('16点幅频特性gr

3、idon;（2）程序:closeall;clearallnl=0:3;xl=nl+l;n2=4:7;x2=8-n2;x=[xl,x2];x8k=fft(x,8);xl6k=fft(x,16);yl=abs(x8k);y2=abs(xl6k);y21=angle(x8k);y22=angle(xl6k);subplot(2,2,l);stemfyl/J);%求系统8位FFT变换%求系统8位FFT变换%8点幅频特性%16点幅频特性%8点相频特性%16点相频特性title('8点相频特性*);subplot(2,2,

4、3);stem(y2；.');title('16点幅频特性*>;subplot(2,2,4);stem(y22/.');titleCl6点和频特性*);gridon;%产生矩形信号R4(n)x8k=fft(xl,8);%求系统8位FFT变换%8点幅频特性%16点幅频特性%8点相频特性%16点相频特性xl6k=fft(xl,16);%求系统8位FFT变换yl=abs(x8k);y2=abs(xl6k);y21=angle(x8k);y22=angle(xl6k);subplot(2,2,l);stem(yl；.'

5、);titleC8点幅频特性*);subplot(2,2,2);stem(y21/.');title('8点相频特性subplot(2,2,3)；stem(y2/.');subplot(2,2,4);stem(y22;?};title('16点相频特性•);gridon;(3)程序:closeall;clearallnl=0:3;xl=4-nl;n2=4:7;%8点幅频特性%16点幅频特性%8点相频特性%16点相频特性%求系统8位FFT变换%求系统8位FFT变换x2=n2-3;x=[xlzx2];x8k=fft

6、(x,8);xl6k=fft(x,16);yl=abs(x8k);y2=abs(xl6k);y21=angle(x8k);y22=angle(xl6k);subplot(2,2,l)；stem(yl/.');title('8点幅频特性')；subplot(2,2,2);stem(y21/.');title('8点相频特性*);subplot(2,2,3);stem(y2/.');title('16点幅频特性•>;subplot(2,2,4);stem(y22;?};title('16点相频特性•)；gridon

7、;%产生妍形信号R4(n)x8k=fft(xl,8);%求系统8位FFT变换xl6k=fft(xl,16);%求系统8位FFT变换yl=abs(x8k);%8点幅频特性y2=abs(xl6k);%16点幅频特性y21=angle(x8k);%8点和频特性y22=angle(xl6k);%16点相频特性subplot(2,2,l);stem(yl/.');subplot(2,2,2);stem(y21/.');title('8点相频特性*);title('16点相频特性gridon;讨论：第一行说明xUn卜R(n

8、＞的8点DFT和16点DFT分别是xjn｝的频谱函数的8点和16点采样；x3n)=x2((n+3))8R8(n),,则8点DFT相同，当N=16时，x2(n)和x3(n)不满足循环移位，则16点DFT不同。TTTTJT1.对以下周期序列进行谱分析：x4(n)=cos—n,x5(n)=cos—n+cos—n448选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。分别打印