892 words
4 minutes
工科数学分析-3.16

5.3.5 高阶偏导数和高阶全微分#

1.高阶偏导数#

z=f(x,y)z=f(x,y)
zx=fx(x,y)\frac{\partial z}{\partial x}=f_{x}(x,y)
x(zx)=2zx2=fxx(x,y)=f11\frac{\partial}{\partial x}(\frac{\partial z}{\partial x})=\frac{\partial^2 z}{\partial x^{2}}=f_{xx}(x,y)=f_{11}
x(zy)=2zxy=fxy(x,y)=f12\frac{\partial}{\partial x}(\frac{\partial z}{\partial y})=\frac{\partial^2 z}{\partial x\partial y}=f_{xy}(x,y)=f_{12}
y(zy)=2zy2=fyy(x,y)=f22\frac{\partial}{\partial y}(\frac{\partial z}{\partial y})=\frac{\partial^2 z}{\partial y^{2}}=f_{yy}(x,y)=f_{22}
y(zx)=2zyx=fyx(x,y)=f21\frac{\partial}{\partial y}(\frac{\partial z}{\partial x})=\frac{\partial^2 z}{\partial y\partial x}=f_{yx}(x,y)=f_{21}

Ex.12 计算z=ex+2yz=e^{x+2y}的所有二阶偏导数#

zx=ex+2y,zy=2ex+2yz_{x}=e^{x+2y},z_{y}=2e^{x+2y}
zxx=ex+2yz_{xx}=e^{x+2y}
zxy=2ex+2yz_{xy}=2e^{x+2y}
zyx=2ex+2yz_{yx}=2e^{x+2y}
zyy=4ex+2yz_{yy}=4e^{x+2y}

Ex.13 计算z=arctan(yx)z=arctan(\frac{y}{x})所有的二阶偏导数#

zx=11+(yx)2(yx2)=yx2+y2,zy=11+(yx)2(1x)=xx2+y2z_{x}=\frac{1}{1+(\frac{y}{x})^{2}}(-\frac{y}{x^{2}})=-\frac{y}{x^{2}+y^{2}},z_{y}=\frac{1}{1+(\frac{y}{x})^{2}}(\frac{1}{x})=\frac{x}{x^{2}+y^{2}}
zxx=2xy(x2+y2)2z_{xx}=-\frac{2xy}{(x^{2}+y^{2})^{2}}
zxy=y2x2(x2+y2)2z_{xy}=\frac{y^{2}-x^{2}}{(x^{2}+y^{2})^{2}}
zyx=y2x2(x2+y2)2z_{yx}=\frac{y^{2}-x^{2}}{(x^{2}+y^{2})^{2}}
zyy=2xy(x2+y2)2z_{yy}=-\frac{2xy}{(x^{2}+y^{2})^{2}}

问:zxy=zyx?z_{xy}=z_{yx}?#

反例:

