用ODE45求解一个二阶微分方程组

oneonly

Pparis 发表于 2015-9-29 10:46
可以考虑直接采用矩阵运算，矩阵运算搞不清楚要展开的话
相关运算规则看一下下面的附件就行了，matlab中 ...

主要是不知道y(1)   y(2)  y(3)  y(4)  ……y(222)  y(223)  y(224)  怎么用矩阵表示，只好挨个运算啦

Pparis

oneonly 发表于 2015-9-29 10:56
主要是不知道y(1)   y(2)  y(3)  y(4)  ……y(222)  y(223)  y(224)  怎么用矩阵表示，只好挨个运算啦

这个表示为列向量的形式就行了

happy

wxl789 发表于 2008-12-5 22:34
能解

具体解法搜索论坛

不建议用ode系列的函数去接，多数情况下无法收敛，即便能够得到解，代价也太大了

xiaoshihanlan

happy 发表于 2015-9-30 09:07
不建议用ode系列的函数去接，多数情况下无法收敛，即便能够得到解，代价也太大了

请问教授，一般自由度多的方程用什么方法求解呢？想要得到分岔图，使用newmark方法可以吗

suffer

xiaoshihanlan 发表于 2017-8-29 15:55
请问教授，一般自由度多的方程用什么方法求解呢？想要得到分岔图，使用newmark方法可以吗

newmark法是一个不错的选择，wilson法也可以考虑