Skip to content
CalcGospel 微積分福音
返回

112 政治大學微積分 第 7 題

考題 / 轉學考微積分 / 政大 / 微積分

112學年度 · 112微積分 · 第 7 題

題目

Problem

3. Let S(a,b)=∑i(yi−a−bxi)2S(a, b) = \sum_i (y_i - a - b x_i)^2, where ∑ixi2=n, ∑ixi=n/2\sum_i x_i^2 = n, \ \sum_i x_i = n/2. (1) Let (a1,b1)(a_1, b_1) be the minimizer of S(a,b)S(a, b). Find (a1,b1)(a_1, b_1). (6%) (2) Under the constraint a+b=ca + b = c, find the minimizer of S(a,b)S(a, b) in terms of (a1,b1)(a_1, b_1) and cc. (6%) Ps., yiy_i is a known number for i=1,2,…,ni = 1, 2, \dots, n.

解答

(1) 求解無約束最小化最優解 (a1,b1)(a_1, b_1)

我們要對偏偏導數設為零:

∂S∂a=−2∑i=1n(yi−a−bxi)=0  ⟹  ∑i=1nyi=na+b∑i=1nxi— (1)\frac{\partial S}{\partial a} = -2\sum_{i=1}^n (y_i - a - b x_i) = 0 \implies \sum_{i=1}^n y_i = na + b\sum_{i=1}^n x_i \quad \text{--- (1)} ∂S∂b=−2∑i=1nxi(yi−a−bxi)=0  ⟹  ∑i=1nxiyi=a∑i=1nxi+b∑i=1nxi2— (2)\frac{\partial S}{\partial b} = -2\sum_{i=1}^n x_i(y_i - a - b x_i) = 0 \implies \sum_{i=1}^n x_i y_i = a\sum_{i=1}^n x_i + b\sum_{i=1}^n x_i^2 \quad \text{--- (2)}

將已知條件 ∑i=1nxi=n2\sum_{i=1}^n x_i = \frac{n}{2} 與 ∑i=1nxi2=n\sum_{i=1}^n x_i^2 = n 代回方程組:

