0%

拉格朗日恒等式:

计算两个向量夹平行四边形的面积

a2b2(ab)2=1i<jn(aibjajbi)2where(ab)2=a2b2cos2θ\| \mathbf{a} \|^2 \| \mathbf{b} \|^2 - (\mathbf{a} \cdot \mathbf{b})^2 = \sum_{1 \le i < j \le n} (a_i b_j - a_j b_i)^2 \\ \quad \text{where} \quad (\mathbf{a} \cdot \mathbf{b})^2 = \| \mathbf{a} \|^2 \| \mathbf{b} \|^2 \cos^2 \theta

即对于任意维度的两个向量a\mathbf{a}b\mathbf{b},可以通过a2b2(ab)2\sqrt{\| \mathbf{a} \|^2 \| \mathbf{b} \|^2 - (\mathbf{a} \cdot \mathbf{b})^2}计算出两个向量夹平行四边形的面积,在面积分求面积微元时很好用。

Extra

nn维空间里,通过mm个参数可以唯一确定一个mm维流形(e.g. n=3,m=2n=3,m=2时,我们通过两个参数u,vu,v在三维空间(x(u,v),y(u,v),z(u,v)x(u,v),y(u,v),z(u,v))中唯一确定一个二维流形,即曲面)。

这一mm维流形的超体积微元的平方等于:所有(nm)\binom{n}{m}种在nn个维度中任选mm个维度,将微元投影到这些维度上的超体积的平方和。具体地,设微分矩阵 J\mathbf{J}(其中 Ji,j=xiuj\mathbf{J}_{i,j}=\frac{\partial x_i}{\partial u_j}),则超体积微元dS=Kdu1du2dum\mathrm{d}S=\sqrt{K}\mathrm{d}u_1\mathrm{d}u_2\cdots\mathrm{d}u_m,其中KKJ\mathbf{J}中所有m×mm\times m子矩阵的行列式的平方和。

