三维旋转的表示 - 版本历史

2022年2月24日 (四) 08:52 Qinz18

2022-02-24T08:52:12Z

←上一版本		2022年2月24日 (四) 08:52的版本
第161行：		第161行：

	</syntaxhighlight>		</syntaxhighlight>

			=== 三维旋转+平移 ===
			描述刚体的空间变换，需要有一个初始刚体和变换后的刚体，选择初始刚体位于世界坐标系中心，刚体自身坐标系与世界坐标系重合，然后来描述末态刚体相对初始刚体的几何变换。这可以通过先旋转后平移来描述，设末态刚体某点坐标为$$\vec{r}$$，则$$\vec{r}=R\vec{r}_0+\vec{d}$$，其中$$\vec{d}$$为末态刚体的坐标系中心相对初始刚体坐标系中心的位矢, $$R$$为将末态刚体平移至使得其自身坐标系中心和初态刚体的坐标系中心重合后，相对于初态刚体的旋转矩阵。

Qinz18：/* 通过欧拉角表示三维旋转 */

2022-02-24T08:41:29Z

通过欧拉角表示三维旋转

←上一版本		2022年2月24日 (四) 08:41的版本
第41行：		第41行：
	[[文件:EulerAngle.png\|缩略图\|欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$]]		[[文件:EulerAngle.png\|缩略图\|欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$]]
	[[文件:EulerAngle ZXZ.gif\|缩略图\|欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$动图]]		[[文件:EulerAngle ZXZ.gif\|缩略图\|欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$动图]]

	<!--		<!--
	{{multiple image		{{multiple image
第54行：		第53行：
	}}		}}
	-->		-->


	考虑刚体的旋转时，我们取初始的刚体轴作为固定的世界轴		考虑刚体的旋转时，我们取初始的刚体轴作为固定的世界轴

Qinz18：/* 通过欧拉角表示三维旋转 */

2022-02-24T08:40:46Z

通过欧拉角表示三维旋转

←上一版本		2022年2月24日 (四) 08:40的版本
第55行：		第55行：
	-->		-->


			考虑刚体的旋转时，我们取初始的刚体轴作为固定的世界轴

	考虑内旋(intrinsic rotations)定义下（即绕刚体轴旋转）的$$Z(\psi)-X'(\theta)-Z''(\varphi)$$欧拉角旋转，依次旋转的角度为$$\psi,\theta,\varphi$$,它与外旋(extrinsic rotations)定义下（即绕世界轴旋转）的欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$的旋转效果相同，旋转矩阵均为''		考虑内旋(intrinsic rotations)定义下（即绕刚体轴旋转）的$$Z(\psi)-X'(\theta)-Z''(\varphi)$$欧拉角旋转，依次旋转的角度为$$\psi,\theta,\varphi$$,它与外旋(extrinsic rotations)定义下（即绕世界轴旋转）的欧拉角$$Z(\varphi)-X(\theta)-Z(\psi)$$的旋转效果相同，旋转矩阵均为''

Qinz18：/* 例子与ROOT代码验证 */

2022-02-24T08:13:57Z

例子与ROOT代码验证

←上一版本		2022年2月24日 (四) 08:13的版本
第109行：		第109行：
	注：在TRotation的文档里有这么一句		注：在TRotation的文档里有这么一句
	<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>		<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>

	这很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。		这很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。

第137行：		第136行：
	//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)		//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)

	//~~通过SetXPhi~~(), SetXTheta(), and SetXPsi()~~来进行欧拉角旋转~~		//通过文档所说的`SetXPhi(), SetXTheta(), and SetXPsi()`来进行欧拉角旋转，会得到错误的结果
	root [13] TRotation R;		root [13] TRotation R;
	root [14] R.SetToIdentity();R.SetXPhi(TMath::Pi()/6);R.SetXTheta(TMath::Pi()/4);R.SetXPsi(TMath::Pi()/3);(R*TVector3(1,0,0)).Print()		root [14] R.SetToIdentity();R.SetXPhi(TMath::Pi()/6);R.SetXTheta(TMath::Pi()/4);R.SetXPsi(TMath::Pi()/3);(R*TVector3(1,0,0)).Print()

Qinz18：/* 例子与ROOT代码验证 */

2022-02-24T08:09:25Z

