没有合适的资源?快使用搜索试试~ 我知道了~
首页Matrix Analysis and Applied Linear Algebra Solutions[Carl D.Meyer]
资源详情
资源评论
资源推荐


Solutions for Chapter 1
Solutions for exercises in section 1. 2
1.2.1. (1, 0, 0)
1.2.2. (1, 2, 3)
1.2.3. (1, 0, −1)
1.2.4. (−1/2, 1/2, 0, 1)
1.2.5.
2
−4 3
4
−7 4
5
−8 4
1.2.6. Every row operation is reversible. In particular the “inverse” of any row operation
is again a row operation of the same type.
1.2.7.
π
2
,π,0
1.2.8. The third equation in the triangularized form is 0x
3
=1, which is impossible
to solve.
1.2.9. The third equation in the triangularized form is 0x
3
=0, and all numbers are
solutions. This means that you can start the back substitution with any value
whatsoever and consequently produce infinitely many solutions for the system.
1.2.10. α = −3,β=
11
2
, and γ = −
3
2
1.2.11. (a) If x
i
= the number initially in chamber #i, then
.4x
1
+0x
2
+0x
3
+ .2x
4
=12
0x
1
+ .4x
2
+ .3x
3
+ .2x
4
=25
0x
1
+ .3x
2
+ .4x
3
+ .2x
4
=26
.6x
1
+ .3x
2
+ .3x
3
+ .4x
4
=37
and the solution is x
1
=10,x
2
=20,x
3
=30, and x
4
=40.
(b) 16, 22, 22, 40
1.2.12. To interchange rows i and j, perform the following sequence of Type II and
Type III operations.
R
j
← R
j
+ R
i
(replace row j by the sum of row j and i)
R
i
← R
i
− R
j
(replace row i by the difference of row i and j)
R
j
← R
j
+ R
i
(replace row j by the sum of row j and i)
R
i
←−R
i
(replace row i by its negative)
1.2.13. (a) This has the effect of interchanging the order of the unknowns— x
j
and
x
k
are permuted. (b) The solution to the new system is the same as the

2 Solutions
solution to the old system except that the solution for the j
th
unknown of the
new system is ˆx
j
=
1
α
x
j
. This has the effect of “changing the units” of the j
th
unknown. (c) The solution to the new system is the same as the solution for
the old system except that the solution for the k
th
unknown in the new system
is ˆx
k
= x
k
− αx
j
.
1.2.14. h
ij
=
1
i+j−1
1.2.16. If x =
x
1
x
2
.
.
.
x
m
and y =
y
1
y
2
.
.
.
y
m
are two different solutions, then
z =
x + y
2
=
x
1
+y
1
2
x
2
+y
2
2
.
.
.
x
m
+y
m
2
is a third solution different from both x and y.
Solutions for exercises in section 1. 3
1.3.1. (1, 0, −1)
1.3.2. (2, −1, 0, 0)
1.3.3.
111
122
123
Solutions for exercises in section 1. 4
1.4.2. Use y
(t
k
)=y
k
≈
y
k+1
− y
k−1
2h
and y
(t
k
)=y
k
≈
y
k−1
− 2y
k
+ y
k+1
h
2
to write
f(t
k
)=f
k
= y
k
−y
k
≈
2y
k−1
− 4y
k
+2y
k+1
2h
2
−
hy
k+1
− hy
k−1
2h
2
,k=1, 2,...,n,
with y
0
= y
n+1
=0. These discrete approximations form the tridiagonal system
−42− h
2+h −42− h
.
.
.
.
.
.
.
.
.
2+h −42− h
2+h −4
y
1
y
2
.
.
.
y
n−1
y
n
=2h
2
f
1
f
2
.
.
.
f
n−1
f
n
.