f(x,y)={xyx2y2x2y2,x2+y200,x2+y2=0f(x,y)= \begin{cases} &xy\frac{x^{2}-y^{2}}{x^{2}-y^{2}}&,x^{2}+y^{2}\neq 0\\ &0&,x^{2}+y^{2}=0 \end{cases}fx(x,y){y(x4y4+4x2y2)(x2+y2)2,x2+y200,x2+y2=0f_{x}(x,y) \begin{cases} &\frac{y(x^{4}-y^{4}+4x^{2}y^{2})}{(x^{2}+y^{2})^{2}}&,x^{2}+y^{2}\neq0\\ &0&,x^{2}+y^{2}=0 \end{cases}fx(x,y){x(x4y4+4x2y2)(x2+y2)2,x2+y200,x2+y2=0f_{x}(x,y) \begin{cases} &\frac{x(x^{4}-y^{4}+4x^{2}y^{2})}{(x^{2}+y^{2})^{2}}&,x^{2}+y^{2}\neq0\\ &0&,x^{2}+y^{2}=0 \end{cases}

fxy(0,0)=1,fyx(0,0)=1f_{xy}(0,0)=-1,f_{yx}(0,0)=1
fxy,fyxf_{xy},f_{yx}在点(x0,y0)(x_{0},y_{0})处连续,则fxy(x0,y0)=fyx(x0,y0)f_{xy}(x_{0},y_{0})=f_{yx}(x_{0},y_{0})

Ex.14 z=xf(2x,y2x)z=xf(2x,\frac{y^{2}}{x}),满足f(u,v)f(u,v)二阶偏导数连续,求zxyz_{xy}#
Method 1:#

zx=f+x(2f1+f2(y2x))=f+2xf1y2xf2z_{x}=f+x(2f_{1}+f_{2}(-\frac{y^{2}}{x}))=f+2xf_{1}-\frac{y^{2}}{x}f_{2}
xxy=f22yx+2xf122yx2yxf2y2xf222yx=4yf122y3x2f22x_{xy}=f_{2}\frac{2y}{x}+2xf_{12}\frac{2y}{x}-\frac{2y}{x}f_{2}-\frac{y^{2}}{x}f_{22}\frac{2y}{x}=4yf_{12}-\frac{2y^{3}}{x^{2}}f_{22}

Method 2:#

zy=xf22yx=2yf2z_{y}=xf_{2}\frac{2y}{x}=2yf_{2}
zyx=2y(2f21+f22(y2x2)=4yf212y3x2f22z_{yx}=2y(2f_{21}+f_{22}(-\frac{y^{2}}{x^{2}})=4yf_{21}-\frac{2y^{3}}{x^{2}}f_{22}

2.高阶全微分#

z=f(x,y)z=f(x,y)的二阶偏导数都连续,则
dz=zxdx+zydydz=z_{x}dx+z_{y}dy
其中,自变量x,yx,y与它们的微分dx,dxdx,dx是互相独立的(从右边来看,变量有四个)
固定dx,dydx,dy,将dzdz作为x,yx,y的函数,定义二阶微分
d2z=d(dz)d^{2}z=d(dz)
对自变量x,yx,y来说,有d(dx)=d(dy)=0d(dx)=d(dy)=0,此时

d2z=d(zxdx+zydy)=d(zx)dx+d(zy)dy=(zxxdx+zxydy)dx+(zyxdx+zyydy)dy=zxx(dx)2+2zxydxdy+zyy(dy)2\begin{align} d^{2}z &=d(z_{x}dx+z_{y}dy)\\ &=d(z_{x})dx+d(z_{y})dy\\ &=(z_{xx}dx+z_{xy}dy)dx+(z_{yx}dx+z_{yy}dy)dy\\ &=z_{xx}(dx)^{2}+2z_{xy}dxdy+z_{yy}(dy)^{2} \end{align}

5.3.6由一个方程确定隐函数的微分法#

F(x,y)=0F(x,y)=0,隐函数存在条件的分析
其一,存在点P0(x0,y0)P_{0}(x_{0},y_{0}),使得F(x0,y0)=0F(x_{0},y_{0})=0,集合Γ=(x,y)R2F(x,y)=0\Gamma={(x,y)\in R^{2}|F(x,y)=0}非空,它是曲面z=F(x,y)z=F(x,y)与平面z=0z=0的交线
Φ(x,y,z)=F(x,y)z\Phi(x,y,z)=F(x,y)-z,则梯度Φ=(Fx,Fy,1)\nabla\Phi=(F_{x},F_{y},-1)与等值面Φ=0\Phi=0正交,设光滑曲线γ(t)=(x(t),y(t),z(t))\gamma(t)=(x(t),y(t),z(t))是曲面z=F(x,y)z=F(x,y)上过点Q0(x0,y0,0)Q_{0}(x_{0},y_{0},0),由Φ(x(t),y(t),z(t))0\Phi(x(t),y(t),z(t))\equiv 0,得
0Φxdxdt+Φydydt+Φzdzdt=Φdγdt0\equiv \Phi_{x}\frac{dx}{dt}+\Phi_{y}\frac{dy}{dt}+\Phi_{z}\frac{dz}{dt}=\nabla\Phi\frac{d\gamma}{dt}
切向量dγdtΦ\frac{d\gamma}{dt}\perp \nabla\Phi,
Q0Q_{0}且垂直于Φ(Q0)\nabla\Phi(Q_{0})的平面称为曲面z=F(x,y)z=F(x,y)在点Q0Q_{0}的切平面,他的方程为
Φ(Q0)(xx0,yy0,zz0)=Fx(P0)(xx0)+Fy(yy0)(zz0)\nabla\Phi(Q_{0})\cdot(x-x_{0},y-y_{0},z-z_{0})=F_{x}(P_{0})(x-x_{0})+F_{y}(y-y_{0})-(z-z_{0})
Q0(x0,y0,0)Q_{0}(x_{0},y_{0},0)Fx(P0)(xx0)+Fy(yy0)z=0F_{x}(P_{0})(x-x_{0})+F_{y}(y-y_{0})-z=0
其二,F(P0)(0,0)\nabla F(P_{0})\neq (0,0),否则曲面z=F(x,y)z=F(x,y)过点Q0Q_{0}的切平面为z=0z=0,此时曲面z=F(x,y)z=F(x,y)与曲面z=0z=0相切于Q0Q_{0},这可能导致交集Γ\Gamma是单点Q0Q_{0},而不是一条曲线
其三,若要求新函数y=y(x)y=y(x)可微,则FF可微时,由F(x,y(x))0F(x,y(x))\equiv 0以及链式法则
Fx+Fydydx=0F_{x}+F_{y}\frac{dy}{dx}=0
Fy(x0,y0)F_{y}(x_{0},y_{0}),则dydxx=x0=Fx(x0,y0)Fy(x0,y0)\frac{dy}{dx}|_{x=x_{0}}=-\frac{F_{x}(x_{0},y_{0})}{F_{y}(x_{0},y_{0})}

定理3.6(隐函数的存在唯一性,可微性)#


(i)F(x,y)F(x,y)在区域DD上联系,P0(x0,y0)P_{0}(x_{0},y_{0})DD的内点
(ii)F(x0,y0)=0F(x_{0},y_{0})=0
(iii)Fy(x0,y0)0F_{y}(x_{0},y_{0})\neq 0
(iv)Fy(x,y)F_{y}(x,y)DD内连续
注:(iii)(iv)表明,存在某个领域U(P0)U(P_{0}),使得FyU(P0)0F_{y}|_{U(P_{0})}\neq 0,不妨设Fy>0F_{y}>0
对某个固定的xx,若F(x,y)=0F(x,y)=0(x,y)U(p0)(x,y)\in U(p_{0})
在点P0P_{0}的某个领域U(B)DU(B)\subset D内,方程F(x,y)=0F(x,y)=0唯一确定一个函数,y=f(x),x(x0δ,x0+δ)y=f(x),x\in (x_{0}-\delta,x_{0}+\delta),使得,
1.f(x0)=y0f(x_{0})=y_{0},当x(x0δ,x0+δ)x\in (x_{0}-\delta,x_{0}+\delta)时,(x,f(x))U(P0),F(x,f(x))=0(x,f(x))\in U(P_{0}),F(x,f(x))=0
2.f(x)f(x)(x0δ,x0+δ)(x_{0}-\delta,x_{0}+\delta)上连续,若Fx(x,y)F_{x}(x,y)也在DD内连续,则f(x)f(x)可导且导数连续,且dydx=Fx(x,y)Fy(x,y)\frac{dy}{dx}=-\frac{F_{x}(x,y)}{F_{y}(x,y)}

Proof:#

只证明Fx,FyF_{x},F_{y}DD内都连续的情形
考虑初值问题,dydx=Fx(x,y)Fy(x,y)\frac{dy}{dx}=-\frac{F_{x}(x,y)}{F_{y}(x,y)}y(x0)=y0y(x_{0})=y_{0}
因为Fx(x,y)Fy(x,y)-\frac{F_{x}(x,y)}{F_{y}(x,y)}P0P_{0}附近连续,由ODE存在性定理,上述问题有解y=y(x)y=y(x),满足
ddxF(x,y(x))=Fx+Fydydx=0\frac{d}{dx}F(x,y(x))=F_{x}+F_{y}\frac{dy}{dx}=0
说明F(x,y(x))=F(x0,y(x0))=F(x0,y0)=0F(x,y(x))=F(x_{0},y(x_{0}))=F(x_{0},y_{0})=0
类似的,由三元方程F(x,y,z)=0F(x,y,z)=0确定一个二元的隐函数,z=f(x,y)z=f(x,y),可微,带入方程得
F(x,y,f(x,y))0F(x,y,f(x,y))\equiv0
取微分
Fxdx+Fydy+Fzdz=0F_{x}dx+F_{y}dy+F_{z}dz=0
Fz0F_{z}\neq0时,dz=FxFzdxFyFzdydz=-\frac{F_{x}}{F_{z}}dx-\frac{F_{y}}{F_{z}}dy

例3.25 设方程ez+xyz=ee^{z}+xyz=e确定的函数z=z(x,y)z=z(x,y),求zx,zy,zxyz_{x},z_{y},z_{xy}#

解:
对原方程取微分得,

0=d(ez)+d(zyz)=ezdz+yzdx+xzdy+xydz\begin{align} 0&=d(e^{z})+d(zyz)\\ &=e^{z}dz+yzdx+xzdy+xydz \end{align}

dz=yzdx+xzdyez+xydz=-\frac{yzdx+xzdy}{e^{z}+xy}
于是 zx=yzez+xyz_{x}=-\frac{yz}{e^{z}+xy}
zy=xzez+xyz_{y}=-\frac{xz}{e^{z}+xy}

zxy=z+yzyez+xy+yzezzy+x(ez+xy)2=zez+xyyez+xyxzez+xy+yzezez+xyxzez+xy+xyz(ez+xy)2=zez+xy+2xyz(ez+xy)2xyz2ez(ez+xy)3=ze2zx2y2z+xyz2ez(ez+xy)3\begin{align} z_{xy}&=-\frac{z+yz_{y}}{e^{z}+xy}+yz\frac{e^{z}z_{y}+x}{(e^{z}+xy)^{2}}\\ &=-\frac{z}{e^{z}+xy}-\frac{y}{e^{z}+xy}\frac{-xz}{e^{z}+xy}+\frac{yze^{z}}{e^{z}+xy}\frac{-xz}{e^{z}+xy}+\frac{xyz}{(e^{z}+xy)^{2}}\\ &=-\frac{z}{e^{z}+xy}+\frac{2xyz}{(e^{z}+xy)^{2}}-\frac{xyz^{2}e^{z}}{(e^{z}+xy)^{3}}\\ &=-\frac{ze^{2z}-x^{2}y^{2}z+xyz^{2}e^{z}}{(e^{z}+xy)^{3}} \end{align}

5.4多元函数的TaylorTaylor公式与极值#

5.4.1 TaylorTaylor公式#

中值定理:#

f(x,y)f(x,y)在凸区域DR2D\subset R^{2}上可微,则对DD中任意两点(x0,y0),(x0+Δx,y+Δy)(x_{0},y_{0}),(x_{0}+\Delta x,y+\Delta y),存在θ(0,1)\theta\in(0,1)使得
f(x0+Δx,y0+Δy)f(x0,y0)=fx(x0+θΔx,y0+θΔy)Δx+fy(x0+θΔx,y0+θΔy)Δyf(x_{0}+\Delta x,y_{0}+\Delta y)-f(x_{0},y_{0})=f_{x}(x_{0}+\theta\Delta x,y_{0}+\theta\Delta y)\Delta x+f_{y}(x_{0}+\theta\Delta x,y_{0}+\theta\Delta y)\Delta y

Proof#

φ(t)=f(x0+tΔx)\varphi(t)=f(x_{0}+t\Delta x),在[0,1][0,1]上可微,用拉格朗日定理,存在θ(0,1)\theta\in(0,1),有
φ(1)φ(0)=φ(θ)=fx(x0+θΔx,y0+θΔy)Δx+fy(x0+θΔx,y0+θΔy)Δy\varphi(1)-\varphi(0)=\varphi'(\theta)=f_{x}(x_{0}+\theta\Delta x,y_{0}+\theta\Delta y)\Delta x+f_{y}(x_{0}+\theta\Delta x,y_{0}+\theta\Delta y)\Delta y

Ex.1#

f(x,y)f(x,y)的二阶偏导数连续,写出一元函数Φ(t)=f(x0+tΔx,y0+tΔy)\Phi(t)=f(x_{0}+t\Delta x,y_{0}+t\Delta y)t=0t=0的二阶泰勒公式(PeanoPeano余项)

Φ(t)=Φ(0)+Φ(0)t+12Φ(0)(t2)+o(t2)\Phi(t)=\Phi(0)+\Phi'(0)t+\frac{1}{2}\Phi''(0)(t^{2})+o(t^{2})
其中Φ(t)=fxΔx+fyΔy\Phi'(t)=f_{x}\Delta x+f_{y}\Delta y

Φ(t)=(fxxΔx+fxyΔy)Δx+(fyxΔx+fyyΔy)Δ=fxx(Δx)2+fxyΔxΔy+fyy(Δy)2=(Δx,Δy)(fxxfxyfxyfyy)(Δx,Δy)T\begin{align} \Phi''(t)&=(f_{xx}\Delta x+f_{xy}\Delta y)\Delta x+(f_{yx}\Delta x+f_{yy}\Delta y)\Delta\\ &=f_{xx}(\Delta x)^{2}+f_{xy}\Delta x\Delta y+f_{yy}(\Delta y)^{2}\\ &=(\Delta x,\Delta y)\begin{pmatrix}f_{xx}&f_{xy}\\ f_{xy}&f_{yy}\end{pmatrix}(\Delta x,\Delta y)^{T} \end{align}

X0=(x0,y0),ΔX=(Δx,Δy)X_{0}=(x_{0},y_{0}),\Delta X=(\Delta x,\Delta y),则
f(X0+ΔX)=f(X0)+f(X0,ΔX)+12ΔXHf(X0)(ΔX)T+o(ΔX2)f(X_{0}+\Delta X)=f(X_{0})+\langle\nabla f(X_{0},\Delta X) \rangle+\frac{1}{2}\Delta X\cdot H_{f}(X_{0})\cdot(\Delta X)^{T}+o(||\Delta X||^{2})
HesseHesse矩阵 Hf(X0)=(fxxfxyfxyfyy)(X0)H_{f}(X_{0})=\begin{pmatrix}f_{xx}&f_{xy}\\f_{xy}&f_{yy}\end{pmatrix}(X_{0})
对自变量x,yx,y的一阶微分dx,dydx,dydx=Δx,dy=Δydx=\Delta x,dy=\Delta y
则有 f(x0+dx,y0+dy)=f(x0,y0)+12d2z+o((dx)2+(dy)2)f(x_{0}+dx,y_{0}+dy)=f(x_{0},y_{0})+\frac{1}{2}d^{2}z+o((dx)^{2}+(dy)^{2})

例4.1#

z=z(z,y)z=z(z,y)是由方程z32xz+y=0z^{3}-2xz+y=0确定的隐函数且z(1,1)=1z(1,1)=1,求z(x,y)z(x,y)(1,1)(1,1)处的二阶TaylorTaylor公式(PeanoPeano余项)
(1,1)(1,1)带入有三个实根为1,1+52,1521,\frac{-1+\sqrt{5}}{2},\frac{-1-\sqrt{5}}{2}

Method 1#

两边取微分
3z2dz2(zdx+xdz)+dy=03z^{2}dz-2(zdx+xdz)+dy=0
dz=2zdxdy3z22xdz=\frac{2zdx-dy}{3z^{2}-2x}
zx=2z3z22xz_{x}=\frac{2z}{3z^{2}-2x}
zy=13z22xz_{y}=\frac{-1}{3z^{2}-2x}
z(1,1)=1,zx(1,1)=2,zy(1,1)=1z(1,1)=1,z_{x}(1,1)=2,z_{y}(1,1)=-1
zxx=2(3z2+2x)zx(3z22x)2z_{xx}=\frac{-2(3z^{2}+2x)z_{x}}{(3z^{2}-2x)^{2}}
zxy=2(3z2+2x)zy(3z22x)2z_{xy}=\frac{-2(3z^{2}+2x)z_{y}}{(3z^{2}-2x)^{2}}
zxx=6zzy(3z22x)2z_{xx}=\frac{6zz_{y}}{(3z^{2}-2x)^{2}}

(x,y,z)=(1,1,1)(x,y,z)=(1,1,1)时,zxx=16,zxy=10,zyy=6z_{xx}=-16,z_{xy}=10,z_{yy}=-6,则
z(1+x,1+y)=1+2xy+12(16x2+10xy6y2)+o(x2+y2)z(1+x,1+y)=1+2x-y+\frac{1}{2}(-16x^{2}+10xy-6y^{2})+o(x^{2}+y^{2})

Method 2#

两边取微分
3z2dz2(zdx+xdz)+dy=03z^{2}dz-2(zdx+xdz)+dy=0
(3z22x)dz2zdx+dy=0()(3z^{2}-2x)dz-2zdx+dy=0(*)
z(1,1)=1z(1,1)=1dz(1,1,1)=2dxdydz|_{(1,1,1)}=2dx-dy
()(*)取微分得,

0=d(3z22x)dz+(3z22x)d2z2dzdx=6z(dz)2+(3z22x)d2z4dzdx\begin{align} 0&=d(3z^{2}-2x)dz+(3z^{2}-2x)d^{2}z-2dzdx\\ &=6z(dz)^{2}+(3z^{2}-2x)d^{2}z-4dzdx \end{align}

于是 d2z(1,1,1)=4dx(2dxdy)6(2dxdy)2=16(dx)2+20dxdy6(dy)2d^{2}z|_{(1,1,1)}=4dx(2dx-dy)-6(2dx-dy)^{2}=-16(dx)^{2}+20dxdy-6(dy)^{2}