看看这篇Matlab 傅里叶变换傅里叶逆变换-FFT-IFFT

uperrua

Matlab 傅里叶变换傅里叶逆变换-FFT-IFFT
+ r. g# m" l3 ?
( Z2 r; e# q% A; d$ F! j/ I# j( I%% 信号经过傅里叶变换然后进行傅里叶逆变换后信号的变化
7 p: a0 r- p! M6 Q& d7 K: i  m4 J/ Rclear   all;   clc;
%------Author&Date------
%Author:
8 A1 L/ R& _# w6 U/ F%Date: 2013/07/31
%==========================================================================
Fs=8e3;                      %采样率
( o0 ?1 P4 b: P$ \0 It=0:1/Fs:1;                       %采样点
( y! d; z" `0 j' P9 ulen=length(t);                       %采样长度
f1=10;                       %频率 1
* V- f* W4 p/ ]% ~! ]* rf2=100;                      %频率 2
, p) l) @& l% _2 R* ~6 x9 {* `( Qf3=1000;                      %频率 3
A1=1;                       %幅度 1
A2=0.8;                      %幅度 2
A3=0.3;                       %幅度 3
MaxS=A1+A2+A3;                      %信号幅度的最大值
signal=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t)+A3*sin(2*pi*f3*t);
0 M6 o, C  d4 l4 ]  e4 [X=fft(signal,len);                       %傅里叶变换
- W6 c- a3 E3 M, y" P; b# [magX=abs(X);                       %信号的幅度
angX=angle(X);                       %信号的相位
& H3 `( Z/ g! ~) R( K  i& m. NY=magX.*exp(1i*angX);                       %信号的频域表示
y=ifft(Y,len);                       %信号进行傅里叶逆变换
y=real(y);
er=signal-y;                       %原始信号和还原信号的误差
, P; T7 n( y( O- M0 z; xsubplot(311);plot(t,signal);axis([0 1 -MaxS MaxS]);xlabel( ' 时间 ' );ylabel( ' 振幅' );title( ' 原始信号 ' );
subplot(312);plot(t,y);axis([0 1 -MaxS MaxS]);xlabel( ' 时 间 ' );ylabel( ' 振 幅 ' );title( ' 还原信号 ' );
* u% D+ o4 B' ^1 p* z' |) ssubplot(313);plot(t,er);xlabel( ' 时间 ' );ylabel( ' 振幅 ' );title( ' 误差 ' );
% End Script


; X9 A. O3 h2 L6 _

误差er的数量级为 ，误差绝对值的最大值（ ）如下图：
( @9 z/ F, [) B; h/ Y8 E
7 f5 g) d- i: v- K

3 O: P1 f% \* z! `7 T6 c

0 ?7 V2 s" H7 Q
/ ^1 z$ }0 o7 q0 Q5 U

WNSKS

谢谢分享