MATLAB中ifft函数用法、性质、特性-以及与fft的组合应用全面深入解析(含程序).docx

_MATLAB中ifft函数用法、性质、特性，以及与fft的组合应用全面深入解析（含程序）前言在我之前的MATLAB中fft函数用法、性质、特性、缺陷全面深入解析（含程序）中，我已经详细阐述了fft的所有性质，大家应该能够正确使用fft来获得合适的信号频谱图。鉴于网上的程序代码基本都只介绍了使用fft实现信号时域到频域的转换，很少有介绍使用ifft将频域信号转换为时域信号，尤其是对同一信号进行连续fft和ifft的控制和变换，更不用说网上的很多程序都有很大局限性，考虑不周全，直接拿来用会出很多问题。所以这里我在前面那篇文章的基础上，利用程序全面讲解fft和ifft的组合用法，来保证信号的正确变换，同时说明如何调整参数来减少失真，降低误差。让大家能够将信号在时域频域间自由轻松变换，同时根据需要修改参数，来获得所需的时域或频域信号。ifftifft是fft的逆运算，也就是将频域信号反变换为时域信号，要使用ifft，就必须先清除fft后的信号频谱的结构，大家应该先看我之前的那篇MATLAB中fft函数用法、性质、特性、缺陷全面深入解析（含程序）就知道了。我们需要先将频域值变换为最初的fft后的结构，再使用ifft，才能经过变换，直接获得正确的信号时域值。这里先说明一下，对信号进行fft后，如果直接进行ifft，那么可以直接还原时域信号，误差可以忽略，这个我就不演示了。我这篇文章要讲解的是在利用fft变换后，又进行了各个运算步骤获得了最后的正确的频域图后，如何利用这个频域图进行反运算，最后利用ifft获得正确的原来的时域图。或者是已知一个信号的准确的频域图，如何进行正确运算，利用ifft获得正确的时域图。ifft的定性解析首先我通过程序给大家展示一下fft和ifft的组合用法，下面的程序演示的是fft后运算获得正确的频域图后，再反变换利用ifft获得原来的时域图。最后测量输出恢复的时域信号与原时域信号的误差。组合程序解释如下：%这个程序用于初步测试ifft的正确性clear all;clc;Fs = 1e3; %这是原始时域信号x的采样频率multiple=0; %x进行fft前的补零倍数，可以增加频率分辨率，也会影响信号频域值time=10;t = 0:1/Fs:time-1/Fs; %这个会影响x的采样点数，当然点数越多，频域越精确越平滑%x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*202.5*t)+2*sin(2*pi*1450*t);x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*200*t);%x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*200*t)+2*cos(2*pi*0.1*t);%频率最小分辨率x=x,zeros(1,length(x)*multiple); freqres = Fs/length(x); %变换运算获得正确频谱图xdft = fft(x); xdft2 = xdft(1:length(x)/2+1);xdft2 = 1/length(x).*xdft2;xdft2(2:end-1) = (1+multiple)*2*xdft2(2:end-1);freq = 0:Fs/length(x):Fs/2;plot(freq,abs(xdft2);%利用频谱图恢复原来的fft图xdft3=xdft2;xdft3(2:end-1)=xdft2(2:end-1)./2;xdft3=xdft3.*length(x);for j=length(x)/2+2:length(x)xdft3(j)=conj(xdft3(length(x)-j+2);enderror2=xdft3-xdft; %原来的fft图与恢复后的fft图的误差subplot(4,1,1);plot(real(xdft);title(最初的fft图实部);subplot(4,1,2);plot(real(xdft3);title(恢复后的fft图实部);subplot(4,1,3);plot(imag(xdft);title(最初的fft图虚部);subplot(4,1,4);plot(imag(xdft3);title(恢复后的fft图虚部);%现在利用恢复的fft图进行ifft反变换back=ifft(xdft3);error=x-back; %测量恢复的时域信号与原来的时域信号的误差min(error) max(error) %输出误差的最大值ans = -2.6645e-015ans = 2.6645e-015ifft的定量详细测试与解析利用程序解释如下：%这个程序用于全面测试ifft的通用程序的正确性%结果完全正确，完整展示了如何使用fft,以及正确利用频域的值进行ifft变换%注意ifft后的时域信号的时间起点可以是任意的，不影响频域；同时必须为实信号时，才能保证此程序的ifft变换方法的正确性%注意在特殊点，即第1、L/2+1点的振幅是实际值的两倍，且虚部为零，不过实际的频域信号可能虚部并不为零，所以会有误差，-%且时域信号有补零时也不能用此程序，除非你实现知道补了多少个零，否则误差很大。clear all;clc;Fs = 1e4; %这是原始时域信号x的采样频率multiple=0; %这里补零的个数是零，才能进行正确的ifft变换time=2;t = 0:1/Fs:time-1/Fs; %这个会影响x的采样点数，当然点数越多，频域越精确越平滑%x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*202.5*t)+2*sin(2*pi*1450*t);x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*200*t);%x = 1.5*cos(2*pi*100*t)+3*sin(2*pi*200*t)+2*cos(2*pi*0.1*t);%频率最小分辨率x=x,zeros(1,length(x)*multiple); freqres = Fs/length(x); %fft变换与横纵坐标的运算伸缩xdft = fft(x); xdft2 = xdft(1:length(x)/2+1);xdft2 = 1/length(x).*xdft2;xdft2(2:end-1) = (1+multiple)*2*xdft2(2:end-1);freq = 0:Fs/length(x):Fs/2;plot(freq,abs(xdft2);%反运算%现在我们只利用fft变换获得的正确的最终频域的值，对最终频域进行反变换与横纵坐标的反伸缩，最终进行ifft%注意这里的频域的值是正的那一半，且必须保证没有频谱混叠xdft3=xdft2;L=(length(xdft3)-1)*2;%L=length(x)resolution=freq(2)-freq(1);%频率分辨率Fs2=resolution*L;%时域的采样率xdft3(2:end-1)=xdft3(2:end-1)./2;%如果fft有补零的话，这里就不是除以2了-%而是除以(1+multiple)*2，但是multiple也就是补零倍速并不知道怎么估算，因此这个-%程序对于时域有补零的情况误差很大xdft3=xdft3.*L;for j=L/2+2:Lxdft3(j)=conj(xdft3(L-j+2);end%开始ifft反变换-back=ifft(xdft3);back=real(back);%信号频域对称，时域应该是实信号，但反变换后有误差虚部存在，需要去掉虚部。error=x-back;max(abs(error),max(real(error),max(imag(error),min(imag(error),min(real(error)%输出误差值subplot(2,1,1);plot(t,x);title(原来的时域信号);subplot(2,1,2);plot(0:1/Fs2:(L/Fs2-1/Fs2),back); title(恢复的时域信号);ans = 3.9968e-015ans = 3.9968e-015ans = 0ans = 0ans = -3.5527e-015ifft变换的特性与问题通过理论分析，和我的其他程序的实验验证，ifft和fft的各项参数设置会产生如下特性：(1). 通过提高频域补零倍数multiple也就是提高频域点数L，能够提高时域的采样率，也就是时域分辨率（影响曲线平滑度），但不会改变时域时间长短。反之同理。(2). 通过降低频域的采样率，也就是频域的分辨率，可以降低时域的时间长短，但不改变时域采样率。反之同理。解释如下：频域补零，能够减小栅栏效应，观察到的时域的点数更多，时域曲线更平滑。频域采样率越高，点数越多越复杂，原始的时域信号就需要更长的时间序列来代表这些多余的频率。最后注意fft和ifft的横轴的单位都要转换为Hz，而不是角速度w。那些实验验证程序我就不贴出来了，太繁杂，上面的程序已经是可以通用的程序，大家可以根据需要，参考上面的几种变换特性，自行修改参数使用，来获得任意精度和形状的时频图。欢迎批评、指正。 2015年5月23日BY XUHAO百度账号：xuhaobestone2邮箱：xuhaobestoneWelcome ToDownload !欢迎您的下载，资料仅供参考！精品资料