数学建模社区-数学中国

标题: 求助:Mathematica 求解微分方程组 [打印本页]

作者: Magic111    时间: 2017-9-27 08:43
标题: 求助:Mathematica 求解微分方程组
萌新求助如何求解一个非线性方程组:) {0 t2 m: X/ i. P

" c& b$ h( T* N2 u( Q求助大神,我要解的非线性方程组不是规范的,无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到,具体的:7 ]  H7 N$ ?, l  r
! N# A0 {) h9 U
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*8 L. A6 n; U+ E
    omega,
; V, C/ G5 D- m1 l1 r, S- Y  segma'[t] == 8 E, d# p% ?! T) _
   shear*str -
& w0 I# k' K2 n: ?: p% S    1/(str*omega)*Df[EIS[t]]*EIS'[t],EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0, $ u8 N0 F+ n) ?6 x
  10}]
  V1 d9 }2 {  X- Y以上是我想要求解的微分方程,形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到:7 K2 }' N0 W2 M4 ]* P
8 x: G: `5 j3 W8 N
f[EIS_] :=
' U' ?3 }9 e! J: T8 \, O# Q  Re[Coef /. ; U2 t. R+ _3 |5 e( ]' o* X3 Y2 C9 E
    FindRoot[: C/ _" ^) F4 C5 }  N, ^( ]" v+ O
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];! ]! H4 U! E; L* S" H8 W
也就是f[EIS]每一点的值都需要解一个上述的方程。
) K. i# c$ X9 w. ?! [' G
$ n. L, d- p/ v
- ^. D- F' S; ]1 v
6 n# O& `" a& a: S; \* X% h2 i8 N# f& p程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢?我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程,但是也无法得到结果。
( S8 c" W, v8 n7 m  F
2 G& U! J  Y: u! m- M2 U' T% a& Y' Q! Y; j2 K* |- i

( |1 R2 c4 r4 H+ u9 X9 K  i十分感谢大神帮助5 x" c3 r" L; L1 \

+ m# p  d5 n  n0 @& x
( d/ j% n" g) t6 ?
作者: 有生之莲    时间: 2017-12-21 21:08
泛函就是没办法。。。尝试重新转化处理问题吧
1 N. _& a( \+ v- `6 J9 z% [. ]! V( A




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