QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4199|回复: 1
打印 上一主题 下一主题

求助:Mathematica 求解微分方程组

[复制链接]
字体大小: 正常 放大
Magic111        

1

主题

4

听众

3

积分

升级  60%

该用户从未签到

自我介绍
Mathematica初学者
跳转到指定楼层
1#
发表于 2017-9-27 08:43 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
萌新求助如何求解一个非线性方程组:
" V* L' O! g8 N( E- b3 M2 v: P1 D0 N- z
求助大神,我要解的非线性方程组不是规范的,无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到,具体的:
: Z# t$ ~" U% n0 k0 }+ o. F, w
% C3 ~8 g$ N2 G- e' JNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*( U" f& Z! |6 E
    omega,
6 q9 T4 F1 |9 J5 {  segma'[t] == 4 @; W. \' i& B2 R0 ]; ]
   shear*str - , F8 [# _% C" F% a3 w, T
    1/(str*omega)*Df[EIS[t]]*EIS'[t],EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0, : c7 J4 e. s; k" F# ^2 `4 }8 T3 A
  10}]
9 i" V; {$ i  f* {7 ?. ~6 s以上是我想要求解的微分方程,形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到:4 {) N; H. g: E. @& w

2 A. t$ Z- f/ ^2 Z8 K0 n$ F6 sf[EIS_] :=
. z4 i- O) f# @/ @: S" u4 ~# w  Re[Coef /.
2 g8 c' B3 q2 m  R    FindRoot[
+ O; P) l" M1 k7 T" K     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];0 J+ ]) f; r" \2 r  M1 W
也就是f[EIS]每一点的值都需要解一个上述的方程。
4 P7 n+ T" ^" d
1 m$ G2 s5 k; y6 b+ T8 s) d$ \; M8 ~1 ~
( ]* c& Y* E) m5 K; Y+ @  R# P( a* K" k1 R* q3 [
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢?我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程,但是也无法得到结果。' L7 K, M9 v0 V+ p& o& J3 y( h( R

3 V' t5 h1 q) A( C% h! W2 s6 z8 Q2 [4 C
" H0 j6 p4 C* [- {# d$ T: k5 P
十分感谢大神帮助
& l' {5 g* A# o' z' M, T3 T) O
" j8 R& G. e8 V" T9 t: N
: o2 Z8 p- r; i* o3 @* j5 K+ d8 X
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

0

主题

1

听众

10

积分

升级  5.26%

  • TA的每日心情
    奋斗
    2017-12-21 21:05
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    物理学本科生
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-20 02:21 , Processed in 0.402292 second(s), 60 queries .

    回顶部