电动力学CAS表示的简单方法.docVIP

电动力学CAS表示的简单方法 摘 要 论述了“场论算符Math化”、“Math函数算符化”。提出一种电动力学CAS表示的简单方法。例示和验证了在电动力学中应用计算机进行公式推导及字符解析运算的可行性。 关键词 算符; 代数; 四维协变量; 相对论 我们曾基于数值运算建立变分的电磁场方程组[1],应用中有许多编程之苦，有应用的局限性。实践表明字符运算与数值运算结合更好。应用字符运算软件不难使教材电子化，印在纸上是本书，输入电脑既是能说话，又能推演解难。用它来克服电动力学中的计算瓶颈不但效率高，而且可使教材突出物理模型和数学模型，并可充分阐明原理和规律。“Mathematica”系统（以下简称Ｍath）由物理学家Stephen Wolfram领导开发的数学软件。它有强大的符号计算和数值计算能力。系统提供了数学、物理所需得全部字符。我们用Math作为工具经过多年的摸索，找到了“场论算符Math化”、“Math函数算符化”的简单而又可行方法，使电动力学的表达式具有代数风格。实现了电动力学的计算机代数系统（CAS）表示，并用计算机进行了公式推导、字符解析运算。 1．场量Math表示 我们主张以Lorentz协变量展开电动力学教材，因此要把电动力学中的物理量用熟悉的字符定义,并采用“表”或矩阵的输入/输出形式，例如： 三维物理量可表为： 等等。又如四维量： uk ：= {1,ux/c,uy/c,uz/c} (* 速度 *) um ：= gnn uk jk ：=ρ c uk; (* 电流密度 *) jm ：= gnn.jk kl ：= {ω/c,k1,k2,k3}; (* 波矢 *) km ：= gnn.kl; pm ：= m0 c uk; (* 动量 *) pk ：= gnn.pm; am ：= {φ/c,-ax,-ay,-az}; (* 四维势 *) ak ：= gnn.am; fmn={{0,e1/c,e2/c,e3/c},{-(e1/c),0,-b3,b2}, {-(e2/c),b3,0,-b1},{-(e3/c),-b2,b1,0}} ; 等等。四维协变量缩并可得Lorentz不变量,在惯性系变换中很有应用天地，现在可由计算机导出。 例如: cs=(xm.xk/.{(-x^2-y^2-z^2)-0})/.t-t0;Print[xm.xk==,cs] cu=um.uk/.{ux-0,uy-0,uz-0,gmu-1} ;Print[um.uk==,cu] cj=jm.jk/.{ux-0,uy-0,uz-0,gmu-1} ;Print[jm.jk==,cj] cp=pm.pk/.{ux-0,uy-0,uz-0,gmu-1} ;Print[pm.pk==,cp] fmnfmn=Sum[fmn[[i,j]] fkl[[i,j]],{i,1,4},{j,1,4}]/2//Simplify; 即得： uiui = 1 2．Math函数算符化 Math系统的数学操作与演算函数等都是函数表达式，表达式的“头”一般是单词、单词组合或它们的缩写组成，以数学符号命名加以定义，使用时酷似一个算符。 例如： ?:=Map[Function[g,D[#,g]],{x,y,z}] ?a?b:=((#1.▽)[#2])+（（#2.▽）[#1]+#1×(▽×#2)+#2×(▽×#1) 等等, 以上各式中左边的字符(串)既是电动力学的算符名,又是Math系统的操作函数的“头”，右边是操作函数。其它的Math函数都可以由此定义、储存、调用。我们把它简称Math算符。应用举例如下： 前式计算梯度，后式是梯度公式展开的操作函数。结果是： 式中“{{}，{}，{} }”用以表示矢量和张量的各个分量的。 3．场论算符Math化 电磁场理论的表达式及其演绎体系简洁、优美，是经典物理宝贵财富，必需沿用发展。我们以算符记号命名并定义场论的算符展开函数(为了易于辩认，这里用字符串表示，)。举例如下： Delcc:=(#1 Del[#2]+#2 Del[#1]) Divcv:=(#2.Del[#1]+#1 Div[#2]) Rotcv:=