{na+n2b=∑yin2a+nb=∑xiyi\begin{cases} n a + \frac{n}{2} b = \sum y_i \\[2mm] \frac{n}{2} a + n b = \sum x_i y_i \end{cases}

我們將兩式同除以 nn 簡化:

{a+12b=yˉ(其中 yˉ=1n∑yi)12a+b=1n∑xiyi\begin{cases} a + \frac{1}{2} b = \bar{y} \quad (\text{其中 } \bar{y} = \frac{1}{n}\sum y_i) \\[2mm] \frac{1}{2} a + b = \frac{1}{n}\sum x_i y_i \end{cases}

我們使用消元法解此線性方程組:

  • 第一式乘以 2 減去第二式: 2a+b−(12a+b)=2yˉ−1n∑xiyi  ⟹  32a=2yˉ−1n∑xiyi2a + b - \left(\frac{1}{2}a + b\right) = 2\bar{y} - \frac{1}{n}\sum x_i y_i \implies \frac{3}{2}a = 2\bar{y} - \frac{1}{n}\sum x_i y_i a1=4yˉ3−23n∑xiyi=4∑yi−2∑xiyi3na_1 = \frac{4\bar{y}}{3} - \frac{2}{3n}\sum x_i y_i = \frac{4\sum y_i - 2\sum x_i y_i}{3n}
  • 第二式乘以 2 減去第一式: a+2b−(a+12b)=2n∑xiyi−yˉ  ⟹  32b=2n∑xiyi−yˉa + 2b - \left(a + \frac{1}{2}b\right) = \frac{2}{n}\sum x_i y_i - \bar{y} \implies \frac{3}{2}b = \frac{2}{n}\sum x_i y_i - \bar{y} b1=4∑xiyi−2∑yi3nb_1 = \frac{4\sum x_i y_i - 2\sum y_i}{3n}

(2) 求解在約束 a+b=ca + b = c 下的最優解

我們可以使用拉格朗日乘子法 (Lagrange Multipliers) 來求解此有約束的最優化問題。 定義拉格朗日函數:

L(a,b,μ)=∑i=1n(yi−a−bxi)2+μ(a+b−c)L(a, b, \mu) = \sum_{i=1}^n (y_i - a - b x_i)^2 + \mu (a + b - c)

對 aa 與 bb 求偏導數並令其為零:

∂L∂a=−2∑i=1n(yi−a−bxi)+μ=0  ⟹  2(∑yi−na−b∑xi)=μ— (3)\frac{\partial L}{\partial a} = -2\sum_{i=1}^n (y_i - a - b x_i) + \mu = 0 \implies 2\left( \sum y_i - na - b\sum x_i \right) = \mu \quad \text{--- (3)} ∂L∂b=−2∑i=1nxi(yi−a−bxi)+μ=0  ⟹  2(∑xiyi−a∑xi−b∑xi2)=μ— (4)\frac{\partial L}{\partial b} = -2\sum_{i=1}^n x_i(y_i - a - b x_i) + \mu = 0 \implies 2\left( \sum x_i y_i - a\sum x_i - b\sum x_i^2 \right) = \mu \quad \text{--- (4)}

將式 (3) 與 (4) 聯立消去乘子 μ\mu:

∑yi−na−b∑xi=∑xiyi−a∑xi−b∑xi2\sum y_i - na - b\sum x_i = \sum x_i y_i - a\sum x_i - b\sum x_i^2

我們代入已知條件 ∑xi=n2\sum x_i = \frac{n}{2} 與 ∑xi2=n\sum x_i^2 = n:

∑yi−na−n2b=∑xiyi−n2a−nb\sum y_i - na - \frac{n}{2}b = \sum x_i y_i - \frac{n}{2}a - nb ∑yi−∑xiyi=n2a−n2b=n2(a−b)\sum y_i - \sum x_i y_i = \frac{n}{2}a - \frac{n}{2}b = \frac{n}{2}(a - b)

我們兩側同除以 nn:

∑yin−∑xiyin=12(a−b)— (5)\frac{\sum y_i}{n} - \frac{\sum x_i y_i}{n} = \frac{1}{2}(a - b) \quad \text{--- (5)}

注意到無約束最優解 (a1,b1)(a_1, b_1) 滿足:

  • ∑yin=a1+12b1\displaystyle \frac{\sum y_i}{n} = a_1 + \frac{1}{2}b_1
  • ∑xiyin=12a1+b1\displaystyle \frac{\sum x_i y_i}{n} = \frac{1}{2}a_1 + b_1

我們將這兩個關係代入式 (5) 的左側:

(a1+12b1)−(12a1+b1)=12(a−b)\left( a_1 + \frac{1}{2}b_1 \right) - \left( \frac{1}{2}a_1 + b_1 \right) = \frac{1}{2}(a - b) 12a1−12b1=12(a−b)  ⟹  a−b=a1−b1— (6)\frac{1}{2}a_1 - \frac{1}{2}b_1 = \frac{1}{2}(a - b) \implies a - b = a_1 - b_1 \quad \text{--- (6)}

現在,我們聯立式 (6) 與約束條件 a+b=ca + b = c:

{a−b=a1−b1a+b=c\begin{cases} a - b = a_1 - b_1 \\ a + b = c \end{cases}

相加可解得 aa:

2a=a1−b1+c  ⟹  a=a1−b1+c22a = a_1 - b_1 + c \implies a = \frac{a_1 - b_1 + c}{2}

相減可解得 bb:

2b=c−(a1−b1)  ⟹  b=−a1+b1+c22b = c - (a_1 - b_1) \implies b = \frac{-a_1 + b_1 + c}{2}

結論:

  1. 無約束最優解為 a1=4∑yi−2∑xiyi3n\displaystyle a_1 = \frac{4\sum y_i - 2\sum x_i y_i}{3n}, b1=4∑xiyi−2∑yi3n\displaystyle b_1 = \frac{4\sum x_i y_i - 2\sum y_i}{3n}。
  2. 有約束最優解為 a=a1−b1+c2\displaystyle a = \frac{a_1 - b_1 + c}{2}, b=−a1+b1+c2\displaystyle b = \frac{-a_1 + b_1 + c}{2}。