Triple product formula of vector and its proof (with Python code)

Wutong snow 2021-11-25 09:20:48
triple product formula vector proof

The triple product formula of a vector is an identity that is often used in vector algebra , Its expression is as follows : a ⃗ × ( b ⃗ × c ⃗ ) = ( a ⃗ ⋅ c ⃗ ) b ⃗ − ( a ⃗ ⋅ b ⃗ ) c ⃗ \vec{a}\times\left(\vec{b}\times\vec{c}\right) = \left(\vec{a}\cdot\vec{c}\right)\vec{b}-\left(\vec{a}\cdot\vec{b}\right)\vec{c} a ×(b ×c )=(a c )b (a b )c

We can use the following python Code to prove ( Let's change the above expression a little and prove this identity ) a ⃗ × ( b ⃗ × c ⃗ ) − ( a ⃗ ⋅ c ⃗ ) b ⃗ + ( a ⃗ ⋅ b ⃗ ) c ⃗ = 0 \vec{a}\times\left(\vec{b}\times\vec{c}\right) - \left(\vec{a}\cdot\vec{c}\right)\vec{b}+\left(\vec{a}\cdot\vec{b}\right)\vec{c}=0 a ×(b ×c )(a c )b +(a b )c =0

import sympy as sym
from sympy import sin,cos,diff
x_a,y_a,z_a,x_b,y_b,z_b,x_c,y_c,z_c = sym.symbols('x_a,y_a,z_a,x_b,y_b,z_b,x_c,y_c,z_c')
a = sym.Matrix([x_a,y_a,z_a])
b = sym.Matrix([x_b,y_b,z_b])
c = sym.Matrix([x_c,y_c,z_c])
A = a.cross(b.cross(c))
B = a.dot(c)*b
C = a.dot(b)*c
print(sym.simplify(A-B+C))

The output result is :
 Insert picture description here
so , This identity holds .
Of course , We can also calculate... Manually , The amount of calculation is a little larger , It takes patience and care .

Using the triple product formula , We can get another interesting identity , namely :
a ⃗ × ( b ⃗ × c ⃗ ) + b ⃗ × ( c ⃗ × a ⃗ ) + c ⃗ × ( a ⃗ × b ⃗ ) = 0 \vec{a}\times\left(\vec{b}\times\vec{c}\right) + \vec{b}\times\left(\vec{c}\times\vec{a}\right) + \vec{c}\times\left(\vec{a}\times\vec{b}\right) = 0 a ×(b ×c )+b ×(c ×a )+c ×(a ×b )=0

We can also use python Procedures to prove :

import sympy as sym
from sympy import sin,cos,diff
x_a,y_a,z_a,x_b,y_b,z_b,x_c,y_c,z_c = sym.symbols('x_a,y_a,z_a,x_b,y_b,z_b,x_c,y_c,z_c')
a = sym.Matrix([x_a,y_a,z_a])
b = sym.Matrix([x_b,y_b,z_b])
c = sym.Matrix([x_c,y_c,z_c])
A = a.cross(b.cross(c))
B = b.cross(c.cross(a))
C = c.cross(a.cross(b))
print(sym.simplify(A+B+C))

The output result is :
 Insert picture description here
so , This identity holds .
Of course , We can also calculate... Manually , Expand directly by using the triple product formula, and then merge the similar terms .

版权声明
本文为[Wutong snow]所创,转载请带上原文链接,感谢
https://pythonmana.com/2021/11/20211109004946211Q.html

  1. The difference between single quotation marks, double quotation marks, three single quotation marks and three double quotation marks in Python
  2. Python exercise (11.7)
  3. python11~20
  4. Matlab, python example exercise -- Day2
  5. Matlab, python example exercise -- Day1
  6. White Python road -- Day5
  7. Python正则表达式选择字符
  8. The white way of Python -- day4 (I can't learn it gradually)
  9. The white way of Python -- Day3
  10. The white Python road -- Day2
  11. The white way of Python -- Day1
  12. Difference between input in python2 and python3
  13. Exception handling in Python
  14. Metaclass and singleton pattern of Python
  15. Python methods and functions
  16. Python regular expression select Characters
  17. Introduction to Python object-oriented reflection and built-in methods
  18. Object oriented composition of Python
  19. Python object oriented polymorphism
  20. Object oriented encapsulation of Python
  21. Object oriented inheritance in Python
  22. Python Software Development Directory
  23. Python package
  24. Python module
  25. Python is process oriented and functional
  26. Python function recursion
  27. Python generator
  28. Python iterator
  29. Python decorator
  30. Python function closure
  31. Python namespace and scope
  32. Python function
  33. Python file processing
  34. Python data type
  35. Process control of Python
  36. Garbage collection mechanism of Python
  37. Python: the field notes of dtype = 'object' can also be arranged in ascending and descending order
  38. . value in pandas_ Usage of counts()
  39. 202111 | some hands-on practice of Python time function
  40. 202111 | in Python, timestamp, time string and real date time are converted to each other
  41. Pandas draws line chart, bar chart and bar chart
  42. Exploration and practice of easy transformation between various data structures in Python
  43. A simple practice of Python list generation
  44. Use Python to guess the number of machines and judge the number of guesses
  45. Question about Python: did you learn Python
  46. 20210928 | Python case: building tax calculation function
  47. Python basic syntax collation
  48. Some small accumulation of writing programs in python (4)
  49. Some small accumulation of writing programs in python (3)
  50. Python leak detection tips (3)
  51. Python leak detection tips (2)
  52. Python leak detection tips (1)
  53. Python foundation and MySQL Foundation
  54. Some small accumulation of writing programs in python (2)
  55. Some small accumulation of writing programs in Python
  56. Python 3.7.3 + cuda9.2 installing Python
  57. Python knowledge used to write programs
  58. Python installation + vscode configuration Python environment
  59. Some small problems during Python installation
  60. Answer and Q & A of Python practice introduction course