完成实验03 快速傅立叶变换.docVIP

实验三 快速傅立叶变换 一、实验目的 1．掌握fft函数的用法。 2．了解FFT的用途。 二、实验内容 1．利用FFT进行信号检测。分析信号频谱所对应频率轴的数字频率和频率轴的频率范围。 2. 对例2，进一步增加截取长度和FFT点数，如N加大到256，观察信号频谱的变化，分析产生这一变化的原因。在截取长度不变的条件下改变采样频率，观察信号频谱的变化，分析产生这一变化的原因。 3. 对例3，加大噪声到2*randn(1,N)和8*randn(1,N)，画出并比较不同噪声下时域波形和频谱。 4．比较DFT和FFT的运算时间。(计时函数 tic， toc) 5．对给定语音信号进行谱分析，写出采样频率，画出语音信号的波形及频谱，并分析语音信号的频率分布特点。 三、实验例程 %例1 求信号功率谱 实验结果与实验流程图如下图示： %例2 模拟信号，以进行取样，求N点DFT的幅值谱 （1）%N分别为 (1) N=45 (2) N=50 (3) N=55 (4) N=60 实验运行结果如下图示： （2）. 对例2，进一步增加截取长度和FFT点数，如N加大到256，观察信号频谱的变 化，分析产生这一变化的原因。 N=256;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); plot(q,abs(y));title(FFT N=256) 实验程序运行结果如右图示： 产生这一变化原因为： （3）.在截取长度不变的条件下改变采样频率，观察信号频谱的变化，分析产生变化的原因。 改编程序如下： N=256;n=0:N-1;t=0.1*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); plot(q,abs(y));title(FFT N=256) 实验运行结果如右图示： 信号频谱产生这一变化的原因: 实验流程图如下图示： %例3 参数同上，N取64，并在信号中加入噪声w(t),比较有无噪声时的信号谱 %由结果可以看出这种噪声不影响信号检测 （1）加大噪声到8*randn(1,N) figure(2) subplot(2,1,1) N=64;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); plot(q,abs(y));title(FFT N=64) subplot(2,1,2) N=64;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t) +0.8*randn(1,N); y=fft(x,N); plot(q,abs(y)); title(FFT N=64(with noise)) 实验运行结果如右图示： （2） 对例3，加大噪声到2*randn(1,N)和8*randn(1,N)，画出并比较不同噪声下时域波形和频谱。 加大噪声到2*randn(1,N) figure(2) subplot(2,1,1) N=64;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); plot(q,abs(y));title(FFT N=64) subplot(2,1,2) N=64;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t) + 2*randn(1,N); y=fft(x,N); plot(q,abs(y)); title(FFT N=64(with noise)) 实验运行结果如右图示： 比较两种不同噪声频谱与有无噪声时的信号谱并分析结果 4．比较DFT和FFT的运算时间。(计时函数 tic， toc) N分别取256，512，1024，2048，4096，… 程序如下： N=256; N=4096; x=randn(1,N); tic y=dft(x,N); toc tic z=fft(x); toc N=256：Elapsed time is 0.172000 seconds. Elapsed time is 0.015000 seconds. N=512：Elapsed time is 0.687000 seconds. Elapsed time is 0.000000 seconds. N=1024：Elapsed time is 3.031000 seconds. Elapsed time is 0.047000 seconds. N=2048：