Orthonormal Basis
标准正交的向量 满足: 相互垂直并且模都为 1

0 \quad i\neq j\quad \text{Orthogonal}\\ \\ 1 \quad i=j\quad \text{unit} \left\lvert \left\lvert \mathbf{q}_{i} \right\rvert \right\rvert=1 \end{cases}$$ 标准正交的向量作为矩阵的列 将矩阵记为 $Q$ $$\begin{align} Q^{T}Q=\begin{bmatrix} \mathbf{q}_{1}^{T}\\\mathbf{q}_{2}^{T}\\ \vdots\\\mathbf{q}_{n}^{T} \end{bmatrix}\begin{bmatrix} \mathbf{q}_{1} &\mathbf{q}_{2} & \cdots & \mathbf{q}_{n} \end{bmatrix}=\begin{bmatrix} 1 & 0 & \cdots & 0 \\ 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 \end{bmatrix}=I=Q^{-1}Q \end{align}$$ 矩阵的转置等于[[逆矩阵\|逆矩阵]] **Transpose=Inverse** $Q^{T}=Q^{-1}$ - $\left\lvert \left\lvert Q\mathbf{x} \right\rvert \right\rvert=\left\lvert \left\lvert \mathbf{x} \right\rvert \right\rvert$ - $(Q\mathbf{x})^{T}(Q\mathbf{y})=\mathbf{x}^{T}Q^{T}Q\mathbf{y}=\mathbf{x}^{T}\mathbf{y}$ 如果使用标准正交基进行[[向量投影\|向量投影]]: $Q^{T}Q \hat{\mathbf{x}}=\hat{\mathbf{x}}=Q^{T}\mathbf{b}$ $\mathbf{p}=Q \hat{\mathbf{x}}$ $$\begin{align} \hat{\mathbf{x}}&=Q^{T}\mathbf{b} \\ \mathbf{p}&=Q \hat{\mathbf{x}}=QQ^{T}\mathbf{b} =\mathbf{b}\\ \mathbf{b}&=\mathbf{q}_{1}(\mathbf{q}_{1}^{T}\mathbf{b})+\cdots+\mathbf{q}_{n}(\mathbf{q}_{n}^{T}\mathbf{b}) \end{align}$$ $\mathbf{b}=QQ^{T}\mathbf{b}$ 本质上也为[[变换\|变换]]的思想 将向量或函数通过**变换**分解为 perpendicular pieces 然后将 pieces 通过**逆变换**重新得到原来的向量或函数 ### Gram -Schmidt Process 利用线性独立的向量组来构建标准正交基 假设有三个独立的向量 $\mathbf{a},\mathbf{b},\mathbf{c}$ 目的是得到标准正交的向量 $\mathbf{p}_{1},\mathbf{p}_{2},\mathbf{p}_{3}$ 先构造三个相互垂直的向量 $\mathbf{A},\mathbf{B},\mathbf{C}$,再标准化 $\mathbf{A}=\mathbf{a}$ $\mathbf{B}=\mathbf{b}- \dfrac{\mathbf{A}^{T}\mathbf{b}}{\mathbf{A}^{T}\mathbf{A}}\mathbf{A}$ $\mathbf{C}=\mathbf{c}- \dfrac{\mathbf{A}^{T}\mathbf{c}}{\mathbf{A}^{T}\mathbf{A}}\mathbf{A}- \dfrac{\mathbf{B}^{T}\mathbf{c}}{\mathbf{B}^{T}\mathbf{B}}\mathbf{B}$ $\cdots$ 则标准正交的向量: $\mathbf{p}_{1}=\dfrac{\mathbf{A}}{\left\lvert \left\lvert \mathbf{A} \right\rvert \right\rvert},\mathbf{p}_{2}=\dfrac{\mathbf{B}}{\left\lvert \left\lvert \mathbf{B} \right\rvert \right\rvert},\mathbf{p}_{3}=\dfrac{\mathbf{C}}{\left\lvert \left\lvert \mathbf{C} \right\rvert \right\rvert},\cdots$ 其实本质思想很简单 就是向量减去在所有正交基上的投影,得到新的正交基 <svg xmlns="http://www.w3.org/2000/svg" version="1.1" viewBox="0 0 478.9786030719383 260.33435122696505" width="478.9786030719383" height="260.33435122696505"> <!-- svg-source:excalidraw --> <defs> <style class="style-fonts"> @font-face { font-family: "Virgil"; src: url("https://excalidraw.com/Virgil.woff2"); } @font-face { font-family: "Cascadia"; src: url("https://excalidraw.com/Cascadia.woff2"); } @font-face { font-family: "Assistant"; src: url("https://excalidraw.com/Assistant-Regular.woff2"); } </style> </defs> <rect x="0" y="0" width="478.9786030719383" height="260.33435122696505" fill="#ffffff"></rect><g stroke-linecap="round"><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M0 0 C-21.38 54.33, -39.68 111.58, -52.13 138.17 M0 0 C-10.4 28.29, -20.71 58.1, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M-51.62 113.17 C-52.57 122.03, -50.32 134.25, -52.13 138.17 M-51.62 113.17 C-52.11 117.59, -51.64 123.77, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M-35.67 119.34 C-42.97 125.88, -47.08 135.64, -52.13 138.17 M-35.67 119.34 C-39.54 122.54, -42.4 127.43, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M0 0 C40.6 12.91, 77.65 27.97, 104.57 39.57 M0 0 C32.85 12, 64.65 23.28, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M79.58 38.87 C89.93 37.52, 96.87 37.87, 104.57 39.57 M79.58 38.87 C88.24 39.14, 95.5 38.43, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M85.87 22.97 C93.9 27.67, 98.45 34.04, 104.57 39.57 M85.87 22.97 C92.45 28.16, 97.77 32.36, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M0 0 C16.54 48.35, 28.67 100.1, 46.16 172.08 M0 0 C12.29 49.74, 25.28 96.74, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M31.66 151.71 C37.84 156.52, 40.87 164.68, 46.16 172.08 M31.66 151.71 C35.84 158.35, 39.85 162.35, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M48.14 147.16 C49.62 153.3, 47.92 162.76, 46.16 172.08 M48.14 147.16 C47.62 155.03, 46.97 160.32, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(199.150630841636 63.41647208551984) rotate(0 -27.47634311348748 67.82735191621965)"><path d="M0 0 C-21.47 51.77, -42 103.47, -54.95 135.65" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(49.05121663940909 162.3313617138553) rotate(0 47.41634024589943 18.21289469777965)"><path d="M0 0 C39.12 17.04, 73.97 28.2, 94.83 36.43" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M0 0 C-19.14 31.12, -37.06 65.96, -53.02 93.76 M0 0 C-21.89 37.7, -42.94 74.45, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M-49.32 69.03 C-51.5 76.38, -52.49 87.52, -53.02 93.76 M-49.32 69.03 C-51.33 79.14, -52.9 88.63, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M-34.29 77.2 C-41.63 81.77, -47.77 90.11, -53.02 93.76 M-34.29 77.2 C-42.19 84.18, -49.7 90.44, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M0 0 C35.48 6.01, 74.13 7.17, 114.88 12.2 M0 0 C41.43 4.13, 83.31 8.45, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M90.58 18.05 C97.64 17.86, 107.8 13.34, 114.88 12.2 M90.58 18.05 C98.97 15.87, 107.73 13.59, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M92.51 1.05 C98.88 6.2, 108.44 7.03, 114.88 12.2 M92.51 1.05 C100.13 5.1, 108.19 9.03, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M0 0 C19 -19.38, 35.28 -36.71, 69.49 -72.85 M0 0 C22.82 -24.34, 47.23 -48.09, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M59.75 -49.82 C64.09 -56.66, 65.86 -61.89, 69.49 -72.85 M59.75 -49.82 C62.47 -57.58, 67.02 -64.78, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M47.23 -61.47 C54.57 -65.57, 59.38 -67.97, 69.49 -72.85 M47.23 -61.47 C54.04 -65.36, 62.76 -68.69, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M0 0 C1.27 -16.21, -0.56 -29.69, 0 -74.92 M0 0 C0.45 -16.93, 0.81 -35.1, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M8.89 -51.55 C8.14 -57.21, 4.48 -60.11, 0 -74.92 M8.89 -51.55 C7.33 -56.13, 5.58 -62.03, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M-8.21 -51.3 C-5.45 -56.95, -5.6 -59.91, 0 -74.92 M-8.21 -51.3 C-5.7 -56.09, -3.38 -62.05, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(386.4832976507976 52.11191262299832) rotate(0 -45.307933309774455 53.83120468560068)"><path d="M0 0 C-27.37 31.56, -51.96 61.87, -90.62 107.66" stroke="#1e1e1e" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(386.4832976507976 53.009158754423595) rotate(0 4.037299860511936 36.78466193394823)"><path d="M0 0 C2.77 16.63, 2.35 31.25, 8.07 73.57" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(320.5399990273022 47.62597350041045) rotate(0 34.09311789603143 1.5700916499963427)"><path d="M0 0 C15.86 -0.57, 34.23 -0.14, 68.19 3.14" stroke="#1e1e1e" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g transform="translate(10 162.80683187129404) rotate(0 8.033189916846382 13.850489672149791)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.08780812547788px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">A</text></g><g transform="translate(183.47343229601904 10) rotate(0 8.161442583831239 14.07479310682448)"><text x="0" y="22.52014705492415" font-family="Helvetica, Segoe UI Emoji" font-size="24.477901055346916px" fill="#e03131" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">B</text></g><g transform="translate(136.99130975692674 211.70922407754875) rotate(0 8.154903443856668 16.86253607635753)"><text x="0" y="26.980630498533333" font-family="Helvetica, Segoe UI Emoji" font-size="29.326149698013097px" fill="#1971c2" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">b</text></g><g transform="translate(255.53551830304207 222.63337188266547) rotate(0 8.03318991684639 13.850489672149791)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.08780812547788px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">A</text></g><g transform="translate(452.9167749114527 131.5689307864937) rotate(0 8.030914080242809 13.850489672149783)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.087808125477885px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">B</text></g><g transform="translate(418.7018463102315 19.99716409420634) rotate(0 7.329805907797294 16.86253607635753)"><text x="0" y="26.98063049853334" font-family="Helvetica, Segoe UI Emoji" font-size="29.3261496980131px" fill="#1971c2" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">c</text></g><g stroke-linecap="round"><g transform="translate(323.6802795054746 128.82138739816088) rotate(0 27.81278368877254 17.495165817365116)"><path d="M0 0 C16.77 10.95, 31.31 17.83, 55.63 34.99" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(301.2505838925355 153.94265944174325) rotate(0 39.027631495242105 4.485922926224545)"><path d="M0 0 C30.49 4.5, 58.13 5.51, 78.06 8.97" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(394.55802694272774 131.0643893521813) rotate(0 -6.504556660117274 15.925057971005513)"><path d="M0 0 C-4.7 12.84, -10.45 24.22, -13.01 31.85" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(385.13757422092914 62.42954669076906) rotate(0 -2.0186823229824427 50.01817586703695)"><path d="M0 0 C-2.54 25.56, -3.3 46.41, -4.04 100.04" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g transform="translate(289.5004159248151 12.915847472591508) rotate(0 8.837454908730663 14.074793106824472)"><text x="0" y="22.52014705492415" font-family="Helvetica, Segoe UI Emoji" font-size="24.477901055346916px" fill="#e03131" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">C</text></g></svg> ### 矩阵分解 $$\begin{align} \large A=QR \end{align}$$ $$\begin{align} \begin{bmatrix} \\ \mathbf{a} & \mathbf{b} & \mathbf{c} \\ \\ \end{bmatrix}=\begin{bmatrix} \\ \mathbf{q}_{1} & \mathbf{q}_{2} & \mathbf{q}_{3} \\ \\ \end{bmatrix}\begin{bmatrix} \mathbf{q}_{1}^{T}\mathbf{a} & \mathbf{q}_{1}^{T}\mathbf{b} & \mathbf{q}_{1}^{T}\mathbf{c} \\ & \mathbf{q}_{2}^{T}\mathbf{b} & \mathbf{q}_{2}^{T}\mathbf{c} \\ & & \mathbf{q}_{3}^{T}\mathbf{c} \end{bmatrix} \end{align}$$ 由[[向量投影\|向量投影]]知:$A^{T}A \hat{\mathbf{x}}=A^{T}\mathbf{b}$ 所以有:$A^{T}A=(QR)^{T}QR=R^{T}R$ $R^{T}R\hat{\mathbf{x}}=R^{T}Q^{T}\mathbf{b}$ [[最小二乘法\|最小二乘法]]的最优近似解则为:$\hat{\mathbf{x}}=R^{-1}Q^{T}\mathbf{b}$ --- ## AI 结构化补充(2026-05-02) ### 定义 **Orthonormal Basis** 标准正交基是一组既两两正交又已经单位化的基向量;它把坐标分解变成互不耦合的点积计算。 向量 $q_1,\dots,q_n$ 标准正交,当且仅当

