谁帮我看看这个分岔图对不对啊

无水1324 · 发表于 2007-6-25 18:31

那个不是重要的问题

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 编辑 ]

无水1324 · 发表于 2007-6-25 18:32

你的采样频率不是1000，应该是你1/dt

sssssxxxxx921 · 发表于 2007-6-25 18:36

现在没时间优化　我现在又要弄清楚分叉庞加莱混沌分形李雅普诺夫指数频谱时域图这些概念和算法及实现思路，导师还要求用ＶＣ编程，ＶＣ很麻烦的，我想弄个ＶＣ和ＭＡＴＬＡＢ混合编程，可看了一些简单实例，按照书上做都做不出来，所以还得弄清楚它们的接口部分，而且时间还比较紧　　这年代没法混了

sssssxxxxx921 · 发表于 2007-6-25 18:37

说了半天忘了
那为什么出不来幅值谱图啊

应该没出错呀
很是郁闷

无水1324 · 发表于 2007-6-25 18:42

可以出来哈

无水1324 · 发表于 2007-6-25 18:43

FFT图

sssssxxxxx921 · 发表于 2007-6-25 18:47

最后出来的是这样一个图　　怎么看着有点不太像啊

sssssxxxxx921 · 发表于 2007-6-25 18:49

sssssxxxxx921 · 发表于 2007-6-25 18:52

我就用的你的代码怎么和你的结果不一样啊

无水1324 · 发表于 2007-6-25 18:52

plot(f(1:512),log(Pyy(1:512)))
刚才写错了，不好意思！

sssssxxxxx921 · 发表于 2007-6-25 18:57

那也不对啊　　出来的是这图啊　　另外对ＰＹＹ取ＬＯＧ是什么意思啊

sssssxxxxx921 · 发表于 2007-6-25 18:58

无水1324 · 发表于 2007-6-25 19:00

取log其实有时候没有什么意义，就是将幅值作一定的转换，使图比较清楚而已
还有我取的yc=q(8000:n,2);

sssssxxxxx921 · 发表于 2007-6-25 19:00

这怎么有点像衰减信号啊

无水1324 · 发表于 2007-6-25 19:02

是有点像，你再好好看一下吧，因为你的信号波动非常微弱，所以直接变换后得到很小的幅值X10-7

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 谁帮我看看这个分岔图对不对啊

回复 #30 无水1324 的帖子

回复 #32 无水1324 的帖子

FFT图

回复 #32 无水1324 的帖子

我画出来的图

回复 #36 无水1324 的帖子

回复 #40 无水1324 的帖子

用ＬＯＧ出来的图

回复 #40 无水1324 的帖子

回复 #44 sssssxxxxx921 的帖子