跳到主要内容

第八章:积分的进一步应用

《斯图尔特微积分》第八章:积分的进一步应用读书笔记。本章继续扩展积分的应用范围,从几何中的弧长和曲面面积,延伸到流体压力、质心、经济学、生物学与概率密度。

本章地图

小节核心问题需要掌握的内容
8.1 弧长的计算曲线长度怎样由折线极限得到?弧长公式、弧长函数
8.2 旋转曲面的面积旋转曲面面积怎样累积?表面积微元、绕坐标轴旋转
8.3 物理学和工程学中的应用积分怎样描述连续物理量?液体压力、矩、质心、帕普斯定理
8.4 经济学和生物学中的应用积分怎样进入跨学科模型?消费者剩余、血流量、心输出量
8.5 概率连续随机变量怎样分配概率?概率密度、均值、中位值、正态分布

8.1 弧长的计算

引文

让一根线和曲线重合,然后用尺子测量线的长度,那么这个长度也就是曲线的长度,最初我们是这么设想的,但是当曲线变得复杂时,这种方法就不那么奏效了,此时,我们需要一种可以求解曲线长度的方法。

曲线的弧长

如果曲线是一条折线,我们很容易就可以求出它的长度,仅仅是几条线段相加就可以了。可以使用类似的方法,要定义一般曲线的长度,首先将其近似为一个折线轨迹,可以回想一下,曾经我们用圆内接多边形的边长之和的极限来求圆的周长。

原讲义插图待补充

原始笔记引用了 8-1.1.png,但当前目录中没有找到对应图片。正文与公式已保留,后续补入图片后可替换此提示。

小结 L=limn+i=1nPi1PiL=\lim_{ n \to +\infty }\sum_{i=1}^n|P_{i-1}-P_{i}| 其中 Pi1Pi|P_{i-1}-P_{i}|Pi1P_{i-1}PiP_{i} 之间的距离。也就是很多折线的长度之和。

我们可以用积分推导出 LL 的一个公式:

Δyi=yiyi1\Delta y_{i}=y_{i}-y_{i-1} ,那么 Pi1Pi=(xixi1)2+(yiyi1)2=(Δx)2+(Δyi)2|P_{i-1}P_{i}|=\sqrt{ (x_{i}-x_{i-1})^2+(y_{i}-y_{i-1})^2 }=\sqrt{ (\Delta x)^2+(\Delta y_{i})^2 }

公式

弧长公式:

如果 ff' 在区间 [a,b][a,b] 上连续,那么曲线 y=f(x),axby=f(x),a\leqslant x\leqslant b 的长度为

