24小时热门版块排行榜    

查看: 1694  |  回复: 6

bluethky

木虫 (著名写手)

[求助] 有界的扩散方程如何通过调和函数分析法来求解 已有1人参与

求解的是一个有界(-L,L)的fick扩散方程,作者说可以通过harmonic analysis ,就得到扩散方程的解,用一个级数的形式表示的,没有看出来怎么来的,貌似是离散的fourier变换,不过我也没有算出来
求大侠指点,谢谢!


附件是图片,已经相关的文献有界的扩散方程如何通过调和函数分析法来求解
01.jpg


有界的扩散方程如何通过调和函数分析法来求解-1
02.jpg
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : supplementary.pdf
  • 2014-05-31 08:57:45, 112.63 K

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

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

xcenxcen

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
这是一个两点边值问题
P=X(x)*T(t)
XT'=DX''T
T'/(D*T)=X''/X=-lamda^2
X''+lamda^2*X=0
X=C1*cos(lamda*x)+C2*sin(lamda*x)
利用边值条件: X'=C1*lamda*(-sin(lamda*x))+C2*cos(lamda*x)*lamda
x=0, X'=0 =====> C2=0
x=L, X'=0 =====> lamda*sin(lamda*L)=0
lamda*L=k*pi
lamda=k*pi/L
Xk(x)=C1*cos(k*pi/L*x)
Tk(t)=C3*exp(-k^2*pi^2*D/L^2*t)
Pk(x,t)=Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*t)
P=sigma(Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*t))
P(x,ts;xs,ts)=delta(x-xs)
sigma(Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*ts))=delta(x-xs)
sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n*pi/L*xs)
Cn*exp(-k^2*pi^2*D/L^2*ts)*integrate(0,pi)dx(1+cos(2*n*x))/2*L/pi=cos(n*pi/L*xs)
Cn*exp(-k^2*pi^2*D/L^2*ts)*L/2=cos(n*pi/L*xs)
Cn=2/L*cos(n*pi/L*xs)*exp(k^2*pi^2*D/L^2*ts)
P(x,t;xs,ts)=sigma(k=0,+infinity)2/L*cos(n*pi/L*xs)*cos(n*pi/L*x)*exp(-k^2*pi^2*D*(t-ts)/L^2).......
待我长发及腰,遮住一身肥膘。纵然虎背熊腰,也要高冷傲娇。
2楼2014-05-31 17:08:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

引用回帖:
2楼: Originally posted by xcenxcen at 2014-05-31 17:08:32
这是一个两点边值问题
P=X(x)*T(t)
XT'=DX''T
T'/(D*T)=X''/X=-lamda^2
X''+lamda^2*X=0
X=C1*cos(lamda*x)+C2*sin(lamda*x)
利用边值条件: X'=C1*lamda*(-sin(lamda*x))+C2*cos(lamda*x)*lamda
x=0, X'=0 = ...

你用的分离变量法吗?先谢谢再说我好好看看的你如何算的
3楼2014-05-31 20:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

送红花一朵
首先先谢谢你的提示,我有几点不太懂,
1)T'/(D*T)=X''/X=-lamda^2     为什么他们的值要设为 -lamda^2,是从题中隐含的吗?

2)sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n*pi/L*xs),这个式子的含义是什么?  不太理解

3)Cn*exp(-k^2*pi^2*D/L^2*ts)*integrate(0,pi)dx(1+cos(2*n*x))/2*L/pi=cos(n*pi/L*xs) ,这个我也不懂,式子中引入了一个(0,pi)的积分,没有明白您的思路
4楼2014-06-01 13:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

昨天晚上我按您的思路计算了一下,自己不懂的地方,再向您请教一下,先谢谢啦
5楼2014-06-01 13:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcenxcen

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
bluethky: 金币+10, ★★★★★最佳答案 2014-06-01 23:44:46
引用回帖:
4楼: Originally posted by bluethky at 2014-06-01 13:16:55
首先先谢谢你的提示,我有几点不太懂,
1)T'/(D*T)=X''/X=-lamda^2     为什么他们的值要设为 -lamda^2,是从题中隐含的吗?

2)sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n ...

1 设为-lamda^2比较好算,左边是纯t的函数 右边是纯x的函数 所以应该和t,x都没有关系,所以就得是常数啊。
2 这个就是两边同乘以cos(n*pi/L*x),然后在[0,L]上积分。
3 左边的级数,除了k=n的时候不是零以外,其他项都是零,右边呃呃呃。
待我长发及腰,遮住一身肥膘。纵然虎背熊腰,也要高冷傲娇。
6楼2014-06-01 17:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

引用回帖:
6楼: Originally posted by xcenxcen at 2014-06-01 17:45:46
1 设为-lamda^2比较好算,左边是纯t的函数 右边是纯x的函数 所以应该和t,x都没有关系,所以就得是常数啊。
2 这个就是两边同乘以cos(n*pi/L*x),然后在上积分。
3 左边的级数,除了k=n的时候不是零以外,其他项 ...

好的,谢谢!我在好好看看,谢谢指点!
7楼2014-06-01 23:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 逃课的小屁孩 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-22 21:54 by OEbVnUOu01ol
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见