数学建模社区-数学中国
标题:
求助：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