24小时热门版块排行榜    

查看: 1164  |  回复: 8
当前主题已经存档。

lizhixiao

银虫 (小有名气)

[交流] 【求助】非线性偏微分方程组用什么算法解

小弟遇到一个非线性的偏微分方程组,不知道什么算法可以解出来,求高人指点一二。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)

★ ★
余泽成(金币+1):辛苦了! 1-25 17:44
lizhixiao(金币+1):谢谢了。 1-27 16:52
lizhixiao(金币+1): 2010-04-09 16:45
你查查matlab帮助。

用过matlab解过偏微方程组,速度快,好像当时给出是数值解。
开心努力一辈子
2楼2010-01-25 17:36:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

★ ★
wangen994(金币+1):谢谢应助 1-26 11:18
lizhixiao(金币+1):谢谢区长 了,区长在这个方面是专家。还请区长指点指点。 1-27 17:03
非线性偏微分方程组的解法往往应方程类型而异,所以最好把具体方程给出来。
3楼2010-01-25 20:06:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

★ ★
wangen994(金币+1):谢谢应助 1-26 11:18
lizhixiao(金币+1):谢谢 1-27 16:53
推荐用COMSOL进行计算,很方便的。
明德厚学、求是创新
4楼2010-01-26 09:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyinsky878

木虫 (正式写手)

★ ★
余泽成(金币+1):谢谢应助 1-27 13:02
lizhixiao(金币+1):谢谢了,我回头查查这种方法 1-27 16:54
gear法
5楼2010-01-27 10:16:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lizhixiao

银虫 (小有名气)

nono2009(金币+0):Senlia版主是这方面的真正专家。建议PM他看一下你这个帖子。 1-27 17:15
nono2009(金币+0):我已经提示senlia版主来看了。请耐心等待一下。谢谢! 1-27 17:17
引用回帖:
Originally posted by nono2009 at 2010-01-25 20:06:56:
非线性偏微分方程组的解法往往应方程类型而异,所以最好把具体方程给出来。

我的那个问题有点复杂,简单一点的
Ut=Uxx+exp(1/V);
Vt=Vxx+U+W;
Wt=Wxx+exp(1/V);
我用一般的隐式解法结果很不稳定,Crank-Nicolson 法不知道怎么样,Newton 迭代法怎么样,就是newton 法不知道怎么样在matlab 上面实现。
6楼2010-01-27 17:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)

★
nono2009(金币+1):辛苦了。 2010-02-01 23:23
引用回帖:
Originally posted by lizhixiao at 2010-01-27 17:02:14:

我的那个问题有点复杂,简单一点的
Ut=Uxx+exp(1/V);
Vt=Vxx+U+W;
Wt=Wxx+exp(1/V);
我用一般的隐式解法结果很不稳定,Crank-Nicolson 法不知道怎么样,Newton 迭代法怎么样,就是newton 法不知道怎么样在m ...

一般的隐式是指newton向后差分么?比CN差一阶精度 CN比Newton向后要稳定些
非线性项exp(1/V)显式处理 就是放在第n层
其它的隐式处理 就是取第n+1层和第n-1层

不知你想用什么方法?如果在差分和有限元里面选一个的话 比较推荐有限元 因为可以通过变分把二阶导数降到一阶导数 这样对解的要求就会降低一些 计算效果会好些
书到用时方恨少, 天涯何处无芳草.
7楼2010-02-01 22:55:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

visage1979

银虫 (初入文坛)

★
nono2009(金币+1):谢谢交流!欢迎常来程序软件版。 2010-02-02 19:49
建议看一下数值分析的书,都有关于偏微分方程计算的内容。不过一般来说,采用梯度下降法之类的常用方法,当然程序需要自己编了,时间步长、迭代步长、差分格式之类需要自己设定(这个需要在仿真中自己观察),此外收敛条件也要设定。总体来说,程序本身并不复杂,MATLAB就行,关键是要把问题理解。

还有一点,如果涉及三阶以上的高阶微分,情况会更复杂,建议最多用二阶
8楼2010-02-02 15:40:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ftp0000

木虫 (小有名气)

★
余泽成(金币+1):谢谢参与应助! 2010-02-05 22:13
如果没有用过有限元方法的话,建议用有限差分法,编程比较简单

看方程的样子,应该用Crank-Nicolson 法
9楼2010-02-05 11:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lizhixiao 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[基金申请] 申请科研项目纯靠实力能成功吗? 200+4 David1982 2026-10-13 7/350 2026-10-14 14:24 by 13108017953
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见