Skip to content
CalcGospel 微積分福音
返回

113 台大微積分(B) 第 12 題

考題 / 轉學考微積分 / 台大 / 微積分(B)

113學年度 · 113微積分(B) · 第 12 題

題目

Problem

2. (a) (10%) Find the maximum value of f(x,y,z)=zf(x, y, z) = z on the curve of the intersection of x+y+z=1x + y + z = 1 and x2+y2+z2=3x^2 + y^2 + z^2 = 3.

(b) (5%) f(x,y,z)f(x, y, z), g(x,y,z)g(x, y, z), h(x,y,z)h(x, y, z) are differentiable functions. Assume that ff obtains a local maximum value at (x0,y0,z0)(x_0, y_0, z_0) when restricted to g(x,y,z)=cg(x, y, z) = c and h(x,y,z)=kh(x, y, z) = k. It is known that f(x0,y0,z0)=λg(x0,y0,z0)+μh(x0,y0,z0)\nabla f(x_0, y_0, z_0) = \lambda \nabla g(x_0, y_0, z_0) + \mu \nabla h(x_0, y_0, z_0) for some constants λ\lambda and μ\mu. Suppose that f(x,y,z)f(x, y, z) obtains new local maximum at (x1,y1,z1)(x_1, y_1, z_1) when restricted to g(x,y,z)=cg(x, y, z) = c and h(x,y,z)=k+ϵh(x, y, z) = k + \epsilon where ϵ|\epsilon| is small and (x1,y1,z1)(x_1, y_1, z_1) is close to (x0,y0,z0)(x_0, y_0, z_0). Show by the linear approximation that we can approximate f(x1,y1,z1)f(x0,y0,z0)f(x_1, y_1, z_1) - f(x_0, y_0, z_0) by μϵ\mu \cdot \epsilon.

(c) (5%) Estimate, by linear approximation, the maximum value of f(x,y,z)=zf(x, y, z) = z on the curve of the intersection of x+y+z=1x + y + z = 1 and x2+y2+z2=3.02x^2 + y^2 + z^2 = 3.02.

解答

解法一

思路

展開

本題是一道結合多元函數拉格朗日乘子法與全微分微擾線性近似的綜合分析題。

(a) 求 f(x,y,z)=zf(x,y,z)=z 在兩約束交線上的最大值

  • 目標函數為 f(x,y,z)=zf(x,y,z) = z,約束條件為 g(x,y,z)=x+y+z=1g(x,y,z) = x+y+z = 1,以及 h(x,y,z)=x2+y2+z2=3h(x,y,z) = x^2+y^2+z^2 = 3
  • 根據拉格朗日乘子法,建立方程式: f=λg+μh    0,0,1=λ1,1,1+μ2x,2y,2z\nabla f = \lambda \nabla g + \mu \nabla h \implies \langle 0, 0, 1 \rangle = \lambda \langle 1, 1, 1 \rangle + \mu \langle 2x, 2y, 2z \rangle
  • 解此方程組可得 x=yx = y。代回約束條件求得兩組候選解,比較 zz 分量大小,最大者即為最大值。

(b) 證明線性近似關係式 Δfμϵ\Delta f \approx \mu \cdot \epsilon

  • 根據多元函數的線性近似(一階泰勒展開式): f(x1,y1,z1)f(x0,y0,z0)f(x0,y0,z0)Δx,Δy,Δzf(x_1, y_1, z_1) - f(x_0, y_0, z_0) \approx \nabla f(x_0, y_0, z_0) \cdot \langle \Delta x, \Delta y, \Delta z \rangle
  • 代入拉格朗日乘子關係式 f=λg+μh\nabla f = \lambda \nabla g + \mu \nabla h,並利用兩個約束條件的微擾差值進行化簡。

