用matlab怎样求解二元二阶非线性微分方程组?

方程1 0.003448645y1''-0.000000571y2''+0.00003183(f1+f2)cos[0.0002411(y2-y1)]=a1-0.0002411(M1+M2)
方程2 -0.000000571y1''+ 0.003448645y2''-0.00003183(f1+f2)cos[0.0002411(y2-y1)]=a2-0.0002411(M1+M2)
其中y1和y2为转角,y1'',y2''即为角加速度,f1、f2是力,a1、a2、M1、M2是力矩,且f1、f2、a1、a2、M1、M2为一系列已知的离散数值点,这些离散的数值点与时间t一一对应。现在求角加速度关于时间t的曲线图。
方程如下图所示。要得到的曲线图也和下面的曲线图类似。
请问用matlab怎样才能解出角加速度和角速度关于时间t的曲线图,具体的离散数值点我有,哪位高手能解下面这个方程,我再将离散的数值点补充进来。感激不尽。

第1个回答  2014-03-11
你用一下ode45试试本回答被提问者采纳
相似回答