快速傅里叶变换与信号频谱分析-赵恒伟.docVIP

快速傅里叶变换与信号频谱分析 实验报告 [数字信号处理实验一] 实验一 快速傅里叶变换与信号频谱分析实验报告 一． 实验目的 1. 在理论学习的基础上，通过本实验加深对离散傅里叶变换的理解。 2. 熟悉并掌握按时间抽取编写快速傅里叶变换（FFT）算法的程序。 3. 了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如频谱混淆、 泄漏、栅栏效应等，以便在实际中正确使用FFT 算法进行信号处理。 二． 实验内容 1. 仔细分析教材第六章“时间抽取法FFT 的FORTRAN 程序”，编写出相应的 使用FFT 进行信号频谱分析的Matlab 程序。 2. 用FFT 程序分析正弦信号，分别在以下情况进行分析，并讨论所得的结果 三． 实验过程记录a) 信号频率F＝50Hz，采样点数N=32，采样间隔T=0.000625s； ，满足：这里k=1FFT采样采到频域谱峰单线谱，恰好位于频谱Hz的位置。FFT之后的频域是的频谱 b) 信号频率F＝50Hz，采样点数N=32，采样间隔T=0.005s； ，满足：这里k=8FFT采样采到频域谱峰单线谱，恰好位于频谱Hz的位置FFT之后的频域是的频谱 c) 信号频率F＝50Hz，采样点数N=32，采样间隔T=0.0046875s； ，不满足：这里k=7.5FFT采样采到频域谱峰一定相位差，FFT频域频谱混。 d) 信号频率F＝50Hz，采样点数N=32，采样间隔T=0.004s； ，不满足：这里k=6.4FFT采样采到频域谱峰一定相位差，FFT频域频谱混。 e) 信号频率F＝50Hz，采样点数N=64，采样间隔T=0.000625s； ，满足：这里k=2FFT采样采到频域谱峰单线谱，恰好位于频谱Hz的位置FFT之后的频域是的频谱 f) 信号频率F＝250Hz，采样点数N=32，采样间隔T=0.005s； ，满足：这里k=40FFT采样采到频域谱峰单线谱，位于频谱Hz的位置FFT之后的频域的频谱频率不满足内奎斯特采样，发生欠采样，单峰不会在频域幅度上重叠，的重叠。显示的是FFT的频域是重叠频率为=50Hz，和=150Hz的频谱图=8,150/6.25=24，表示单峰的位置与图相符）。 g) 将c)中信号后补32 个0，做64 点FFT 将c）直接采样64 个点做FFT ，满足：这里k=15FFT采样采到频域谱峰单线谱 补零之后做FFT，补零之前的有较大不同补零之后原周期，，但是包络一致。补零FFT分辨率一样，FFT补零之后循环移位与FFT有差别，频谱不同。. 实验总结： 【一】a)、b)与c)，d)，FFT得到的频谱有明显不同，a)b)出现单峰，主要是因为峰值，而c)d)受栅栏效应影响，在其他频率也有非零值。可以是否单峰决定于是否满足 【二】信号与采样及FFT变换参数，可以发现频率F决定信号本周期，决定采样频率，采样频率需满足内奎斯特采样定律。采样点之间的频谱宽度，决定频谱单峰的位置1) 在实验a)、b)、c)和d)中，正弦信号的初始相位对频谱图中的幅度特性 是否有影响？为什么？ 答：不影响，频率相位变化率，与初始值无关。2) 信号补零后做FFT 是否可以提高信号频谱的分辨率？为什么？ 答：通过补零可以观察到更多，但是不意味着补零能真正提高分辨率，因为xn]只是x（）主值序列，而n]末尾补零得到x’n],周期延拓之后，更不是（）采样有意义的数据并未增加的真正提高只能在时域采样时增加采样长度和数量来实现。 五． 实验代码 Matlab实现 clc;clear; f=input(正弦频率) N=input(采样点数) T=input(采样间隔) choise=input(是否补0 yes=1/ no=0 ) n=0:T:(N-1)*T; %采样点 k=0:N-1; x=sin(2*pi*f*n); %抽样后的离散信号 if choise==1 e=input(请输入要补入0的个数) x=[x zeros(1,e)] %将所加0存入x矩阵 N=N+e; else end %加零 k=0:N-1; %给k重新赋值，因为有可能出现加零状况 bianzhi=bi2de(fliplr(de2bi(k,leng