13.8

Метод наименьших квадратов и интерполяция функций

[13/54%]
Показать
LaTeX
Задача 1395

Рассмотрим некоторые вектор xx и подпространство LL в евклидовом пространстве Rn\mathbb {R}^{n}. Доказать, что среди векторов подпространства LL ортогональная проекция x∥x^{\| }вектора xx на LL «наименее уклоняется» от данного вектора xx, т. е.

∣x−x∥∣⩽∣x−l∣∀l∈L \left|x-x^{\| }\right| \leqslant \left|x-l\right| \quad \forall l \in L
?
Задача 1396

Доказать, что псевдорешение можно найти как решение новой неоднородной системы уравнений, составленной при помощи матрицы Грама G(a1,a2,…,an)G\left(a_{1}, a_{2}, \ldots , a_{n}\right) :

[(a1,a1)(a1,a2)…(a1,an)⋮⋮⋮(an,a1)(an,a2)…(an,an)][x1⋮xn]=[(b,a1)⋮(b,an)] \left[\begin{smallmatrix} \left(a_{1}, a_{1}\right) & \left(a_{1}, a_{2}\right) & \ldots & \left(a_{1}, a_{n}\right) \\ \vdots & \vdots & & \vdots \\ \left(a_{n}, a_{1}\right) & \left(a_{n}, a_{2}\right) & \ldots & \left(a_{n}, a_{n}\right) \end{smallmatrix}\right]\left[\begin{smallmatrix} x_{1} \\ \vdots \\ x_{n} \end{smallmatrix}\right] =\left[\begin{smallmatrix} \left(b, a_{1}\right) \\ \vdots \\ \left(b, a_{n}\right) \end{smallmatrix}\right]

Каков геометрический смысл найденных констант x1,x2,…,xnx_{1}, x_{2}, \ldots , x_{n}?

?
Задача 1397

Методом наименьших квадратов найти псевдорешение следующих несовместных систем линейных уравнений:

?
(1)