L=ab1+[f(x)]2dxL=\int_{a}^{b} \sqrt{ 1+[f'(x)]^2 } \, dx 也可以写为 L=ab1+(dydx)2dxL=\int_{a}^{b} \sqrt{ 1+\left( \frac{dy}{dx} \right)^2 } \, dx

例题 例1:求半立方抛物线 y2=x3y^2=x^3 在点 (1,1)(1,1) 和点 (4,8)(4,8) 之间的弧长 8-1.2 解: y=x32    y=32x12y=x^{\frac{3}{2}}\implies y'=\frac{3}{2}x^{\frac{1}{2}} L=141+(32x12)2dx=141+94xdxL=\int_{1}^{4} \sqrt{ 1+\left( \frac{3}{2}x^{\frac{1}{2}} \right)^2 } \, dx =\int_{1}^{4} \sqrt{ 1+\frac{9}{4}x } \, dx u=1+94xu=1+\dfrac{9}{4}x,则 du=94dxdu=\dfrac{9}{4}\,dx。当 x=1x=1 时,u=134u=\dfrac{13}{4};当 x=4x=4 时,u=10u=10L=4913410udu=49×23u3213410=8010131327L=\frac{4}{9}\int_{\frac{13}{4}}^{10}\sqrt{ u } \, du=\frac{4}{9}\times \frac{2}{3}u^{\frac{3}{2}}\bigg|_{\frac{13}{4}}^{10}=\frac{80\sqrt{ 10 }-13\sqrt{ 13 }}{27}

对于上述问题,如果曲线方程是 x=g(y),cydx=g(y),c\leqslant y\leqslant d ,我们可以把 xxyy 互换,得到如下变形公式:

公式

弧长公式变形

L=cd1+[g(y)]2dy=cd1+(dxdy)2dyL=\int_{c}^{d} \sqrt{ 1+[g'(y)]^2 } \, dy=\int_{c}^{d} \sqrt{ 1+\left( \frac{dx}{dy} \right)^2 } \, dy

例题 例2:求抛物线 y2=xy^2=x 在点 (0,0)(0,0) 和点 (1,1)(1,1) 之间的弧长。 解: 由 x=y2x=y^2 可得

dxdy=2y\frac{dx}{dy}=2y

L=011+(2y)2dy=011+4y2dyL=\int_{0}^{1} \sqrt{ 1+(2y)^2 } \, dy=\int_{0}^{1} \sqrt{ 1+4y^2 } \, dy 作三角换元

y=12tanθ,dy=12sec2θdθy=\frac{1}{2}\tan\theta, \qquad dy=\frac{1}{2}\sec^2\theta\,d\theta

1+4y2=1+tan2θ=secθ\sqrt{ 1+4y^2 }=\sqrt{ 1+\tan^2\theta }=\sec\thetay=0y=0 时,θ=0\theta=0;当 y=1y=1 时,tanθ=2\tan\theta=2,记 A=arctan2A=\arctan 2L=0Asecθ12sec2θdθ=120Asec3dθL=\int_{0}^{A} \sec\theta\cdot \frac{1}{2}\sec^2\theta \, d\theta=\frac{1}{2}\int_{0}^{A} \sec^3 \, d\theta L=14(secAtanA+lnsecA+tanA)L=\frac{1}{4}(\sec A\tan A+\ln |\sec A+\tan A|)

tanA=2,sec2A=1+tan2A=5\tan A=2, \qquad \sec^2 A=1+\tan^2 A=5

L=52+ln(5+2)4L= \frac{\sqrt{ 5 }}{2}+\frac{\ln(\sqrt{ 5 }+2)}{4}

弧长函数

8-1.3 如果用 s(x)s(x) 表示曲线 CC 从起点 P0(a,f(a))P_{0}(a,f(a)) 到终点 Q(x,f(x))Q(x,f(x)) 的距离,则 ss 为一个函数,称为弧长函数: s(x)=ax1+[f(t)]2dts(x)=\int_{a}^{x} \sqrt{ 1+[f'(t)]^2 } \, dt dsdx=1+[f(t)]=1+(dydx)2\frac{ds}{dx}=\sqrt{ 1+[f'(t)] }=\sqrt{ 1+\left( \frac{dy}{dx} \right)^2 } ds=1+(dydx)dxds=\sqrt{ 1+\left( \frac{dy}{dx} \right) }dx 也可以推导出: ds=1+(dxdy)dyds=\sqrt{ 1+\left( \frac{dx}{dy} \right) }dy

例题 例3:求曲线 y=x218lnxy=x^2-\frac{1}{8}\ln xP0(1,1)P_{0}(1,1) 为起点的弧长函数。 解: y=2x18xy'=2x-\frac{1}{8x} 1+[f(x)]2=2x+18x\sqrt{ 1+[f'(x)]^2 }=2x+\frac{1}{8x} s(x)=1x(2t+18t)dt=[t2+18lnt]1xs(x)=\int_{1}^{x} \left( 2t+\frac{1}{8t} \right) \, dt =\left[ t^2+\frac{1}{8}\ln t \right]_{1}^x s(x)=x2+18lnx1s(x)=x^2+\frac{1}{8\ln x}-1


8.2 旋转曲面的面积

小结 圆锥体侧面积:

A=πrlA=\pi rl 底面圆半径为 rr ,母线长为 ll

8-2.2 小结 圆台侧面积:

A=2πrlA=2\pi rl 其中 rr 为条带平均半径。

8-2.1 观察上图,它是由曲线 y=f(x),axby=f(x),a\leqslant x\leqslant bxx 轴旋转而成。 我们计算其中一个条带的面积: 母线长和平均半径分别为

l=Pi1Pi,r=12(yi1+yi)l=|P_{i-1}P_i|, \qquad r=\frac{1}{2}(y_{i-1}+y_i)

Pi1Pi=1+[f(x)]2Δx |P_{i-1}P_{i}|=\sqrt{ 1+[f'(x^*)]^2 }\Delta x r=f(xi)r=f(x_{i}^*) A=2πPi1Pi12(yi1+yi)A=2\pi |P_{i-1}P_{i}|\frac{1}{2}(y_{i-1}+y_{i}) 那么,整个旋转曲面的面积为: S=i=1n2πf(xi)1+f(xi)2ΔxS=\sum_{i=1}^n2\pi f(x_{i}^*)\sqrt{ 1+f'(x_{i}^*)^2 }\Delta x

小结 旋转曲面面积公式:

由曲线 y=f(x),axby=f(x),a\leqslant x\leqslant bxx 轴旋转而成, S=ab2πf(x)1+[f(x)]2dxS=\int_{a}^{b} 2\pi f(x)\sqrt{ 1+[f'(x)]^2 } \, dx S=ab2πy1+(dydx)2dxS=\int_{a}^{b} 2\pi y\sqrt{ 1+\left( \frac{dy}{dx} \right)^2 } \, dx xx 轴变形为 S=cd2πy1+(dxdy)dyS=\int_{c}^{d}2\pi y\sqrt{ 1+\left( \frac{dx}{dy} \right) } \, dy

例题 例1:曲线 y=4x2,1x1y=\sqrt{ 4-x^2 },-1\leqslant x\leqslant1 是一段弧,求这段弧绕 xx 轴旋转所得曲面的面积。 解: dydx=x4x2\frac{dy}{dx}=-\frac{x}{\sqrt{ 4-x^2 }} S=112πy1+(dydx)2dxS=\int_{-1}^{1} 2\pi y\sqrt{ 1+\left( \frac{dy}{dx} \right)^2 } \, dx S=2π114x21+x24x2dxS=2\pi\int_{-1}^{1} \sqrt{ 4-x^2 }\sqrt{ 1+ \frac{x^2}{4-x^2} } \, dx S=4π111dx=4π×2=8πS=4\pi \int_{-1}^{1} 1 \, dx=4\pi \times 2=8\pi

例2:求曲线 x=23y32x=\frac{2}{3}y^{\frac{3}{2}}y=0y=0y=3y=3 之间的弧绕 xx 轴旋转所得曲面的面积。 8-2.3 解: S=032πy1+(dxdy)dy=2π03y1+(y12)2dyS=\int_{0}^{3} 2\pi y\sqrt{ 1+\left( \frac{dx}{dy} \right) } \, dy=2\pi \int_{0}^{3} y\sqrt{ 1+\left( y^{\frac{1}{2}} \right)^2 } \, dy S=2π03y1+ydyS=2\pi \int_{0}^{3} y\sqrt{ 1+y } \, dy u=1+y,du=dyu=1+y,du=dy S=2π14(u1)udu=23215πS=2\pi \int_{1}^{4} (u-1)\sqrt{ u } \, du=\frac{232}{15}\pi


8.3 物理学和工程学中的应用

液体的压强和压力

引文 潜水员经验:越下潜,水的压力就越大,这是因为他们上方的水越来越重。

液体对平板的压力: F=mg=ρgAdF=mg=\rho gAd 单位面积上所受压力: P=FA=ρgdP=\frac{F}{A}=\rho gd 例题 例1:梯形大坝高20m,顶部宽50m,底部宽30m,如果水位到坝顶的距离为4m,求流体静压强造成的大坝所受的力。 解: 8-3.1 建立坐标系,令 xx 轴垂直,向下为正方向,原点在水平面处,水深16m,将区间 [0,16][0,16] 分成等宽子区间,端点为 xix_{i} ,样本点为 xi[xi1,xi]x_{i}^*\in [x_{i-1},x_{i}] ,那么,大坝中的第 ii 个水平条形可以用高为 Δx\Delta x 、宽为 wiw_{i} 的矩形来近似, a16xi=1020    a=8xi2\frac{a}{16-x_{i}^*}=\frac{10}{20}\implies a=8-\frac{x_{i}^*}{2} wi=2(15+a)=46xiw_{i}=2(15+a)=46-x_{i}^*ii 个条形的面积为 AiA_{i} AiwiΔx=(46xi)ΔxA_{i} \approx w_{i}\cdot\Delta x=(46-x_{i}^*)\Delta x 因为 Δx\Delta x 很小,那么第 ii 个条形上的压强 PiP_{i} 几乎恒定的 Pi=1000gxiP_{i}=1000gx_{i}^* Fi=PiAi1000gxi(46xi)ΔxF_{i}=P_{i}A_{i}\approx1000gx_{i}^*(46-x_{i}^*)\Delta x F=limn+i=1n1000gxi(46xi)ΔxF=\lim_{ n \to +\infty }\sum_{i=1}^n1000gx_{i}^*(46-x_{i}^*)\Delta x =0161000gx(46x)dx=9800[23x2x33]016=\int_{0}^{16} 1000gx(46-x) \, dx=9800\left[ 23x^2-\frac{x^3}{3} \right]_{0}^{16} 4.43×107N\approx 4.43\times 10^7N

物体的矩和质心

引文 在任意形状的薄板上找到一点 PP ,使得以这一点为支点,薄板可以维持水平平衡,这一点称为薄板的质心(重心)。

8-3.2

质心所在处: xˉ=m1x1+m2x2m1+m2\bar{x}=\frac{m_{1}x_{1}+m_{2}x_{2}}{m_{1}+m_{2}} xˉ=i=1nmixii=1nmi=i=1nmixim\bar{x}=\frac{\displaystyle\sum_{i=1}^nm_{i}x_{i}}{\displaystyle\sum_{i=1}^nm_{i}}=\frac{\displaystyle\sum_{i=1}^nm_{i}x_{i}}{m} 系统关于原点的矩: M=i=1nmixiM=\sum_{i=1}^nm_{i}x_{i} 系统关于 yy 轴的矩: My=i=1nmixiM_{y}=\sum_{i=1}^nm_{i}x_{i} My=limn+i=1nρxˉif(xˉi)Δx=ρabxf(x)dxM_{y}=\lim_{ n \to +\infty }\sum_{i=1}^n\rho \bar{x}_{i}f(\bar{x}_{i})\Delta x=\rho \int_{a}^{b} xf(x) \, dx 系统关于 xx 轴的矩: Mx=i=1nmiyiM_{x}=\sum_{i=1}^nm_{i}y_{i} Mx=limn+i=1nρ12[f(xˉi)]2Δx=ρab12[f(x)]2dxM_{x}=\lim_{ n \to +\infty }\sum_{i=1}^n\rho \frac{1}{2}[f(\bar{x}_{i})]^2\Delta x=\rho \int_{a}^{b} \frac{1}{2}[f(x)]^2 \, dx

质心位于点 (xˉ,yˉ)(\bar{x},\bar{y}) 处: xˉ=1Aabxf(x)dx,yˉ=1Aab12[f(x)]2dx\bar{x}=\frac{1}{A}\int_{a}^{b} xf(x) \, dx ,\quad \bar{y}=\frac{1}{A}\int_{a}^{b} \frac{1}{2}[f(x)]^2 \, dx

例题 例:求密度均匀的半径为 rr 的半圆形平板的质心。 解: 8-3.3 f(x)=r2x2,a=r,b=rf(x)=\sqrt{ r^2-x^2 },a=-r,b=r 由对称性可知 xˉ=0\bar{x}=0A=12πr2A=\frac{1}{2}\pi r^2 yˉ=1Arr[f(x)]2dx=112πr212rr(r2x2)2dx\bar{y}=\frac{1}{A}\int_{-r}^{r} [f(x)]^2 \, dx =\dfrac{1}{\dfrac{1}{2}\pi r^2} \frac{1}{2}\int_{-r}^{r}(\sqrt{ r^2-x^2 })^2 \, dx =2πr20r(r2x2)dx=4r3π=\frac{2}{\pi r^2}\int_{0}^{r} (r^2-x^2) \, dx=\frac{4r}{3\pi} 质心为

(0,4r3π)\left(0,\frac{4r}{3\pi}\right)

如果区域 RR 位于曲线 y=f(x)y=f(x)y=g(x)y=g(x) 之间,其中 f(x)g(x)f(x)\geqslant g(x) ,那么: 8-3.4

公式

形心

xˉ=1Aabx[f(x)g(x)]dx\bar{x}=\frac{1}{A}\int_{a}^{b} x[f(x)-g(x)] \, dx yˉ=12Aab([f(x)]2[g(x)]2)dx\bar{y}=\frac{1}{2A}\int_{a}^{b} ([f(x)]^2-[g(x)]^2) \, dx

帕普斯定理

RR 为一个平面区域,它完全位于平面上直线 ll 的一侧,那么 RRll 旋转所得立体的体积是 RR 的面积 AARR 的形心所经过 的距离 dd 的乘积。

V=A×dV=A\times d


8.4 经济学和生物学中的应用

消费者剩余

消费者剩余 对于给定的价格,一些购买商品的消费者本来愿意支付更多的钱,但他们因不必这样做而受益,消费者愿意支付的钱与消费者实际支付的钱的差额称为消费者剩余

通过求出某商品所有购买者的总消费者剩余,经济学家可以评估市场对社会的总体效益。

在每个商品节约的金额 ×\times 商品数量= [p(xi)P]Δx[p(x_{i})-P]\Delta x

商品的总消费者剩余: 0X[p(x)P]dx\int_{0}^{X} [p(x)-P] \, dx 几何解释:需求曲线下方和直线 p=Pp=P 上方区域的面积。

例题 例:某商品的需求函数 p=12000.2x0.0001x2p=1200-0.2x-0.0001x^2 求销量为 500500 时的消费者剩余。

解: px=500=12000.2×5000.0001×5002=1075p|_{x=500}=1200-0.2\times 500-0.0001\times 500^2=1075 0500(12000.2x0.0001x21075)dx\int_{0}^{500} (1200-0.2x-0.0001x^2-1075) \, dx 33333.33\approx33333.33

血流量

公式 层流定律: v(r)=P4ηl(R2r2)v(r)=\frac{P}{4\eta l}(R^2-r^2) 泊肃叶定律: F=πPR48ηlF=\frac{\pi PR^4}{8\eta l}

心输出量

经过测量点的燃料量近似为

浓度×体积=c(ti)(FΔt)\text{浓度}\times\text{体积} =c(t_i)(F\Delta t)

公式 染料的总量 A=F0Tc(t)dtA=F\int_{0}^{T} c(t) \, dt 心输出量 F=A0Tc(t)dtF=\frac{A}{\displaystyle\int_{0}^{T}c(t) \, dt }


8.5 概率

概率密度函数

每个连续随机变量 XX 都有一个概率密度函数 ff ,这表示 XXaabb 之间的概率为 ffaabb 的积分: P(aXb)=abf(x)dxP(a\leqslant X\leqslant b) =\int_{a}^{b} f(x) \, dx

均值和中位值

一般来说,我们将任何概率密度函数 ff 的均值定义为 μ=+xf(x)dx\mu=\int_{-\infty}^{+\infty} xf(x) \, dx

正态分布

正态分布是指随机变量 XX 的概率密度函数是一下函数族中的一个成员函数: f(x)=1σ2πe(xμ)22σ2f(x)=\frac{1}{\sigma \sqrt{ 2\pi }}e^{\frac{-(x-\mu)^2}{2\sigma^2}}


本章总结

  • 弧长来自折线长度之和的极限,其微元包含水平与竖直变化。
  • 旋转曲面面积由周长与弧长微元的乘积累积得到。
  • 压力、质心等问题需要先建立正确的连续分布模型。
  • 积分可以表示消费者剩余、血流量与心输出量等跨学科量。
  • 概率密度曲线下的面积给出连续随机变量落在区间内的概率。