0%

以二元函数f(x,y)f(x,y)为例,设Δx0,Δy0\Delta x\to0,\Delta y\to 0为增量
偏导数为只沿着一个方向的导数,即可以将其他变量视为常数,原函数变为这个位置上的单元函数。
可微可以简单理解成在该点往任何方向进行小的位移,函数变化量接近于偏导数的对应变化量。
f(x+Δx,y+Δy)f(x,y)=fxΔx+fyΔy+o(ρ)f(x+\Delta x,y+\Delta y)-f(x,y)=f'_x\Delta x+f'_y\Delta y+o(\rho)

考虑函数f(x,y)f(x,y)只在xxyy至少有一个为00时为11,其他时候是00。其在(0,0)(0,0)处存在偏导数,但不连续。以绝对值函数为代表的连续函数在0处不可导,因此存在偏导数和连续无任何关系。
f(x+Δx,y+Δy)f(x,y)=f(x+Δx,y+Δy)f(x,y+Δy)+[f(x,y+Δy)f(x,y)]f(x+\Delta x,y+\Delta y)-f(x,y)=f(x+\Delta x,y+\Delta y)-f(x,y+\Delta y)+[f(x,y+\Delta y)-f(x,y)]
即当且仅当f(x+Δx,y+Δy)f(x,y+Δy)=f(x+Δx,y)f(x,y)+o(ρ)f(x+\Delta x,y+\Delta y)-f(x,y+\Delta y)=f(x+\Delta x,y)-f(x,y)+o(\rho)时原函数在(x,y)(x,y)可微。
如果偏导数连续,则上述显然正确。
即:偏导数连续可以推出可微,可微可以推出连续+偏导数存在,剩下两个互不相干。

extra:
一个显然的,连续+偏导数存在,但不可微的例子
f(x,y)=xyf(x,y)=\sqrt{|xy|}(0,0)(0,0)的偏导数均为0且连续,但在不同方向上的斜率不同。

大概是韦达定理
对于多项式方程a0x0+a1x1+a2x2++anxn=0a_0x^0+a_1x^1+a_2x^2+\cdots+a_nx^n=0
若方程的解为 x1,,xnx_1,\ldots,x_n,则有an(xx1)(xx2)(xxn)=0a_n(x-x_1)(x-x_2)\cdots(x-x_n)=0,
此时原方程任取 kk 个解的乘积之和等于(1)kankan(-1)^k\frac{a_{n-k}}{a_n}

举例:当 n=3n=3 时:

  • x1+x2+x3=a2a3x_1+x_2+x_3=-\dfrac{a_2}{a_3}
  • x1x2+x2x3+x1x3=a1a3x_1x_2+x_2x_3+x_1x_3=\dfrac{a_1}{a_3}
  • x1x2x3=a0a3x_1x_2x_3=-\dfrac{a_0}{a_3}

拉格朗日恒等式:

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

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(虽然是打星)…
彻底蜕变为旅游选手啦~

阅读全文 »