Solutions 3
Solutions for exercises in section 1. 5
1.5.1. (a) (0, −1) (c) (1, −1) (e)
1
1.001
,
−1
1.001
1.5.2. (a) (0, 1) (b) (2, 1) (c) (2, 1) (d)
2
1.0001
,
1.0003
1.0001
1.5.3. Without PP: (1.01, 1.03) With PP: (1, 1) Exact: (1, 1)
1.5.4. (a)
1 .500 .333
.333
.500 .333 .250
.333
.333 .250 .200
.200
−→
1 .500 .333
.333
0 .083 .083
.166
0 .083 .089
.089
−→
1 .500 .333 .333
0 .083 .083
.166
00.006
−.077
z = −.077/.006 = −12.8,
y =(.166 −.083z)/.083 = 14.8,x= .333 − (.5y + .333z)=−2.81
(b)
1 .500 .333
.333
.500 .333 .250
.333
.333 .250 .200
.200
−→
1 .500 .333
.333
1 .666 .500
.666
1 .751 .601
.601
−→
1 .500 .333
.333
0 .166 .167
.333
0 .251 .268
.268
−→
1 .500 .333
.333
0 .251 .268
.268
0 .166 .167
.333
−→
1 .500 .333
.333
0 .251 .268
.268
00−.01
.156
z = −.156/.01 = −15.6,
y =(.268 −.268z)/.251 = 17.7,x= .333 − (.5y + .333z)=−3.33
(c)
1 .500 .333
.333
.500 .333 .250
.333
.333 .250 .200
.200
−→
1 .500 .333
.333
1 .666 .500
.666
1 .751 .601
.601
−→
1 .500 .333
.333
0 .166 .167
.333
0 .251 .268
.268
−→
1 .500 .333
.333
0 .994 1
1.99
0 .937 1
1
−→
1 .500 .333
.333
0 .994 1
1.99
00.057
−.880
z = −.88/.057 = −15.4,
y =(1.99 −z)/.994 = 17.5,x= .333 − (.5y + .333z)=−3.29
(d) x = −3,y=16,z= −14
1.5.5. (a)
.0055x + .095y + 960z = 5000
.0011x + . 01y + 112z = 600
.0093x + .025y + 560z = 3000

4 Solutions
(b) 3-digit solution = (55, 900 lbs. silica, 8, 600 lbs. iron, 4.04 lbs. gold).
Exact solution (to 10 digits) = (56, 753.68899, 8, 626.560726, 4.029511918). The
relative error (rounded to 3 digits) is e
r
=1.49 ×10
−2
.
(c) Let u = x/2000,v= y/1000, and w =12z to obtain the system
11u +95v +80w = 5000
2.2u +10v +9.33w = 600
18.6u +25v +46.7w = 3000.
(d) 3-digit solution = (28.5 tons silica, 8.85 half-tons iron, 48.1 troy oz. gold).
Exact solution (to 10 digits) = (28.82648317, 8.859282804, 48.01596023). The
relative error (rounded to 3 digits) is e
r
=5.95 × 10
−3
. So, partial pivoting
applied to the column-scaled system yields higher relative accuracy than partial
pivoting applied to the unscaled system.
1.5.6. (a) (−8.1, −6.09) = 3-digit solution with partial pivoting but no scaling.
(b) No! Scaled partial pivoting produces the exact solution—the same as with
complete pivoting.
1.5.7. (a) 2
n−1
(b) 2
(c) This is a famous example that shows that there are indeed cases where par-
tial pivoting will fail due to the large growth of some elements during elimination,
but complete pivoting will be successful because all elements remain relatively
small and of the same order of magnitude.
1.5.8. Use the fact that with partial pivoting no multiplier can exceed 1 together with
the triangle inequality |α + β|≤|α|+ |β| and proceed inductively.
Solutions for exercises in section 1. 6
1.6.1. (a) There are no 5-digit solutions. (b) This doesn’t help—there are now infinitely
many 5-digit solutions. (c) 6-digit solution = (1.23964, −1.3) and exact solution
=(1, −1) (d) r
1
= r
2
= 0 (e) r
1
= −10
−6
and r
2
=10
−7
(f) Even if computed
residuals are 0, you can’t be sure you have the exact solution.
1.6.2. (a) (1, −1.0015) (b) Ill-conditioning guarantees that the solution will be very
sensitive to some small perturbation but not necessarily to every small perturba-
tion. It is usually difficult to determine beforehand those perturbations for which
an ill-conditioned system will not be sensitive, so one is forced to be pessimistic
whenever ill-conditioning is suspected.
1.6.3. (a) m
1
(5) = m
2
(5) = −1.2519,m
1
(6) = −1.25187, and m
2
(6) = −1.25188
(c) An optimally well-conditioned system represents orthogonal (i.e., perpen-
dicular) lines, planes, etc.
1.6.4. They rank as (b) = Almost optimally well-conditioned. (a) = Moderately well-
conditioned. (c) = Badly ill-conditioned.
1.6.5. Original solution = (1, 1, 1). Perturbed solution = (−238, 490, −266). System
is ill-conditioned.
剩余171页未读,继续阅读















安全验证
文档复制为VIP权益,开通VIP直接复制

评论0