QQ登录

只需要一步,快速开始

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

求助:Mathematica 求解微分方程组

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

1

主题

4

听众

3

积分

升级  60%

该用户从未签到

自我介绍
Mathematica初学者
跳转到指定楼层
1#
发表于 2017-9-27 08:43 |只看该作者 |正序浏览
|招呼Ta 关注Ta
萌新求助如何求解一个非线性方程组:; p' w1 j1 @& A' K: Z

. J7 c* ]% j( k/ ?, o求助大神,我要解的非线性方程组不是规范的,无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到,具体的:8 p9 D* I9 L) A0 l+ [6 k9 K+ S
- S$ |8 f3 P5 v0 Z% x8 K
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*/ D& F: i$ q; t5 Z2 ]3 l4 m
    omega,
% c" Z" w; ~' R5 ?3 i9 \  segma'[t] == ) e1 |7 N8 s9 T: Z8 e
   shear*str -
# S7 I7 J5 [+ m; V    1/(str*omega)*Df[EIS[t]]*EIS'[t],EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
. M0 J# j5 W" G; G# f  10}]
/ w, [! R, w3 f; v3 l以上是我想要求解的微分方程,形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到:+ R2 a  \+ X3 T* G# |

% t/ n. d9 O* u. d8 [! vf[EIS_] :=
1 J/ w: R( L8 t8 E5 w  Re[Coef /.
; E$ ?* }' Q4 `" |5 I2 ~  s0 {    FindRoot[
  {& N* O$ u; q& d     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];/ m  A5 Y0 P2 s# ]
也就是f[EIS]每一点的值都需要解一个上述的方程。
$ q7 c1 s1 }0 p
) H# r% r  }* v% K& r, d. T' L& M0 b, ?2 Y

2 |" V, j6 f; G9 C+ E程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢?我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程,但是也无法得到结果。! C" `1 K' K) G. T8 z9 d
5 F  ^, g: u/ F  ~! C' g

3 F% X- B8 ]  K( t6 J6 X
6 J9 [  {9 G5 o) a) G& L; Q十分感谢大神帮助
/ I/ E5 M4 J' ^, v/ O  R% L1 j4 o- [, {% b& L# n
9 T9 x/ `- V; T+ L7 j) B9 n3 t, z- B  ?
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 00:14 , Processed in 0.468885 second(s), 60 queries .

    回顶部