| 查看: 1531 | 回复: 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 向步长数组（与梯度函数定义对应）。计算结果某一平面的电荷密度分布如下图所示：  图中黄色的两个凸起不符合实际，实际上对应位置是空气应该没有电荷分布，而中间接近圆形的凸起处才应该有电荷分布，这里比较符合实际。 我想求助下大家以上求电荷密度的过程有没有错误？如果有错误的话那应该怎么改正或者说该怎么求取电荷密度? 十分感谢！ |
回复此楼» 猜你喜欢
不应该看fileCode
已经有7人回复
-
FileCode能看出啥？
已经有25人回复
-
关于Filecode分析方法
已经有11人回复
-
重要来源：本周末出结果
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
好奇怪的filecode
已经有6人回复
-
filecode
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab 3D 等值面图
已经有3人回复
- 【求助】在三维空间中有五个点，这五个点能否相互之间距离相等？非常感谢
已经有8人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 物理EPI: 1
- 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2015-05-14 00:27:20
茶壶bxy
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5793.6
- 散金: 13
- 红花: 2
- 帖子: 327
- 在线: 207.4小时
- 虫号: 1985317
- 注册: 2012-09-08
- 性别: GG
- 专业: 电力电子学
3楼2015-05-14 14:14:17
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 物理EPI: 1
- 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-05-14 22:18:36
茶壶bxy
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5793.6
- 散金: 13
- 红花: 2
- 帖子: 327
- 在线: 207.4小时
- 虫号: 1985317
- 注册: 2012-09-08
- 性别: GG
- 专业: 电力电子学
5楼2015-05-15 21:38:45
