| 查看: 1690 | 回复: 6 | ||
bluethky木虫 (著名写手)
|
[求助]
有界的扩散方程如何通过调和函数分析法来求解 已有1人参与
|
|
求解的是一个有界（-L,L）的fick扩散方程，作者说可以通过harmonic analysis ，就得到扩散方程的解，用一个级数的形式表示的，没有看出来怎么来的，貌似是离散的fourier变换，不过我也没有算出来 求大侠指点，谢谢！ 附件是图片，已经相关的文献  01.jpg  02.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : supplementary.pdf
2014-05-31 08:57:45, 112.63 K
» 本帖已获得的红花（最新10朵）
» 猜你喜欢
时间戳又变了
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有4人回复
-
filecode，4个jtjc了
已经有14人回复
-
看来今天不会放榜了？
已经有9人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 对流扩散方程迭代求解发散
已经有11人回复
- 求助用Maple解扩散方程-误差函数
已经有3人回复
- 求解扩散方程问题
已经有9人回复
【答案】应助回帖
感谢参与，应助指数 +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
bluethky
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4937.8
- 散金: 849
- 红花: 11
- 帖子: 1451
- 在线: 228.1小时
- 虫号: 2609601
- 注册: 2013-08-21
- 专业: 植物学研究的新技术、新方
3楼2014-05-31 20:47:31
bluethky
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4937.8
- 散金: 849
- 红花: 11
- 帖子: 1451
- 在线: 228.1小时
- 虫号: 2609601
- 注册: 2013-08-21
- 专业: 植物学研究的新技术、新方
送红花一朵 |
首先先谢谢你的提示，我有几点不太懂， 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
bluethky
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4937.8
- 散金: 849
- 红花: 11
- 帖子: 1451
- 在线: 228.1小时
- 虫号: 2609601
- 注册: 2013-08-21
- 专业: 植物学研究的新技术、新方
5楼2014-06-01 13:17:49
6楼2014-06-01 17:45:46
bluethky
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4937.8
- 散金: 849
- 红花: 11
- 帖子: 1451
- 在线: 228.1小时
- 虫号: 2609601
- 注册: 2013-08-21
- 专业: 植物学研究的新技术、新方
7楼2014-06-01 23:44:40
40