例如,在n=3,m=2n=3,m=2时,J=(xuxvyuyvzuzv)\mathbf{J}=\begin{pmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \\ \frac{\partial z}{\partial u} & \frac{\partial z}{\partial v} \end{pmatrix},则K=xuxvyuyv2+xuxvzuzv2+yuyvzuzv2K=\begin{vmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \end{vmatrix}^2+\begin{vmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial z}{\partial u} & \frac{\partial z}{\partial v} \end{vmatrix}^2+\begin{vmatrix} \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \\ \frac{\partial z}{\partial u} & \frac{\partial z}{\partial v} \end{vmatrix}^2

拉格朗日恒等式的面积项1i<jn(aibjajbi)2\sum\limits_{1 \le i < j \le n} (a_i b_j - a_j b_i)^2等同于在m=2m=2时累加所有子行列式的平方和。

Here's something encrypted, password is required to continue reading.
阅读全文 »

曲率是衡量曲线在某一点弯曲程度的值。具体地,曲率的定义为某一点处,单位位移距离,切线角度的变化率,即κ=dθds\kappa = \frac{d\theta}{ds}

曲率的倒数为曲率半径,即R=1κR = \frac{1}{\kappa},代表在该点和曲线相切的圆的半径。

对于参数方程曲线r(t)=(x(t),y(t))\mathbf{r}(t) = (x(t), y(t)),其曲率的计算公式为κ=r(t)×r(t)r(t)3\kappa = \frac{|\mathbf{r}'(t) \times \mathbf{r}''(t)|}{|\mathbf{r}'(t)|^3},以下是推导:

令速度向量为v(t)=r(t)\mathbf{v}(t) = \mathbf{r}'(t),则加速度向量为a(t)=r(t)\mathbf{a}(t) = \mathbf{r}''(t)

速度向量=速率*切线单位向量,即v(t)=VT(t)\mathbf{v}(t) = \mathbf{V}\mathbf{T}(t),其中T(t)\mathbf{T}(t)为切线单位向量,V\mathbf{V}为速率。

对速度向量求导,得到加速度向量a(t)=dvdt=dVdtT+VdTdt\mathbf{a}(t) = \frac{d\mathbf{v}}{dt} = \frac{d\mathbf{V}}{dt}\mathbf{T} + \mathbf{V}\frac{d\mathbf{T}}{dt},其中dVdt\frac{d\mathbf{V}}{dt}为切向加速度,VdTdt\mathbf{V}\frac{d\mathbf{T}}{dt}为法向加速度。切向加速度影响速率变化的大小,法向加速度影响速率变化的方向。

根据曲率的定义,有dT=Ndθ=κNdsd\mathbf{T} = \mathbf{N}d\theta = \kappa \mathbf{N} ds,其中N\mathbf{N}为法线单位向量。这是因为单位向量转动时,方向为法线方向,位移长度为转动的角度。

基于上述,有a(t)=dvdt=dVdtT+VdTdt=dVdtT+VκNdsdt=dVdtT+V2κN\mathbf{a}(t) = \frac{d\mathbf{v}}{dt} = \frac{d\mathbf{V}}{dt}\mathbf{T} + \mathbf{V}\frac{d\mathbf{T}}{dt} = \frac{d\mathbf{V}}{dt}\mathbf{T} + \mathbf{V}\kappa \mathbf{N} \frac{ds}{dt} = \frac{d\mathbf{V}}{dt}\mathbf{T} + \mathbf{V}^2\kappa \mathbf{N}

由于切向和法向加速度向量正交,因此我们可以对a(t)\mathbf{a}(t)叉乘T\mathbf{T}

a(t)×T=dVdtT×T+VκNdsdt×T=V2κN×T\mathbf{a}(t) \times \mathbf{T} = \frac{d\mathbf{V}}{dt}\mathbf{T} \times \mathbf{T} + \mathbf{V}\kappa \mathbf{N} \frac{ds}{dt} \times \mathbf{T} = \mathbf{V}^2\kappa \mathbf{N} \times \mathbf{T}

由于相互垂直的单位向量叉乘为1,平行的向量叉乘为0,因此有:
κ=a(t)×TV2=a(t)×v(t)V3=r(t)×r(t)r(t)3\kappa = \frac{\mathbf{a}(t) \times \mathbf{T}}{\mathbf{V}^2}=\frac{\mathbf{a}(t) \times \mathbf{v}(t)}{\mathbf{V}^3}=\frac{|\mathbf{r}'(t) \times \mathbf{r}''(t)|}{|\mathbf{r}'(t)|^3}

留下些文字,作为我曾经存在过的证明。

阅读全文 »

随想

Rf(x,y)dxdy=Sf(aRcos(θ),bRsin(θ))(x,y)(R,θ)dRdθ=Sf(aRcos(θ),bRsin(θ))xRxθyRyθdRdθ=Sf(aRcos(θ),bRsin(θ))acos(θ)aRsin(θ)bsin(θ)bRcos(θ)dRdθ=Sf(aRcos(θ),bRsin(θ))(acos(θ))(bRcos(θ))(aRsin(θ))(bsin(θ))dRdθ=Sf(aRcos(θ),bRsin(θ))abRcos2(θ)+abRsin2(θ)dRdθ=Sf(aRcos(θ),bRsin(θ))abR(cos2(θ)+sin2(θ))dRdθ=Sf(aRcos(θ),bRsin(θ))abRdRdθ(假设 a,b,R0)\begin{aligned} \iint_R f(x,y)\,\mathrm{d} x\,\mathrm{d} y &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left|\frac{\partial(x,y)}{\partial(R,\theta)}\right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left| \begin{vmatrix} \frac{\partial x}{\partial R} & \frac{\partial x}{\partial \theta} \\ \frac{\partial y}{\partial R} & \frac{\partial y}{\partial \theta} \end{vmatrix} \right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left| \begin{vmatrix} a\cos(\theta) & -aR\sin(\theta) \\ b\sin(\theta) & bR\cos(\theta) \end{vmatrix} \right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left| (a\cos(\theta))(bR\cos(\theta)) - (-aR\sin(\theta))(b\sin(\theta)) \right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left| abR\cos^2(\theta) + abR\sin^2(\theta) \right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta))\left| abR(\cos^2(\theta) + \sin^2(\theta)) \right|\,\mathrm{d} R\,\mathrm{d} \theta \\ &= \iint_S f(aR\cos(\theta),bR\sin(\theta)) abR \, \mathrm{d} R\,\mathrm{d} \theta \quad (\text{假设 } a,b,R \ge 0) \end{aligned}

