关于非线性振动的数值解法

wanyeqing2003

按gghhjj的方法应该是可以的。
或者将其表示为下式的形式，再用线性加速度方法来分析，是不是也行呢？

咕噜噜

是吗？看来我matlab学的不行啊，我怎么记得ode45是有局限性的，例如二自由度系统，不能解方程中同时含有两个未知变量的二阶导数的方程组那？记错了？！！！！还是编程序错了？！！:@L

咕噜噜

万老师，你没有看清楚题目，那个左边应该是二阶导数和一阶导数，而不是零阶，零阶就简单多了，正则化就可以了

无水1324

原帖由 咕噜噜 于 2007-5-31 10:46 发表

                               
登录/注册后可看大图

是吗？看来我matlab学的不行啊，我怎么记得ode45是有局限性的，例如二自由度系统，不能解方程中同时含有两个未知变量的二阶导数的方程组那？记错了？！！！！还是编程序错了？！！:@L

是可以的！

咕噜噜

:@L 晕，我再找书去，我以前明明看书上这么说得，所以我同学问我我直接告诉他ode45不能解

无水1324

原帖由 咕噜噜 于 2007-5-31 11:50 发表

                               
登录/注册后可看大图

:@L 晕，我再找书去，我以前明明看书上这么说得，所以我同学问我我直接告诉他ode45不能解

不要晕,晕了，编程序是会调试不出来的，
好好看书，有什么收获发出来看看

wanyeqing2003

原帖由 咕噜噜 于 2007-5-31 10:48 发表

                               
登录/注册后可看大图

万老师，你没有看清楚题目，那个左边应该是二阶导数和一阶导数，而不是零阶，零阶就简单多了，正则化就可以了

是的，d太多了。搞糊涂了。:@L

咕噜噜

看来是我记错了，实际自己编一下还是挺简单的上面哪个微分方程组
只是换了一下z(1)=x  z(2)=y  z(3)=x'  z(4)=y'
可是真的记得在那里看过说ode45有局限，右边不能怎么着了
^_^，不好意思，见笑了
看来先入为主的主观看法要不得，记忆力不行了啊

[ 本帖最后由 咕噜噜 于 2007-5-31 15:50 编辑 ]

无水1324

原帖由 咕噜噜 于 2007-5-31 15:49 发表

                               
登录/注册后可看大图

看来是我记错了，实际自己编一下还是挺简单的上面哪个微分方程组
只是换了一下z(1)=x  z(2)=y  z(3)=x'  z(4)=y'
可是真的记得在那里看过说ode45有局限，右边不能怎么着了
^_^，不好意思，见笑了
看来先入为 ...

做出来了就好！

wanyeqing2003

是的。

咕噜噜

数值解法来说的话肯定有其局限性，在特定数值下的解，很难看出系统长期特性
靠大量的数值解来分析系统特性有时候也不见得就全面或者说正确，有没有比较好得数值解法，能够比较完善的来表征系统特性？
^_^，似乎有点强求了

wanyeqing2003

原帖由 咕噜噜 于 2007-6-1 15:10 发表

                               
登录/注册后可看大图

数值解法来说的话肯定有其局限性，在特定数值下的解，很难看出系统长期特性
靠大量的数值解来分析系统特性有时候也不见得就全面或者说正确，有没有比较好得数值解法，能够比较完善的来表征系统特性？
^_^，似 ...

任何方法都有它的局限性。我们应该根据自己的问题选用合适方法。

无水1324

原帖由 咕噜噜 于 2007-6-1 15:10 发表

                               
登录/注册后可看大图

数值解法来说的话肯定有其局限性，在特定数值下的解，很难看出系统长期特性
靠大量的数值解来分析系统特性有时候也不见得就全面或者说正确，有没有比较好得数值解法，能够比较完善的来表征系统特性？
^_^，似 ...

  这个的话，可能要用点映射或者胞映射的方法来全局分析了

咕噜噜

好的，相关内容我去看看

无水1324

原帖由 咕噜噜 于 2007-6-1 16:36 发表

                               
登录/注册后可看大图

好的，相关内容我去看看

  

  比较难哈，做好了教我！