| 查看: 1533 | 回复: 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 向步长数组（与梯度函数定义对应）。计算结果某一平面的电荷密度分布如下图所示：  图中黄色的两个凸起不符合实际，实际上对应位置是空气应该没有电荷分布，而中间接近圆形的凸起处才应该有电荷分布，这里比较符合实际。 我想求助下大家以上求电荷密度的过程有没有错误？如果有错误的话那应该怎么改正或者说该怎么求取电荷密度? 十分感谢！ |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有12人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab 3D 等值面图
已经有3人回复
- 【求助】在三维空间中有五个点，这五个点能否相互之间距离相等？非常感谢
已经有8人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 物理EPI: 1
- 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-05-14 22:18:36
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
茶壶bxy
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5793.6
- 散金: 13
- 红花: 2
- 帖子: 327
- 在线: 207.4小时
- 虫号: 1985317
- 注册: 2012-09-08
- 性别: GG
- 专业: 电力电子学
5楼2015-05-15 21:38:45
