《电动力学数学基础)资料》由会员分享，可在线阅读，更多相关《电动力学数学基础)资料（7页珍藏版）》请在金锄头文库上搜索。
1、电动力学习题课（一）Feb 20th, 20091特特特殊殊殊函函函数数数与与与爱爱爱因因因斯斯斯坦坦坦求求求和和和约约约定定定1)Kronecker delta函数ij:ij=1i = j0i 6= j(1.1)2)Levi-civita张量ijk:ijk=1ijk是偶置换，即ijk = 123,231,3121ijk是奇置换，即ijk = 213,132,3210otherwise(1.2)3)爱因斯坦求和约定：公式中重复指标自动求和，略去求和号。3Xi=1AiBi AiBi(1.3)注注注意意意：以上的定义只适用于平直空间，而不适用于弯曲空间（广义相对论考虑的空间）。Levi-civit
2、a张量的性质： Levi-civita张量对于下标反称：ijk= jik 下标重复，Levi-civita张量为0：iik= 0 单重求和：ijkmnk= imjn injmflflflfliminjmjnflflflfl 两重求和：ijkmjk= imjj ijjm= 3im im= 2im 三重求和：ijkijk= 2ii= 62标标标量量量，矢矢矢量量量和和和张张张量量量的的的引引引入入入2.1定定定义义义讨论物理量是标量还是矢量的大前提是物理规律的协协协变变变性性性，即描述物理规律的方程在惯性系变换下形式不变。由于运动总是在四维时空中进行，所以不妨设四维空间惯性系的变换为：x0= ax
3、(2.1)那么物理量可按如下分类：1 标量C：不同惯性系变换下的不变量。 矢量A：由四个分量组成，且不同惯性系变换下按如下方式变换的量：A0= aA 二阶张量T：由十六个分量组成，不同惯性系变换下按如下方式变换的量：T0= aaT 类似可以知道n阶张量的定义。注注注意意意：这里所提到的惯性系之间的变换并没有特指某种坐标变换，比如说伽利略变换或者洛伦兹变换。2.2从从从伽伽伽利利利略略略变变变换换换到到到洛洛洛伦伦伦兹兹兹变变变换换换在19世纪末，很多人都认为经典力学是完美的，这是因为经典力学的基本方程Lagranges Equa-tion在伽利略变换下是协变的。我们首先来验证这一点。为了简单，
4、我们考虑下面的变换：x0=x + vtt0=t(2.2)从Eq.(2.2)可以得到：x0=x(2.3)t0=vx+t(2.4)进一步，可得：x0t0=xt(2.5)注注注意意意：Eq.(2.5)并不表示不同参考系下的速度相等，因为速度是位移对时间的全微分。又经典力学的时间观是绝对的，其数学表示为：ddt0=ddt(2.6)综合Eq.(2.3,2.5,2.6)可知：ddt0(L0 x0) L0x0=ddt(L x) Lx= 0(2.7)Eq.(2.7)说明Lagranges Equation在变换Eq.(2.2)下是协变的。然而由麦克斯韦方程组推出的达朗贝尔方程：2 002t2= 0(2.8)显
5、然在伽利略变换下是不协变的，那么现在的问题就是如果承认伽利略变换，麦克斯韦方程组就是错误的；如果承认麦克斯韦方程组，伽利略变换就是不正确的。Einstein选择了后者，用洛伦兹变换取代伽利略变换作为惯性系之间变换，进而建立了狭义相对论。也许你会问：四维时空除了上面两类变换，还有别的变换么？幸运的是数学上可以严格证明R4只有两种微分结构。23正正正交交交曲曲曲线线线坐坐坐标标标系系系3.1基基基本本本概概概念念念定定定义义义：三维空间R3，如果p R3，一组独立、连续、单值函数：u1= f1(x,y,z),u2= f2(x,y,z),u3= f3(x,y,z)(3.1)并且其反函数x = 1(u
6、1,u2,u3), y = 2(u1,u2,u3), z = 3(u1,u2,u3)(3.2)也独立连续单值，则称(u1,u2,u3)为p点的曲曲曲线线线坐坐坐标标标(Curvilinear Coordinates)，(u1,u2,u3)为一一一般般般曲曲曲线线线坐坐坐标标标系系系。在曲线坐标系中，位置矢量为 r(u1,u2,u3)，那么微分线元为：d d r = a1du1+ a2du2+ a3du3(3.3)如果 r， a1, a2, a3两两垂直，则称此曲线坐标系为正正正交交交曲曲曲线线线坐坐坐标标标系系系(Orthogonal CurvilinearCoordinates)。下面考虑正
7、交曲线坐标系的微分线元和基矢：d exdx + eydy + ezdz(3.4)= ex(11d1+12d2+13d3) + ey(21d1+22d2+23d3) + ez(31d1+32d2+33d3)(3.5)=( ex1i+ ey2i+ ez3i)di gi eidi(3.6)其中度量因子gigi= (1i)2+ (2i)2+ (3i)212(3.7)基矢 ei ei=1gi( ex1i+ ey2i+ ez3i)(3.8)3.2基基基矢矢矢 ei ej=ij(3.9) ei ej=ijk ek(3.10)有了基矢的运算公式，便可以考虑矢量A = Ai ei的运算：1)点乘(dot pro
8、duct):A B = AiBj ei ej= AiBi(3.11)2)叉乘(cross product):A B = AiBj ei ej= ijkAiBj ek(3.12)3)三重标积(scalar triple product)：A (B C) = Am em ijkBiCj ek= ijmAmBiCj(3.13)3用Eq.(3.13)不难证明：A (B C) =B (C A) =C (A B)(3.14)4)三重矢积(vector triple product):A (B C)=Am em ijkBiCj ek=AmBiCjijkmkn en=AmBiCj(injm imjn) en=
9、AmBiCm ei AmBmCj ej=(A C)B (A B)C(3.15)利用Eq.(3.15)不难证明：A (B C) +B (C A) +C (A B) = 0(3.16)5)并积(dyadic product):T AB = AiBj ei ej(3.17)我们一般称AB为并并并矢矢矢(dyad)，一共9个分量，其中6个独立。两个或两个以上的并矢之和称为并并并矢矢矢式式式(dyadic)，有9个独立分量，也称二阶张量。并并并矢矢矢相相相关关关运运运算算算:a)点乘： ei ej ek=ij ek(3.18) em ei ej ek=ij em ek(3.19)b)叉乘： ei ej
10、ek=ijn en ek(3.20) em ei ej ek=ij em e ek(3.21)事实上，在张量代数中只要知道了基矢的运算公式就可以计算任何矢量的运算结果了。作为练习，大家可以验证下面的公式： a (b c)=( a b) c(3.22)(b c) a=b( c a)(3.23) a (b c)=( a b) c(3.24)(b c) a=b( c a)(3.25)c)双点积(double dot product)：( ei ej) : ( ek e) ( ej ek)( ei e)(3.26)那么A :B = (Aij ei ej) : (Bk ek e) = AijBji(3.
11、27)很明显，双点积的作用相当于矩阵相乘再求迹(trace)。43.3梯梯梯度度度根据我们已知的直角坐标下的运算可得：dT=Tidi (T) (d)(3.28)T=1giTi ei(3.29)= ei1gii(3.30)虽然这里的推导是对于标量T而言的，但是实际上Eq.(3.30)确是恒成立的。综合Eq.(3.9,3.10,3.29)可得到下面一些有用的公式： ej=gjj(3.31) em=12ijmgigj(i) (j)(3.32)(i) (j)=ijkgigj ek(3.33)细细细论论论：原则上来说，有了Eq.(3.30)以及前面张量分析的基矢运算法则，便可以计算正交曲线坐标系下任何形
12、式的微分运算，但是这个过程中会涉及到联络的概念，比较难以计算，所以下面来介绍一种较为简单的算法。 = ei1gii(3.34)很明显，既是矢量又是线性算符，所以计算时可以分为两步：1)忽略的矢量特征，仅把其当做算符作用于函数或矢量，但要保持等式的运算顺序；2)再考虑是矢量，运用矢量公式进行计算，但要保证求导顺序的正确。下面来看两个例子：(A B)=A(A B) + B(A B)(3.35)=B (AA) + (B A)A +A (BB) + (A B)B(3.36)=B ( A) + (B )A +A ( B) + (A )B(3.37)Eq.(3.35)利用了的算符特征，Eq.(3.36)利
13、用了的矢量特性。 ( A)= ( A) + A ( A)(3.38)=(AA) ( A)A(3.39)=( A) ( )A( A) 2A(3.40)Eq.(3.38)利用了的算符特征，Eq.(3.39)利用了的矢量特性。一般我们用Eq.(3.40)来定义矢量的Laplacian。因为电动力学经常碰到这类运算，所以应该熟悉这种计算方法，作为练习，大家可以尝试推导书后附录一中的公式。3.4散散散度度度和和和旋旋旋度度度散度： A= (Am em)=1g1g2g3(g2g3A1)1+(g3g1A2)2+(g1g2A3)3(3.41)5旋度： A= (Am em)=1g1g2g3flflflflflf
14、lflflflflg1 e1g2 e2g3 e3123g1A1g2A2g3A3flflflflflflflflflfl(3.42)Eq.(3.41,3.42)的证明留给大家。(提示：需要用到Eq.(3.35,3.36,3.37)。)3.5积积积分分分曲线积分：ZbaP(T) d = T(b) T(a)(3.43)高斯定理：ZV( v)d =IS v d a(3.44)斯托克斯定理：ZS( v) d a =IP v d(3.45)3.6举举举例例例：球球球坐坐坐标标标系系系坐标变换：x=rsincosy=rsinsinz=rcos(3.46)易求得：g1=1g2=rg3=rsin(3.47)将E
15、q.(3.47)代入Eq.(3.29,3.41,3.42)可得：T =Tr er+1rT e+1rsinT e(3.48) A =1r2sinhsin(r2Ar)r+ r(sinA)+ rAi(3.49) A=1rsinh(Asin)Ai er+1rh1sinAr(rA)ri e+1rh(rA)rAri e(3.50)下面我们来计算 v =14 rr3的散度：a)当 r 6= 0时， v=14(1r3) r +1r3( r)(3.51)=14(3r3+3r3) = 0(3.52)6b)当 r = 0时，取球心为原点，半径为的球形邻域V，则 v14limVOHS n vdV(3.53)=14li
16、m0+R20dR0sind(r2 rr3 er)433=14lim0+33= +(3.54)若考虑球形邻域V内的积分：ZV vd=IS v nd=IS14 rr3 erd=1(3.55)因此我们定义：( r) v =14 rr3(3.56)称之为Dirac Delta函函函数数数。满足：( r) =0 r 6= 0+ r = 0(3.57)Z( r)d = 1(3.58)类似我们可以求出柱坐标系下的算符的Eq.(3.48,3.49,3.50)。4势势势定定定理理理4.1标标标量量量势势势存存存在在在定定定理理理对于无无无旋旋旋场场场F(Curl-less or irrotational fie
17、lds)，下面说法等价：(a) 场的旋度处处为零，即 F = 0；(b)RbaPF d的值与积分路径P无关；(c) 对任意闭合路径HbaPF d = 0;(d)F可以表示为某标量函数的梯度，即F = V ，其中V称为标（量）势。注注注意意意：标量势不唯一，可以差一常数。在物理上，比如说漩涡、电场。4.2矢矢矢量量量势势势存存存在在在定定定理理理对于无无无散散散场场场F(Divergence-less or solenoidal fields)，下面说法等价：(a) 场的散度处处为零，即 F = 0；(b) 给定积分表面S的边界线，面积分值RSF nd与积分表面S无关；(c) 对任意闭合曲面HF nd = 0;(d)F可以表示为某矢量函数的旋度，即F = A，其中A称为矢（量）势。注注注意意意：矢量势不唯一，可以差一标量函数的梯度V 。在物理上，比如说水杯中旋转的水、磁场。7
