从全局坐标结构到局部坐标结构

计算科学 局部坐标
2021-12-16 20:00:47

我也将把这个问题放在数学交换和堆栈交换上(因为它并不真正适合任何特定领域,我不知道它应该放在哪个领域),但问题是:

我正在寻求评估给定点的球谐函数,但要根据局部坐标系。让实际点被调用(xg,yg,zg)(“全球”点)。现在,我有一个以原子为中心的原子,比如说,(xat,yat,zat)它有一个以下列方式定义的局部坐标系:

(1) z 轴由从给定原子中心到选定原子的向量定义(例如,atom1 的 z 轴定义为从 atom1 到 atom5 的向量)

(2) y 轴由位于第二个原子定义的平面中的正交向量定义(例如,atom1 的 y 轴定义为位于第一个向量和从 atom1 到 atom3 的向量定义的平面内) . 这是通过简单的克施密特正交化来完成的。

(3) x 轴由前两个向量的简单叉积定义。

(请注意,所有向量都是标准化的)。

现在,查找此信息将为您留下一个 3x3 矩阵,该矩阵定义了定义给定轴的“端点”的点。所以局部坐标矩阵(对于给定的原子)可以定义为:

(xxxyxzyxyyyzzxzyzz)

假设“零”点(或局部坐标系的中心)是(xat,yat,zat).

现在,给定这些信息,在全局坐标系中某个点的局部坐标系中评估球谐函数(或任何函数)的最佳方法是什么所以,例如,我想评估一个点(XG,是的G,zG)但“x”方向定义为(XX,X是的,Xz)等等。

2个回答

“全局”和“局部”系统之间的坐标转换在Geodesy领域是相当标准的新不伦瑞克大学有关于坐标系和最小二乘法主题的精彩讲义。我推荐大地测量学笔记中的坐标系。在第 52 页左右,他们解释了如何在本地天文系统(有点像您的本地系统)到平均(或常规)陆地系统(全球系统)之间转换

(XG是的GzG)=R3(180-λ)R2(90-φ)2(XXX是的Xz)

在哪里λ是经度(X一个,是的一个,z一个)θ是纬度。角度 180 和 90 您需要适应您的问题,但逻辑是相同的。矩阵2用于反映是的轴,因为全球系统是右旋的,而本地天文系统是左旋的。如果这不是您的情况,那么只需将矩阵排除在外。

之后,您可以使用Deathbreath 的答案来计算球谐函数。

您必须将每个全局点转换为局部坐标系。根据您的描述,您实际上只有平移和旋转。

XlC一个l=ü(XGlb一个l-)
在哪里ü是单一的并且是一个平移向量。现在你有,llC一个l(θ(XlC一个l),η(XlC一个l)).