(c) 估計在微擾約束 x2+y2+z2=3.02x^2+y^2+z^2 = 3.02 下的新最大值

  • 對比 (a) 與 (b) 小題,此時常數 c=1c = 1 未變(Δc=0\Delta c = 0),而常數 kk33.023 \to 3.02(微擾 ϵ=0.02\epsilon = 0.02)。
  • 我們首先計算出 (a) 小題極值點處的乘子 μ\mu 值。
  • 代入 (b) 小題證得的公式 Δfμϵ\Delta f \approx \mu \epsilon 進行線性估算。

答題過程

展開

(a) 計算 f(x,y,z)=zf(x, y, z) = z 的最大值

我們設定約束方程式為:

g(x,y,z)=x+y+z=1g(x, y, z) = x + y + z = 1 h(x,y,z)=x2+y2+z2=3h(x, y, z) = x^2 + y^2 + z^2 = 3

建立拉格朗日乘子關係式 f=λg+μh\nabla f = \lambda \nabla g + \mu \nabla h

0,0,1=λ1,1,1+μ2x,2y,2z\langle 0,\, 0,\, 1 \rangle = \lambda \langle 1,\, 1,\, 1 \rangle + \mu \langle 2x,\, 2y,\, 2z \rangle

這給出分量方程組:

