§B1.1 向量与内积§B1.2 叉积与混合积(三维)§B1.3 矩阵与二次型自测题 本节只作简要回顾。正文用到的地方:刚体转动中的叉积恒等式(§P1.5 )、能量均分定理中的二次型(§5.5 )、熵产生与昂萨格系数的半正定性(§P5.4 )、晶格的原胞体积与倒格子(§10.1 、§10.2 )。
§B1.1 向量与内积
n n n 维实向量 a = ( a 1 , … , a n ) \mathbf a = (a_1,\dots,a_n) a = ( a 1 , … , a n ) 。内积 (点积)a ⋅ b = ∑ i a i b i = a T b \mathbf a\cdot\mathbf b = \sum_ia_ib_i = \mathbf a^{\mathsf T}\mathbf b a ⋅ b = ∑ i a i b i = a T b ,长度 ∣ a ∣ = a ⋅ a \lvert\mathbf a\rvert = \sqrt{\mathbf a\cdot\mathbf a} ∣ a ∣ = a ⋅ a 。二维、三维中 a ⋅ b = ∣ a ∣ ∣ b ∣ cos θ \mathbf a\cdot\mathbf b = \lvert\mathbf a\rvert\lvert\mathbf b\rvert\cos\theta a ⋅ b = ∣ a ∣ ∣ b ∣ cos θ ;一般维数中,柯西–施瓦茨不等式 ∣ a ⋅ b ∣ ≤ ∣ a ∣ ∣ b ∣ \lvert\mathbf a\cdot\mathbf b\rvert\le\lvert\mathbf a\rvert\lvert\mathbf b\rvert ∣ a ⋅ b ∣ ≤ ∣ a ∣ ∣ b ∣ 保证了"夹角"有定义。(证明:对一切实数 t t t ,∣ a + t b ∣ 2 = ∣ a ∣ 2 + 2 t a ⋅ b + t 2 ∣ b ∣ 2 ≥ 0 \lvert\mathbf a + t\mathbf b\rvert^2 = \lvert\mathbf a\rvert^2 + 2t\,\mathbf a\cdot\mathbf b + t^2\lvert\mathbf b\rvert^2\ge0 ∣ a + t b ∣ 2 = ∣ a ∣ 2 + 2 t a ⋅ b + t 2 ∣ b ∣ 2 ≥ 0 ,这个关于 t t t 的二次式的判别式必须非正。)
内积为零称为正交 。一组两两正交的单位向量 e 1 , … , e n \mathbf e_1,\dots,\mathbf e_n e 1 , … , e n 构成正交归一基 ,e i ⋅ e j = δ i j \mathbf e_i\cdot\mathbf e_j = \delta_{ij} e i ⋅ e j = δ ij ;任意向量的分量就是它在基向量上的投影:a = ∑ i a i e i \mathbf a = \sum_ia_i\mathbf e_i a = ∑ i a i e i ,a i = e i ⋅ a a_i = \mathbf e_i\cdot\mathbf a a i = e i ⋅ a 。
正交补 (附录 A8 用到):给定向量 g 1 , … , g m \mathbf g_1,\dots,\mathbf g_m g 1 , … , g m ,与它们都正交的向量构成一个子空间 W W W ;而与 W W W 中所有向量都正交的向量,恰好是 g 1 , … , g m \mathbf g_1,\dots,\mathbf g_m g 1 , … , g m 的线性组合。(理由:把任一向量分解为"在 g k \mathbf g_k g k 张成的子空间中的分量"加"与之正交的分量";后者属于 W W W ,若它不为零,它就与自身不正交。)
§B1.2 叉积与混合积(三维)
a × b = ( a y b z − a z b y , a z b x − a x b z , a x b y − a y b x ) (B1.1) \mathbf a\times\mathbf b = (a_yb_z - a_zb_y,\ a_zb_x - a_xb_z,\ a_xb_y - a_yb_x) \tag{B1.1} a × b = ( a y b z − a z b y , a z b x − a x b z , a x b y − a y b x ) ( B1.1 )
它垂直于 a \mathbf a a 与 b \mathbf b b (按右手定则),长度 ∣ a ∣ ∣ b ∣ sin θ \lvert\mathbf a\rvert\lvert\mathbf b\rvert\sin\theta ∣ a ∣ ∣ b ∣ sin θ 等于以 a \mathbf a a 、b \mathbf b b 为边的平行四边形面积;b × a = − a × b \mathbf b\times\mathbf a = -\mathbf a\times\mathbf b b × a = − a × b 。混合积 a ⋅ ( b × c ) \mathbf a\cdot(\mathbf b\times\mathbf c) a ⋅ ( b × c ) 等于以三个向量为棱的平行六面体的有向体积,也等于以它们为行的 3 × 3 3\times3 3 × 3 行列式(附录 B3 )。例如布拉维格子的原胞体积 v c = a 1 ⋅ ( a 2 × a 3 ) v_{\mathrm c} = \mathbf a_1\cdot(\mathbf a_2\times\mathbf a_3) v c = a 1 ⋅ ( a 2 × a 3 ) ,倒格子基矢 b 1 = 2 π a 2 × a 3 v c \mathbf b_1 = 2\pi\frac{\mathbf a_2\times\mathbf a_3}{v_{\mathrm c}} b 1 = 2 π v c a 2 × a 3 等(§10.1 、§10.2 )。
两个常用恒等式:
a × ( b × c ) = b ( a ⋅ c ) − c ( a ⋅ b ) , ∣ a × b ∣ 2 = ∣ a ∣ 2 ∣ b ∣ 2 − ( a ⋅ b ) 2 (B1.2) \mathbf a\times(\mathbf b\times\mathbf c) = \mathbf b(\mathbf a\cdot\mathbf c) - \mathbf c(\mathbf a\cdot\mathbf b),\qquad \lvert\mathbf a\times\mathbf b\rvert^2 = \lvert\mathbf a\rvert^2\lvert\mathbf b\rvert^2 - (\mathbf a\cdot\mathbf b)^2 \tag{B1.2} a × ( b × c ) = b ( a ⋅ c ) − c ( a ⋅ b ) , ∣ a × b ∣ 2 = ∣ a ∣ 2 ∣ b ∣ 2 − ( a ⋅ b ) 2 ( B1.2 )
(第一个俗称"BAC–CAB"规则,可以逐个分量验证,例如 x x x 分量:a y ( b x c y − b y c x ) − a z ( b z c x − b x c z ) = b x ( a y c y + a z c z ) − c x ( a y b y + a z b z ) a_y(b_xc_y - b_yc_x) - a_z(b_zc_x - b_xc_z) = b_x(a_yc_y + a_zc_z) - c_x(a_yb_y + a_zb_z) a y ( b x c y − b y c x ) − a z ( b z c x − b x c z ) = b x ( a y c y + a z c z ) − c x ( a y b y + a z b z ) ,在第一项中加上、第二项中减去 a x b x c x a_xb_xc_x a x b x c x ,即得 b x ( a ⋅ c ) − c x ( a ⋅ b ) b_x(\mathbf a\cdot\mathbf c) - c_x(\mathbf a\cdot\mathbf b) b x ( a ⋅ c ) − c x ( a ⋅ b ) ;第二个由 sin 2 θ = 1 − cos 2 θ \sin^2\theta = 1 - \cos^2\theta sin 2 θ = 1 − cos 2 θ 得到。)§P1.5 中刚体的动能 1 2 ∑ m ∣ ω × r ∣ 2 \frac12\sum m\lvert\boldsymbol\omega\times\mathbf r\rvert^2 2 1 ∑ m ∣ ω × r ∣ 2 与角动量 ∑ m r × ( ω × r ) \sum m\,\mathbf r\times(\boldsymbol\omega\times\mathbf r) ∑ m r × ( ω × r ) 正是用这两个恒等式化成转动惯量张量的。
§B1.3 矩阵与二次型
n × n n\times n n × n 矩阵 A A A 作用于向量:( A x ) i = ∑ j A i j x j (A\mathbf x)_i = \sum_jA_{ij}x_j ( A x ) i = ∑ j A ij x j 。转置 ( A T ) i j = A j i (A^{\mathsf T})_{ij} = A_{ji} ( A T ) ij = A ji ,( A B ) T = B T A T (AB)^{\mathsf T} = B^{\mathsf T}A^{\mathsf T} ( A B ) T = B T A T 。
二次型 Q ( x ) = x T A x = ∑ i , j A i j x i x j Q(\mathbf x) = \mathbf x^{\mathsf T}A\mathbf x = \sum_{i,j}A_{ij}x_ix_j Q ( x ) = x T A x = ∑ i , j A ij x i x j 。由于 x i x j = x j x i x_ix_j = x_jx_i x i x j = x j x i ,只有 A A A 的对称部分 1 2 ( A + A T ) \frac12(A + A^{\mathsf T}) 2 1 ( A + A T ) 起作用,所以总可以假定 A A A 是对称矩阵 。例:x 2 + 4 x y + 3 y 2 x^2 + 4xy + 3y^2 x 2 + 4 x y + 3 y 2 对应 A = ( 1 2 2 3 ) A = \begin{pmatrix}1&2\\2&3\end{pmatrix} A = ( 1 2 2 3 ) (交叉项系数平分到两个非对角元上)。
二次型的梯度:∂ Q / ∂ x k = 2 ∑ j A k j x j \partial Q/\partial x_k = 2\sum_jA_{kj}x_j ∂ Q / ∂ x k = 2 ∑ j A kj x j ,即 ∇ Q = 2 A x \nabla Q = 2A\mathbf x ∇ Q = 2 A x ;由此 x ⋅ ∇ Q = 2 Q \mathbf x\cdot\nabla Q = 2Q x ⋅ ∇ Q = 2 Q ,这正是二次齐次函数的欧拉定理(附录 A9 ),§5.5 由它得到能量均分。
正定性 :若对一切 x ≠ 0 \mathbf x\ne0 x = 0 都有 Q > 0 Q>0 Q > 0 ,称 A A A 为正定 的;若 Q ≥ 0 Q\ge0 Q ≥ 0 ,称为半正定 的。判断方法见附录 B2 (所有本征值为正)。两个简单的必要条件:
正定(半正定)矩阵的对角元 A i i > 0 A_{ii}>0 A ii > 0 (≥ 0 \ge0 ≥ 0 ):取 x = e i \mathbf x = \mathbf e_i x = e i 。
半正定时,任意两个指标 i , j i,j i , j 满足 A i i A j j ≥ A i j 2 A_{ii}A_{jj}\ge A_{ij}^2 A ii A jj ≥ A ij 2 :只让 x i x_i x i 、x j x_j x j 不为零,A i i x i 2 + 2 A i j x i x j + A j j x j 2 ≥ 0 A_{ii}x_i^2 + 2A_{ij}x_ix_j + A_{jj}x_j^2\ge0 A ii x i 2 + 2 A ij x i x j + A jj x j 2 ≥ 0 对一切 x i , x j x_i,x_j x i , x j 成立,要求判别式 A i j 2 − A i i A j j ≤ 0 A_{ij}^2 - A_{ii}A_{jj}\le0 A ij 2 − A ii A jj ≤ 0 。
物理中的正定(或半正定)二次型 :动能 1 2 ∑ m i v i 2 \frac12\sum m_iv_i^2 2 1 ∑ m i v i 2 ;稳定平衡附近的势能 1 2 u T K u \frac12\mathbf u^{\mathsf T}K\mathbf u 2 1 u T K u (§6.5 );熵产生 σ s = ∑ i j X i L i j X j ≥ 0 \sigma_s = \sum_{ij}X_iL_{ij}X_j\ge0 σ s = ∑ ij X i L ij X j ≥ 0 (§P5.4 )——由条件 1,热导率、扩散系数、电导率都不能为负;由条件 2,交叉效应(如热电效应)的系数受到 L 12 2 ≤ L 11 L 22 L_{12}^2\le L_{11}L_{22} L 12 2 ≤ L 11 L 22 的限制(这里 L L L 指昂萨格系数矩阵的对称部分)。
自测题
验证 a × ( b × c ) \mathbf a\times(\mathbf b\times\mathbf c) a × ( b × c ) 的 y y y 分量满足 (B1.2) 。
写出 2 x 2 − 2 x y + 5 y 2 2x^2 - 2xy + 5y^2 2 x 2 − 2 x y + 5 y 2 对应的对称矩阵,并判断它是否正定。[答:( 2 − 1 − 1 5 ) \begin{pmatrix}2&-1\\-1&5\end{pmatrix} ( 2 − 1 − 1 5 ) ;2 > 0 2>0 2 > 0 ,2 × 5 − 1 = 9 > 0 2\times5 - 1 = 9>0 2 × 5 − 1 = 9 > 0 ,正定(判据见 §B2.4 )。]
设 a 1 = a 2 ( 0 , 1 , 1 ) \mathbf a_1 = \frac a2(0,1,1) a 1 = 2 a ( 0 , 1 , 1 ) ,a 2 = a 2 ( 1 , 0 , 1 ) \mathbf a_2 = \frac a2(1,0,1) a 2 = 2 a ( 1 , 0 , 1 ) ,a 3 = a 2 ( 1 , 1 , 0 ) \mathbf a_3 = \frac a2(1,1,0) a 3 = 2 a ( 1 , 1 , 0 ) (面心立方的原胞基矢)。求原胞体积。[答:a 3 / 4 a^3/4 a 3 /4 。]
这一篇已记为读完。标为未读