| 查看: 1821 | 回复: 0 | |||
茶壶bxy银虫 (正式写手)
|
[求助]
100金币求助_MATLAB计算球坐标系下的拉普拉斯算子
|
|
大家好！我有这样一个问题想请教大家： 我通过计算得到一个三维矩阵V1，它是球坐标系下对应坐标位置 (r,θ,φ) 的电势，其中 r,θ,φ 坐标上取值分别为数组r，theta 和 phi。现在我想对V1求拉普拉斯算子。球坐标系下的拉普拉斯算符如下式：  因为我求过直角坐标系下的拉普拉斯算符，我想到三种思路。 球坐标系下的拉普拉斯算符如下式：  思路一是利用MATLAB自带的离散拉普拉斯算符函数 del2。但是该函数所求的是直角坐标系下的拉式算符，无法直接在球坐标系下使用。将球坐标转化成直角坐标，则又不满足 del2 函数对输入坐标 vectors的要求（转化的直角坐标为三维矩阵）。 思路二是根据拉式算符的表达式用 diff 函数求偏微分计算。这种思路我遇到的问题是由于我的r轴坐标不是均匀分布的，无法用如下所示的二阶中心差分公式，而二阶后向差分公式在非均匀步长下的具体形式我没有推导出来。  思路三是根据拉式算符的定义求梯度的散度，利用 gradient 函数和 divergence 函数计算。但这两个函数默认按照直角坐标系进行计算，似乎不能求球坐标系下对应的梯度和散度。 我个人认为思路二的方法较为可行，想请教下大家，这个问题该怎么做？十分感谢！ |
回复此楼» 猜你喜欢
放榜前的不淡定
已经有18人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有4人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有7人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
100金币求助拉普拉斯方程球坐标下的形式推导问题
已经有16人回复
- 帮忙：求解一组微分代数方程，奖励100金币
已经有7人回复
- 【求助】拉普拉斯算符如何读
已经有7人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20
