24小时热门版块排行榜    

查看: 2470  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ammazn

铁虫 (初入文坛)

[求助] 关于matlab中fft之后变换值的单位问题

研究连续信号  f(t)=2exp(-3t)  (t>=0) 傅立叶变换。
1,解析解                   F(ω)=2/(3+jω)
2,利用FFT,程序如下


%%%模拟采样%%%
N=128;%取2的幂次方个点
t=linspace(0,3,N);%t在[0,3]间取N个点
f=2*exp(-3*t);%人为模拟N个采样点
Fk=fft(f);%FFT变换
Ts=t(2)-t(1);%计算t的间隔
Fwfft=Fk(1:N/2+1)*Ts;%仅从Fk中取正频率分量,并乘采样间隔计算F(w)
%%%解析式%%%
Ws=2*pi/Ts;%时间t等差值作周期
W=Ws*(0:N/2)/N;%建立频率轴,始于0,终于Nyquist频率ws/2
Fw=2./(3+j*W);%解析结果
%%%绘图%%%
plot(W,abs(Fwfft),‘:b’);%abs()复数模
hold on
plot(W,abs(Fw),'-r');%解析图
xlabel('频率(rad/s)');ylabel('|F(w)|');



标记红字体的那行,为了跟解析结果比较,为什么对f进行傅里叶变换后,要乘以采样间隔时间呢???对DFT变换不是很了解, 万永革主编 《数字信号处理的MATLAB实现》之中提到“为了与真实振幅对应,需要将变换后结果乘以2除以N。 ”这个程序中乘以采样间隔时间和“乘以2除以N ”两者怎么联系起来呢?新手没接触过这个,还请大家不吝赐教啊,非常感谢!

[ Last edited by ammazn on 2013-7-18 at 21:03 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ammazn

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by justyet at 2013-07-19 13:24:21
"标记红字体的那行,为了跟解析结果比较,为什么对f进行傅里叶变换后,要乘以采样间隔时间呢???"
其实任何一列数都可以做fft,无论时间间隔是多少。事实上,fft也不关心时间间隔。
如果没有时间间隔 ...

非常感谢!
5楼2013-07-20 13:16:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

justyet

金虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2013-07-24 21:44:51
"标记红字体的那行,为了跟解析结果比较,为什么对f进行傅里叶变换后,要乘以采样间隔时间呢???"
其实任何一列数都可以做fft,无论时间间隔是多少。事实上,fft也不关心时间间隔。
如果没有时间间隔,相当于默认两个点之间的间隔是一秒。

“乘以2除以N“  说的是另外一回事,跟做fft变换的数据的个数有关系,跟采样间隔没关系。
2楼2013-07-19 13:24:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_强鹏

银虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2013-07-24 21:44:59
你可以去百度上搜下fft的物理意思;TS的倒数是相当于基频,W=Ws*(0:N/2)/N;我感觉不要再除以N了,fft之后Fwfft=Fk(1:N/2+1)*Ts不用再乘以TS;Fwfft算出来的是频率,而W算出来的是弧度,2个相差一个2*pi
3楼2013-07-19 22:54:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

simazairu

新虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2013-07-24 21:45:08
标记红字体的那行,为了跟解析结果比较,为什么对f进行傅里叶变换后,要乘以采样间隔时间呢??
     请比较在matlab语言傅里叶变换的离散表达式与高等代数中的傅里叶表达式中的异同。仔细点。
     你会发现:前者表达式中没有乘上dt.为与定义相同,必需乘上该因子。
4楼2013-07-19 23:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见