24小时热门版块排行榜    

查看: 3908  |  回复: 14

zjppzhu

铜虫 (初入文坛)

[求助] 抛物型偏微分方程的数值解求解问题,谢谢!

需要求解如下抛物型偏微分方程:
du/dt=-1/2(du2/dx2)+f(x,t)
加上一定的初始条件和边界条件。
打算采用有限差分法求解,疑惑在于:采用Crank-Nicolson格式进行离散化时,由于系数a=-1/2(看到所有的文献资料中的例子,均设a>0),导致网格比r=a*delta_t/(delta_x*delta_x)是负数,在这种情况下,Crank-Nicolson格式还是恒稳定的吗?谢谢!

另外,如果这个偏微分方程只有初始条件,而不能确定边界条件(两个边界条件都没有),应该属于反边界问题吗?能求解吗?谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

archdevil

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zjppzhu: 金币+2, ★★★很有帮助, 回答了部分疑惑,谢谢! 2013-01-11 21:46:52
这是一个线性微分方程,如果边界条件不复杂,可以求出解析解的,这样不是更好。
如果一定要用数值解,这里有个问题,一般而言,方程的解中会包含因子exp(-at),即若a>0,则解会随着时间的增长而衰减,此时我们称解是稳定的(注意:是数学上的稳定,不是数值意义上的稳定)。如你所说,若a<0,则解在数学上都不稳定,那么通过数值方法也一定得不到稳定的解。
鲜花因生之绚烂而凋谢,珊瑚因死之静默而永恒。
3楼2013-01-11 00:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zjppzhu

铜虫 (初入文坛)

谢谢!可惜没有我这种情况的:(
我的问题特殊性在于:那个系数a=-1/2,是个负数,而所有看到的资料上,这个a>0?就是搞不清楚,在这种情况下进行有限差分法求数值解行不行?
我对于偏微分方程的了解得实在太少了,不好意思呀。
2楼2013-01-10 22:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjppzhu

铜虫 (初入文坛)

感谢您的答复!我根据问题推导得到的偏微分方程如上面所述,如果初始条件不能用数学表达式表达出来(如随机信号、脑电信号等),应该得不到解析解了吧?还有,根据问题的先验知识,还定不出边界条件(也可能是我的功底太弱,不知道如何定,只好先用Dirichlet、Neumann、Robin边界条件试试)。但如你上面所述,若a<0会导致解不稳定,那整个模型就废掉了,呜呜.......
4楼2013-01-11 21:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjppzhu

铜虫 (初入文坛)

我想了一下,我的模型即使在a<0时,应该还是有希望有解的,因为du/dt=-1/2(d2u/dx2)+f(x,t)中最右边的非齐次项是-detla×u,而delta是征得、可调的微小量,这样的话,在解析解中,会贡献一个exp(-1/(detla)^2*t)因子,应该有希望可以抵消掉发散的exp(-at)项吧?还请上面的大侠和其他各位们不吝赐教,谢谢!
5楼2013-01-11 22:18:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archdevil

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
zjppzhu: 金币+8, ★★★很有帮助, 感谢您的回复! 2013-01-12 12:16:08
引用回帖:
4楼: Originally posted by zjppzhu at 2013-01-11 21:44:23
感谢您的答复!我根据问题推导得到的偏微分方程如上面所述,如果初始条件不能用数学表达式表达出来(如随机信号、脑电信号等),应该得不到解析解了吧?还有,根据问题的先验知识,还定不出边界条件(也可能是我的功 ...

从方程上看,如果边界条件不是很复杂(如你所说,是Dirichlet、Neumann边界条件),这个问题可以得到解析解。用分离变量法中关于解非齐方程非齐边界的本征函数展开法可以试试,这里初始条件不需要是一个连续函数,只要能够做Fourier级数展开即可。
另一方面,即便是真的找不到解析解,那么由于a<0导致解不稳定,通常的数值方法不能再用,不过在数值方法里还有一种称为相对稳定的稳定性分析,楼主可能要找找这方面的文献了(我一般只做绝对稳定的数值方法),看看有没有新路子。
鲜花因生之绚烂而凋谢,珊瑚因死之静默而永恒。
6楼2013-01-11 22:18:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjppzhu

铜虫 (初入文坛)

非常感谢楼上的答复!关于分离变量法求解非齐方程非齐边界的本征函数展开法、数值方法关于“相对稳定的稳定性分析”有没有好的文献或者书籍推荐?万分感谢!
7楼2013-01-12 12:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archdevil

铁杆木虫 (正式写手)

引用回帖:
7楼: Originally posted by zjppzhu at 2013-01-12 12:14:10
非常感谢楼上的答复!关于分离变量法求解非齐方程非齐边界的本征函数展开法、数值方法关于“相对稳定的稳定性分析”有没有好的文献或者书籍推荐?万分感谢!

本征函数展开法(也有称特征函数展开法,Fourier级数展开法)是数学物理方程或偏微分方程求解的一个基本方法,基本在各类数学物理方程或偏微分方程教材中都有提及,比如中科大严镇军的《数学物理方程》,陈祖曦的《偏微分方程》,复旦姜礼尚的《数学物理方程讲义》等等。
至于相对稳定的稳定性分析,在李荣华的《微分方程数值解法》和胡建伟、汤怀民的《微分方程数值方法》中略有提及,但这方面详细的讲解我就不清楚哪本书里有了。
鲜花因生之绚烂而凋谢,珊瑚因死之静默而永恒。
8楼2013-01-14 21:48:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjppzhu

铜虫 (初入文坛)

非常感谢!
9楼2013-01-16 20:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuqiang_sdu

铜虫 (初入文坛)

Cracknicolson格式不再无条件稳定,而是条件稳定的。需要对网格比加以限制。换句话说dt不能太大。具体我没推导。仅供参考。如果没有边界条件(连周期性或者紧支集的那种都没有的话)那就考虑解析的方法吧。Fourier transform试试
10楼2013-03-14 16:02:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zjppzhu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 03:31 by xxniao123
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见