sssssxxxxx921 发表于 2007-8-14 12:19

帮我看看这个高维的分岔程序

这个高维的微分方程在求解时总的到10的几十次幂，大家帮我看看这个程序是怎么回事啊

chuandong418 发表于 2007-8-14 12:44

呵呵，sssssxxxxx921好象也问题和偶一样的多呢,加油呢!

无水1324 发表于 2007-8-14 12:48

回复 #1 sssssxxxxx921 的帖子

最怕的就是你这种一计算就达到非常大的数的情况，太麻烦了，检查你的参数吧

sssssxxxxx921 发表于 2007-8-14 13:47

我也没办法不过我正在查发上来是想这来个艳遇了   

说不定有个人太牛了非这种问题不解决那我不占大便宜了   哈哈

咕噜噜 发表于 2007-8-14 14:22

回复 #4 sssssxxxxx921 的帖子

可能是参数问题，换一个解法试试看，顺便换换参数

sssssxxxxx921 发表于 2007-8-14 15:44

参数我试了好多了ODE的解法换完了呵呵，其他的方法老师强烈建议不换

说是这是最苦的时候撑一撑就过去了   但愿如此吧

sssssxxxxx921 发表于 2007-8-14 18:59

回复 #2 chuandong418 的帖子

是啊   对于对称模型解是求出来了   但分岔图出来似乎都是混沌   
现在怀疑是周期选取的问题    周期激励是个多频的无法通约的   以前就没解决了   本来选其中之一估计肯定会出来结果（即使有点偏差）   现在看来不是那么一回事    周期选不对   庞加莱截面就选取就有问题，    结果就可能是一直振荡   没有一条稳定曲线
    对于非对称模型来说    解根本无法获得   一会就发散到10的几十次    所以最近一直很是烦恼

咕噜噜 发表于 2007-8-14 19:09

回复 #7 sssssxxxxx921 的帖子

我没仔细看你的程序和问题，你的问题是自由振动还是强迫振动？
周期怎么取得？

sssssxxxxx921 发表于 2007-8-14 22:40

回复 #8 咕噜噜 的帖子

支承运动而导致的位移激励是强迫振动
周期选取见我另外一个帖子，就是从三个周期中选了一个模型耦合处的周期
见http://forum.vibunion.com/forum/viewthread.php?tid=47614
这样选择实在不合适   但多频激励不好选如果按照传统的选三者的最小公倍数不太可能做到
查了一些资料   遇到多频激励的都是用近似解析法来解的
页: [1]
查看完整版本: 帮我看看这个高维的分岔程序