24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Z_Wang

银虫 (正式写手)

应助: 3 (幼儿园)
金币: 7890.3
散金: 2
红花: 1
帖子: 921
在线: 406.6小时
虫号: 1842712
注册: 2012-06-01
专业: 数学物理

[求助] 关于matlab，哪位大侠帮我看下，频谱为什么画不出来？哪里有问题？已有1人参与

c=299792458;
lambda1=1.0*10^(-6);
lambda2=1.3*10^(-6);
omega1=c/lambda1;
omega2=c/lambda2;
dt=1*10^(-18);
t=0:dt:8*10^(-15);
n=size(t,2);
f=sin(2*pi*omega1*t)+2*sin(2*pi*omega2*t);
subplot(2,1,1)
plot(t,f)
y=fft(f,n); %计算频谱
subplot(2,1,2)
plot(1/dt*(-n/2+1:n/2)/n,abs(y)) %画不出来，为什么？

答案可以发到我邮箱：wykszkd@163.com

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有7人回复
三无产品还有机会吗已经有5人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请教一下大侠，两种不同半导体材料接触后构成异质结的能带图用什么软件画图比较好啊？已经有6人回复
请教图中第五列的值是怎么算出来的？已经有8人回复
哪位大神来指导下数字信号处理怎么学，感激不尽已经有30人回复

1楼2014-06-20 09:54:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sxf2012

木虫 (正式写手)

应助: 23 (小学生)
金币: 1133.8
红花: 3
帖子: 313
在线: 148.9小时
虫号: 900812
注册: 2009-11-12
专业: 光学

引用回帖:

3楼: Originally posted by Z_Wang at 2014-06-21 10:31:13
对这个信号做频谱分析，代码怎么写？
我写的对吗，总感觉有点问题。
你能不能帮我重写一下，答案发我邮箱吧，谢谢了
可以加金币

捕获.PNG
...

CODE:

你把总的时间加长，比如到200e-15，然后，由于对称性，频谱只要一半的数据c=299792458;

lambda1=1.0*10^(-6);

lambda2=1.3*10^(-6);

omega1=c/lambda1;

omega2=c/lambda2;

tmax = 200*10^(-15);

dt = tmax/8191;

t = 0:dt:tmax;

n=length(t);

f=cos(2*pi*omega1*t)+2*sin(2*pi*omega2*t);

wmax = 0.5/dt;

dw = 2*wmax/(n-1);

w = -wmax:dw:wmax;

w = fftshift(w);

subplot(2,1,1)

plot(t,f)

y=fft(f);    %计算频谱

subplot(2,1,2)

plot(w(1:n/2),abs(y(1:n/2)))

 axis([0 1e15 0 1e4])