24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1416  |  回复: 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

专家顾问 (正式写手)

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

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的回帖

茶壶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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科西工大 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-13 23:08 by pies112
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[教师之家] 转长聘了 +3 简单化xn 2026-04-13 3/150 2026-04-13 14:18 by wwwkkk83
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[基金申请] 山东省基金2026 +5 jerry681 2026-04-08 6/300 2026-04-12 08:33 by kudofaye
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 电子信息279求调剂,有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
信息提示
请填处理意见