数学建模社区-数学中国

标题: 求教——一个非线性微分方程的解法 [打印本页]
作者: 神七    时间: 2009-5-2 11:01
标题: 求教——一个非线性微分方程的解法
本帖最后由 神七 于 2009-5-2 11:10 编辑
6 D: f! Y' h& M, W" C4 Q) V2 ]# J& w# p" ]5 Y6 @  H7 a
我是学气象的，特来这里向各位高手请教用数值方法求解一个非线性微分方程，好久了都没有算出来这个例子，心急如焚！还请各位高手帮帮忙啊！谢谢了先！！！8 ?) k$ t  W& W! R0 e  C
, G' e$ C& ]" e( \- Q3 M
方程和参数请见附件。- W' W3 r9 p# u7 r0 a$ j+ U
[attach]6388[/attach]
- q' \9 m' V$ V3 m/ b! I7 q
' K% b2 `% k6 X" p0 D; ^2 G
8 Q: N! j: a( A( L, C! M' ^我把最后的方程做二阶差分，求非线性方程组的梯度法和拟牛顿法都不能求解这一方程。" A* _  ]$ A4 ~' B
- u! y: d3 h9 n$ O0 G. b* [0 U7 ]
采用直接迭代算出来的fai的量级在10E7左右，基本正确。但是算出来的水平风速u在10E4左右，明显偏大。实际的u的大小为-50-50之间。不知道是哪里出了问题，是不是有别的方法可以求解这一方程？还请各位高手指点，不胜感激！！！

图片2.gif (42.44 KB, 下载次数: 342)

图片2.gif

作者: 神七    时间: 2009-5-4 16:52
顶起来，让各位高手看到，呵呵！
作者: artin    时间: 2009-8-30 10:26
试试用经典runge_kutta方法
作者: willzhang1982    时间: 2009-11-11 12:55
帮忙顶了，貌似很麻烦啊！！！！！！！！！！！！！！！
作者: justbetoo    时间: 2010-8-8 21:40
不会啊，貌似比较难啊
作者: yangyuhonglirui    时间: 2013-9-22 09:56
比较厉害了，都已经知道这么多算法了。值得学习。这些算法我还一个都不会用呢。



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5