|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
比如 时域信号 exp(-(t.^2)),波形图是:
. Y. d7 J% K9 g S* h( O W" }- s; G. J7 J! k
对它做fft变换,求得频谱是:# X& j2 h' ~, J
5 z2 K" k6 V( c L; k. j另外,我知道exp(-(t.^2)) 的傅里叶变换的公式是sqrt(pi).*exp(-(pi.*freq).^2),直接仿真该公式:
, Q% H O( a/ l+ C8 n' e9 J3 T! r9 ]) V* k$ r
再对其 ifft得到时域波形是:
0 P- F/ v' p4 U
# G0 w( h; E+ O) y3 n7 y8 T3 v不管时域波形和频域波形,都不一致啊,为什么呢?* \# z* O, L0 k& ?9 E
, {; b; |; L+ T; d4 Vclear all;close all;clc;6 Y1 G" e- N, W: q4 D/ a
t=linspace(-10,10,1024);6 W: @$ G& M Y+ |% s. M
y=exp(-(t.^2));
7 }$ ?. @4 u2 I% y4 wfigure(1);. F0 y. V' M4 Q2 s3 Y" j
plot(t,y);
; X7 e0 ?% j3 C/ s$ yfigure(2);
5 I1 X9 E2 m) W% b. @y_spec=fft(y);
$ F7 m$ W0 _1 R1 p9 `freq=linspace(-1/10,1/10,1024);
: d" j7 V4 J% ?; ~2 d5 Zplot(freq,fftshift(abs(y_spec)));# I% H: k# | K2 e( Y( p
figure(3);9 s6 O' n' n$ L) Z
z_spec=sqrt(pi).*exp(-(pi.*freq).^2);& A( v$ x8 d5 L
plot(freq,z_spec);
* m9 s6 l, p2 w" O4 Ufigure(4);/ u9 m0 w: s2 b+ o( D7 Q) u* \8 t" u
plot(t,fftshift(ifft(z_spec)));
& i- q& S7 Y& \& m* T |
|