QQ登录

只需要一步,快速开始

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

求助:Mathematica 求解微分方程组

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

1

主题

4

听众

3

积分

升级  60%

该用户从未签到

自我介绍
Mathematica初学者
跳转到指定楼层
1#
发表于 2017-9-27 08:43 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
萌新求助如何求解一个非线性方程组:/ X  _+ ~* @6 n( G1 N! x
" S8 ^/ ]3 H5 S5 x- B7 P
求助大神,我要解的非线性方程组不是规范的,无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到,具体的:
6 n9 X" S4 k% C0 Y2 N6 v/ q& ^! m* `2 ^7 z* F
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*- i. R; O% u4 d9 d* y# k
    omega,
7 w- m% w; B6 Y4 m2 D) Z- w* S  segma'[t] ==
# v7 x8 h% q$ Z- B   shear*str - : i5 P: A$ Y0 [8 A3 Z% d. M: }2 u- F' H
    1/(str*omega)*Df[EIS[t]]*EIS'[t],EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0, 4 a/ k5 ]) J0 E
  10}]0 y  N, L# U3 J0 I0 E  w7 y$ i3 Q
以上是我想要求解的微分方程,形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到:, M# }% T" ]: t6 X( X* z6 N

6 N5 v0 I( }* V$ \1 ef[EIS_] := + I& }' {+ [8 ?+ k  G: j9 b4 T" U# |
  Re[Coef /.
- X' `- s1 p3 F: y5 M* [    FindRoot[
! h1 H6 p0 ^& d) _     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];( q1 x% j- d" [( d" a
也就是f[EIS]每一点的值都需要解一个上述的方程。
, Q/ x9 H( o3 C: d3 Z. u. k3 h
3 s& a7 r- [  x7 t- o. c2 t5 }. l! a9 [
, W. \. o- o" M8 s2 s
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢?我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程,但是也无法得到结果。
* r  i$ C& K4 I; |; n  P1 S) t: B/ U& ?+ N( e

/ f/ I$ m/ k; d9 D) ~
, C- j' C! ?* R; N; b$ a十分感谢大神帮助8 ?8 D4 @) V" L) h% f5 i3 ^1 M* x2 G

- Z% x9 A% A6 t8 X, e% H- h* [
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-6-12 12:37 , Processed in 0.468182 second(s), 59 queries .

    回顶部