24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2215  |  回复: 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的回帖

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的回帖

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的回帖

ammazn

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by _强鹏 at 2013-07-19 22:54:03
你可以去百度上搜下fft的物理意思;TS的倒数是相当于基频,W=Ws*(0:N/2)/N;我感觉不要再除以N了,fft之后Fwfft=Fk(1:N/2+1)*Ts不用再乘以TS;Fwfft算出来的是频率,而W算出来的是弧度,2个相差一个2*pi

非常感谢你的帮助!
6楼2013-07-20 13:16:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ammazn

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by simazairu at 2013-07-19 23:27:25
标记红字体的那行,为了跟解析结果比较,为什么对f进行傅里叶变换后,要乘以采样间隔时间呢??
     请比较在matlab语言傅里叶变换的离散表达式与高等代数中的傅里叶表达式中的异同。仔细点。
     你会发现:前 ...

多谢多谢,我得好好看基础的理论!
7楼2013-07-20 13:17:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ammazn 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械专硕274求调剂,不挑专业学校 +4 泛泛2333 2026-04-05 5/250 2026-04-05 17:13 by lbsjt
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊!
[考研] 290求调剂085701 +10 1314捧花 2026-04-02 10/500 2026-04-05 10:19 by Sealedwind
[考研] 一志愿郑州大学材料与化工085600,求调剂 +24 吃的不少 2026-04-02 24/1200 2026-04-04 23:20 by 永字号
[考研] 283求调剂 +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
[考研] 083000环境科学与工程调剂,总分281 +4 橙子(胜意) 2026-03-30 4/200 2026-03-31 00:44 by Linzejun
信息提示
请填处理意见