24小时热门版块排行榜    

查看: 2839  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 材料调剂 +6 匹克i 2026-03-23 6/300 2026-03-24 21:09 by greychen00
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-24 3/150 2026-03-24 19:17 by pswait
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 0854 考研调剂 招生了!AI 方向 +5 pk3725069 2026-03-19 17/850 2026-03-24 17:30 by zhouxuan..
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 求材料,环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 0854电子信息求调剂 324 +3 Promise-jyl 2026-03-23 3/150 2026-03-23 13:43 by wangkm
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 289材料与化工(085600)B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
信息提示
请填处理意见