您所在位置：主题：[讨论]解微分方程组遇到这种问题该如何处理？请各位高手指点
frankcjn
[专家分：0] 发布于 2007-08-19 01:23:00
用ode15s解微分方程组时，其中有一方程为：
dy12=P-T=((Ay6-B(y12+y14))/(D-y5)+C/(C-Ay5-B(y11+y13)))
-(Ey5-2y7+y11)/(Ey5+2y7-y11)
其中，P和T为力，后边为P和T的具体表达式，y5、y6、y7、y11、y12、y13和y14均为状态变量，A、B、C、D、E为常数。且当 P-T≤0时，y12(位移）由于限位不能为负，也即无论P、T关系如何状态变量y12≥0。请问该问题如何处理？请各位高手给出部分程序段，谢谢！！
dy12=P-T=((Ay6-B(y12+y14))/(D-y5)+C/(C-Ay5-B(y11+y13)))
-(Ey5-2y7+y11)/(Ey5+2y7-y11)
其中，P和T为力，后边为P和T的具体表达式，y5、y6、y7、y11、y12、y13和y14均为状态变量，A、B、C、D、E为常数。且当 P-T≤0时，y12(位移）由于限位不能为负，也即无论P、T关系如何状态变量y12≥0。请问该问题如何处理？请各位高手给出部分程序段，谢谢！！