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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 软件使用与操作 » Mathematica偏微分方程求解编程求助
汕头大学海洋科学接受调剂
61/1返回列表
查看: 1661  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 计算强帖 ，点击这里进行查看

北-极-星

木虫 (正式写手)

[求助] Mathematica偏微分方程求解编程求助 已有1人参与

如图，我想在Mathematica中编程求出此方程的含边界条件的解，自己编程如下，不能运行，求大神帮忙指点。。
DSolve[{(\[PartialD]/\[PartialD]t)c[z, t] -
    D*(\[PartialD]/\[PartialD]z)(\[PartialD]/\[PartialD]z)c[z, t] ==
   0, c[z, t == 0] == Subscript[c, 0],
  D (\[PartialD]/\[PartialD]z)c[z == 1/2*Subscript[h, s], t] == 0,
  D (\[PartialD]/\[PartialD]z)
     c[z == 1/2*Subscript[h, s] + Subscript[h, f], t] == Subscript[J,
   0]}, c[z, t], {z, t},
Constants -> {D, Subscript[c, 0], Subscript[h, s], Subscript[h, f],
   Subscript[J, 0]}]

Mathematica偏微分方程求解编程求助
1.jpg


Mathematica偏微分方程求解编程求助-1
2.jpg


Mathematica偏微分方程求解编程求助-2
程序.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-07-25 20:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
dbb627: 金币+5, 计算强帖+1, 应助指数+1, 谢谢指导 2015-01-06 23:45:02
首先要指出的是，楼主的代码里语法错误非常多，并且都是很明显的错误。Mathematica（啊，确切地说是Mathematica内所使用的Wolfram）终究是编程语言，它有它明确的语法规范。就算这个方程DSolve可解，那正确的语法也应该是
CODE:
With[{c = c[z, t]}, eqn = D[c, t] - d D[c, z, z] == 0];
bc = {c[z, 0] == c0, Derivative[1, 0][c][hs/2, t] == 0, Derivative[1, 0][c][hs/2 + hf, t] == j0/d};
DSolve[{eqn, bc}, c[z, t], t]

然后，上面这段代码依旧会将输入原样返回，并提示软件无法求解这个问题。这个方程之所以会无法求解，个人觉得和Mathematica的DSolve对偏微分方程符号解的求解力弱有一定关系，但是更多的应该是因为Mathematica在总体上排斥使用级数来表示解（顺便Maple似乎也一样），理由嘛，个人猜测和级数的收敛性判定的困难有关，当然，这只是猜测。

题外话好像说的多了一点。既然Mathematica无法直接求解这个问题，那么我们有没有迂回策略呢？答案是有的。我们利用Mathematica来复现一般偏微分方程里面的相应求解方法就好了。

先拉普拉斯变换并求解所得的常微分方程：
CODE:
sol = FullSimplify[
  DSolve[LaplaceTransform[{eqn, Rest@bc}, t, s] /.
      bcrule[[1]] /. {l_[c[z, t], t, s] :> c[z],
      l_[d_[a_, 0][c][z_, t], t, s] :> d[a][c][z]}, c[z], z][[1, 1,
    2]], d > 0 && hf != 0 && Abs[hs/2] <= Abs[z] <= Abs[hs/2 + hf]]

然后，因为我们已经知道这个东西反变换没有有限形式的解，而Mathematica不能直接给出级数解，那么，我们就把上面的结果在{hs/2, hs/2+hf}上做傅立叶余弦展开：
CODE:
(* 平移到{0, hf} *)
midsol = sol /. z -> z + hs/2 // Simplify
(* term1为余弦展开的第0项。这项没包含在通项里所以要单独求。 *)
term1 = FourierCosSeries[midsol, z, 0, FourierParameters -> {1, Pi/hf}, Assumptions -> hf > 0]
(* Mathematica虽然不乐意使用无穷级数表示解，却不缺求解级数系数通项的能力 *)
term2 = FourierCosCoefficient[midsol, z, n,
  FourierParameters -> {1, Pi/hf}, Assumptions -> hf > 0]

因为拉普拉斯变换是线性变换，我们只要把上面的term2拉普拉斯反变换后再（在形式上）将它加合再把term1的反变换也加进来（别忘了把坐标系换回去）就获得了这个方程的级数解了：
CODE:
# + HoldForm[
    Sum[#2, {n, 1,
      Infinity}]] & @@ ((Simplify[InverseLaplaceTransform[#, s, t],
       d > 0] /. z -> z - hs/2) & /@ {term1, term2 Cos[Pi n/hf z]})

Mathematica偏微分方程求解编程求助-3

嗯……形式上和你的有一点区别，主要原因是Sum内部的两项其实可以展开，而其中的一部分是可以求和的——这一求和Sum也确实可以做到，但是，Mathematica默认情况下给出的形式和你的文献所给的不太一样，并且我暂时没找到将该式子化简成你的那个形式的方法。不过至此为止你的问题也基本解决了，而且我回答这个半年前的问题也纯粹是饭前血糖低导致的脑抽。所以，到此为止吧。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

北-极-星 北-极-星
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
2楼2015-01-06 19:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

北-极-星

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by xzczd at 2015-01-06 19:02:51
首先要指出的是，楼主的代码里语法错误非常多，并且都是很明显的错误。Mathematica（啊，确切地说是Mathematica内所使用的Wolfram）终究是编程语言，它有它明确的语法规范。就算这个方程DSolve可解，那正确的语法也 ...

多谢指导  还要自己在研究研究
赞一下 回复此楼
3楼2015-01-06 22:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北-极-星

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by xzczd at 2015-01-06 19:02:51
首先要指出的是，楼主的代码里语法错误非常多，并且都是很明显的错误。Mathematica（啊，确切地说是Mathematica内所使用的Wolfram）终究是编程语言，它有它明确的语法规范。就算这个方程DSolve可解，那正确的语法也 ...

多谢指导   我自个再研究研究
赞一下 回复此楼
4楼2015-01-06 22:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)
引用回帖:
3楼: Originally posted by 北-极-星 at 2015-01-06 22:04:03
多谢指导  还要自己在研究研究...

……仔细一看发现我漏掉了bcrule的定义。需要在bc的定义后面增加一行：
CODE:
bcrule = Rule [at=315114]@@[/at] bc

或者直接把bcrule换成
CODE:
c[z, 0] -> c0

也可以。
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
5楼2015-01-07 13:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)
引用回帖:
2楼: Originally posted by xzczd at 2015-01-06 19:02:51
首先要指出的是，楼主的代码里语法错误非常多，并且都是很明显的错误。Mathematica（啊，确切地说是Mathematica内所使用的Wolfram）终究是编程语言，它有它明确的语法规范。就算这个方程DSolve可解，那正确的语法也 ...

请收下小的这个月的膝盖。
赞一下 回复此楼
Be water, my friend.
6楼2015-01-07 19:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 北-极-星 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +8 超人不会飞@ 2026-04-08 8/400 2026-04-14 16:06 by zs92450
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 求调剂 +3 我爱高数高数爱� 2026-04-12 3/150 2026-04-14 01:00 by 王珺璞
[考研] 310求调剂 +15 666真好 2026-04-11 17/850 2026-04-13 22:27 by pies112
[考研] 297工科，求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 考研二轮调剂 +10 故人?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
信息提示
关闭
请填处理意见
关闭