QQ登录

只需要一步,快速开始

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

求助:Mathematica 求解微分方程组

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

1

主题

4

听众

3

积分

升级  60%

该用户从未签到

自我介绍
Mathematica初学者
跳转到指定楼层
1#
发表于 2017-9-27 08:43 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
萌新求助如何求解一个非线性方程组:
8 u$ c2 @$ U) d# [
" D; E2 ^2 D9 Y" `& }- K5 @7 S8 y& S) D求助大神,我要解的非线性方程组不是规范的,无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到,具体的:# A& p4 u- i/ Y: U1 R8 y$ W+ O

8 Y; _- b4 h  R3 E+ S0 DNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*8 R) b2 I9 z  I( V) X( v( B. ?
    omega,1 L+ b  S" t5 ^* _- a% k4 i$ M  ?
  segma'[t] ==
: E( h0 x( F0 v6 _- p   shear*str - ! F$ \& g- U" P9 E9 j2 S- K& z
    1/(str*omega)*Df[EIS[t]]*EIS'[t],EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
5 V7 S  Z( d$ M! Y, z! A) D! A  10}]
- \% K* v, p6 g4 W, d以上是我想要求解的微分方程,形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到:0 }7 ~% w$ ]' K3 k8 r, ]- w+ V

9 s/ n) B$ B2 y# n5 j" c( of[EIS_] := 1 L  \3 J$ z6 T( D$ z
  Re[Coef /.
% I$ U5 j) `! }% Z    FindRoot[. u3 F% y- D& i4 x0 |# `
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
# O3 h( `7 ]; G8 ?6 q; y也就是f[EIS]每一点的值都需要解一个上述的方程。
' x8 ]; A" E2 Y2 `! [+ s. M1 ?
* g9 x& f. l+ b' a: c0 |% X( @2 C* p; B
- _( j  F, @( \" q. a( s$ b9 h2 t
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢?我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程,但是也无法得到结果。) f# _3 r) Z# R" t/ Z

6 m7 T. K* a: _/ r. t
9 m, g* C  u2 q! p: Z- E4 i0 _0 d) h; e- c
十分感谢大神帮助8 j) _! P6 l- v- A; z! r& Q/ u1 [1 m

* f8 {# ?/ c& E/ B: [: x0 Y! E, l1 l+ D
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-8-9 17:50 , Processed in 0.448683 second(s), 60 queries .

    回顶部