数学建模社区-数学中国

标题: 求助：Mathematica 求解微分方程组 [打印本页]
作者: Magic111    时间: 2017-9-27 08:43
标题: 求助：Mathematica 求解微分方程组
萌新求助如何求解一个非线性方程组：
9 H! e+ D) d0 w4 G( W# Q& k3 b+ m" L) O* u8 u' b+ C9 `5 X: d
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
$ Y7 @/ J. Y8 A- P5 f% i+ M& T! I! W! P8 I4 V& C
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*/ s! o2 o1 n; [* Y& i
    omega,
0 x4 p0 K) a- V: K  segma'[t] == 9 V, T7 r( F* ?6 G7 W, }
   shear*str -
0 F' y% V5 |" Q1 l    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0, . Q+ d$ @: Q# |& `5 Y
  10}]1 Y7 [9 A6 x0 }* x$ a& c# g
以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
* Y" c6 |/ G" o6 c
9 z) b* M/ x% I/ q/ ]9 k1 c, cf[EIS_] :=
3 H4 h( a$ w5 M8 _  Re[Coef /. ; |  i* J4 ]8 g1 i( R& i/ k
    FindRoot[
( C' e% Y; O1 b( S     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];3 y6 R3 F4 T. l2 i
也就是f[EIS]每一点的值都需要解一个上述的方程。/ l4 `8 O8 N9 P* G' [
2 u2 ~! X* ~5 P. }3 f7 x
! S5 l9 x+ x& X6 T( |

9 Z6 e' ~. J+ u3 P0 `5 v程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。1 ]- A2 H$ O/ \
  Z8 J& \4 \" a1 {2 Z6 T

3 X' }' @4 u) O. O% \7 C* s6 {, d0 x- ^
十分感谢大神帮助
( D  Q' s$ L  [$ x- p# v5 [& i) k0 t  q; F8 ~0 v
0 L* g, ^7 L$ Y
作者: 有生之莲    时间: 2017-12-21 21:08
泛函就是没办法。。。尝试重新转化处理问题吧, i# R- K) i1 Z& T# B




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5