2楼: Originally posted by lxm3515 at 2014-06-21 10:00:24
里面有一句话有问题：y=fft(f,n);    %计算频谱
改为：y=fftshift(fft(f,n)); 即可

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

捕获.PNG
...

4楼: Originally posted by lxm3515 at 2014-06-21 18:39:00
这两个频率中心不一样，需要分开计算，再在频域相加。...

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

捕获.PNG
...

你把总的时间加长，比如到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])