找回密码
 注册
关于网站域名变更的通知
查看: 285|回复: 3
打印 上一主题 下一主题

[仿真讨论] 为什么利用fft画出的频谱和直接利用公式画出的频谱不一致?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2022-7-21 11:25 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

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

该用户从未签到

推荐
发表于 2022-7-21 13:17 | 只看该作者
fft是离散的,你给的那个傅里叶变换公式是连续信号的变换,两个信号都不一样,所以变换不一样
! ]4 w( f+ h+ n5 j) |( Q# \" o

该用户从未签到

3#
发表于 2022-7-21 14:47 | 只看该作者
同问,求结果* ?6 ^% Z# R: B% N

该用户从未签到

4#
发表于 2022-7-21 15:37 | 只看该作者
再看看别人是怎么说的$ R/ N% k2 k5 H
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-3 14:13 , Processed in 0.156250 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表