q_i^Tq_j=\delta_{ij}

\begin{cases} 1,&i=j,\ 0,&i\ne j. \end{cases}

这里的 $\delta_{ij}$ 是 Kronecker delta。等式同时表达两个条件:$q_i^Tq_i=\|q_i\|^2=1$ 表示单位长度,$q_i^Tq_j=0\ (i\ne j)$ 表示正交。 把这些向量放成矩阵列

Q=[q_1\ q_2\ \cdots\ q_n],

Q^TQ= \begin{bmatrix} q_1^T\ q_2^T\ \vdots\ q_n^T \end{bmatrix} \begin{bmatrix} q_1&q_2&\cdots&q_n \end{bmatrix} =I_n.

这条式子不要求 $Q$ 是方阵;当 $Q\in\mathbb R^{m\times n}$ 且 $m>n$ 时,$Q^T$ 只是左逆,不能写成普通[[逆矩阵\|逆矩阵]]。 ### 与正交矩阵的边界 若 $Q$ 同时是方阵,则 $Q^TQ=I$ 推出

QQ^T=I,\qquad Q^T=Q

这时 $Q$ 是[[正交矩阵\|正交矩阵]],它的行向量和列向量都构成全空间的标准正交基。正交矩阵保持长度与角度:

|Qx|^2=(Qx)^T(Qx)=x^TQ^TQx=|x|