这里可以理解为坐标轴改变,即假设以(θ,R\theta,R)为新坐标轴,原先的x,yx,y就变成了x,y\partial x,\partial y的新坐标,新的面积即为他们的叉乘(即二阶行列式)


当函数f(x)f(x)满足axb,f(x)=f(a+bx)\forall a\leq x \leq b,f(x)=f(a+b-x)(即积分区间内函数对称)时,有

abxf(x)dx=ab(a+bx)f(x)dx=a+b2abf(x)dxe.g.0πxsin(x)dx=π20πsin(x)dx\begin{aligned} \int_a^b xf(x)\,\mathrm{d} x &=\int_a^b(a+b-x)f(x)\,\mathrm{d} x\\ &=\frac{a+b}{2}\int_a^b f(x)\,\mathrm{d} x \\ \text{e.g.} \quad \int_{0}^{\pi} x\sin(x)\,\mathrm{d} x &=\frac{\pi}{2}\int_{0}^\pi \sin(x)\,\mathrm{d} x \end{aligned}

没有包袱的一场XCPC(虽然是打星)…
彻底蜕变为旅游选手啦~

阅读全文 »

数学公式渲染测试

ex=n=0xnn!=10!+x11!+x22!+x33!+e^x = \sum_{n=0}^{\infty} \frac{x^n}{n!} = \frac{1}{0!} + \frac{x^1}{1!} + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots

sinx=n=0(1)nx2n+1(2n+1)!=x11!x33!+x55!\sin x = \sum_{n=0}^{\infty} (-1)^n \frac{x^{2n+1}}{(2n+1)!} = \frac{x^1}{1!} - \frac{x^3}{3!} + \frac{x^5}{5!} - \cdots

cosx=n=0(1)nx2n(2n)!=10!x22!+x44!\cos x = \sum_{n=0}^{\infty} (-1)^n \frac{x^{2n}}{(2n)!} = \frac{1}{0!} - \frac{x^2}{2!} + \frac{x^4}{4!} - \cdots

ln(1+x)=n=1(1)n+1xnn=x11x22+x33\ln(1+x)=\sum_{n=1}^{\infty}(-1)^{n+1}\frac{x^n}{n}=\frac{x^1}{1}-\frac{x^2}{2}+\frac{x^3}{3}-\cdots

11x=n=0xn=x0+x1+x2+\frac{1}{1-x}=\sum_{n=0}^{\infty}x^n=x^0+x^1+x^2+\cdots

arctanx=n=0(1)nx2n+12n+1=x11x33+x55\arctan x=\sum_{n=0}^{\infty}(-1)^n\frac{x^{2n+1}}{2n+1}=\frac{x^1}{1}-\frac{x^3}{3}+\frac{x^5}{5}-\cdots

泰勒展开公式

f(x)k处展开:f(x)=n=0f(n)(k)n!(xk)nf(x)在k处展开:f(x)=\sum_{n=0}^{\infty}\frac{f^{(n)}(k)}{n!}(x-k)^n

e.g.sinx=n=0sin(n)0n!xne.g. \sin x=\sum_{n=0}^{\infty}\frac{\sin^{(n)}0}{n!}x^n

4n+1阶导为1,4n+3阶导为-1,2n阶导为0.