跳到主要内容

九州大学 工学府 機械系専攻 2024年度 数学 問4

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

(1)

mm 個の点 (xk,yk)(x_k,y_k) を二次多項式

f(x)=a0+a1x+a2x2f(x)=a_0+a_1x+a_2x^2

で最小二乗近似する。

F=k=1m{ykf(xk)}2,a=(a0a1a2),b=(ykykxkykxk2)F=\sum_{k=1}^m\{y_k-f(x_k)\}^2,\qquad \boldsymbol a=\begin{pmatrix}a_0\\a_1\\a_2\end{pmatrix},\qquad \boldsymbol b= \begin{pmatrix} \sum y_k\\ \sum y_kx_k\\ \sum y_kx_k^2 \end{pmatrix}

とする。F/a0=F/a1=F/a2=0\partial F/\partial a_0=\partial F/\partial a_1=\partial F/\partial a_2=0Ma=bM\boldsymbol a=\boldsymbol b と表すとき、MM を求めよ。

(2)

A(1,1),B(0,3),C(1,6),D(2,0)A(-1,1),B(0,3),C(1,6),D(2,0) について、次の問に答えよ。

  1. (1) の行列は

    M=(4262686818)M=\begin{pmatrix}4&2&6\\2&6&8\\6&8&18\end{pmatrix}

    である。M=LLTM=LL^{\mathsf T} となる下三角行列 LL を求めよ。

  2. 4 点を近似する二次多項式を求めよ。

Kai

(1)

j=0,1,2j=0,1,2 に対して

Faj=2k=1mxkj(yki=02aixki)=0.\frac{\partial F}{\partial a_j} =-2\sum_{k=1}^m x_k^j \left(y_k-\sum_{i=0}^2a_ix_k^i\right)=0.

したがって

M=(mxkxk2xkxk2xk3xk2xk3xk4).\boxed{ M= \begin{pmatrix} m&\sum x_k&\sum x_k^2\\ \sum x_k&\sum x_k^2&\sum x_k^3\\ \sum x_k^2&\sum x_k^3&\sum x_k^4 \end{pmatrix}}.

(2)

(2-1)

対角成分を正に取ると

L=(200150352)\boxed{ L=\begin{pmatrix} 2&0&0\\ 1&\sqrt5&0\\ 3&\sqrt5&2 \end{pmatrix}}

であり、実際に LLT=MLL^{\mathsf T}=M となる。

(2-2)

b=(1057).\boldsymbol b= \begin{pmatrix}10\\5\\7\end{pmatrix}.

Lz=bL\boldsymbol z=\boldsymbol bLTa=zL^{\mathsf T}\boldsymbol a=\boldsymbol z の順に解くと

z=(504),a=(9/222).\boldsymbol z= \begin{pmatrix}5\\0\\-4\end{pmatrix},\qquad \boldsymbol a= \begin{pmatrix}9/2\\2\\-2\end{pmatrix}.

よって

f(x)=92+2x2x2.\boxed{f(x)=\frac92+2x-2x^2}.

Reference