{0=λ+2μx— (1)0=λ+2μy— (2)1=λ+2μz— (3)\begin{align*} \begin{cases} 0 = \lambda + 2\mu x & \text{--- (1)} \\[2mm] 0 = \lambda + 2\mu y & \text{--- (2)} \\[2mm] 1 = \lambda + 2\mu z & \text{--- (3)} \end{cases} \end{align*}

由式 (1) 與 (2) 可得 λ=2μx=2μy\lambda = -2\mu x = -2\mu y。 若 μ=0\mu = 0,則會導致 λ=0\lambda = 0,這與式 (3) 矛盾(得 1=01 = 0)。因此 μ0\mu \neq 0,這直接說明:

x=y— (4)x = y \quad \text{--- (4)}

我們將式 (4) 代回兩個約束條件中:

  • g=1g = 1 得: 2y+z=1    z=12y2y + z = 1 \implies z = 1 - 2y
  • h=3h = 3 得: 2y2+z2=32y^2 + z^2 = 3

zz 代入:

2y2+(12y)2=3    2y2+14y+4y2=3    6y24y2=0    3y22y1=0    (3y+1)(y1)=0\begin{align*} 2y^2 + (1 - 2y)^2 = 3 \implies&\, 2y^2 + 1 - 4y + 4y^2 = 3 \\[2mm] \implies&\, 6y^2 - 4y - 2 = 0 \\[2mm] \implies&\, 3y^2 - 2y - 1 = 0 \implies (3y + 1)(y - 1) = 0 \end{align*}

解得:

y=13y=1y = -\frac{1}{3} \quad \text{或} \quad y = 1

我們對應求出兩個候選點:

  1. P1(1,1,1)P_1(1, 1, -1): 此時, z=1z = -1
  2. P2(13,13,53)P_2\left(-\frac{1}{3}, -\frac{1}{3}, \frac{5}{3}\right): 此時, z=53z = \frac{5}{3}

比較兩點的 zz 值大小,最大值為:

zmax=53(發生在 (13,13,53) 處)z_{\max} = \frac{5}{3} \quad \text{(發生在 } \left(-\frac{1}{3}, -\frac{1}{3}, \frac{5}{3}\right) \text{ 處)}

(b) 證明線性近似公式

Δx=x1x0,Δy=y1y0,Δz=z1z0\Delta x = x_1 - x_0, \Delta y = y_1 - y_0, \Delta z = z_1 - z_0。 根據多元函數的一階線性近似(全微分關係):

f(x1,y1,z1)f(x0,y0,z0)f(x0,y0,z0)Δx,Δy,Δz=(λg(x0,y0,z0)+μh(x0,y0,z0))Δx,Δy,Δz=λ(g(x0,y0,z0)Δx,Δy,Δz)+μ(h(x0,y0,z0)Δx,Δy,Δz)— (5)\begin{align*} f(x_1, y_1, z_1) - f(x_0, y_0, z_0) \approx&\, \nabla f(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \\[2mm] =&\, \Big( \lambda \nabla g(x_0, y_0, z_0) + \mu \nabla h(x_0, y_0, z_0) \Big) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \\[2mm] =&\, \lambda \Big( \nabla g(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \Big) \\ &\, + \mu \Big( \nabla h(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \Big) \quad \text{--- (5)} \end{align*}

另一方面,對約束函數 gghh(x0,y0,z0)(x_0, y_0, z_0) 處作一階線性近似:

  • 對於 ggg(x1,y1,z1)g(x0,y0,z0)g(x0,y0,z0)Δx,Δy,Δzg(x_1, y_1, z_1) - g(x_0, y_0, z_0) \approx \nabla g(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle 由於 g(x1,y1,z1)=cg(x_1, y_1, z_1) = cg(x0,y0,z0)=cg(x_0, y_0, z_0) = c,故: g(x0,y0,z0)Δx,Δy,Δzcc=0— (6)\nabla g(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \approx c - c = 0 \quad \text{--- (6)}
  • 對於 hhh(x1,y1,z1)h(x0,y0,z0)h(x0,y0,z0)Δx,Δy,Δzh(x_1, y_1, z_1) - h(x_0, y_0, z_0) \approx \nabla h(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle 由於 h(x1,y1,z1)=k+ϵh(x_1, y_1, z_1) = k + \epsilonh(x0,y0,z0)=kh(x_0, y_0, z_0) = k,故: h(x0,y0,z0)Δx,Δy,Δz(k+ϵ)k=ϵ— (7)\nabla h(x_0, y_0, z_0) \cdot \langle \Delta x,\, \Delta y,\, \Delta z \rangle \approx (k + \epsilon) - k = \epsilon \quad \text{--- (7)}

將式 (6) 與 (7) 代回式 (5) 中:

f(x1,y1,z1)f(x0,y0,z0)λ(0)+μ(ϵ)=μϵf(x_1, y_1, z_1) - f(x_0, y_0, z_0) \approx \lambda(0) + \mu(\epsilon) = \mu \cdot \epsilon

公式得證。


(c) 估計微擾後的最大值

此時,約束為 g=1g = 1h=3.02h = 3.02。這對應了 (b) 小題中的條件:

c=1,k=3,ϵ=3.023=0.02c = 1, \quad k = 3, \quad \epsilon = 3.02 - 3 = 0.02

我們需要求出 (a) 小題最大值點 P2(13,13,53)P_2\left(-\frac{1}{3}, -\frac{1}{3}, \frac{5}{3}\right) 處對應的拉格朗日乘子 μ\mu。 將 x0=13,z0=53x_0 = -\frac{1}{3}, z_0 = \frac{5}{3} 代回 (a) 小題的方程組中:

  1. 由式 (1) 得: λ=2μx0=23μ\lambda = -2\mu x_0 = \frac{2}{3}\mu
  2. 將其代入式 (3): 1=λ+2μz0    1=23μ+2μ(53)=123μ=4μ1 = \lambda + 2\mu z_0 \implies 1 = \frac{2}{3}\mu + 2\mu\left(\frac{5}{3}\right) = \frac{12}{3}\mu = 4\mu 由此解得乘子: μ=14\mu = \frac{1}{4}

根據 (b) 小題已證的線性近似公式:

Δfμϵ=14(0.02)=0.005\Delta f \approx \mu \cdot \epsilon = \frac{1}{4} \cdot (0.02) = 0.005

因此,新的最大值估計為:

fnewfold+Δf=53+0.005=53+1200=10036001.6717f_{\text{new}} \approx f_{\text{old}} + \Delta f = \frac{5}{3} + 0.005 = \frac{5}{3} + \frac{1}{200} = \frac{1003}{600} \approx 1.6717

結論:

  • (a) 最大值為 53\displaystyle \frac{5}{3}
  • (b) 證明如上。
  • (c) 新最大值估計為 1003600\displaystyle \frac{1003}{600}(約 1.67171.6717)。