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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】偏微分方程解法
41/1返回列表
查看: 2734  |  回复: 19
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yuehua8500

禁虫 (小有名气)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-03-31 11:09:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
9楼: Originally posted by flyapple5288 at 2010-04-01 09:31:23
我媳妇会这个问题

你这是借用你媳妇的虫号吧？
赞一下 回复此楼
17楼2013-11-13 13:37:37
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
写的潦草模糊，看来楼主有的是时间，用不着别人帮忙的，慢慢自己弄吧。
赞一下 回复此楼
18楼2013-11-13 13:39:20
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
理论上应该可以。因为若用拉氏积分变换，条件已足够，可以顺利解出T1(z,t)的拉氏变换。但因为其太复杂，反变换的解析式求不出，还是不行。倒是可以考虑将拉氏变换展开成复平面的罗伦级数，然后逐项求反变换，应是可行的。
将T2(z,t)关于t的拉氏反变换记为U(z,s),对方程两边求拉氏变换，经整理得到：
d^2U/dz^--s/a*U=0
U(z,s)=C1(s)*exp[sqrt(s/a)*z]+C2(s)*exp[-sqrt(s/a)*z]
由边界条件可推出：dU/dz(z=L）=0，由此C2(s)=C1(s)*exp[2*sqrt(s/a)*L]                                            (1)
故：  U(z,s)=C1(s)*exp[sqrt(s/a)*z]+C1(s)*exp[2*sqrt(s/a)*L]*exp[-sqrt(s/a)*z]                                  （2）
再由hg*[T2(0,t)-Tm]=-k*PT2(z,t)/Pz(z=0)得到：
hg*[U(0,s)-Tm/s]=-k*[dU(z,s)/dz(z=0)]                                                                                           （3）
（2）代入（3），整理后得到：
C1(s)=(hg*Tm/s)/{[hg-k*sqrt(s/a)]*exp[2*sqrt(s/a)]+hg+k*sqrt(s/a)]}                                                  (4)
(4)代入（2）
U(z,s)={exp[sqrt(s/a)*z]+exp[2*sqrt(s/a)*L]*exp[-sqrt(s/a)*z] }* (hg*Tm/s)/{[hg-k*sqrt(s/a)]*exp[2*sqrt(s/a)]+hg+k*sqrt(s/a)]}     (5)
理论上，对（5）求拉氏反变换即可求出T2(z,t),但它太复杂了，估计解不出反变换的解析解了，但可尝试将U(z,s)在复平面展开成s的罗伦级数，然后反变换后就可得到t的幂级数形式了。
赞一下 回复此楼
19楼2013-11-13 21:33:46
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuehua8500 的主题更新
41/1返回列表
信息提示
关闭
请填处理意见
关闭