应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab中利用FFT进行Fourier变换后真实频率f值计算问题
111/2返回列表12下一页
查看: 2810  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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范围内已经可以看到完整频谱了。请各位同学帮我看看是不是我的程序有什么问题，到底怎样才能让横轴对应真实频率啊，跪谢了~！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-26 16:27:41
已阅   回复此楼   关注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 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 272求调剂 +3 材紫有化 2026-02-28 3/150 2026-02-28 22:52 by ms629
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考博] 26申博 +4 想申博！ 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 284求调剂 +4 天下熯 2026-02-28 4/200 2026-02-28 21:13 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 5/250 2026-02-28 20:11 by iwuli
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂；
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