24小时热门版块排行榜    

查看: 1525  |  回复: 4

茶壶bxy

银虫 (正式写手)

[求助] 由三维空间离散点电势求电荷密度分布问题 已有1人参与

大家好!我现在遇到这样一个求取电荷密度分布的问题:
我利用有限差分法通过MATLAB编程计算得到一三维空间电势分布三维矩阵V,想根据V求取此三维空间内电荷密度分布。根据维基百科上关于电荷密度的词条(http://zh.wikipedia.org/wiki/%E9%9B%BB%E8%8D%B7%E5%AF%86%E5%BA%A6 ),电荷密度可以通过电势的二阶梯度求取,如下公式所示:
由三维空间离散点电势求电荷密度分布问题
式中场位置r的电势为φ(r),电荷密度为ρ(r)。

我在MATLAB编程中运用了如下语句求取电荷密度:

e=8.854*10^(-15);   %介电常数
Rho=e*gradient(gradient( gradient(V, step_y, step_x, step_z) ),step_y, step_x, step_z);

其中step_y, step_x, step_z 分别为y、x、z 向步长数组(与梯度函数定义对应)。计算结果某一平面的电荷密度分布如下图所示:
由三维空间离散点电势求电荷密度分布问题-1

图中黄色的两个凸起不符合实际,实际上对应位置是空气应该没有电荷分布,而中间接近圆形的凸起处才应该有电荷分布,这里比较符合实际。

我想求助下大家以上求电荷密度的过程有没有错误?如果有错误的话那应该怎么改正或者说该怎么求取电荷密度?
十分感谢!
回复此楼
加油
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
茶壶bxy: 金币+100, ★★★很有帮助, 能不能请你再给我写计算建议?十分感谢! 2015-05-14 14:15:03
gradient(gradient( gradient
这里不对吧, 拉普拉斯算符对应的是div(grad(
2楼2015-05-14 00:27:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶壶bxy

银虫 (正式写手)

引用回帖:
2楼: Originally posted by jerkwin at 2015-05-14 00:27:20
gradient(gradient( gradient
这里不对吧, 拉普拉斯算符对应的是div(grad(

谢谢回复!你的意见我研究了一下,的确是我之前搞错了。求拉普拉斯算符,我找到三种思路,一是用divergence(gradient())函数嵌套,二是diff(V1,x,2)+diff(V1,y,2)+diff(V1,z,2),三是用del2 函数,但是目前不知道为什么,这三种方法算出来的电荷密度在z=0平面都是NaN。不知道你对MATLAB熟不熟,能不能请你给我些计算的建议?
十分感谢!
加油
3楼2015-05-14 14:14:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

你先看看你的势函数V对z=0是不是有奇点, 然后再考虑结果
4楼2015-05-14 22:18:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶壶bxy

银虫 (正式写手)

引用回帖:
4楼: Originally posted by jerkwin at 2015-05-14 22:18:36
你先看看你的势函数V对z=0是不是有奇点, 然后再考虑结果

十分感谢!解决了
加油
5楼2015-05-15 21:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 茶壶bxy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 11:31 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
信息提示
请填处理意见