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

uperrua

Matlab 傅里叶变换傅里叶逆变换-FFT-IFFT

' P8 [; w( W7 A* j9 X- x%% 信号经过傅里叶变换然后进行傅里叶逆变换后信号的变化
clear   all;   clc;
, _" |" e. A3 @( R%------Author&Date------
%Author:
3 x2 D9 }- n+ f4 c6 G%Date: 2013/07/31
# I$ c- w4 z( L5 \( I6 g6 h( W%==========================================================================
Fs=8e3;                      %采样率
9 ^- ?9 {1 d- {$ @& \5 b) ?t=0:1/Fs:1;                       %采样点
9 d: y( L+ f& T2 g7 Y  h- m0 N, Slen=length(t);                       %采样长度
# @& @" [" I( R+ X7 I" t. N' af1=10;                       %频率 1
f2=100;                      %频率 2
f3=1000;                      %频率 3
A1=1;                       %幅度 1
- q& C, S8 v) _, JA2=0.8;                      %幅度 2
/ h  c( _- V6 M# tA3=0.3;                       %幅度 3
MaxS=A1+A2+A3;                      %信号幅度的最大值
! S& P4 u+ Z( I* usignal=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t)+A3*sin(2*pi*f3*t);
6 ^# X* P: [) _/ m) L1 _X=fft(signal,len);                       %傅里叶变换
6 |  d% M* A7 Z, R+ fmagX=abs(X);                       %信号的幅度
angX=angle(X);                       %信号的相位
Y=magX.*exp(1i*angX);                       %信号的频域表示
( O  Y  v* K& G( D: p: Cy=ifft(Y,len);                       %信号进行傅里叶逆变换
y=real(y);
* o6 G* _- |8 G6 l8 r& Her=signal-y;                       %原始信号和还原信号的误差
, {. S$ T4 B' l' e4 r1 Jsubplot(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( ' 还原信号 ' );
subplot(313);plot(t,er);xlabel( ' 时间 ' );ylabel( ' 振幅 ' );title( ' 误差 ' );
5 ]8 ^- a/ D; \/ ?! g% End Script
$ N7 h5 i6 ^0 M2 I6 x* t, C
2 X: t5 b7 x- E+ V

0 k1 L& T$ D  [0 w. @( C4 J+ G( X! g3 c
误差er的数量级为 ，误差绝对值的最大值（ ）如下图：

. z. |( p5 `! L* Y




% M3 H. [$ Z! G7 w8 G! n7 W" P

WNSKS

谢谢分享