| 查看: 1920 | 回复: 6 | ||||
[求助]
单个数列的傅里叶变换（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 |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有203人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有10人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个傅立叶变换的问题
已经有14人回复
- 傅立叶变换-图像处理
已经有4人回复
- 傅里叶变换问题求助
已经有7人回复
- 【紧急求助】请问哪位大神会FFT变换的？急需靠谱FFT代码?非常感谢！
已经有30人回复
- 求问FFT变换是怎么整的？
已经有9人回复
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
感谢参与，应助指数 +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个等间距分布的强烈的峰值吧 |
2楼2013-06-07 09:47:20
3楼2013-06-07 09:59:40
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
4楼2013-06-07 10:30:11
5楼2013-06-07 10:39:44
|
先查看一下 y 函数， 转换到频域的话，大的峰值 的原因是因为你的 y 函数 主要分为三部分 ： y(401:500)=exp(1i*0.63*2*pi); y(701:800)=exp(1i*0.55*2*pi); y(1 :100)=exp(1i*0.85*2*pi); y(101:200)=exp(1i*0.91*2*pi); y(401:400)=exp(1i*0.92*2*pi); y(801:900)=exp(1i*0.96*2*pi); y(901:1000)=exp(1i*0.97*2*pi); y(201:300)=exp(1i*0.12*2*pi); y(501:600)=exp(1i*0.1*2*pi); y(601:700)=exp(1i*0.28*2*pi); 也就不难理解转换到频域，为什么只有左边的峰值大，且总体只有三个值比较大 |
6楼2013-06-08 06:48:00
7楼2013-06-08 06:52:12