{x+2y=12x+4y=1x+y=0\begin{cases} x+2 y = 1 \\ 2 x+4 y = 1 \\ x+y = 0 \end{cases}

(2)

{x=22x−y=1,x−2y=−1,x+y=2\begin{cases} x = 2 \\ 2 x-y = 1, \\ x-2 y = -1, \\ x+y = 2 \end{cases}

(3)

{2x1−x3=1,x2+x3=−1,x1−x2+x3=0,x1−x3=−1;\begin{cases} 2 x^{1}-x^{3} = 1, \\ x^{2}+x^{3} = -1, \\ x^{1}-x^{2}+x^{3} = 0, \\ x^{1}-x^{3} = -1 ; \end{cases}

(4)

{x1−x3=1,x2+x3−x4=−1,x1−x2=0,−x1+x3+2x4=1,x2+x3−x4=1.\begin{cases} x^{1}-x^{3} = 1, \\ x^{2}+x^{3}-x^{4} = -1, \\ x^{1}-x^{2} = 0, \\ -x^{1}+x^{3}+2 x^{4} = 1, \\ x^{2}+x^{3}-x^{4} = 1. \end{cases}

Задача 1398

Доказать, что система многочленов φk(x),k=0,…,n\varphi_{k}(x), k = 0, \ldots , n, из (13.1) образует ортонормированный базис пространства многочленов степени не выше nn относительно скалярного произведения

(f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(xn)g(xn)−∞<x0<x1<…<xn<+∞ \begin{aligned} (f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right) & +\ldots +f\left(x_{n}\right) g\left(x_{n}\right) \\ & -\infty < x_{0} < x_{1} < \ldots < x_{n} < +\infty \end{aligned}
?
Задача 1399

Даны 2n+12 n+1 попарно различных точек xi,i=0,…,2nx_{i}, i = 0, \ldots , 2 n, отрезка [0,2π][0,2 \pi ]. В пространстве V2n+1V_{2 n+1} тригонометрических многочленов

V2n+1={c0+c1cos⁡x+c2sin⁡x+…+c2n−1cos⁡nx+c2nsin⁡nx} V_{2 n+1} = \left\{ c_{0}+c_{1} \cos x+c_{2} \sin x+\ldots +c_{2 n-1} \cos n x+c_{2 n} \sin n x\right\}

найти такие тригонометрические многочлены tkt_{k}, что

tk(xj)={1,j=k0,j≠k t_{k}\left(x_{j}\right) = \begin{cases} 1, & j = k \\ 0, & j \neq k\end{cases}

Показать, что tkt_{k} образуют ортонормированный базис пространства V2n+1V_{2 n+1} относительно скалярного произведения

(f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(x2n)g(x2n) (f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right)+\ldots +f\left(x_{2 n}\right) g\left(x_{2 n}\right)
?
Задача 1400

Доказать, что в пространстве V2n+1V_{2 n+1} тригонометрических многочленов существует, и притом единственное, решение T2n+1\mathscr {T}_{2 n+1} системы уравнений

T2n+1(xi)=f(xi),i=0,…,2n \mathscr {T}_{2 n+1}\left(x_{i}\right) = f\left(x_{i}\right), \quad i = 0, \ldots , 2 n

где x0,…,x2n−2n+1x_{0}, \ldots , x_{2 n}-2 n+1 попарно различных точек отрезка [0,2π][0,2 \pi ], a f:[a,b]→Rf:[a, b] \rightarrow \mathbb {R} — некоторая функция.

?
Задача 1401

Найти многочлен hh минимальной степени такой, что в данных попарно различных точках x0,x1,…,xnx_{0}, x_{1}, \ldots , x_{n} его значения и производные совпадают со значениями и производными данной функции:

h(xi)=f(xi),h′(xi)=f′(xi),0⩽i⩽n h\left(x_{i}\right) = f\left(x_{i}\right), \quad h^{\prime }\left(x_{i}\right) = f^{\prime }\left(x_{i}\right), \quad 0 \leqslant i \leqslant n

С помощью интерполяционного многочлена Лагранжа LnfL_{n}^{f} можно интерполировать произвольную функцию ff на числовой прямой по таблице ее значений f(x0),f(x1),…,f(xn)f\left(x_{0}\right), f\left(x_{1}\right), \ldots , f\left(x_{n}\right), где −∞<x0<……<xn<+∞-\infty < x_{0} < \ldots \ldots < x_{n} < +\infty.

?
Задача 1402

Доказать, что если множество {xi,i=0,…,n}\left\{ x_{i}, i = 0, \ldots , n\right\} узлов интерполяции симметрично относительно нуля, то из четности (нечетности) функции f(x)f(x) следует четность (нечетность) функции Lnf(x)L_{n}^{f}(x).

?
Задача 1403

Введем пространство Vn+1V_{n+1} дискретных функций f:Mn+1→Rf: M_{n+1} \rightarrow \mathbb {R} на конечном множестве Mn+1={x0,…,xn}M_{n+1} = \left\{ x_{0}, \ldots , x_{n}\right\}. Показать, что Vn+1V_{n+1} является евклидовым пространством относительно скалярного произведения (f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(xn)g(xn)(f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right)+\ldots +f\left(x_{n}\right) g\left(x_{n}\right). Указать какой-нибудь ортонормированный базис пространства Vn+1V_{n+1}.

?
Задача 1404

Доказать, что решение Lmf(x)L_{m}^{f}(x) задачи интерполяции по методу наименьших квадратов в пространстве Vn+1V_{n+1} существует и единственно, причем коэффициенты ai,i=0,…,ma_{i}, i = 0, \ldots , m, находятся из системы неоднородных линейных уравнений

G(φ0,φ1,…,φm)[a0⋮am]=[(f,φ0)⋮(f,φm)] G\left(\varphi _{0}, \varphi _{1}, \ldots , \varphi _{m}\right)\left[\begin{smallmatrix} a_{0} \\ \vdots \\ a_{m} \end{smallmatrix}\right] =\left[\begin{smallmatrix} \left(f, \varphi _{0}\right) \\ \vdots \\ \left(f, \varphi _{m}\right) \end{smallmatrix}\right]

где G(φ0,φ1,…,φm)G\left(\varphi_{0}, \varphi_{1}, \ldots , \varphi_{m}\right) — матрица Грама системы функций φ0(x),……,φm(x)\varphi_{0}(x), \ldots \ldots , \varphi_{m}(x) в евклидовом пространстве Vn+1V_{n+1}.

?
Задача 1405

Методом наименьших квадратов найти интерполяцию (наилучшее среднеквадратичное приближение) функции ff, заданной значениями f(0)=1,f(1)=2,f(2)=3,f(3)=5f(0) = 1, f(1) = 2, f(2) = 3, f(3) = 5 :

?
(1)

линейным многочленом L1f(x)=b1x+b0L_{1}^{f}(x) = b_{1} x+b_{0} (т. е. интерполировать по системе функций φ0(x)=1,φ1(x)=x\varphi_{0}(x) = 1, \varphi_{1}(x) = x);

(2)

квадратичным многочленом L2f(x)=b2x2+b1x+b0L_{2}^{f}(x) = b_{2} x^{2}+b_{1} x+b_{0};

(3)

многочленом 3 -й степени L3f(x)=b3x3+b2x2+b1x+b0L_{3}^{f}(x) = b_{3} x^{3}+b_{2} x^{2}+b_{1} x+b_{0}.

Найти квадратичное отклонение ∥Lif−f∥4\left\| L_{i}^{f}-f\right\|_{4} в каждом из этих случаев.

Задача 1406

Доказать, что в случае ортонормированной системы функций φ0(x),…,φn(x)\varphi_{0}(x), \ldots , \varphi_{n}(x) коэффициенты ai,i=0,…,na_{i}, i = 0, \ldots , n, равняются коэффициентам Фурье функции ff, т. е. линейная комбинация с коэффициентами Фурье функции ff есть наилучшая интерполяция функции ff.

?
Задача 1407

Методом наименьших квадратов найти в евклидовом пространстве непрерывно дифференцируемых функций со скалярным произведением

(f,g)=∫02πf(x)g(x)dx (f, g) = \int _{0}^{2 \pi } f(x) g(x) d x

интерполяцию (наилучшее среднеквадратичное приближение) непрерывно дифференцируемой функции ff при помощи системы функций φ0(x)=1,φ1(x)=cos⁡x,φ2(x)=sin⁡x\varphi_{0}(x) = 1, \varphi_{1}(x) = \cos x, \varphi_{2}(x) = \sin x.

?