6.1
-
5, 8, 8/5 3. 3 / 35 -1 / 35 -1 / 7 5. 8 / 13 12 / 13
-
√(35) 9. [ c -0.6 0.8 ] 11. [ c7 / √(69) 2 / √(69) 4 / √(69) ]
-
5√(5) 15.不正交
-
正交
-
写出答案之后参考“学习指导”.
-
提示:用 2.1 节的定理3 和 2.
-
u · v = 0, |u|² = 30, |v|² = 101, |u + v|² = (-5)² + (-9)² + 5² = 131 = 30 + 101.
-
所有 -b a 的倍数的集合(当 v ≠ 0 ).
-
提示:用正交性的定义.
-
提示:考虑一个属于 W 的典型向量 w = c₁v₁ + … + cₚvₚ .
-
提示:如果 x 属于 W^(perp) ,那么 x 与 W 中任一向量正交.
-
[M]说明你的猜测并给出代数验证.
6.2
-
不正交 3. 不正交 5. 正交
-
证明 ±bu₁ · ±bu₂ = 0 ,注意到定理4且观察到 ℝ² 中两个线性无关向量构成一个基,那么得到 ±bx = (39)/(13) 2 -3 + (26)/(52) 6 4 = 3 2 -3 + 1/2 6 4
-
证明 ±bu₁ · ±bu₂ = 0, ±bu₁ · ±bu₃ = 0 且 ±bu₂ · ±bu₃ = 0.注意到定理4且观察到 ℝ³ 中线性无关向量构成一个基,那么得到 x = 5/2 u₁ - (27)/(18) u₂ + (18)/(9) u₃ = 5/2 u₁ - 3/2 u₂ + 2u₃ .
-
[ c-2 1 ] 13. y = [ c-4 / 5 7 / 5 ] + [ c14 / 5 8 / 5 ]
-
y - haty = 0.6 -0.8 ,距离是1.
-
单位正交 21. 单位正交
-
见“学习指导”.
-
提示: |Ux|²=(Ux)^(T)(Ux) ,同样,(a)和(c)可从(b)得到.
-
提示:你需要两个定理,其中一个仅应用于方阵.
-
提示:如果你有一个候补逆,则可以检查这个候补是否可行.
-
假设 haty = (y · u)/(u · u)u.对 c ≠ 0 用 cu 代替 u ; 那么 (y · (cu))/((cu) · (cu))(cu) = (c(y · u))/(c²u · u)(c)u = haty
-
设 L = operatornameSpan{±bu} ,其中 ±bu 非零,记 T(±bx) = proj_L±bx ,根据定义, T(±bx) = frac±bx · ±bu±bu · ±bu±bu = (±bx · ±bu)(±bu · ±bu)⁻1±bu 对 ℝⁿ 中的 ±bx 和 ±by 以及任意数 c 和 d ,内积的性质(定理1)说明 T(cx + dy) = [(cx + dy) · u](u · u)⁻1u = [c(x · u) + d(y · u)](u · u)⁻1u = c(x · u)(u · u)⁻1u + d(y · u)(u · u)⁻1u = cT(x) + dT(y)
因此 T 是线性的.
6.3
-
x = -8/9u₁ - 2/9u₂ + 2/3u₃ + 2u₄; x = 0 -2 4 -2 + 10 -6 -2 2 3. -1 4 0 5. -1 2 6=y
-
y=10/32/38/3+-7/37/37/3 9. y=2400+2-13-1
-
√(40)
-
a. U^(T)U=1 & 0 0 & 1, UU^(T)=8/9 & -2/9 & 2/9 -2/9 & 5/9 & 4/9 2/9 & 4/9 & 5/9
b. operatornameproj_(w)y=6u₁+3u₂=2 4 5 ,且 (UU^(T))y=2 4 5
-
0 2 / 5 1 / 5 的任何倍数,如 0 2 1 .
-
在查阅“学习指导”之前,写出答案.
-
提示:用定理3 和正交分解定理.对唯一性,假设 Ap = b 和 Ap₁ = b ,考虑方程 p = p₁ + (p - p₁) 和 p = p + 0 .
-
[M]由6.2节定理6,因为 U^(T)U = Iₖ ,所以 U 有单位正交列。Col U 中离 ±by 最近的点是 ±by 在Col U 。由定理10, haty 。由定理10, haty = UU^(T)y = (1.2,0.4,1.2,1.2,0.4,1.2,0.4,0.4)
6.4
- R=6 & 12 0 & 6
-
见“学习指导”.
-
若 x 满足 Rx = 0 ,那么 QRx = Q0 = 0 ,且 Ax = 0 。由于 A 的列是线性无关的,因此 x 必为零。这个事实反过来说明 R 的列是线性无关的。因为 R 是方阵,由可逆矩阵定理可知它是可逆的。
-
记 Q 的列为 q₁, …, qₙ ,注意 n ≤slant m ,因为 A 是 m × n 且有线性无关的列。利用以下事实:Q 的列可以扩充为 R^(m) 的一个标准正交基,如 {q₁, …, qₘ} (“学习指导”描述了一个方法)。取 Q₀ = [qₙ+1 … qₘ] 那么 Q₁ = [Q Q₀] 。利用分块矩阵乘积, Q₁ R 0 = QR = A 。
-
提示:将 R 作为 2 × 2 分块矩阵.
-
[M] R 的对角元素是 20, 6, 10.3923 和 7.0711, 精确到 4 位小数.
6.5
-
a. 6 & -11 -11 & 22 x₁ x₂ = -4 11 b. hatx = 3 2
-
a. 6 & 6 6 & 42 x₁ x₂ = 6 -6 b. hatx = 4 / 3 -1 / 3
-
hatx=5 -3 0+x₃-1 1 1
-
2√(5)
-
a. hat±bb = 1 1 0 b. hat±bx = 2 / 7 1 / 7
-
a.
-
A±b u = [ c11 -11 11 ], A±b v = [ c7 -12 7 ], ±b b - A±b u = [ c0 2 -6 ], ±bb - A±bv = 4 3 -2.否,±bu 不可能是 Ax = ±bb 的一个最小二乘解.为什么?
-
hatx=4 -1 17. 见“学习指导”.
-
a.如果 Ax=0,那么 A^(T)Ax=A^(T)0=0 ,这说明 Nul A 包含在 Nul A^(T)A 中. b. 如果 A^(T)Ax=0 ,那么 x^(T)A^(T)Ax=x^(T)0=0 ,所以 (Ax)^(T)(Ax)=0 (这说明 |Ax|²=0 ),因此 Ax=0。这说明 Nul A^(T)A 包含在 Nul A 中。
-
提示:对(a),利用第2章的一个重要定理.
-
由定理14, hat±bb = Ahat±bx = A(A^(T)A)⁻1A^(T)±bb.矩阵 A(A^(T)A)⁻1A^(T) 在统计中经常出现,通常称为帽子矩阵.
-
法方程是 2 & 2 2 & 2 x y = 6 6 ,它的解是 (x, y) 的集合使得 x + y = 3 。解对应于位于直线 x + y = 2 和 x + y = 4 中间的点集。
6.6
-
y = 0.9 + 0.4x 3. y = 1.1 + 1.3x
-
如果两个点集有不同的 x 坐标,那么设计矩阵 X 的两列不可能是相互的倍数,因而它们线性无关。由 6.5 节的定理14,法方程有唯一解。
-
a. y = Xβ + ε ,其中 y = 1.8 2.7 3.4 3.8 3.9 , X = 1 & 1 2 & 4 3 & 9 4 & 16 5 & 25 , β=β₁ β₂, ε=ε₁ ε₂ ε₃ ε₄ ε₅ b. [M] y = 1.76x - 0.20x² 9. y = Xβ + ε ,其中 y = 7.9 5.4 -0.9 , X = cos 1 & sin 1 cos 2 & sin 2 cos 3 & sin 3
-
[M] β = 1.45 和 e = 0.811 ;轨迹一个椭圆。对于方程 r = β / (1 - e · cos θ) ,当 θ = 4.6 时, r = 1.33 。
-
[M] a. y = -0.8558 + 4.7025t + 5.5554t² - 0.0274t³ b. 速度函数是 v(t)=4.7025+11.1108t-0.0822t² v(4.5)=53.0 英尺/秒
-
提示:写出在方程(1)中的 X 和 y,且计算 X^(T)X 和 X^(T)y .
-
a. x 数据的平均是 barx=5.5.数据的平均偏差形式是 (-3.5,1),(-0.5,2),(1.5,3),(2.5,3) . X 的列是正交的,因为第 2 列中元素之和为零. b. 4&0 0&21β₀ β₁=9 7.5,y=9/4+(5)/(14)x^(*)=9/4+(5)/(14)(x-5.5)
-
提示:方程有一个好的几何解释.
6.7
-
a. 3, √(105) , 225 b. 所有 1 4 的倍数 3.28 5. 5√(2) , 3√(3) 7. (56)/(25)+(14)/(25)t
-
a.常数多项式, p(t)=5 b. t²-5 正交于 p₀ 和 p₁ ;值: (4,-4,-4,4) ;答案: q(t)=1/4(t²-5) 11. (17)/(5)t
-
验证 4 个公理中的每一个,例如
-
langle u,vrangle=(Au)·(Av) 定义 =(Av)·(Au) 点积的性质 =langle v,urangle 定义
-
langle u, cv rangle = langle cv, u rangle 公理1 公理3 公理1= clangle v, u rangle 公理3 = clangle u, v rangle 公理1
-
提示:计算 4 乘右手边.
-
langle ±bu,±b vrangle = √(a)√(b) +√(b)√(a) = 2√(ab) , | ±b u| ² = (√(a))²+ (√(b))² = a + b.因为 a,b 非负,所以 | ±b u| = √(a + b) 类似地, | ±b v| = √(b + a).由柯西-施瓦茨不等式, 2√(ab)≤slant √(a + b)√(b + a) = a + b ,因而 √(ab)≤slant (a + b)/(2)
21.0 23. 2 / √(5) 25.1, t , 3t² - 1
- [M] 新的正交多项式是 -17t + 5t³ 和 72 - 155t² + 35t⁴ 的倍数.重新度量这些多项式使得它们在 -2, -1, 0, 1 和 2 的值是小整数.
6.8
-
y = 2 + 3/2t
-
p(t) = 4p₀ - 0.1p₁ - 0.5p₂ + 0.2p₃ = 4 - 0.1t - 0.5(t² - 2) + 0.2(5/6 t³ - (17)/(6) t) (这个多项式精确地拟合了数据.)
-
利用等式 sin mt sin nt = 1/2 [cos(mt - nt) - cos(mt + nt)]
-
利用等式 cos²kt=(1+cos2kt)/(2) .
-
π + 2sin t + sin 2t + 2/3sin 3t (提示:利用例4 的结果可节约时间).
-
1/2 -1/2cos 2t (为什么?)
-
提示:选取 C[0,2π] 中的函数 f 和 g ,且固定一个整数 m ≥slant 0 。写出 f + g 的包含 cos mt 的傅里叶系数且写出包含 sin mt 的傅里叶系数 (m > 0) 。
-
[M]立方曲线是下列函数的图形: g(t) = -0.2685 + 3.6095t + 5.8576t² - 0.0477t³.在 t = 4.5 秒的速度是 g'(4.5) = 53.4 英尺/秒.这比 6.6 节习题 13 中估计得到的结果快 0.7%.
a. F b. T c. T d. F e. F f. T g. T h. T i. F j. T k. T l. F m. T n. F o. F p. T q. T r. F s. F
- 提示:若 {v₁, v₂} 是单位正交集合且 x = c₁v₁ + c₂v₂ ,则向量 c₁v₁, c₂v₂ 是正交的且 |x|²=|c₁v₁+c₂v₂|²=|c₁v₁|²+|c₂v₂|² =(|c₁||v₁|)²+(|c₂||v₂|)²=|c₁|²+|c₂|²
(解释为什么). 故对于 p=2 状态方程成立.假设等式对 p=k ( k≥slant2 ) 成立,令 {v₁,…,vₖ+1} 是单位正交集,考虑 x=c₁v₁+…+cₖvₖ+cₖ+1vₖ+1=uₖ+cₖ+1vₖ+1 ,其中 uₖ=c₁v₁+…+cₖvₖ .
-
给定 x 和 Rⁿ 中一个单位正交集 {v₁,…,vₚ} ,设 hatx 是 x 在由 v₁,…,vₚ 生成的子空间上的正交投影。由 6.3 节的定理10, hatx=(x· v₁)v₁+…+(x· vₚ)vₚ 。由习题 2, |hatx|²=|x· v₁|²+…+|x· vₚ|² 。贝塞尔不等式可从 |hatx|²≤slant|x|² 得到,该不等式在 6.7 节柯西-施瓦茨不等式的证明前面给出。
-
假设对任意 x, y 属于 Rⁿ 有 (Ux) · (Uy) = x · y ,且设 e₁, …, eₙ 是 Rⁿ 的标准基。对 j = 1, …, n , Ueⱼ 是 U 的第 j 列。由于 |Ueⱼ|² = (Ueⱼ) · (Ueⱼ) = eⱼ · eⱼ = 1 ,故 U 的列是单位向量;由于对 j ≠ k , (Ueⱼ) · (Ueₖ) = eⱼ · eₖ = 0 ,故列是两两正交的。
-
提示: 计算 Q^(T)Q ,利用事实 (uu^(T))^(T)=u^(TT)u^(T)=uu^(T)
-
设 W = operatornameSpan{±bu,±bv}.给定 ℝⁿ 中的 z,设 hatz = operatornameproj_W z,那么 hatz 属于 operatornameCol A,其中 A = [±bu ±bv],例如,对 ℝ² 中的一些 hatx,有 hatz = Ahatx.所以 hatx 是 Ax = z 的一个最小二乘解,法方程可以解出 hatx,那么 hatz 可通过计算 Ahatx 得到.
-
提示:设 x=x y z, b=a b c, v=1 -2 5 和 A=v^(T) v^(T) v^(T)= 1 & -2 & 5 1 & -2 & 5 1 & -2 & 5.给定的方程组是 Ax=b,且所有
最小二乘解的集合.与 A^(T)Ax=A^(T)b 的解集相一致(6.5节定理13). 这可利用事实 (vv^(T))x=v(x^(T)x)=(v^(T)x)v ,其中 v^(T)x 是数.
- a. Au 的行-列计算表明 A 的每一行是与 operatornameNul A 中的每一个 ±bu 正交的。所以 A 的每一行属于 (operatornameNul A)^(perp) 。由于 (operatornameNul A)^(perp) 是一个子空间,故它必须包含 A 的行的所有线性组合;因此 (operatornameNul A)^(perp) 包含 operatornameRow A 。
b.如果 rank A = r ,那么由秩定理, dim Nul A = n - r. 由6.3节的习题24(c),有 dimNul A+dim(Nul A)^(perp)=n 所以 dim(Nul A)^(perp) 一定是 r.但由秩定理和(a)可知 Row A 是一个 (Nul A)^(perp) 的 r 维子空间,所以 Row A 必和 (Nul A)^(perp) 一致.
c.在(b)中用 A^(T) 代替 A 且可得 operatornameRow A^(T) 和 (operatornameNul A^(T))^(perp) 一致.又由于 operatornameRow A^(T) = operatornameCol A 这就证明了(c).
-
如果 A = URU^(T) 且 U 是正交的,那么 A 与 R 相似(因为 U 是可逆的且 U^(T) = U⁻1 ),且 A 与 R 具有相同特征值(5.2节定理4),即 n 个在 R 对角线上的实数.
-
[M] (|Δ x|)/(|x|) = 0.4618, cond(A) × (|Δ b|)/(|b|) = 3363 × (1.548 × 10⁻4) = 0.5206. 注意 |Δ x| / |x| 约等于 cond(A) 乘 |Δ b| / |b| .
-
[M] (|Δ x|)/(|x|) = 7.178 × 10⁻8, (|Δ b|)/(|b|) = 2.832 × 10⁻4,注意 x 的相对改变比 b 的相对改变小很多.事实上,由于 operatornamecond(A) × (|Δ b|)/(|b|) = 23683 × (2.832 × 10⁻4) = 6.707 故 x 相对改变的理论上界是 6.707 (到 4 位有效数字). 这个习题表明,既使一个条件数很大,一个解的相对误差也并不像你想象的那么大.