24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2141  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ldlckp

银虫 (正式写手)

[求助] 求助二阶偏微分方程求解 已有2人参与

万能的虫虫们,求助这道题咋求解,本人非数学专业,望能解释的稍微详细一些,谢谢!

求助二阶偏微分方程求解
IMG_1222.JPG
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

若使用分离变量法求解:
先令T(x,t)=F(x,t)-1/3*x^3+1/3*L^2*x
代入原方程得到:
PF(x,t)/Pt=P^2F(x,t)/Px^2
F(0,t)=0 ,F(L,t)=0;
F(x,0)=1/3*x^3-x^2+(1-L^2/3)*x
令F(x,t)=X(x)*T(t) ,代入方程中得到:
X(x)*dT(t)/dt=T(t)*d^2X(x)/dx^2
整理后:1/T(t)*dT(t)/dt=1/X(x)*d^2X(x)/dx^2
由于等式左边为t的函数,右边为x的函数,它们若想等必然只有一种情况,就是都为一常数C。即:
1/T(t)*dT(t)/dt=1/X(x)*d^2X(x)/dx^2≡C
由:1/T(t)*dT(t)/dt=C  得到:
dT(t)/dt-C*T(t)=0
T(t)=D*e^(C*t)
由于t≥0,并且当t-->∞时定解问题应该有界,即e^(C*t)<∞,因此C<0,故而可记为C=-α^2,此处α为正实数。
因此T(t)=D*e^(-α^2*t)
因此,对于X(x),有:
d^2X(x)/dx^2+α^2*X(x)=0
X(x)=A*Cos(α*x)+B*Sin(α*x)
由边界条件,再加上T(t)≠0得到: X(0)=0 , X(L)=0
由此得到A=0,B*Sin(α*L)=0
由于B不能为零,否则只有零解,不符合实际了,故而只能有Sin(α*L)=0
即α*L=k*π ,k=0,±1,±2,±3,±4,......。
即 α=k*π/L
因此:Fk(x,t)=Xk(x)*Tk(t)=Dk*e^[-(k*π/L)^2*t]*Sin[k*π/L*x]
故:F(x,t)=SUM{Fk(x,t),k=-∞~∞}
           =SUM{Dk*e^[-(k*π/L)^2*t]*Sin(k*π/L*x) , k=-∞~∞}
   =SUM{Dk*e^[-(k*π/L)^2*t]*Sin(k*π/L*x) , k=-∞~0} +
      +SUM{Dk*e^[-(k*π/L)^2*t]*Sin(k*π/L*x) , k=0~∞}
    =SUM{Jk*e^[-(k*π/L)^2*t]*Sin(k*π/L*x) , k=0~∞}
由F(x,0)=1/3*x^3-x^2+(1-L^2/3)*x得到:
1/3*x^3-x^2+(1-L^2/3)*x=SUM{Jk*Sin(k*π/L*x) , k=0~∞}
这是一个傅立叶正弦级数展开的问题,由此由傅立叶正弦级数系数公式可求得Jk,从而得到F(x,t)的表达式,再代入T(x,t)=F(x,t)-1/3*x^3+1/3*L^2*x
而得到定解问题的最终解答。


βγ
10楼2015-03-30 11:21:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

P.F.Yin

木虫 (著名写手)

用分离变量方法试试呀

[ 发自小木虫客户端 ]
我想所以我能!
2楼2015-03-30 05:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
目前本人觉得可用分离变量法和拉氏变换法求解。不过用分离变量法求解前须先用一个代换F=U-x^3/3,如此就可将非齐次方程变为其次方程,方便求解。而拉氏变换法最后涉及到复杂式子的反变换问题,若查不到直接可用的公式,就要按照基本定义做复变函数的积分了。

[ 发自手机版 http://muchong.com/3g ]
3楼2015-03-30 06:50:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldlckp

银虫 (正式写手)

引用回帖:
3楼: Originally posted by peterflyer at 2015-03-30 06:50:16
目前本人觉得可用分离变量法和拉氏变换法求解。不过用分离变量法求解前须先用一个代换F=U-x^3/3,如此就可将非齐次方程变为其次方程,方便求解。而拉氏变换法最后涉及到复杂式子的反变换问题,若查不到直接可用的公 ...

谢谢!利用你所说的代换后,能换到两边分别对t和对x的齐次偏微分方程,接下来应该怎么做?是要进行拉氏变换吗?分离变量不是很懂,数学好几年没学了,刚刚接触拉氏变换,基本属于只了解一丁点皮毛,能麻烦你再解释一下吗?十分感激!
4楼2015-03-30 09:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 求调剂 +9 小聂爱学习 2026-04-16 11/550 2026-04-17 22:34 by chixmc
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +4 qut2026 2026-04-11 8/400 2026-04-17 16:10 by 会飞的猪157
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 22408 312求调剂 +23 门路摸摸 2026-04-14 25/1250 2026-04-16 21:21 by Art1977
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
信息提示
请填处理意见