1027 words
5 minutes
工科数学分析-3.30
方法一:#

曲线
Γ:r(t)=(x(t),y(t),z(t))\Gamma:r(t)=(x(t),y(t),z(t))
切向量
r′(t)=(x′(t),y′(t),z′(t))r'(t)=(x'(t),y'(t),z'(t))
过P0P_0且垂直于r′(t0)r'(t_0)的平面称为Γ\Gamma在P0P_0的法平面
x′(t0)(x−x0)+y′(t0)(y−y0)+z′(t0)(z−z0)=0x'(t_{0})(x-x_0)+y'(t_{0})(y-y_0)+z'(t_{0})(z-z_0)=0
曲线Γ\Gamma是两个曲面的交线,方程为
F(x,y,z)=0,G(x,y,z)=0(∗)F(x,y,z)=0,G(x,y,z)=0(*)
设(∗)(*)在点p0(x0,y0,z0)p_0(x_0,y_0,z_0)满足隐函数组定理的条件,其中,
J(P0)=∂(F,G)∂(x,y)∣P0≠0J(P_0)=\frac{\partial(F,G)}{\partial(x,y)}|_{P_0}\neq0
则(∗)(*)在P0P_0附近唯一确定隐函数组x=x(z),y=y(z)x=x(z),y=y(z),且

