数学建模社区-数学中国

标题: 求教——一个非线性微分方程的解法 [打印本页]

作者: 神七    时间: 2009-5-2 11:01
标题: 求教——一个非线性微分方程的解法
本帖最后由 神七 于 2009-5-2 11:10 编辑
) E7 U0 v) a* f! `6 _, A* p; _% o1 k9 U2 s- n
我是学气象的,特来这里向各位高手请教用数值方法求解一个非线性微分方程,好久了都没有算出来这个例子,心急如焚!还请各位高手帮帮忙啊!谢谢了先!!!2 w& G8 E6 P, f

0 U1 W8 b1 p, j/ |" T. n方程和参数请见附件。
# H% t7 H9 M; ]) j[attach]6388[/attach]9 n3 B2 S! ?# @5 G" d2 h2 @- v3 P
% q. T3 O+ [8 s  L- ^0 G% j
2 X' V3 h. E' Z+ e* X5 P' L
我把最后的方程做二阶差分,求非线性方程组的梯度法和拟牛顿法都不能求解这一方程。8 s% Y1 i) F( G+ g
: c$ Y: w3 p. X- t" f# y
采用直接迭代算出来的fai的量级在10E7左右,基本正确。但是算出来的水平风速u在10E4左右,明显偏大。实际的u的大小为-50-50之间。不知道是哪里出了问题,是不是有别的方法可以求解这一方程?还请各位高手指点,不胜感激!!!

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

图片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