Skip to content
CalcGospel 微積分福音
返回

112 政治大學微積分 第 7 題

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

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

題目

Problem

3. Let S(a,b)=i(yiabxi)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)

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

Sa=2i=1n(yiabxi)=0    i=1nyi=na+bi=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)} Sb=2i=1nxi(yiabxi)=0    i=1nxiyi=ai=1nxi+bi=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ˉ=1nyi)12a+b=1nxiyi\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ˉ1nxiyi    32a=2yˉ1nxiyi2a + 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ˉ323nxiyi=4yi2xiyi3na_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)=2nxiyiyˉ    32b=2nxiyiyˉ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=4xiyi2yi3nb_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(yiabxi)2+μ(a+bc)L(a, b, \mu) = \sum_{i=1}^n (y_i - a - b x_i)^2 + \mu (a + b - c)

aabb 求偏導數並令其為零:

La=2i=1n(yiabxi)+μ=0    2(yinabxi)=μ— (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)} Lb=2i=1nxi(yiabxi)+μ=0    2(xiyiaxibxi2)=μ— (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

yinabxi=xiyiaxibxi2\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

yinan2b=xiyin2anb\sum y_i - na - \frac{n}{2}b = \sum x_i y_i - \frac{n}{2}a - nb yixiyi=n2an2b=n2(ab)\sum y_i - \sum x_i y_i = \frac{n}{2}a - \frac{n}{2}b = \frac{n}{2}(a - b)

我們兩側同除以 nn

yinxiyin=12(ab)— (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(ab)\left( a_1 + \frac{1}{2}b_1 \right) - \left( \frac{1}{2}a_1 + b_1 \right) = \frac{1}{2}(a - b) 12a112b1=12(ab)    ab=a1b1— (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

{ab=a1b1a+b=c\begin{cases} a - b = a_1 - b_1 \\ a + b = c \end{cases}

相加可解得 aa

2a=a1b1+c    a=a1b1+c22a = a_1 - b_1 + c \implies a = \frac{a_1 - b_1 + c}{2}

相減可解得 bb

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

結論:

  1. 無約束最優解為 a1=4yi2xiyi3n\displaystyle a_1 = \frac{4\sum y_i - 2\sum x_i y_i}{3n}b1=4xiyi2yi3n\displaystyle b_1 = \frac{4\sum x_i y_i - 2\sum y_i}{3n}
  2. 有約束最優解為 a=a1b1+c2\displaystyle a = \frac{a_1 - b_1 + c}{2}b=a1+b1+c2\displaystyle b = \frac{-a_1 + b_1 + c}{2}