资源描述：
《实验一FFT与谱分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《数字信号处理实验报告》姓名：朱道萍班级：通信1401学号：201403090129实验一FFT与谱分析1•实验口的⑴加深对DFT算法原理和物理意义的理解。⑵熟悉FFT算法，增强对FFT结果的分析能力。⑶拿握用FFT对连续时间信号和离散吋间信号进行谱分析的方法，了解误差及其产生的原因。2.实验原理⑴复习DFT的定义、性质和用DFT作谱分析的有关内容。⑵复习FFT算法原理与编程思想。(3)查询Matlab中fft函数的用法,学会通过改变fft的参数计算不同DFT。⑷编制信号产生子程序，产生一下信号:[I刃=0,】23u.w.n-t-18-/r0n-30ow
2、30・W・xM=cosxjwjssinnnTxj/）=cos（S^/）4-cos（l6刃）+cos（2C/rZ）巧同=兀同+巧同心仏卜{1A3,4QO,OQOQOQO,0,0,0}x10H={1.23・4・1.23.4儿234丄234}{】・OQ020・0・0・3.0QO・4・000}⑸编写实验程序:参考图14所示的程序流程图，编写实验程序。读入长度N调用信号产生子程序产生实鲨信号3.实验内容对2中所给出的11个信号逐一进行谱分析。各信号的FFT变换区间N、连续时间信号x6(t)的采样频率fs等参数的参考取值如下：a)xl[n],x2[n],x3[n],
3、x4[n],x5[n],x7[n],x8[n]:N=8,16b)x6(t):fs=64Hz,N=16,32,64c)x9[n],xlO[n],xll[n]:N=16xl[n]:N=8n=0:N-l;xl=[l1110000];m=0:2*N-1;x2=[l11100000000fl=fft(xlzN);f2=fft(x2z2*N);subplot(2,2Z1)stem{n,xl)subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2)subplot(2,2,4)stem(mzabs(f2))x2[n]:N
4、=8n=0:N-l;xl=[l234432m=0:2*N-1;x2=[l234432fl=fft(xlzN);f2=fft(x2,2*N);subplot(2,2,1)stem(n,xl)subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)1]；100020Q20(3—15?1051510°0GoOggqGgQo5108stem(mzx2)subplot(2,2,4)stem(m,abs(f2))x3[n]:N=8n=0:N-l;15105°°-xl=[43211234];m=0:2*N-1;x2=[432112340f
5、l=fft(xl,N);f2=fft(x2,2*N)subplot(2,2,1)stem(n,xl)0526
6、<>subplot(2,2,2)stem(n,abs(f1))subplot(2,2,3)stem(m,x2)subplot(2#2#4)stem(mzabs(f2))x4[n]:N=8;n=0:N-l;xl=cos(0.25*pi*n);m=0:2*N-1;x2=cos(0.25;fl=fft(xlzN);f2=fft(x2,2*N);subplot(2,2,1)stem(n,xl)subplot(2,2,2)stem(n,abs(f1))su
7、bplot(2,2,3)stem(mfx2)subplot(2,2,4)stem(m,abs(f2))x5[n]:N=8;n=0:N-l;xl=sin(0.125*pi*n);m=0:2*N-1;x2=sin(0・125*pi*m);fl=fft(xl,N);f2=fft(x2,2*N);subplot(2,2,1)stem(n,xl)subplot(2,2,2)stem(nfAbs(f1))subplot(2,2,3)IV05005-1r—<>Cj><>>024681<>)I5101542QCj)0丨1^
8、0009102468stem(m,x2)subplot(2#2#4)stem(mzabs(f2))x6[n]:N=16;fs=64;n=0:N-l;xl=cos(8*pi*n/fs)+cos(16*pi*n/fs)+cos(20*pi*n/fs);m=0:2*N-1;x2=cos(8*pi*m/fs)+cos(16*pi*m/fs)+cos(20*pi*m/fs);1=0:4*N-1;x3=cos(8*pi*l/fs)+cos(16*pi*l/fs)+cos(20*pi*l/fs);fl=fft(xl,N);f2=fft(x2,2*N);f3=fft(x
9、3z4*N);subplot(3,2,1)stem(nfxl)subplot(3
