24小时热门版块排行榜    

查看: 2784  |  回复: 5

liulbcd

铁虫 (小有名气)

[求助] 如何用Matlab对函数sin(pi*t)/(pi*t)进行FFT变换 已有2人参与

如题,本人现在想利用Matlab的FFT函数对sin(pi*t)/(pi*t)进行傅里叶变换,通过书本得到其傅里叶变换应为-pi到pi为1,其余频率成分为0. 但本人利用如下程序总是得不到正确结果,烦请大家指正。谢谢。

clc
clear
Fs=1000;
T=1/Fs;
L=1000;
t=(1:L-1)*T; %这一步不从0取值是为了y=sin(pi*t)./pi./t在整个t范围内有值
L=size(t);
L=L(1,2);
y=sin(pi*t)./pi./t;
NFFT=2^nextpow2(L);
Y=fft(y,NFFT)/L;
f=Fs/2*linspace(0,1,NFFT/2+1);
plot(f,2*abs(Y(1:NFFT/2+1)))
回复此楼

» 猜你喜欢

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

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

fanlet

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
先画出y=sin(pi*t)./pi./t;
看看是否正确?
野花不栽年年有烦恼无根日日生
2楼2014-09-27 13:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulbcd

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by fanlet at 2014-09-27 13:20:54
先画出y=sin(pi*t)./pi./t;
看看是否正确?

我尝试画了下,采样长度L取1000时间有点短,不能准确反映信号随时间的变化趋势,将采样长度改为10000后作出的时域图如附图所示,是没有什么问题的
3楼2014-09-27 13:28:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulbcd

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by fanlet at 2014-09-27 13:20:54
先画出y=sin(pi*t)./pi./t;
看看是否正确?

附图
如何用Matlab对函数sin(pi*t)/(pi*t)进行FFT变换
1111.png

4楼2014-09-27 13:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulbcd

铁虫 (小有名气)

Matlab做FFT要对时域函数有什么要求吗?
5楼2014-09-27 13:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by liulbcd at 2014-09-27 13:32:17
Matlab做FFT要对时域函数有什么要求吗?

不是MATLAB有要求,要做DFT,采样必须等间隔,且满足香农——奈奎斯特采样定理,才能较准确重建傅里叶谱。你这数据没衰减完,变换后矩形底部是不是有锯齿?

[ 发自手机版 http://muchong.com/3g ]
6楼2014-09-30 13:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liulbcd 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2rKfVQg2jQww 2026-11-06 4/200 2026-11-07 10:03 by byEJ59EAM6hR
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见