24小时热门版块排行榜    

查看: 2820  |  回复: 10

yudiandida

新虫 (初入文坛)

[求助] matlab中利用FFT进行Fourier变换后真实频率f值计算问题 已有4人参与

CODE:
clear;
clc;
f1 = 3e8/1550e-9;
f2 = 10e9;
N = 2000;
fs = 2e12;
t = (0:N-1)/fs;
f = ((0:N-1)*fs/N-fs/2).*10^3;
A = 20;
Y1 = cos(pi/4*cos(2*pi*f2*t)-pi/4).^2.*cos(2*pi*f1*t+pi/2+A*(cos(2*pi*f2*t)-1/16*cos(2*2*pi*f2*t)));
Y2 = Y1.*exp(i*2*pi*f1*t);
Z1 = fft(Y2);
Z2 = fftshift(Z1);
plot(f,abs(Z2))
axis([-0.3e15,0.3e15,0,120])

函数Y有些复杂,不用细看,简单来说这个程序实现的功能是波长为1550nm的连续光,被频率为10GHz的信号光调制后观察输出光的频谱,函数Y就是一个调制过程。我的问题在于,按理说fft再shift之后,显示频谱的横轴应该是实际的频率,而1550nm的光频率在193THz=193e12Hz附近(根据公式f=c/入,其中f是光频率,c是真空中的光速,入是光的波长),经过调制后频谱的中心频率应该也在这个频率附近才对,但是我的坐标横轴对应的频率值明显数量级都不对啊,想把fs加长,但完全没有作用,程序中的fs范围内已经可以看到完整频谱了。请各位同学帮我看看是不是我的程序有什么问题,到底怎样才能让横轴对应真实频率啊,跪谢了~!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与,应助指数 +1
频率要除以2N
字数不够15个。。。。。。。。。。。
showmethemoney
2楼2014-03-26 17:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yudiandida

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by csgt0 at 2014-03-26 17:44:28
频率要除以2N
字数不够15个。。。。。。。。。。。

能不能说的具体点啊。。。。哪里频率除以2N呢?
3楼2014-03-26 18:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
补充下,N为奇数还是偶数x轴不完全一样,略有差别。

[ 发自手机版 http://muchong.com/3g ]
4楼2014-03-27 00:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yudiandida

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 戴钢盔的猪头 at 2014-03-27 00:32:29
补充下,N为奇数还是偶数x轴不完全一样,略有差别。

那会差到两个数量级么....我差了两个数量级啊
5楼2014-03-27 10:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

引用回帖:
3楼: Originally posted by yudiandida at 2014-03-26 18:45:46
能不能说的具体点啊。。。。哪里频率除以2N呢?...

横坐标啊,你看看fft的介绍吧
showmethemoney
6楼2014-03-27 11:24:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yudiandida

新虫 (初入文坛)

啊~算对了 谢谢大家的 帮助了!
7楼2014-03-27 12:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

野望洋叹

新虫 (初入文坛)

挺有用的啊
8楼2014-03-27 17:43:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)

木虫博士

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
7楼: Originally posted by yudiandida at 2014-03-27 12:20:28
啊~算对了 谢谢大家的 帮助了!

9楼2014-03-27 18:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

兜兜不爱哭

铜虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
(0:N-1)/N*2*pi............................................................................
我要触手可及的阳光!
10楼2014-03-27 20:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yudiandida 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
信息提示
请填处理意见