tc=t_span(8000:n);
xc=q(8000:n,1);
yc=q(8000:n,2);
xcc=q(8000:n,3)
ycc=q(8000:n,4)
Y = fft(yc,2000);
Pyy = Y.* conj(Y) / 2000;
f = 1.1849e+005*(0:1000)/2000;
plot(f(1:512),log(Pyy(1:512)))
把这个放在后面
[ 本帖最后由 无水1324 于 2007-6-25 19:01 编辑 ] 你的采样频率不是1000，应该是你1/dt
回复 #30 无水1324 的帖子
现在没时间优化 我现在又要弄清楚分叉庞加莱混沌分形李雅普诺夫指数频谱时域图这些概念和算法及实现思路，导师还要求用ＶＣ编程，ＶＣ很麻烦的，我想弄个ＶＣ和ＭＡＴＬＡＢ混合编程，可看了一些简单实例，按照书上做都做不出来，所以还得弄清楚它们的接口部分，而且时间还比较紧 这年代没法混了回复 #32 无水1324 的帖子
说了半天忘了那为什么出不来幅值谱图啊
应该没出错呀
很是郁闷 可以出来哈
FFT图
FFT图回复 #32 无水1324 的帖子
最后出来的是这样一个图 怎么看着有点不太像啊我画出来的图
回复 #36 无水1324 的帖子
我就用的你的代码怎么和你的结果不一样啊 plot(f(1:512),log(Pyy(1:512)))刚才写错了，不好意思！
回复 #40 无水1324 的帖子
那也不对啊 出来的是这图啊 另外对ＰＹＹ取ＬＯＧ是什么意思啊用ＬＯＧ出来的图
取log其实有时候没有什么意义，就是将幅值作一定的转换，使图比较清楚而已还有我取的yc=q(8000:n,2);