例子与ROOT代码验证

←上一版本		2022年2月24日 (四) 08:09的版本
第109行：		第109行：
	注：在TRotation的文档里有这么一句		注：在TRotation的文档里有这么一句
	<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>		<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>

	这很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。		这很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。

Qinz18：/* 例子与ROOT代码验证 */

2022-02-24T08:09:02Z

例子与ROOT代码验证

←上一版本		2022年2月24日 (四) 08:09的版本
第107行：		第107行：

	使用ROOT代码验证（两种方法+一种错误使用）		使用ROOT代码验证（两种方法+一种错误使用）
	注：在TRotation的文档里有这么一句<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>~~很容易让人困惑，因为SetXEulerAngles~~(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。		注：在TRotation的文档里有这么一句
			<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>
			这很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。

	<syntaxhighlight lang="cpp" line="1">		<syntaxhighlight lang="cpp" line="1">
第154行：		第156行：
	root [7] R.GetXPhi()		root [7] R.GetXPhi()
	(double) 0.26179939		(double) 0.26179939
			//这也许是SetXPhi(), SetXTheta(), and SetXPsi()不能正常使用的原因



	</syntaxhighlight>		</syntaxhighlight>

Qinz18：/* 例子与ROOT代码验证 */

2022-02-24T08:07:29Z

例子与ROOT代码验证

←上一版本		2022年2月24日 (四) 08:07的版本
第105行：		第105行：

	几何计算可得坐标为$$(\frac{\sqrt{3}}{4}-\frac{\sqrt{3}}{4\sqrt{2}},\frac{1}{4}+\frac{3}{4\sqrt{2}},\frac{\sqrt{3}}{2\sqrt{2}})\approx(0.13,0.78,0.61)$$		几何计算可得坐标为$$(\frac{\sqrt{3}}{4}-\frac{\sqrt{3}}{4\sqrt{2}},\frac{1}{4}+\frac{3}{4\sqrt{2}},\frac{\sqrt{3}}{2\sqrt{2}})\approx(0.13,0.78,0.61)$$

			使用ROOT代码验证（两种方法+一种错误使用）
			注：在TRotation的文档里有这么一句<code>With a minor caveat, the Euler angles of the rotation may be set using SetXEulerAngles() or individually set with SetXPhi(), SetXTheta(), and SetXPsi(). These routines set the Euler angles using the X-convention which is defined by a rotation about the Z-axis, about the new X-axis, and about the new Z-axis. This is the convention used in Landau and Lifshitz, Goldstein and other common physics texts. </code>很容易让人困惑，因为SetXEulerAngles(phi,theta,psi)并不是如文档所说的按新的轴一次转phi,theta,psi,而是按固定轴依次转phi,theta,psi；另外，通过依次SetXPhi(), SetXTheta(), and SetXPsi()并不能正确的实现欧拉角旋转。下面验证这两件事。

	<syntaxhighlight lang="cpp" line="1">		<syntaxhighlight lang="cpp" line="1">
第130行：		第133行：
	r.Print()		r.Print()
	//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)		//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)

			//通过SetXPhi(), SetXTheta(), and SetXPsi()来进行欧拉角旋转
			root [13] TRotation R;
			root [14] R.SetToIdentity();R.SetXPhi(TMath::Pi()/6);R.SetXTheta(TMath::Pi()/4);R.SetXPsi(TMath::Pi()/3);(R*TVector3(1,0,0)).Print()
			TVector3 A 3D physics vector (x,y,z)=(0.324469,0.928023,0.183013) (rho,theta,phi)=(1.000000,79.454709,70.728583)
			root [17] R.SetToIdentity();R.SetXPhi(TMath::Pi()/3);R.SetXTheta(TMath::Pi()/4);R.SetXPsi(TMath::Pi()/6);(R*TVector3(1,0,0)).Print()
			TVector3 A 3D physics vector (x,y,z)=(0.573223,0.739199,0.353553) (rho,theta,phi)=(1.000000,69.295189,52.207654)
			//均不能得到满意的结果，检查SetXPhi()的源代码可以发现，它调用了GetXPhi()等函数，但GetXPhi()并不是一个定义良好的函数，比如：
			root [0] TRotation R
			(TRotation &) Name: TRotation Title: Rotations of TVector3 objects
			root [1] R.SetXPhi(TMath::Pi()/3)
			root [2] R.GetXPhi()
			(double) 0.52359878
			root [3] R.SetXTheta(0)
			root [4] R.SetXPsi(0)
			root [5] R.GetXPhi()
			(double) 0.26179939
			root [6] R.GetXTheta()
			(double) 0.0000000
			root [7] R.GetXPhi()
			(double) 0.26179939



	</syntaxhighlight>		</syntaxhighlight>

Qinz18：/* 通过欧拉角表示三维旋转 */

2022-02-24T07:25:24Z

通过欧拉角表示三维旋转

←上一版本		2022年2月24日 (四) 07:25的版本
第98行：		第98行：
	r.SetXEulerAngles(phi,theta,psi);		r.SetXEulerAngles(phi,theta,psi);

			</syntaxhighlight>


			====例子与ROOT代码验证====
			我们考虑沿`X`轴的单位矢量$$\hat{\vec{x}}$$，即<code>TVector3(1,0,0)</code>，经过欧拉角$$(Z(30^{\circ})-X'(45^{\circ})-Z''(60^{\circ}))$$旋转后的坐标,即上文中令$$\psi=\frac{\pi}{6},\theta=\frac{\pi}{4},\varphi=\frac{\pi}{3}$$。

			几何计算可得坐标为$$(\frac{\sqrt{3}}{4}-\frac{\sqrt{3}}{4\sqrt{2}},\frac{1}{4}+\frac{3}{4\sqrt{2}},\frac{\sqrt{3}}{2\sqrt{2}})\approx(0.13,0.78,0.61)$$

			<syntaxhighlight lang="cpp" line="1">
			TRotation R;
			R.SetXEulerAngles(TMath::Pi()/3,TMath::Pi()/4,TMath::Pi()/6);
			r = TVector3(1,0,0);
			(R*r).Print()//并不会改变r的值
			//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)

			//错误的示范
			TRotation R2;
			R2.SetXEulerAngles(TMath::Pi()/6,TMath::Pi()/4,TMath::Pi()/3);
			(R2*r).Print()
			//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.926777,0.353553)

			//另一种方式
			TVector3 Z(0,0,1);
			TVector3 X(1,0,0);
			r = TVector3(1,0,0);
			r.Rotate(TMath::Pi()/6,Z);//会改变r的值
			X.Rotate(TMath::Pi()/6,Z);
			r.Rotate(TMath::Pi()/4,X);
			Z.Rotate(TMath::Pi()/4,X);
			r.Rotate(TMath::Pi()/3,Z);
			r.Print()
			//输出TVector3 A 3D physics vector (x,y,z)=(0.126826,0.780330,0.612372)
	</syntaxhighlight>		</syntaxhighlight>

Qinz18：/* 三维旋转矩阵 */

2022-02-24T06:26:26Z

三维旋转矩阵

←上一版本		2022年2月24日 (四) 06:26的版本
第19行：		第19行：


	在ROOT中，TRotation类用来描述主动转动的旋转矩阵<code>(TRotation class specifies the rotation of the object in the original system (an active rotation))</code>~~，使用方法如下~~		在ROOT中，TRotation类用来描述主动转动的旋转矩阵<code>(TRotation class specifies the rotation of the object in the original system (an active rotation))</code>，与上文的定义一致，使用方法如下
	<syntaxhighlight lang="cpp" line="1">		<syntaxhighlight lang="cpp" line="1">
	TRotation r;		TRotation r;
第36行：		第36行：

	</syntaxhighlight>		</syntaxhighlight>


	===通过欧拉角表示三维旋转===		===通过欧拉角表示三维旋转===

2022年2月24日 (四) 06:25 Qinz18

2022-02-24T06:25:24Z

←上一版本		2022年2月24日 (四) 06:25的版本
第19行：		第19行：


	~~在ROOT中，TRotation即为旋转矩阵，使用方法如下~~		在ROOT中，TRotation类用来描述主动转动的旋转矩阵<code>(TRotation class specifies the rotation of the object in the original system (an active rotation))</code>，使用方法如下
	<syntaxhighlight lang="cpp" line="1">		<syntaxhighlight lang="cpp" line="1">
	TRotation r;		TRotation r;