漫步微积分十七——最大最小值问题（续）

我们用其他的例子继续讨论上一篇文章的基本方法。

例1：圆柱形汤罐头的制造商接了一笔大订单，订单要求罐头的体积为V0。哪种尺寸可以最小化罐头的表面积，也就是所需的金属最少？

解：r,h分别表示圆柱底的半径和高(图1，左)，那么体积为

V0=πr2h(1)

总的表面积为

A=2πr2+2πrh(2)

我们必须最小化A，它是两个变量的函数，利用等式(1)得到h=V0/πr2，代入A得

A=2πr2+2πr⋅V0πr2=2πr2+2V0r(3)

图1

函数图像(图1，右)表明，r小或大时，A都大，那么在中间的某个位置存在最小值。跟之前一样，为了找到最小值的精确位置，我们对(3)求导，并让它等于零得

dAdr=4πr−2V0r2,4πr−2V0r2=0,4πr3=2V0,

2πr3=V0(4)

如果需要更多参数的解，可以根据等式(4)解出r，然后用它来计算h

r=V02π−−−√3,h=V0πr2=V0π(2πV0)2/3=2V02π−−−√3

从中我们可以看到h=2r。如果主要对形状感兴趣，那么将V0=πr2h代入(4)得

2πr3=πr2hor2r=h.

从降低原材料的角度考虑(最制造商非常重要)，这个结果告诉我们圆柱形罐头最好的形状是高等于底的直径。

例2：找出使圆柱体体积最大时的高与直径比。圆柱可以内接于半径为R的球内。

解：画一个外接球，如图2左，那么我们看到

V=2πx2y(5)

其中

x2+y2=R2(6)

图2

可视化极端的情况(图2，右)，当x趋近0,R时，V都很小，所以在极端情况之间存在最大体积的尺寸。为了找到它，将(6)代入(5)的

V=2πy(R2−y2)=2π(R2y−y3)

求导得

dVdy=2π(R2−3y2)

令其等于零得出y，然后利用(6)得出x

y=R3√x=R2−13R2−−−−−−−−−√=2√3√R

最大圆柱的高与直径比为

2y2x=yx=12√=frac122√≅0.707

这个结果利用隐函数求导会效率更高。x作为独立变量，y 作为函数的因变量，那么根据(6)得到

2x+2ydydx=0ordydx=−xy

根据(5)可得

dVdx==2π(x2dydx+2xy)=2π[x2(−xy)+2xy]2π(−x3+2xy2y)=2πxy(2y2−x2)

令上式等于零得到

2y2=x2oryx=12√=122√

例3：一束光从点A入射到平面镜上的一点P，反射到达点B，如图3，准确测量得出入射光线和出射光线与镜子的夹角相当α=β。假设光线从A经过镜子的反射到B的路径是最短的，证明当α=β时，路径就是最短的。

图3

解：将点P看做镜面上可变化的点，每个位置用x 表示，那么我们希望路径L可以表示成x的函数。从图中可以明显得出

L==a2+x2−−−−−−√+b2+(c−x)2−−−−−−−−−−√(a2+x2)1/2+[b2+(c−x)2]1/2

求导得

dLdx=12(a2+x2)−1/2⋅(2x)+12[b2+(c−x)2]−1/2⋅2(c−x)⋅(−1)=xa2+x2−−−−−−√−c−xb2+(c−x)2−−−−−−−−−−√(7)

为了最小化L，令上式等于零得

xa2+x2−−−−−−√=c−xb2+(c−x)2−−−−−−−−−−√(8)

变换一下等式的形式

a2+x2−−−−−−√x=b2+(c−x)2−−−−−−−−−−√c−x,(ax)2+1−−−−−−−−−√=(bc−x)2+1−−−−−−−−−−−√ax=bc−x

最后的等式可以轻松的解出x。然而，没必要这么做，因为最后比值的形式已经告诉了我们想要的：对于两个直角三角形中的角α,β，对边比邻边相等，所以两角相等。

直观上看，很明显，我们可以最小化L。如果想用二阶导来证实，利用(7)式计算二阶导

d2Ldx2=a2(a2+x2)3/2+b2[b2+(c−x)2]3/2

