跳转至

最佳平方逼近

希尔伯特矩阵

若取 \(\varphi_k(x)=x^k\),则 \((\varphi_i,\varphi_j)=\frac{1}{i+j+1}\)。此时记 \(G_n\)\(H_n\),则:

定义

具有存在唯一性

求解方法

核心:解法方程

\[ \sum^n_{j=0}(\varphi_j(x),\varphi_k(x))a_j=(f(x),\varphi_k(x))\quad k=0,1,\dots,n \]

常见一般情况:

\[ \begin{aligned} \varphi_k(x)&=x^k,f(x)\in C[a,b]\\ \rho(x)&=1,\Phi=\text{span}\{1,x,x^2,……,x^n\} \end{aligned} \]

此时,要求解的 \(P^*_n(x)=a^*_0+a^*_1x+…+a^*_nx^n\)。设:

\[ \begin{aligned} d_k&=\int_a^bf(x)x^k\mathrm dx=(f,\varphi_k) \\h_{ij}&=\int_a^b\varphi_i(x)\varphi_j(x)\mathrm dx=(\varphi_i,\varphi_j) \end{aligned} \]

则只需求解方程:

\[ \pmb{H_na=d} \]

解得的 \(\pmb a\) 即为 \(\pmb a^*=(a_1^*,a^*_2,\dots,a^*_n)\), 由此得到 \(P^*_n(x)\)

特殊情况一:\({\Phi=\text{span}\{\phi_0,\dots,\phi_n\}}\) 为正交基。此时:

\[ a^*_k=\frac{(f,\varphi_k)}{(\varphi_k,\varphi_k)}=\frac{(f,\varphi_k)}{\left\Vert \varphi_k\right\Vert ^2_2} \]

由此得到\(P^*_n(x)\)

特殊情况二:在离散意义下的最佳平方逼近

已知节点\((x_i,y_i), \quad i=1,2,\dots,n\)。矩阵方程变为:

其余步骤类似。

例题9:【2024真题回忆版】 求解\(e^x\) 的最佳平方逼近多项式,空间是 \(\text{span}\{1,x\}\)

解:(因为是回忆版,条件不完整,缺少一个区间,这里以[0,1]为例)

适用情况一

\[ \begin{aligned} d_0&=\int_0^1e^x\mathrm dx=e-1 \\ d_1&=\int_0^1 e^x\cdot x\mathrm dx=1 \\ \pmb d&=\begin{pmatrix} e-1\\ 1 \end{pmatrix} \\ h_{11}&=\int_0^1\mathrm dx=1 \\ h_{12}&=h_{21}=\int_0^1x\mathrm dx=\frac{1}{2} \\ h_{22}&=\int_0^1x^2\mathrm dx=\frac{1}{3} \\ H_2&=\begin{pmatrix} 1&\frac{1}{2}\\ \frac{1}{2}&\frac{1}{3} \end{pmatrix} \end{aligned} \]

\(\pmb{H_2a=d}\) 得:

\[ a_0^*=4e-10,a_1^*=6(3-e) \]

所以 \(P_n^*(x)=a_0^*+a_1^*x=4e-10+6(3-e)x\)

例题10:【课件原题】