这个大概的意思是：周期运动的话，会出现相同的最大值，其实道理上与 闪频差不多，这是这种方法在系统的解出现跳跃的时候，是不准确的。个人的想法
个人认为这种方法是不可取的
1. 正如无水所说，这个方法只能用于个别系统，对于未知系统根本无法判断结果是否正确
2. 这种方法效率十分低 确实是有一定的问题，今天我也用一下无水上面的方法试试，呵呵 :@L :@L :@L :@L :@L
无水，我用你前面写的程序画出来的分岔图怎么是这样的？？？
为避免你翻页看麻烦，把程序贴出来了
function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=;
clear;clc;close all;
omega=1:0.01:2;
for h=1:length(omega)
T=2*pi/omega(h);
=ode45('dddd',,,[],omega(h));
plot(omega(h),x(10000:100:end,2),'k.');hold on
end 是这样的也很正常，你再去看一下这类系统出现混沌是在那个区域。将omega调整到那个区域就可以了 呵呵，方法是可以算的，和分岔参数的设置有关系，我再算算吧！ 楼主看看这个分岔图吧，算了一个小时左右，看看效果如何～
回复 #63 octopussheng 的帖子
前辈你好！你能关心这个问题我很高兴！！谢谢！就你画出的分岔图来说代表什摸？其实我的问题是a的前面有负号．如果去掉的话系统应该是周期运动，也就是说你画出的分岔图应该代表周期运动，不会发生分岔．你的图怎摸解释？麻烦再把它换成－a试试．谢谢！:@)
回复 #67 mechanic05 的帖子
a= -48.704,他的这个方程很奇怪的回复 #68 无水1324 的帖子
前辈好！我的问题是一个梁的大挠度问题．其中a的值在大于零时是稳定的，不会屈曲．a在小于零时出现屈曲．
前辈，我画出的分岔图瞬态响应太大．就原来的程序我把 plot(omega(h),x(10000:100:end,2),'k.');hold on 修改成 plot(omega(h),x(100000:100:end,2),'k.');hold on为什摸老出现矩阵长度不符的问题？:@L
回复 #66 octopussheng 的帖子
哈，这个图好像人脸啊不过这个方程很奇怪 omega与x矩阵大小看一看是不是一致，主要是后面画图的时候
回复 #70 咕噜噜 的帖子
前辈好！我把１００００改成１０００００，应该可以减小瞬态响应吧？:@) 是的
但是你计算的总点数都没有这么多哈
回复 #73 无水1324 的帖子
:@L :@L :@L我刚装上７．０，怎么同样的程序６．５好象更快一些？
我用７．０已经算了一个小时了，还在算．:lol
[ 本帖最后由 无水1324 于 2007-7-16 11:52 编辑 ] 原帖由 mechanic05 于 2007-7-16 11:34 发表 http://www.chinavib.com/forum/images/common/back.gif
:@L :@L :@L
我刚装上７．０，怎么同样的程序６．５好象更快一些？
我用７．０已经算了一个小时了，还在算．:lol
其实6.5的计算功能还要好一些，只是一些功能、工具箱有限