(我们跳过了具体计算细节)，注意到这个值是正的。

注解1：回顾一下一个锐角A余弦的定义。如果我们将A看作直角三角形(图4)锐角中的一个，那么根据定理

cosA=bc=adjacent sidehypotenuse

利用此定义，最小化(8)可以重新写为

cosα=cosβ

所以α=β。在回顾一下正弦的定义

sinA=ac=opposite sidehypotenuse

图4}

注解2：例3讨论的反射定律在古希腊时期就存在了。然而，反射光线遵从最短路径这个事实却是很晚之后才被发现，发现者是公元一世纪亚历山大时期的Heron，他的几何证明简单而巧妙。论述如下：A,B两点的位置跟之前一样(如图5)，B′是B的镜像，这样的话镜面就是BB′的垂直平分线。AB′与镜面交于点P，也就是光线入射的位置，因为α=γ,γ=β，所以α=β。路径的长度为AP+PB=AP+PB′=AB′。对于镜面上的其他位置P′，总长度为AP′+P′B=AP′+P′B′，它比AB′大，因为三角形任意两边之和大于第三边。这说明光的反射路径是最短的。


图5

例4：前面讨论的反射光线都是以恒定的速度在单一介质中传播。然而，在不同的介质(水，空气，玻璃)中，光的传播速度不一样。如果光线是从空气进入到水中，如图6，它朝与水面垂直的方向折射(弯曲)。路径APB明显不是最短路径。那么它由什么定律确定呢？在1621 年，荷兰科学家Snell发现光线的路径满足

sinαsinβ=a constant(9)

常数与A,B的位置无关。这个事实叫做Snell折射定律。可以通过最小化光线的传播时间来证明它。


图6

解：如果空气中光的速度为va，水中速度为vw，那么总时间T为

T==a2+x2−−−−−−√va+b2+(c−x)2−−−−−−−−−−√vw1va(a2+x2)1/2+1vw[b2+(c−x)2]1/2

如果计算这个函数的导数，那么

dTdx=1vaxa2+x2−−−−−−√−1vwc−xb2+(c−x)2−−−−−−−−−−√=sinαva−sinβvw0(1)

现在最小化T，结果为

sinαva=sinβvworsinαsinβ=vavw1(1)

这是Snell定律的显示表达，因为它告诉我们等式(9)右边常数的物理含义：光在空气中的速度与水中速度之比。这个常数叫做折射率。如果水用其他介质代替，例如酒精、甘油或玻璃，那么会得到不同的数值。

在例3中，我们可以计算二阶导为正值来证实(11)存在最小值T：

d2Tdx2=1vaa2(x2+x2)3/2+1vwb2[b2+(c−x)2]3/2>0

但是有必要提一下另一种方法。首先观察可得(10)式给出的dT/dx有两项。当x从0增加到c时，第一项(sinα/va)从0增加到某个正值。第二项(sinβ/vw)从某个正值减小到0。这表明dT/dx在x=0处为负值，增加到x=c的一个正值。所以T的最小值只存在一个，精确的值可由(11)得到。

注解3：例4的想法是费马在1657年发现的，为此光在光学系统中总是沿最短时间路径传播这个命题叫做费马最短时间原则。(应该注意到，当光在同一均匀介质中传播时，最短路径等价于最短时间)在他之后的两个世纪，费马的想法催生出了许多理论，最大最小原理，欧拉微积分变分法的创立，哈密尔顿的最小行动原则，它是物理学中最统一(deepest unifying)的原则之一。欧拉有句令人难忘的话：

Since the fabric of the world is the most perfect and was established by the wisest Creator,nothing happens in this world in which some reason of maximum or minimum would not come to light.

注解4：Snell的正弦定律(9)由笛卡尔在1637年发表(其中没有提及Snell)，他本意是想用等式

sinαsinβ=vwva

证明它。笛卡尔将他的论点建立在奇特的模型上，根据他的观点，光在高密度的介质中传播速度更快。费马拒绝他的观点(违背常识)。多年来笛卡尔只得被动接受费马的看法但是一直持怀疑态度，直到1657年他证明他的结论是对的，并且伴随着创造出了微积分方法。