11b

11.3 QR-разложение: 11.3.1 Вычисление регрессионных величин, 11.3.2 Алгоритмы QR- и WU-разложений

[7/86%]
Показать
LaTeX
Задача 11b.1

Пусть q1,…,qp\mathbf{q}_{1}, \ldots , \mathbf{q}_{p} — векторы, построенные в алгоритме 11.3. Покажите, что q1,…,qj\mathbf{q}_{1}, \ldots , \mathbf{q}_{j} — ортонормированный базис для C(a(1),…,a(j))\mathcal{C}\left(\mathbf{a}^{(1)}, \ldots , \mathbf{a}^{(j)}\right) при j=1,2,…,pj = 1,2, \ldots , p.

?
Задача 11b.2

Покажите, что верхнетреугольная матрица с положительными диагональными элементами невырождена. Также докажите, что обратная к верхнетреугольной матрице с единичными диагональными элементами также является верхнетреугольной с единичными диагональными элементами.

?
Задача 11b.3

Выведите формулы (11.23), (11.24) и (11.25).

?
Задача 11b.4

Докажите, что в начале этапа jj алгоритма MGSA столбцы a(1),…,a(j)\mathbf{a}^{(1)}, \ldots , \mathbf{a}^{(j)} ортогональны, ортогональны a(j+1),…,a(p)\mathbf{a}^{(j+1)}, \ldots , \mathbf{a}^{(p)} и порождают первые jj столбцов исходной матрицы.

?
Задача 11b.5

Покажите, что произведение ортогональных матриц ортогонально.

?
Задача 11b.6

Покажите, что если Qp,Qn−p\mathbf{Q}_{p}, \mathbf{Q}_{n-p} и r1\mathbf{r}_{1} определены как в (11.36), то предсказанные значения Y^\widehat{\mathbf{Y}} задаются формулой

Y^=(Qp,Qn−p)(r10). \widehat{\mathbf{Y}} = \left(\mathbf{Q}_{p}, \mathbf{Q}_{n-p}\right)\binom {\mathbf{r}_{1}}{\mathbf{0}} .
?
Задача 11b.7

Пусть Gih\mathbf{G}_{i h} задано формулой (11.40)(11.40) и D~=diag⁡(d~1,…,d~n)\widetilde{\mathbf{D}} = \operatorname {diag}\left(\tilde{d}_{1}, \ldots , \tilde{d}_{n}\right). Покажите, что существует диагональная матрица D=diag⁡(d1,…,dn)\mathbf{D} = \operatorname {diag}\left(d_{1}, \ldots , d_{n}\right) и быстрый ротатор вида (11.45), такие что Gih=DMD~−1\mathbf{G}_{i h} = \mathbf{D M} \tilde{\mathbf{D}}^{-1}. Покажите, что

dl=d~l(l≠i,l≠h)di=sd~hdh=sd~iα=r/tη=1/rt \begin{aligned} d_{l} & = \tilde{d}_{l} \quad (l \neq i, l \neq h) \\ d_{i} & = s \tilde{d}_{h} \\ d_{h} & = s \tilde{d}_{i} \\ \alpha & = r / t \\ \eta & = 1 / r t \end{aligned}

где t=s/ct = s / c и r=d~h/d~ir = \tilde{d}_{h} / \tilde{d}_{i}.

?