向量叉乘算子、点乘算子与矩阵运算的关系

最新推荐文章于 2025-03-09 12:39:26 发布

罗伯特祥

最新推荐文章于 2025-03-09 12:39:26 发布

阅读量6.9k

点赞数 6

分类专栏： # ▶ MATLAB/Mathematica/Maple 文章标签：向量叉乘矩阵叉乘算子点乘算子

本文链接：https://blog.csdn.net/weixin_43455581/article/details/105949286

版权

▶ MATLAB/Mathematica/Maple 专栏收录该内容

92 篇文章

订阅专栏

本文详细探讨了向量叉乘的性质及其与矩阵运算的关系，通过MATLAB测试验证了向量叉乘公式及矩阵与向量点乘的一致性。展示了如何使用符号计算来简化复杂的向量和矩阵运算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

文章目录

向量叉乘

$(a \times b) \times c = b (a \cdot c) - a (b \cdot c)$

$a \times (b \times c) = b (a \cdot c) - c (a \cdot b)$

$\times b=\left[\begin{array}{ccc} 0 & -a_{3} & a_{2} \\ a_{3} & 0 & -a_{1} \\ -a_{2} & a_{1} & 0 \end{array}\right]\left[\begin{array}{l} b_{1} \\ b_{2} \\ b_{3} \end{array}\right]$

$a \times b = - [b \times] a$

运算测试

>> syms x1 y1 z1 x2 y2 z2 x3 y3 z3 real
>> clear
>> syms a1 a2 a3 b1 b2 b3 x y z real
>> clear
>> syms a1 a2 a3 b1 b2 b3 c1 c2 c3 real
>> A_=[a1;a2;a3];
>> B_ = [b1; b2; b3];
>> C_ = [c1;c2;c3];
>> test1 = cross(A_,cross(B_,C_))
test1 =
   a2*(b1*c2 - b2*c1) + a3*(b1*c3 - b3*c1)
   a3*(b2*c3 - b3*c2) - a1*(b1*c2 - b2*c1)
 - a1*(b1*c3 - b3*c1) - a2*(b2*c3 - b3*c2)
>> A = [0 -a3 a2;a3 0 -a1;-a2 a1 0]
A =
[   0, -a3,  a2]
[  a3,   0, -a1]
[ -a2,  a1,   0]
>>  B = [0 -b3 b2;b3 0 -b1;-b2 b1 0]
B =
[   0, -b3,  b2]
[  b3,   0, -b1]
[ -b2,  b1,   0]
>> C=C_
C =
 c1
 c2
 c3
>> test2=A*B*C
test2 =
 a2*b1*c2 - c1*(a2*b2 + a3*b3) + a3*b1*c3
 a1*b2*c1 - c2*(a1*b1 + a3*b3) + a3*b2*c3
 a1*b3*c1 - c3*(a1*b1 + a2*b2) + a2*b3*c2
>> simplify(test1-test2)
ans =
 0
 0
 0
 >> test3=cross(cross(A_,B_),C_)
test3 =
 - c2*(a1*b2 - a2*b1) - c3*(a1*b3 - a3*b1)
   c1*(a1*b2 - a2*b1) - c3*(a2*b3 - a3*b2)
   c1*(a1*b3 - a3*b1) + c2*(a2*b3 - a3*b2)
>> simplify(test2-test3)
ans =
 a1*b2*c2 - a2*b2*c1 + a1*b3*c3 - a3*b3*c1
 a2*b1*c1 - a1*b1*c2 + a2*b3*c3 - a3*b3*c2
 a3*b1*c1 - a1*b1*c3 - a2*b2*c3 + a3*b2*c2

############################################################
############              a×b=-[b×]a         ###############
############################################################

>> syms a1 a2 a3 b1 b2 b3 real
>> A=[a1;a2;a3]
A =
 a1
 a2
 a3
>> B=[b1;b2;b3]
B =
 b1
 b2
 b3
>> B_=[0 -b3 b2;
b3 0 -b1;
-b2 b1 0]
B_ =
[   0, -b3,  b2]
[  b3,   0, -b1]
[ -b2,  b1,   0]
>> simplify(cross(A,B)-(-B_*A))
ans =
 0
 0
 0

结论

$\times(b \times c)=\left[\begin{array}{ccc} 0 & -a_{3} & a_{2} \\ a_{3} & 0 & -a_{1} \\ -a_{2} & a_{1} & 0 \end{array}\right]\left[\begin{array}{ccc} 0 & -b_{3} & b_{2} \\ b_{3} & 0 & -b_{1} \\ -b_{2} & b_{1} & 0 \end{array}\right]\left[\begin{array}{c} c_{1} \\ c_{2} \\ c_{3} \end{array}\right]$

矩阵与向量点乘

变量

在这里插入图片描述

算子

在这里插入图片描述

测试

>> I=[Ixx Ixy Ixz;Ixy Iyy Iyz;Ixz Iyz Izz]
I =
[ Ixx, Ixy, Ixz]
[ Ixy, Iyy, Iyz]
[ Ixz, Iyz, Izz]

>> I_=[Ixx Ixy Ixz Iyy Iyz Izz]'
I_ =
 Ixx
 Ixy
 Ixz
 Iyy
 Iyz
 Izz
 
>> w_
w_ =
[ c1, c2, c3,  0,  0,  0]
[  0, c1,  0, c2, c3,  0]
[  0,  0, c1,  0, c2, c3]

>> C_ = [0 -c3 c2;
c3 0 -c1;
-c2 c1 0]
C_ =
[   0, -c3,  c2]
[  c3,   0, -c1]
[ -c2,  c1,   0]

>> simplify(C_*w_*I_-cross(C,I*C))
ans =
 0
 0
 0

结论

在这里插入图片描述