{Fxdx+Fydy=−FzdzGxdx+Gydy=−Gzdz\begin{align} \begin{cases} F_xdx+F_ydy=-F_zdz\\ G_xdx+G_ydy=-G_zdz\\ \end{cases}\\ \end{align}

解出

dxdz=1J∣−FzFy−GzGy∣=1J∂(F,G)∂(y,z)dydz=1J∣Fx−FzGx−Gz∣=1J∂(F,G)∂(z,x)\begin{align} \frac{dx}{dz}&=\frac{1}{J} \begin{vmatrix} -F_z&F_y\\ -G_z&G_y\\ \end{vmatrix} &=\frac{1}{J}\frac{\partial(F,G)}{\partial(y,z)}\\ \frac{dy}{dz}&=\frac{1}{J} \begin{vmatrix} F_x&-F_z\\ G_x&-G_z\\ \end{vmatrix} &=\frac{1}{J}\frac{\partial(F,G)}{\partial(z,x)}\\ \end{align}

曲线Γ\Gamma在点P0P_0附近有参数表示 r(z)=(x(z),y(z),z),r′(z)=(x′(z),y′(z),1)r(z)=(x(z),y(z),z),r'(z)=(x'(z),y'(z),1)
过点P0P_0的切线方程为
x−x0dxdz(z0)=y−y0dydz(z0)=z−z01\frac{x-x_0}{\frac{dx}{dz}(z_0)}=\frac{y-y_0}{\frac{dy}{dz}(z_0)}=\frac{z-z_0}{1}
即
x−x0A=y−y0B=z−z0C\frac{x-x_0}{A}=\frac{y-y_0}{B}=\frac{z-z_0}{C}
其中 A=∂(F,G)∂(y,z)∣P0,B=∂(F,G)∂(z,x)∣P0,C=∂(F,G)∂(x,y)∣P0A=\frac{\partial(F,G)}{\partial(y,z)}|_{P_0},B=\frac{\partial(F,G)}{\partial(z,x)}|_{P_0},C=\frac{\partial(F,G)}{\partial(x,y)}|_{P_0}
法平面方程 A(x−x0)+B(y−y0)+C(z−z0)=0A(x-x_0)+B(y-y_0)+C(z-z_0)=0

方法二:#

曲面S1:F(x,y,z)=0S_1:F(x,y,z)=0
曲面S2:G(x,y,z)=0S_2:G(x,y,z)=0
Γ=S1∩S2\Gamma=S_1\cap S_2
对上述两式子微分得

{Fxdx+Fydy+Fzdz=0=∇F⋅dr,∇F⊥drGxdx+Gydy+Gzdz=0∇G⋅dr,∇G⊥dr\begin{align} \begin{cases} F_xdx+F_ydy+F_zdz=0=\nabla F\cdot dr,\nabla F\perp dr\\ G_xdx+G_ydy+G_zdz=0\nabla G\cdot dr,\nabla G\perp dr\\ \end{cases}\\ \end{align}

于是切向量

dr//∇F×∇G=∣ijkFxFyFzGxGyGz∣=(A,B,C)\begin{align} dr//\nabla F\times \nabla G= \begin{vmatrix} i&j&k\\ F_x&F_y&F_z\\ G_x&G_y&G_z\\ \end{vmatrix} =(A,B,C) \end{align}

例6.2#

求曲线

{2x2+y2+z2=45x2+2y2=z\begin{align} \begin{cases} 2x^2+y^2+z^2=45\\ x^2+2y^2=z \end{cases} \end{align}

在点P0(−2,1,6)P_0(-2,1,6)处的切线与法平面
解:

∇F(P0)=(4x,2y,2z)∣P0=(−8,2,12)∇G(P0)=(2x,4y,−1)∣P0=(−4,4,−1)\begin{align} \nabla F(P_0)&=(4x,2y,2z)|_{P_0}=(-8,2,12)\\ \nabla G(P_0)&=(2x,4y,-1)|_{P_0}=(-4,4,-1)\\ \end{align}

所以切向量

dr//∇F(P0)×∇G(P0)=∣ijk−8212−44−1∣=(−50,−56,−24)\begin{align} dr//\nabla F(P_0)\times \nabla G(P_0)= \begin{vmatrix} i&j&k\\ -8&2&12\\ -4&4&-1\\ \end{vmatrix} =(-50,-56,-24) \end{align}

切线方程
x+225=y−128=z−612\frac{x+2}{25}=\frac{y-1}{28}=\frac{z-6}{12}
法平面方程为
25(x+2)+28(y−1)+12(z−6)=025(x+2)+28(y-1)+12(z-6)=0
即 25x+28y+12z−50=025x+28y+12z-50=0

5.6.2弧长#

1.定义与计算#

定义6.1(弧长)#

设简单曲线Γ=\Gamma=弧ABAB的参数方程为
r(t)=(x(t),y(t),z(t)),α≤t≤βr(t)=(x(t),y(t),z(t)),\alpha\le t\le\beta
其中A=r(α),B=r(β)A=r(\alpha),B=r(\beta),在Γ\Gamma上从AA到BB引入分点
折现长度为sn=∑i=1n∣∣Pi−1,Pi∣∣s_{n}=\sum_{i=1}^{n}||P_{i-1},P_i||
称T={Pi}T=\{P_{i}\}为曲线Γ\Gamma的一个分割,d=max⁡1≤i≤n∣∣Pi−1Pi∣∣d=\max_{1\le i\le n}||P_{i-1}P_i||称为细度
若对任何分割,和式极限
lim⁡d→0sn=lim⁡d→0∑i=1n∣∣Pi−1Pi∣∣=s\lim_{d\to0}s_{n}=\lim_{d\to 0}\sum_{i=1}^{n}||P_{i-1}P_i||=s
存在,且与分割TT的选取无关,则称上述极限值ss为曲线Γ\Gamma的弧长
此时称Γ\Gamma是可求长曲线
注:某些分形为不可求长曲线

定理6.1#

设在[α,β][\alpha,\beta]上r′(t)r'(t)连续且r′(t)≠0r'(t)\neq0,则曲线r=r(t)r=r(t)是可求长曲线,且长度为
s=∫αβ∣∣r′(t)∣∣dt=∫αβx′(t)2+y′(t)2+z′(t)2dts=\int_{\alpha}^{\beta}||r'(t)||dt=\int_{\alpha}^{\beta}\sqrt{x'(t)^{2}+y'(t)^{2}+z'(t)^{2}}dt

proof#

对1≤i≤n1\le i\le n,设Pi=r(ti)P_i=r(t_i)
∣∣Pi−1Pi∣∣2=(x(ti)−x(ti−1))2+(y(ti)−y(ti−1))2+(z(ti)−z(ti−1))2||P_{i-1}P_i||^{2}=(x(t_{i})-x(t_{i-1}))^{2}+(y(t_{i})-y(t_{i-1}))^{2}+(z(t_{i})-z(t_{i-1}))^{2}
依条件,由拉格朗日中值定理,分别存在ξi,ηi,ζi∈[ti−1,ti]\xi_{i},\eta_{i},\zeta_{i}\in[t_{i-1},t_{i}],使得

x(ti)−x(ti−1)=x′(ξi)Δtiy(ti)−y(ti−1)=y′(ηi)Δtiz(ti)−z(ti−1)=z′(ζi)Δti\begin{align} x(t_i)-x(t_{i-1})=x'(\xi_{i})\Delta t_{i}\\ y(t_i)-y(t_{i-1})=y'(\eta_{i})\Delta t_{i}\\ z(t_i)-z(t_{i-1})=z'(\zeta_{i})\Delta t_{i}\\ \end{align}

于是 ∣∣Pi−1Pi∣∣=x′(ξi)2+y′(ηi)2+z′(ζi)2Δti||P_{i-1}P_i||=\sqrt{x'(\xi_{i})^{2}+y'(\eta_{i})^{2}+z'(\zeta_{i})^{2}}\Delta t_{i}

sn=∑i=1nx′(ξi)2+y′(ηi)2+z′(ζi)2Δti=∑i=1n∣∣r′(τi)∣∣Δti+σ,τi∈(ti−1,ti)σ=∑i=1n(x′(ξi)2+y′(ηi)2+z′(ζi)2−x′(τi)2+y′(τi)2+z′(τi)2)Δti≤∑i=1n(x′(ξi)−x′(τi))2+(y′(ηi)−y′(τi))2+(z′(ζi)−z′(τi))2Δti\begin{align} s_{n}&=\sum_{i=1}^{n}\sqrt{x'(\xi_{i})^{2}+y'(\eta_{i})^{2}+z'(\zeta_{i})^{2}}\Delta t_{i}\\ &=\sum_{i=1}^{n}||r'(\tau_{i})||\Delta t_{i}+\sigma,\tau_i\in(t_{i-1},t_i)\\ \sigma&=\sum_{i=1}^{n}(\sqrt{x'(\xi_{i})^{2}+y'(\eta_{i})^{2}+z'(\zeta_{i})^{2}}-\sqrt{x'(\tau_{i})^{2}+y'(\tau_{i})^{2}+z'(\tau_{i})^{2}})\Delta t_{i}\\ &\le\sum_{i=1}^{n}\sqrt{(x'(\xi_{i})-x'(\tau_{i}))^{2}+(y'(\eta_{i})-y'(\tau_{i}))^{2}+(z'(\zeta_{i})-z'(\tau_{i}))^{2}}\Delta t_{i} \end{align}

下面证明lim⁡d→0σ=0\lim_{d\to0}\sigma=0
以为r′(t)r'(t)在[α,β][\alpha,\beta]上连续,且r′(t)≠0r'(t)\neq0,所以存在常数m>0m>0,使
x′(ξi)2+y′(ηi)2+z′(ζi)2≥m,∣∣Pi−1Pi∣∣≥mΔti\sqrt{x'(\xi_{i})^{2}+y'(\eta_{i})^{2}+z'(\zeta_{i})^{2}}\ge m,||P_{i-1}P_{i}||\ge m\Delta t_{i}
令Δt=max⁡1≤i≤nΔtt\Delta t=\max_{1\le i\le n}\Delta t_{t},则ΔT≤dm\Delta T\le \frac{d}{m},当d→0d\to 0时,Δt→0\Delta t\to 0
ξi,ηi,ζi,τi∈[ti−1,ti]\xi_{i},\eta_{i},\zeta_{i},\tau_{i}\in[t_{i-1},t_{i}]
依条件中x′(t),y′(t),z′(t)x'(t),y'(t),z'(t)在[α,β][\alpha,\beta]上连续,从而一致连续
对任何ε>0\varepsilon>0,存在δ>0\delta>0,当t,t′∈[α,β]t,t'\in[\alpha,\beta]且∣t−t′∣<δ|t-t'|<\delta时
有∣x′(t)−x′(t′)∣<ε,∣y′(t)−y′(t′)∣<ε,∣z′(t)−z′(t′)∣<ε|x'(t)-x'(t')|<\varepsilon,|y'(t)-y'(t')|<\varepsilon,|z'(t)-z'(t')|<\varepsilon
当d<mδd<m\delta时,Δt<dm<δ\Delta t<\frac{d}{m}<\delta,此时ξi,ηi,ζi,τi∈[ti−1,ti]\xi_{i},\eta_{i},\zeta_{i},\tau_{i}\in[t_{i-1},t_{i}]
所以 σ≤∑i=1n3ε2Δti=3ε∑i=1nΔti=3ε(β−α)\sigma\le\sum_{i=1}^{n}\sqrt{3\varepsilon^{2}}\Delta t_{i}=\sqrt{3}\varepsilon\sum_{i=1}^{n}\Delta t_{i}=\sqrt{3}\varepsilon(\beta-\alpha)
因此

lim⁡d→0sn=lim⁡d→0∑i=1n∣∣r′(τi)∣∣Δti+0,τi∈(ti−1,ti)=lim⁡Δt→0∑i=1n∣∣r′(τi)∣∣Δti=∫αβ∣∣r′(t)∣∣dt\begin{align} \lim_{d\to 0}s_{n}&=\lim_{d\to0}\sum_{i=1}^{n}||r'(\tau_{i})||\Delta t_{i}+0,\tau_i\in(t_{i-1},t_i)\\ &=\lim_{\Delta t\to0}\sum_{i=1}^{n}||r'(\tau_{i})||\Delta t_{i}\\ &=\int_{\alpha}^{\beta}||r'(t)||dt \end{align}

2.弧微分与自然参数#

因为ds=(dx)2+(dy)2+(dz)2ds=\sqrt{(dx)^{2}+(dy)^{2}+(dz)^{2}},所以用弧长ss做参数时,
drds=(dxds,dyds,dzds)\frac{dr}{ds}=(\frac{dx}{ds},\frac{dy}{ds},\frac{dz}{ds})
是单位切向量,设它与x,y,zx,y,z舟正方向的夹角分别为α,β,γ\alpha,\beta,\gamma则

drds=(cos⁡α,cos⁡β,cos⁡γ)\frac{dr}{ds}=(\cos{\alpha},\cos\beta,\cos\gamma)

dx=cos⁡αds,dy=cos⁡βds,dz=cos⁡γdsdx=\cos{\alpha}ds,dy=\cos{\beta}ds,dz=\cos{\gamma}ds

例6.5#

求螺旋新x=acos⁡θ,y=asin⁡θ,z=kθx=a\cos\theta,y=a\sin\theta,z=k\theta,一个螺距之间的长度
解:
r′(θ)=(−asin⁡θ,acos⁡θ,k)r'(\theta)=(-a\sin\theta,a\cos\theta,k)
∣∣r′(θ)∣∣=a2+k2||r'(\theta)||=\sqrt{a^{2}+k^2}
s=∫02π∣∣r′(θ)∣∣dθ=2πa2+k2s=\int_{0}^{2\pi}||r'(\theta)||d\theta=2\pi\sqrt{a^{2}+k^{2}}

5.6.3 曲面的切平面与法线#

1.曲线的参数方程#

球面的参数方程 设P(x,y,z)P(x,y,z)为球面一点,他在xOyxOy平面内的投影为M(x,y,0)M(x,y,0)
则
z=±∣∣PM∣∣=∣∣OP∣∣cos⁡φ,x=∣∣OM∣∣cos⁡θ,y=∣∣OM∣∣sin⁡θ,∣∣OM∣∣=∣∣OP∣∣sin⁡φz=\pm||PM||=||OP||\cos\varphi,x=||OM||\cos\theta,y=||OM||\sin\theta,||OM||=||OP||\sin\varphi
所以
z=Rcos⁡φ,x=Rsin⁡φcos⁡θ,y=Rsin⁡φsin⁡θz=R\cos\varphi,x=R\sin\varphi\cos\theta,y=R\sin\varphi\sin\theta
其中 0≤φ≤π,0≤θ≤2π0\le\varphi\le\pi,0\le\theta\le2\pi
一般地,曲面SS的参数方程
r(u,v)=(x(u,v),y(u,v),z(u,v)),(u,v)∈△,r:△⊂R2→R3r(u,v)=(x(u,v),y(u,v),z(u,v)),(u,v)\in \triangle,r:\triangle\subset R^2\to R^3

2.曲面上的曲线表示#

在曲面SS上,考虑两类曲面
(1)v=v0(1)v=v_0,称为uu曲线
(2)u=u0(2)u=u_0,称为vv曲线
他们的切向量分别为
ru=(xu,yu,zu)r_{u}=(x_u,y_u,z_u)
rv=(xv.yv,zv)r_v=(x_v.y_v,z_v)

3.曲面的切平面与法线#

设曲面SS的参数方程为
r(u,v)=(x(u,v),y(u,v),z(u,v)),(u,v)∈△,r:△⊂R2→R3r(u,v)=(x(u,v),y(u,v),z(u,v)),(u,v)\in \triangle,r:\triangle\subset R^2\to R^3
其中r(u,v)r(u,v)在△\triangle上连续,在点(u0,v0)∈△(u_0,v_0)\in\triangle处可微,过曲面SS上任意一点r0=r(u0,v0)r_0=r(u_0,v_0)的任意光滑曲线Γ\Gamma的参数方程为
r(t)=r(u(t),v(t)),t∈Ir(t)=r(u(t),v(t)),t\in I
其中u(t0)=u0,v(t0)=v0u(t_0)=u_0,v(t_0)=v_0在上式中对tt求导得
r′(t0)=ru(u0,v0)u′(t0)+rv(u0,v0)v′(t0)r'(t_0)=r_{u}(u_0,v_0)u'(t_0)+r_{v}(u_0,v_0)v'(t_0)
右端是切向量ru(u0,v0)r_{u}(u_0,v_0)与rv(u0,v0)r_{v}(u_0,v_0)的线性组合
若ru×rv∣u0,v0≠0r_u\times r_v|_{u_0,v_0}\neq0,说明ru(u0,v0)r_u(u_0,v_0)与rv(u0,v0)r_v(u_0,v_0)不共线
从而可张成一个平面Π\Pi,称之为曲面SS在正则点r0r_0的切平面
他的法向量为

ru×rv=∣ijkxuyzzuxvyvzv∣=(A,B,C)\begin{align} r_u\times r_v= \begin{vmatrix} i&j&k\\ x_u&y_z&z_u\\ x_v&y_v&z_v\\ \end{vmatrix} =(A,B,C) \end{align}

其中
A=∂(y,z)∂(u,v)∣(u0,v0),B=∂(y,z)∂(uz,x)∣(u0,v0),C=∂(x,y)∂(u,v)∣(u0,v0)A=\frac{\partial(y,z)}{\partial(u,v)}|_{(u_0,v_0)},B=\frac{\partial(y,z)}{\partial(uz,x)}|_{(u_0,v_0)},C=\frac{\partial(x,y)}{\partial(u,v)}|_{(u_0,v_0)}
法平面方程为
A(x−x0)+B(y−y0)+C(z−z0)=0A(x-x_0)+B(y-y_0)+C(z-z_0)=0
法线方程为
x−x0A=y−y0B=z−z0C\frac{x-x_0}{A}=\frac{y-y_0}{B}=\frac{z-z_0}{C}
若ru−(xu,yu,xu),rv=(xv,yv,zv)r_u-(x_u,y_u,x_u),r_v=(x_v,y_v,z_v)在△\triangle内连续,则称曲线SS为光滑曲线
对曲面F(x,y,z)=0F(x,y,z)=0,设在点P0(x0,y0,z0)P_0(x_0,y_0,z_0)处满足一函数定理中的条件,Fz(P0)≠0F_z(P_0)\neq0
隐函数z=z(z,y)z=z(z,y)
曲面S:r(x,y)=(x,y,z(x,y))S:r(x,y)=(x,y,z(x,y))
rx=(1,0,zx)=(1,0,−FxFz),ry=(0,1,zy)=(0,1,−FyFz)r_x=(1,0,z_x)=(1,0,-\frac{F_x}{F_z}),r_y=(0,1,z_y)=(0,1,-\frac{F_y}{F_z})

rx×ry=∣ijk10−FxFz01−FyFz∣=(FxFz,FyFz,1)\begin{align} r_x\times r_y= \begin{vmatrix} i&j&k\\ 1&0&-\frac{F_x}{F_z}\\ 0&1&-\frac{F_y}{F_z}\\ \end{vmatrix} =(\frac{F_x}{F_z},\frac{F_y}{F_z},1) \end{align}

法向量n⃗=(Fx,Fy,Fz)=∇F\vec{n}=(F_x,F_y,F_z)=\nabla F
切平面的方程
Fx(P0)(x−x0)+Fy(P0)(y−y0)+Fz(P0)(z−z0)=0F_{x}(P_0)(x-x_0)+F_{y}(P_0)(y-y_0)+F_{z}(P_0)(z-z_0)=0
法线方程
x−x0Fx(P0)=y−y0Fy(P0)=z−z0Fz(P0)\frac{x-x_0}{F_{x}(P_0)}=\frac{y-y_0}{F_{y}(P_0)}=\frac{z-z_0}{F_{z}(P_0)}

例6.8#

求曲面x=ucos⁡v,y=usin⁡v,z=avx=u\cos v,y=u\sin v,z=av在u=2,v=π4u=\sqrt{2},v=\frac{\pi}{4}处的切平面与法线方程
解:
r=(ucos⁡v,usin⁡v,av)r=(u\cos v,u\sin v,av)
ru=(cos⁡v,sin⁡v,0)→ru=(12,12,0)r_u=(\cos v,\sin v,0)\to r_u=(\frac{1}{\sqrt{2}},\frac{1}{\sqrt{2}},0)
rv=(−usin⁡v,ucos⁡v,a)→rv=(−1,1,a)r_v=(-u\sin v,u\cos v,a)\to r_v=(-1,1,a)

ru×rv∣(2,π4)=∣ijk12120−11a∣=12(a,−a,2)\begin{align} r_u\times r_{v}|_{(\sqrt{2},\frac{\pi}{4})}= \begin{vmatrix} i&j&k\\ \frac{1}{\sqrt{2}}&\frac{1}{\sqrt{2}}&0\\ -1&1&a \end{vmatrix} =\frac{1}{\sqrt{2}}(a,-a,2) \end{align}

取法向量n⃗=(a,−a,2)\vec{n}=(a,-a,2)
切平面
a(x−1)−a(y−1)+2(z−π4a)=0a(x-1)-a(y-1)+2(z-\frac{\pi}{4}a)=0
ax−ay+2z=π2aax-ay+2z=\frac{\pi}{2}a

Ex.1 环面的参数方程:#

考虑yOzyOz平面上的圆周C:(y−a)2+z2=b2C:(y-a)^2+z^2=b^2,其中a>b>0a>b>0
改曲线绕zz轴一周形成环面.求曲面方程