Matlab中的FFT使用说明.docVIP

下载本文档

10
0
约1.42千字
约 4页
2019-12-09 发布于河北
举报
版权申诉

Matlab中的FFT使用说明.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

FFT是Fast Fourier Transform(快速傅里叶变换)的简称，FFT算法在MATLAB中实现的函数是Y=fft(x,n)。刚接触频谱分析用到FFT时，几乎都会对MATLAB的fft函数产生一些疑惑，下面以看一个例子（根据MATLAB帮助修改）。 Fs = 2000; % 设置采样频率 T = 1/Fs; % 得到采用时间 L = 1000; % 设置信号点数，长度1秒 t = (0:L-1)*T; % 计算离散时间， % 两个正弦波叠加 f1 = 80; A1 = 0.5; % 第一个正弦波100Hz，幅度0.5 f2 = 150; A2 = 1.0 ; % 第2个正弦波150Hz，幅度1.0 A3 = 0.5; % 白噪声幅度； x = A1*sin(2*pi*f1*t) + A2*sin(2*pi*f2*t); %产生离散时间信号； y = x + A3*randn(size(t)); % 叠加噪声； % 时域波形图 subplot(2,1,1) plot(Fs*t(1:50),x(1:50)) title(Sinusoids Signal) xlabel(time (milliseconds)) subplot(2,1,2) plot(Fs*t(1:50),y(1:50)) title(Signal Corrupted with Zero-Mean Random Noise) xlabel(time (milliseconds)) NFFT = 2^nextpow2(L); % 设置FFT点数，一般为2的N次方，如1024，512等 Y = fft(y,NFFT)/L; % 计算频域信号， f = Fs/2*linspace(0,1,NFFT/2+1); % 频率离散化，fft后对应的频率是-Fs/2到Fs/2,由NFFT个离散频点表示 % 这里只画出正频率； % Plot single-sided amplitude spectrum. figure; plot(f,2*abs(Y(1:NFFT/2+1))); % fft后含幅度和相位，一般观察幅度谱，并把负频率加上去， title(Single-Sided Amplitude Spectrum of y(t)) xlabel(Frequency (Hz)) ylabel(|Y(f)|) 运行结果时域波形图如图所示：幅度谱如下：由图可见，80Hz的信号幅度为0.4762，频率为80.08，150Hz的信号频率为150.4，幅度0.9348，存在误差。去掉白噪声后，结果为：可见幅度精度提高了，但频率任然有误差。如果修改Fs为2048，L为4096，此时没有补0产生的复查，也没有噪声污染，结果如下：可见结果与信号完全一致。