(Qx)^T(Qy)=x^TQ^TQy=x

若 $Q$ 不是方阵,上述保持长度的结论对输入坐标 $x\in\mathbb R^n$ 仍成立,即 $\|Qx\|=\|x\|$;但 $QQ^T$ 不再是恒等矩阵,而是投影到 $Q$ 的列空间。 ### 投影和坐标 设 $Q$ 的列张成子空间 $C(Q)$。把向量 $b$ 投影到 $C(Q)$ 时,标准正交性使正规方程直接退化为

\hat{x}=Q^Tb,\qquad p=Q\hat{x}=QQ

其中 $\hat{x}$ 的第 $i$ 个分量就是 $q_i^Tb$,也就是 $b$ 沿 $q_i$ 的一维投影系数。投影矩阵为

P=QQ^T,\qquad P^2=P,\qquad P

只有当 $Q$ 是覆盖全空间的方阵标准正交基时,才有 $QQ^T=I$,于是 $p=b$,并且

b=q_1(q_1^Tb)+q_2(q_2^Tb)+\cdots+q_n(q_n

若 $Q$ 只给出真子空间的标准正交基,$QQ^Tb$ 只是 $b$ 在该子空间中的最佳近似。 例如

Q=\frac13 \begin{bmatrix} -1&2&2\ 2&-1&2\ 2&2&-1 \end{bmatrix}

的三列是 $\mathbb R^3$ 的标准正交基。对 $b=(0,0,1)^T$,

q_1^Tb=\frac23,\qquad q_2^Tb=\frac23,\qquad q_3

\frac23q_1+\frac23q_2

是 $b$ 到 $\operatorname{span}\{q_1,q_2\}$ 的投影;三个分量全部相加时,因为 $Q$ 是方阵标准正交矩阵,

\frac23q_1+\frac23q_2-\frac13q_3=b.

这说明标准正交基把一个向量拆成互相垂直的一维投影,再由这些分量相加重构。 ### 构造方式 [[格拉姆-施密特正交化\|格拉姆-施密特正交化]]从线性无关向量 $a,b,c$ 出发,先构造正交向量 $A,B,C$,再单位化:

A=a,

B=b-\frac{A^Tb}{A

C=c-\frac{A^Tc}{A^TA}A-\frac{B^Tc}{B

q_1=\frac{A}{|A|},\qquad q_2=\frac{B}{|B|},\qquad q_3=\frac{C}{|C|}.

a=\begin{bmatrix}1\-1\0\end{bmatrix},\qquad b=\begin{bmatrix}2\0\-2\end{bmatrix},\qquad c=\begin{bmatrix}3\-3\3\end{bmatrix}.

A=a,\qquad B=b-\frac{A^Tb}{A

C=c-\frac{A^Tc}{A^TA}A-\frac{B^Tc}{B

此时 $A^TB=A^TC=B^TC=0$,长度分别为 $\sqrt2,\sqrt6,\sqrt3$,所以除以这些长度后得到标准正交列 $q_1,q_2,q_3$。 <svg xmlns="http://www.w3.org/2000/svg" version="1.1" viewBox="0 0 478.9786030719383 260.33435122696505" width="478.9786030719383" height="260.33435122696505"> <!-- svg-source:excalidraw --> <defs> <style class="style-fonts"> @font-face { font-family: "Virgil"; src: url("https://excalidraw.com/Virgil.woff2"); } @font-face { font-family: "Cascadia"; src: url("https://excalidraw.com/Cascadia.woff2"); } @font-face { font-family: "Assistant"; src: url("https://excalidraw.com/Assistant-Regular.woff2"); } </style> </defs> <rect x="0" y="0" width="478.9786030719383" height="260.33435122696505" fill="#ffffff"></rect><g stroke-linecap="round"><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M0 0 C-21.38 54.33, -39.68 111.58, -52.13 138.17 M0 0 C-10.4 28.29, -20.71 58.1, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M-51.62 113.17 C-52.57 122.03, -50.32 134.25, -52.13 138.17 M-51.62 113.17 C-52.11 117.59, -51.64 123.77, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(97.09558382748908 23.85051396669479) rotate(0 -26.063276824854036 69.08341146930798)"><path d="M-35.67 119.34 C-42.97 125.88, -47.08 135.64, -52.13 138.17 M-35.67 119.34 C-39.54 122.54, -42.4 127.43, -52.13 138.17" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M0 0 C40.6 12.91, 77.65 27.97, 104.57 39.57 M0 0 C32.85 12, 64.65 23.28, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M79.58 38.87 C89.93 37.52, 96.87 37.87, 104.57 39.57 M79.58 38.87 C88.24 39.14, 95.5 38.43, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(97.4095859611251 23.850479954331973) rotate(0 52.283577391434626 19.782990396866794)"><path d="M85.87 22.97 C93.9 27.67, 98.45 34.04, 104.57 39.57 M85.87 22.97 C92.45 28.16, 97.77 32.36, 104.57 39.57" stroke="#e03131" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M0 0 C16.54 48.35, 28.67 100.1, 46.16 172.08 M0 0 C12.29 49.74, 25.28 96.74, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M31.66 151.71 C37.84 156.52, 40.87 164.68, 46.16 172.08 M31.66 151.71 C35.84 158.35, 39.85 162.35, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(98.35168306166736 25.420581322146262) rotate(0 23.080131843314852 86.04026928890787)"><path d="M48.14 147.16 C49.62 153.3, 47.92 162.76, 46.16 172.08 M48.14 147.16 C47.62 155.03, 46.97 160.32, 46.16 172.08" stroke="#1971c2" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(199.150630841636 63.41647208551984) rotate(0 -27.47634311348748 67.82735191621965)"><path d="M0 0 C-21.47 51.77, -42 103.47, -54.95 135.65" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(49.05121663940909 162.3313617138553) rotate(0 47.41634024589943 18.21289469777965)"><path d="M0 0 C39.12 17.04, 73.97 28.2, 94.83 36.43" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M0 0 C-19.14 31.12, -37.06 65.96, -53.02 93.76 M0 0 C-21.89 37.7, -42.94 74.45, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M-49.32 69.03 C-51.5 76.38, -52.49 87.52, -53.02 93.76 M-49.32 69.03 C-51.33 79.14, -52.9 88.63, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(317.93831079307057 120.61215702724274) rotate(0 -26.511883694203448 46.87800309468034)"><path d="M-34.29 77.2 C-41.63 81.77, -47.77 90.11, -53.02 93.76 M-34.29 77.2 C-42.19 84.18, -49.7 90.44, -53.02 93.76" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M0 0 C35.48 6.01, 74.13 7.17, 114.88 12.2 M0 0 C41.43 4.13, 83.31 8.45, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M90.58 18.05 C97.64 17.86, 107.8 13.34, 114.88 12.2 M90.58 18.05 C98.97 15.87, 107.73 13.59, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g><g transform="translate(318.25231292670657 121.06069749150274) rotate(0 57.44237822895681 6.100885790791937)"><path d="M92.51 1.05 C98.88 6.2, 108.44 7.03, 114.88 12.2 M92.51 1.05 C100.13 5.1, 108.19 9.03, 114.88 12.2" stroke="#1e1e1e" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M0 0 C19 -19.38, 35.28 -36.71, 69.49 -72.85 M0 0 C22.82 -24.34, 47.23 -48.09, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M59.75 -49.82 C64.09 -56.66, 65.86 -61.89, 69.49 -72.85 M59.75 -49.82 C62.47 -57.58, 67.02 -64.78, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g><g transform="translate(318.297261074003 122.63079885931703) rotate(0 34.74342455550115 -36.425778058105024)"><path d="M47.23 -61.47 C54.57 -65.57, 59.38 -67.97, 69.49 -72.85 M47.23 -61.47 C54.04 -65.36, 62.76 -68.69, 69.49 -72.85" stroke="#1971c2" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M0 0 C1.27 -16.21, -0.56 -29.69, 0 -74.92 M0 0 C0.45 -16.93, 0.81 -35.1, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M8.89 -51.55 C8.14 -57.21, 4.48 -60.11, 0 -74.92 M8.89 -51.55 C7.33 -56.13, 5.58 -62.03, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g><g transform="translate(317.84855216756523 121.19536215376014) rotate(0 -0.00001619636321947837 -37.45756413979066)"><path d="M-8.21 -51.3 C-5.45 -56.95, -5.6 -59.91, 0 -74.92 M-8.21 -51.3 C-5.7 -56.09, -3.38 -62.05, 0 -74.92" stroke="#e03131" stroke-width="2" fill="none"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(386.4832976507976 52.11191262299832) rotate(0 -45.307933309774455 53.83120468560068)"><path d="M0 0 C-27.37 31.56, -51.96 61.87, -90.62 107.66" stroke="#1e1e1e" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(386.4832976507976 53.009158754423595) rotate(0 4.037299860511936 36.78466193394823)"><path d="M0 0 C2.77 16.63, 2.35 31.25, 8.07 73.57" stroke="#1e1e1e" stroke-width="2.5" fill="none" stroke-dasharray="8 10"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(320.5399990273022 47.62597350041045) rotate(0 34.09311789603143 1.5700916499963427)"><path d="M0 0 C15.86 -0.57, 34.23 -0.14, 68.19 3.14" stroke="#1e1e1e" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g transform="translate(10 162.80683187129404) rotate(0 8.033189916846382 13.850489672149791)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.08780812547788px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">A</text></g><g transform="translate(183.47343229601904 10) rotate(0 8.161442583831239 14.07479310682448)"><text x="0" y="22.52014705492415" font-family="Helvetica, Segoe UI Emoji" font-size="24.477901055346916px" fill="#e03131" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">B</text></g><g transform="translate(136.99130975692674 211.70922407754875) rotate(0 8.154903443856668 16.86253607635753)"><text x="0" y="26.980630498533333" font-family="Helvetica, Segoe UI Emoji" font-size="29.326149698013097px" fill="#1971c2" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">b</text></g><g transform="translate(255.53551830304207 222.63337188266547) rotate(0 8.03318991684639 13.850489672149791)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.08780812547788px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">A</text></g><g transform="translate(452.9167749114527 131.5689307864937) rotate(0 8.030914080242809 13.850489672149783)"><text x="0" y="22.1612539404421" font-family="Helvetica, Segoe UI Emoji" font-size="24.087808125477885px" fill="#1e1e1e" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">B</text></g><g transform="translate(418.7018463102315 19.99716409420634) rotate(0 7.329805907797294 16.86253607635753)"><text x="0" y="26.98063049853334" font-family="Helvetica, Segoe UI Emoji" font-size="29.3261496980131px" fill="#1971c2" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">c</text></g><g stroke-linecap="round"><g transform="translate(323.6802795054746 128.82138739816088) rotate(0 27.81278368877254 17.495165817365116)"><path d="M0 0 C16.77 10.95, 31.31 17.83, 55.63 34.99" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(301.2505838925355 153.94265944174325) rotate(0 39.027631495242105 4.485922926224545)"><path d="M0 0 C30.49 4.5, 58.13 5.51, 78.06 8.97" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(394.55802694272774 131.0643893521813) rotate(0 -6.504556660117274 15.925057971005513)"><path d="M0 0 C-4.7 12.84, -10.45 24.22, -13.01 31.85" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g stroke-linecap="round"><g transform="translate(385.13757422092914 62.42954669076906) rotate(0 -2.0186823229824427 50.01817586703695)"><path d="M0 0 C-2.54 25.56, -3.3 46.41, -4.04 100.04" stroke="#1971c2" stroke-width="1.5" fill="none" stroke-dasharray="8 9"></path></g></g><mask></mask><g transform="translate(289.5004159248151 12.915847472591508) rotate(0 8.837454908730663 14.074793106824472)"><text x="0" y="22.52014705492415" font-family="Helvetica, Segoe UI Emoji" font-size="24.477901055346916px" fill="#e03131" text-anchor="start" style="white-space: pre;" direction="ltr" dominant-baseline="alphabetic">C</text></g></svg> ### 与 QR 分解 若 $A=[a_1\ \cdots\ a_n]$ 的列线性无关,Gram-Schmidt 产生的标准正交列组成 $Q=[q_1\ \cdots\ q_n]$。每个 $a_j$ 都只需要前 $j$ 个 $q_i$ 表示:

a_j=r_{1j}q_1+\cdots+r_{jj}q_j.

A=QR,\qquad R=Q

R= \begin{bmatrix} q_1^Ta_1&q_1^Ta_2&\cdots&q_1^Ta_n\ 0&q_2^Ta_2&\cdots&q_2^Ta_n\ \vdots&\vdots&\ddots&\vdots\ 0&0&\cdots&q_n^Ta_n \end{bmatrix}

是上三角矩阵。最小二乘问题 $\min_x\|Ax-b\|$ 因 $A=QR$ 化为

A^TA=R^TQ^TQR=R

R\hat{x}=Q

这比直接处理 $A^TA\hat{x}=A^Tb$ 更适合计算,因为核心只剩下标准正交投影和上三角回代。