24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

like0508

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2.5
帖子: 43
在线: 18.5小时
虫号: 1327248
注册: 2011-06-20
专业: 凝聚态物性I:结构、力学和

[求助] 单个数列的傅里叶变换（FFT）

自己写的程序如下，是一个数列的傅里叶变换，运行后的图谱不是自己预期的，理论上在频域上，有12个等间距分布的强烈的峰值
x=-0.5:0.001:0.5;
y=zeros(1,1001);
y(1:100)=exp(i*0.85*2*pi);
y(101:200)=exp(i*0.91*2*pi);
y(201:300)=exp(i*0.12*2*pi);
y(401:400)=exp(i*0.92*2*pi);
y(401:500)=exp(i*0.63*2*pi);
y(501:600)=exp(i*0.1*2*pi);
y(601:700)=exp(i*0.28*2*pi);
y(701:800)=exp(i*0.55*2*pi);
y(801:900)=exp(i*0.96*2*pi);
y(901:1000)=exp(i*0.97*2*pi);
Y=fft(y);
Y(1)=[ ];
n=fix(length(Y)/2);
freq=[1:n]./length(Y);
plot(freq,abs(Y(1:n)))

傅里叶变换图.jpg

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有190人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

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

一个傅立叶变换的问题已经有14人回复
傅立叶变换-图像处理已经有4人回复
傅里叶变换问题求助已经有7人回复
【紧急求助】请问哪位大神会FFT变换的？急需靠谱FFT代码?非常感谢！已经有30人回复
求问FFT变换是怎么整的？已经有9人回复

1楼2013-06-07 08:36:23

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

like0508

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2.5
帖子: 43
在线: 18.5小时
虫号: 1327248
注册: 2011-06-20
专业: 凝聚态物性I:结构、力学和

引用回帖:

2楼: Originally posted by dbb627 at 2013-06-07 09:47:20
y=zeros(1,1001);
y(1:100)=exp(i*0.85*2*pi);
y(101:200)=exp(i*0.91*2*pi);
y(201:300)=exp(i*0.12*2*pi);
y(401:400)=exp(i*0.92*2*pi);
y(401:500)=exp(i*0.63*2*pi);
y(501:600)=exp(i*0.1*2*pi);
y(60 ...

你好，我刚运行了下你应助的程序，x轴左端有个很强的峰值其他部分的峰值特别小，这个很大的峰值是什么原因产生的？是不是数据处理上有小问题？

赞一下

回复此楼

3楼2013-06-07 09:59:40

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

查看全部 7 个回答

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
like0508: 金币+10, ★★★★★最佳答案, 谢谢 2013-06-07 10:38:56
csgt0: 金币+2, 谢谢 2013-06-09 10:55:49

y=zeros(1,1001);
y(1:100)=exp(i*0.85*2*pi);
y(101:200)=exp(i*0.91*2*pi);
y(201:300)=exp(i*0.12*2*pi);
y(401:400)=exp(i*0.92*2*pi);
y(401:500)=exp(i*0.63*2*pi);
y(501:600)=exp(i*0.1*2*pi);
y(601:700)=exp(i*0.28*2*pi);
y(701:800)=exp(i*0.55*2*pi);
y(801:900)=exp(i*0.96*2*pi);
y(901:1000)=exp(i*0.97*2*pi);
Y=fft(y);
Pyy = Y.* conj(Y) /length(Y); %功率谱转换
n=fix(length(Y)/2);
f = 1000*(0:n)/length(Y);
plot(f,Pyy(1:n+1));axis([0 100 0 300])
理论上在频域上，有10个等间距分布的强烈的峰值吧

赞一下

回复此楼

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.