向量公式MATLAB符号证明和应用.docVIP

向量公式MATLAB符号证明和应用 　　作者简介：邹正冬（1991―），男，江苏常州人，盐城工学院 纺织服装学院 轻化专业本科三年级学生。研究方向:轻化工程。 　　季红蕾（1961―），女，江苏靖江人，盐城工学院 基础教学部，副教授。研究方向:大学数学。 　　基金项目：国家大学生创新训练项目（盐城工学院2013019）。向量公式的MATLAB符号证明和应用 　　摘要：用MATLAB的符号计算功能证明三向量的混合积、矢量积和四向量的混合积等向量公式，用这些公式分析圆周运动的加速度、由功率公式引入力矩概念以及导出圆周运动质点的转动定律。 　　关键词：向量 ；MATLAB；符号计算 　　中图分类号：O13文献标识码：A文章编号：1671―1580（2014）06―0151―02 　　 　　在高等数学课程中，由于种种原因，一些常用的向量公式没有介绍，如三向量的矢量积和四向量的混合积等。MATLAB是一种具有强大数值计算、可视化和符号计算功能的软件。本文介绍用MATLAB的符号计算功能证明三向量的混合积、矢量积和四向量的混合积等向量公式和在圆周运动的加速度分析和由圆周运动的功率公式引入力矩概念的应用。 　　一、向量公式的证明 　　命题1向量a 、 b和c的混合积 　　(a×b)?c=(b×c)?a=(c×a)?b 　　证明用MATLAB符号计算证明，将命题中的表达式写成 　　(a×b)?c-(b×c)?a=0 　　(b×c)?a-(c×a)?b=0 　　证明的MATLAB程序为 　　 symsa1 a2 a3 b1 b2 b3 c1 c2 c3real 　　 a=［a1,a2,a3］; b=［b1,b2,b3］;c=［c1,c2,c3］; 　　result_1=simple(dot(cross(a,b),c) 　　- dot(cross(b,c),a)) 　　result_2=simple(dot(cross(b,c),a) 　　- dot(cross(c,a),b)) 　　程序执行结果为 　　 result_1= 　　 0 　　 result_2= 　　 0 　　命题1成立。 　　命题2向量a 、 b和 c的二重向量积 　　(a×b)×c=(a?c)b-(b?c)a 　　a×(b×c)=(a?c)b-(a?b)c 　　 证明用MATLAB符号计算证明，将命题中的表达式写成 　　 （a×b)×c-((a?c)b-(b?c)a)=0 　　 证明的MATLAB程序为 　　symsa1 a2 a3 b1 b2 b3 c1 c2 c3real 　　a=［a1,a2,a3］; b=［b1,b2,b3］;c=［c1,c2,c3］; 　　 result_1= simple(cross(cross (a,b),c)- 　　( dot (a,c)*b -dot (b,c)*a)) 　　 result_2= simple(cross(a,cross (b,c))- 　　( dot (a,c)*b -dot (a,b)*c)) 　　 程序执行结果为 　　result_1= 　　［0, 0, 0］ 　　result_2= 　　［0, 0, 0］ 　　命题2成立。 　　注意一般地，（a× b）×c≠a×(b×c).因此，二重向量积不遵循结合律，表达式a×b×c没有意义。 　　命题3向量a、b、c和d的混合积 　　(a×b)?(c×d)=(a?c)(b?d)-(b?c)(a?d) 　　 证明用MATLAB符号计算证明，将命题中的表达式写成 　　(a×b)?(c×d)-((a?c)(b?d)-(b?c) 　　(a?d))=0 　　 证明的MATLAB程序为 　　symsa1 a2 a3 b1 b2 b3 c1 c2 c3 d1 d2 d3real 　　 a=［a1,a2,a3］; b=［b1,b2,b3］; 　　c=［c1,c2,c3］; d=［d1,d2,d3］; 　　 result= simple(dot(cross( a,b), cross( c,d))- 　　( dot( a,c)*dot( b,d) -dot( b,c)*dot( a,d))) 　　程序执行结果为 　　 result= 　　 0 　　命题3成立。 　　二、向量公式的应用 　　应用1圆周运动的质点相对于圆心的位矢为r，则运动速度为:r=w×r 　　加速度为:a=a×r+w×(w×r) 　　式中w为角速度，a为角加速度，且w×a=0，w?r=0，a?r=0. 　　计算可得，v=|w×r|=wr及|a×r|=ar. 　　 由命题1及w×a=0可得 　　(a×r