24小时热门版块排行榜    

查看: 2840  |  回复: 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的回帖

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的回帖
查看全部 11 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 招08考数学 +7 laoshidan 2026-03-20 16/800 2026-03-24 17:52 by 乌拉儿山脉
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 284求调剂 +10 Zhao anqi 2026-03-22 10/500 2026-03-24 00:08 by Equinoxhua
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 0703化学求调剂 +4 奶油草莓. 2026-03-22 5/250 2026-03-23 19:37 by pswait
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
信息提示
请填处理意见