第一章部分习题:
1:
试采用直角坐标中
∇⋅A=∂x∂Ax+∂y∂Ay+∂z∂Az相似的方法推出圆柱坐标下的公式
∇⋅A=ρ1∂ρ∂(ρAρ)+ρ∂ϕ∂Aϕ+∂z∂Az
在直角坐标中选用的是正方体,前后、左右、上下六个面相互抵消后都有一个维度的微元量多出来,将这三个维度的矢量相加,就得到
∇⋅A=∂x∂Ax+∂y∂Ay+∂z∂Az。类似地,柱坐标系中应用角度、半径、高等作为微元分别分析。如下。
固定张角和高度、将半径移动一个微元量,有
ϕr=∫ϕϕ+Δϕ∫zz+ΔzAr∣r+δr(r+δr)drdϕ−∫ϕϕ+Δϕ∫zz+ΔzAr∣rdrdϕ
≈[(r+Δr)Ar(r+Δr,ϕ,z)−rAr(r,ϕ,z)]ΔϕΔz=r1∂r∂(rAr)Δτ
其中,Δτ表示体积元。得到的是半径向外伸展Δr时,通量的变化量。
类似地,
Φϕ=∫rr+Δr∫ϕϕ+ΔϕAz∣z+Δzrdrdϕ−∫rr+Δr∫ϕϕ+ΔϕAz∣zrdrdϕ
≈[Aϕ(r,ϕ+Δϕ,z)−Aϕ(r,ϕ,z)]ΔrΔz≈r1∂ϕ∂(Aϕ)Δτ同理,
ϕz=∫rr+Δr∫ϕϕ+ΔϕAz∣z+Δzrdrdϕ−∫rr+Δr∫ϕϕ+ΔϕAz∣zrdrdϕ
≈[Az(r,ϕ,z+Δz)−Az(r,ϕ,z)]rΔrΔϕΔz≈∂z∂AzrΔrΔϕΔz=∂z∂AzΔτ
所以矢量场穿出六面体表面的通量是
Φ=Φr+Φϕ+Φz≈[r1∂r∂(rAr)+r∂ϕ∂Aϕ+∂z∂Az]Δτ所以有
∇⋅A=ρ1∂ρ∂(ρAρ)+ρ∂ϕ∂Aϕ+∂z∂Az
2:
现有三个矢量ABC分别为
A=ersinθcosϕ+eθcosθcosϕ−eϕsinϕ
B=eρz2sinϕ+eϕz2cosϕ+ez2ρzsinϕ
C=ex(3y2−2x)+eyx2+ez2z
哪些矢量可以由一个标量函数的梯度表示?哪些可以用矢量函数的旋度表示?求出这些矢量的源分布。
此题主要是计算和代入,难度不大,但容易算乱。先根据变量形式,确定ABC分别应该用球坐标系、柱坐标系、直角坐标系公式验证。
在球坐标系中套用公式(要记住柱、球、直角坐标系的旋度公式、散度公式,题目中常用),还要记住无旋场可表示标量函数的梯度、无散场可表示(另一个)矢量函数的旋度。
A⋅∇=r21∂r∂(r2Ar)+rsinθ1∂θ∂(sinθAθ)+rsinθ1∂ϕ∂Aϕ
=r21∂r∂(r2sinθcosϕ)+rsinθ1∂θ∂(sinθcosθcosϕ)+rsinθ1∂ϕ∂(−sinϕ)
=r2sinθcosϕ+rsinθcosϕ−r2sinθcosϕ−rsinθcosϕ=0注意上式中中间一项的求导,是将sin(theta)cos(theta)合并后求导成2cos(2θ),然后展成1-2sin(theta)^2的结果。
Δ×A=r2sinθ1∣∣∣∣∣∣er∂r∂Arreθ∂θ∂rAθrsinθeϕ∂ϕ∂rsinθAϕ∣∣∣∣∣∣至于怎么记住这个三阶行列式,可以认为(非官方)第一行和第三行每个元素表达的都应是边;第一个变量r显然本来就是边,不需要加东西;而θ需要乘个r以表明它是水平方向弧线宽度;φ需要乘个rsinθ表明它是纵向弧线宽度。这样对应起来就是上面的行列式。
=r2sinθ1∣∣∣∣∣∣er∂r∂sinθcosϕreθ∂θ∂rcosθcosϕrsinθeϕ∂ϕ∂−rsinθsinϕ∣∣∣∣∣∣=0
上式代入计算时个人感觉用最直接的三阶行列式行列式展开法简单些。
矢量函数A算得的旋度和散度结果是0,所以A可以用标量函数的梯度来表示,也可以用一个矢量函数的旋度来表示。再看看B,显然给的是柱坐标形式。仍然先求散度。散度的式子好算但不好记,可以与上方的球坐标公式类比来记。
Δ⋅B=r1∂r∂(rBr)+r1∂ϕ∂Bϕ+∂z∂Bz=r1∂r∂(rz2sinϕ)+r1∂ϕ∂(z2cosϕ)+∂z∂(2rzsinϕ)=rz2sinϕ−rz2sinϕ+2rsinϕ=2rsinϕ
再求旋度。对于柱坐标系的这个三阶行列式,(非官方)采取与球坐标系相似的记忆方法记住。显然,只有θ代表的是角度而不是边,要乘r后才能代表它的边(水平弧线)。计算仍然是直接展开简单。
∇×B=r1∣∣∣∣∣∣er∂r∂Brreθ∂ϕ∂rBθez∂z∂Bz∣∣∣∣∣∣=r1∣∣∣∣∣∣er∂r∂z2sinϕreθ∂ϕ∂rz2cosϕez∂z∂2rzsinϕ∣∣∣∣∣∣=0
旋度为0,故可表成一标量函数的梯度。
完全相同地,对C求散度和旋度。直角坐标系相对来说是最简单的。
∇⋅C=∂x∂Cx+∂y∂Cy+∂z∂Cz=∂x∂(3y2−2x)+∂y∂(x2)+∂z∂(2z)=0
∇×C=∣∣∣∣∣∣ex∂x∂3y2−2xey∂y∂x2ez∂z∂2z∣∣∣∣∣∣=ez(2x−6y)
散度为0,故可表示成另一矢量函数的旋度。
至于题目中的源分布的表示,分别将各自上述散度旋度(计算结果)摆在一起就可以了(因为不考虑边界条件时,矢量场可以由它的旋度和散度唯一确定。)
3:
利用直角坐标,证明
∇×(fG)=f∇×G+∇f×G
拿到证明题,先初步理解它。看左边,f没有加粗说明是个标量函数;fG是将这个标量函数作用到矢量场(矢量函数)G上,得到新的矢量函数,然后再对它求梯度。看右边,先对标量函数f求梯度、得到矢量函数、再求与G的叉积。看起来,这是标量函数先作用在函数上再求叉积,还是先对标量函数经过nabla算子转成梯度式再与那个函数相乘的比较问题,感觉差不多是相等的。
接下来具体证一证。带着nabla显然无法证明。先把nabla符号去掉(就是展开),然后看看左右能不能等在一起。例如,
f∇×G展开简单,就是f乘上旋度的展开式罢了,如下。
f∇×G=f[ex(∂y∂Gz−∂z∂Gy)+ey(∂z∂Gx−∂x∂Gz)+ez(∂x∂Gy−∂y∂Gx)]
等号后面的展开相对麻烦一点。
先利用梯度公式
grad u=(ex∂x∂+ey∂y∂+ez∂z∂)
写成
∇f=(ex∂x∂f+ey∂y∂f+ez∂z∂f)
再利用叉乘公式
或者
展开。比如说这里a1指的是左向量的第一个分量,
∂x∂f;
b1指的是右向量第一个(维度的)分量,
Gz
按照上式乘起来,就得到下式:
∇f×G=[ex(Gz∂y∂f−Gy∂z∂f)+ey(Gx∂z∂f−Gz∂x∂f)+ez(Gy∂x∂f−Gx∂y∂f)]
暂时左右还看不出哪里等于了,将右边两者加和看看能不能约去或合并一些项。
f∇×G+∇f×G
=ex[(Gz∂y∂f+f∂y∂Gz)−(Gy∂z∂f+f∂z∂Gy)]+
ey[(Gx∂z∂f+f∂z∂Gx)−(Gz∂x∂f+f∂x∂Gz)]+
ez[(Gy∂x∂f+f∂x∂Gy)−(Gx∂y∂f+f∂y∂Gx)]
将同一方向的向量联立在一起。发现正好是
=ex[∂y∂(fGz)−∂z∂(fGx)]+ey[∂z∂(fGx)−∂x∂(fGz)]+ez[∂x∂(fG)−∂y∂(fG)]
=∇×(fG)
自此完成证明。
4.
用球坐标表示的场
E=err225
求在直角坐标中点(-3,4,-5)处的|E|和Ex
求在直角坐标中,点(-3,4,-5)处E与矢量B=ex2-ey2+ez构成的夹角。
第一问直接出答案,分别是
21和
−2032
。
在直角坐标中点(-3,4,-5)处,
r=−ex3+ey4−ez5
所以
E=r225=r325r=102
ex3+ey4−ez5
所以E与B构成的夹角
θEB=cos−1∣E∣∣B∣(E⋅B)=153.6°