应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 应用数学与力学 » 抛物型偏微分方程的数值解求解问题，谢谢！
152/2返回列表上一页12
查看: 3983  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zsmj12

新虫 (初入文坛)
存在非线性抛物型方程吗
赞一下 回复此楼
11楼2013-07-12 11:19:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
楼主的表达好混乱！两个未知数u和u2，一个方程怎么求解？
赞一下 回复此楼
12楼2013-11-12 14:18:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
如果这个偏微分方程只有如下初始条件：u(x,0)=φ(x),u(x，t0)=ψ(x)。而不能确定边界条件（两个边界条件都没有），是能够能求解的。综合起来，定解问题的描述如下：
    微分方程：Pu/Pt=-1/2*P^2u/Px^2+f(x,t)                     （1）
      初始条件：u(x,0)=φ(x),u(x，t0)=ψ(x)。               （2）
    其中f(x,t)、φ(x)、ψ(x)均为已知函数。
    求解过程：
     将u=u(x,t)关于变量t的拉氏变换记为U（x,s)；已知函数f(x,t)的相应拉氏变换记为F(x,s)。此处U（x,s）是待求的函数，而F(x,s)因为f(x,t)已知而可求的，因此是已知的。
       方程（1）两边同求拉氏变换：
    s*U-u(x,0)=-1/2*d^2U/dx^2      
       d^2U/dx^2+2*s*U=2*[F(x,s)+u(x,0)]                        （3）
       将s看作常数，则上式（3）为一个二阶常系数非齐次方程，解之获得U=U（x，s)的表达式，再对其求关于s的拉氏反变换，即可求出u(x,t)。由于f(x,t)不定，无法求出F，因此下面进行不下去了。但思路已明确，一旦f确定即可求出u(x,t)。

解题完毕。
赞一下 回复此楼
13楼2013-11-12 14:51:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
求解过程：
  将u=u(x,t)关于变量t的拉氏变换记为U（x,s)。  方程（1）两边同求拉氏变换：
  s*U-u(x,0)=-1/2*d^2U/dx^2+ x/(s+1)    ,  再将u(x,0)=0代入得：
  d^2U/dx^2+2*s*U=2* x/(s+1)                                 （4）
  将s看作常数，则上式（3）为一个二阶常系数非齐次方程，解之获得U=U（x，s)的表达式：
      U(x,s）= C1(s)*Cos[sqrt(2*s)*x] +C2(s)*Sin[sqrt(2*s)*x] +x/[s*(s+1)]
     根据边界条件：u(0,t)=0和Pu(x，t)/Px(当x=0)=a，得到U(0,s)=0和PU(x，s)/Px(当x=0)=a/s,代入上式，求得C1(s)=0, C2(s)=a/sqrt(2*s)-1/[sqrt(2*s)*s*(s+1)]。
    因此  U(x,s）= {a/sqrt(2*s)-1/[sqrt(2*s)*s*(s+1)]}*Sin[sqrt(2*s)*x] +x/[s*(s+1)] 。
对U(x,s）求拉氏反变换，求得u(x,t):  
因为1/[s*(s+1)]的反变换为1-exp(-t)，因此x/[s*(s+1)] 的反变换即为x*[1-exp(-t)] ；而1/sqrt(s)的拉氏反变换为sqrt(t/π)，故a/sqrt(2*s)的拉氏反变换为a*sqrt(t/(2*π)]；
根据拉氏变换的卷积公式可得a/sqrt(2*s)-1/[sqrt(2*s)*s*(s+1)]的Laplace反变换为：
f1(x,t)= a*sqrt(t/(2*π)]-1/sqrt(2*)*Integral{sqrt(u/)*{1-exp[-(t-u)]}*du,0,t}
= a*sqrt(t/(2*π)]-1/sqrt(2)*{2/3*t^(2/3)-exp(-t)/sqrt(π)*Integral{sqrt(u)*exp(u)]}*du, 0, t}
而Sin[sqrt(2*s)*x]的拉氏反变换求解如下：
若f(t)的拉氏变换为F(s)，依据拉氏变换的基本性质：
（1）Integral[f(t),0, t]的拉氏变换为F(s)/s；
（2）由欧拉公式有Sint=[exp(i*t)-exp(-i*t)]/(2*i) ，也就是说求Sin[a*sqrt（s）]的拉氏反变换最后归结于求exp[a*sqrt(s) ]的拉氏反变换问题；
（3）查拉氏变换与反变换表知：误差余弦函数2/sqrt(π)*Integral{exp(-ξ^2)*dξ, a/[2*sqrt(t)], ∞}的拉氏变换为exp[-a*sqrt(s)]/s，因此，exp[-a*sqrt(s)]的拉氏反变换应为上述误差余弦函数对t的一阶导数，即为：a/4*exp[-a^2/(4*t)]*1/t^(3/2).
令a=±sqrt(2)*i*x ，则得Sin[sqrt(2*s)*x]的拉氏反变换为：
f2(x,t)=sqrt(2)/4*x*exp[-x^2/(4*t)]/t^(3/2)
再次利用卷积公式，得：
u(x,t)= x*[1-exp(-t)]+Integral{f1(x,u)*f2(x,t-u)*du, 0, t}        

解题完毕
赞一下 回复此楼
14楼2013-11-13 10:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

终之太刀—晓

铁杆木虫 (著名写手)

数学爱好者
比较好奇lz的Pde描述的具体模型，看起来更像是反应扩散方程之类的，但是为何a<0?
赞一下 回复此楼
PreferenceforMathematics
15楼2015-03-12 01:52:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zjppzhu 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +4 马方园 2026-09-10 5/250 2026-09-17 16:14 by @-孟豪
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-17 16:03 by 0ngsZVSCTS7z
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 15:27 by 0ngsZVSCTS7z
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-13 5/250 2026-09-17 12:05 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-17 04:16 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 01:52 by Tql5LQhh5rLK
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 17:12 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