/*! This file is auto-generated */ .wp-block-button__link{color:#fff;background-color:#32373c;border-radius:9999px;box-shadow:none;text-decoration:none;padding:calc(.667em + 2px) calc(1.333em + 2px);font-size:1.125em}.wp-block-file__button{background:#32373c;color:#fff;text-decoration:none} Q2.5-32Q The band matrix A shown below ... [FREE SOLUTION] | 91Ó°ÊÓ

91Ó°ÊÓ

The band matrixAshown below can be used to estimate
the unsteady conduction of heat in a rod when the temperatures at points \({p_{\bf{1}}},...,{p_{\bf{5}}}\) on the rod change with time.

The constantCin the matrix depends on the physical nature
of the rod, the distance \(\Delta x\) between the points on the rod,
and the length of time \(\Delta t\) between successive temperature
measurements. Suppose that for \(k = 0,1,2,...,\) a vector \({{\bf{t}}_k}\)
in \({\mathbb{R}^5}\) lists the temperatures at time \(k\Delta t\). If the two ends of the
rod are maintained at \(0^\circ \), then the temperature vectors satisfy
the equation \(A{t_{k + 1}} = {t_k}\)\(\left( {k = 0,1,...} \right)\), where \[A = \left[ {\begin{array}{*{20}{c}}{\left( {{\bf{1}} + {\bf{2}}C} \right)}&{ - C}&{}&{}&{}\\{ - C}&{\left( {{\bf{1}} + {\bf{2}}C} \right)}&{ - C}&{}&{}\\{}&{ - C}&{\left( {{\bf{1}} + {\bf{2}}C} \right)}&{ - C}&{}\\{}&{}&{ - C}&{\left( {{\bf{1}} + {\bf{2}}C} \right)}&{ - C}\\{}&{}&{}&{ - C}&{\left( {{\bf{1}} + {\bf{2}}C} \right)}\end{array}} \right]\]

  1. Find the LU factorization ofAwhen\[C = 1\]. A matrix
    such as Awith three nonzero diagonals is called a tridiagonal matrix. The Land Ufactors are bidiagonal matrices.
  1. Suppose \[C = 1\] and\({t_0} = \left( {10,12,12,12,10} \right)\). Use the
    LU factorization of Ato find the temperature distributions \({t_1}\), \({t_2}\), \({t_3}\), and \({t_4}\).

Short Answer

Expert verified
  1. The LU factorizations are

\(L = \left[ {\begin{array}{*{20}{c}}1&0&0&0&0\\{ - 1/3}&1&0&0&0\\0&{ - 3/8}&1&0&0\\0&0&{ - 8/21}&1&0\\0&0&0&{ - 21/55}&1\end{array}} \right]\),\(U = \left[ {\begin{array}{*{20}{c}}3&{ - 1}&0&0&0\\0&{8/3}&{ - 1}&0&0\\0&0&{21/8}&{ - 1}&0\\0&0&0&{55/21}&{ - 1}\\0&0&0&0&{144/55}\end{array}} \right]\).

2. The temperature distributions are \({{\bf{t}}_1} = \left[ {\begin{array}{*{20}{c}}{6.5556}\\{9.6667}\\{10.4444}\\{9.6667}\\{6.5556}\end{array}} \right]\), \({{\bf{t}}_2} = \left[ {\begin{array}{*{20}{c}}{4.7407}\\{7.6667}\\{8.5926}\\{7.6667}\\{4.7407}\end{array}} \right]\), \({{\bf{t}}_3} = \left[ {\begin{array}{*{20}{c}}{3.5988}\\{6.0556}\\{6.9012}\\{6.0556}\\{3.5988}\end{array}} \right]\), \({{\bf{t}}_4} = \left[ {\begin{array}{*{20}{c}}{2.7922}\\{4.7778}\\{5.4856}\\{4.7778}\\{2.7922}\end{array}} \right]\).

Step by step solution

01

(a) Step 1: State the LU factorization of matrix A

Consider the matrix at\(C = 1\)as shown below:

\(A = \left[ {\begin{array}{*{20}{c}}3&{ - 1}&0&0&0\\{ - 1}&3&{ - 1}&0&0\\0&{ - 1}&3&{ - 1}&0\\0&0&{ - 1}&3&{ - 1}\\0&0&0&{ - 1}&3\end{array}} \right]\)

Use the following MATLAB command to obtain theLU factorization of the matrix:

\[ > > \left[ {\begin{array}{*{20}{c}}{\rm{L}}&{\rm{U}}\end{array}} \right] = {\rm{lu}}\left( {\rm{A}} \right)\]

The LU factorization is shown below:

\(L = \left[ {\begin{array}{*{20}{c}}1&0&0&0&0\\{ - 1/3}&1&0&0&0\\0&{ - 3/8}&1&0&0\\0&0&{ - 8/21}&1&0\\0&0&0&{ - 21/55}&1\end{array}} \right]\)

\(U = \left[ {\begin{array}{*{20}{c}}3&{ - 1}&0&0&0\\0&{8/3}&{ - 1}&0&0\\0&0&{21/8}&{ - 1}&0\\0&0&0&{55/21}&{ - 1}\\0&0&0&0&{144/55}\end{array}} \right]\)

Use the following MATLAB command to obtain the\(LU - A\)factorization of the matrix:

\[ > > LU - A\]

The output is shown below:

\(\begin{array}{c}LU - A = \left[ {\begin{array}{*{20}{c}}0&0&0&0&0\\0&0&0&0&0\\0&0&0&0&0\\0&0&0&0&0\\0&0&0&0&0\end{array}} \right]\\ = {\bf{0}}\end{array}\)

Thus, it is verified that \(LU = A\).

02

(b) Step 2: Solve the temperature distribution

Suppose the solution of\(L{{\bf{s}}_{k + 1}} = {{\bf{t}}_k}\)is\({{\bf{s}}_{k + 1}}\), where \(\left( {k = 0,1,...} \right)\). So, the solution of\(U{{\bf{t}}_{k + 1}} = {{\bf{s}}_{k + 1}}\)is\({{\bf{t}}_{k + 1}}\).

Use the following MATLAB command to solve the system using LU factorization:

\[\begin{array}{l} > > {\rm{b}} = \left[ {\begin{array}{*{20}{c}}{10}&{12}&{12}&{12}&{10}\end{array}} \right]\\ > > {\rm{X}} = {\rm{A}}\backslash {\rm{b}}'\end{array}\]

The solutions are shown below:

\({{\bf{s}}_1} = \left[ {\begin{array}{*{20}{c}}{10.0000}\\{15.3333}\\{17.7500}\\{18.7619}\\{17.1636}\end{array}} \right]\),\({{\bf{t}}_1} = \left[ {\begin{array}{*{20}{c}}{6.5556}\\{9.6667}\\{10.4444}\\{9.6667}\\{6.5556}\end{array}} \right]\), \({{\bf{s}}_2} = \left[ {\begin{array}{*{20}{c}}{6.5556}\\{11.8519}\\{14.8889}\\{15.3386}\\{12.4121}\end{array}} \right]\), \({{\bf{t}}_2} = \left[ {\begin{array}{*{20}{c}}{4.7407}\\{7.6667}\\{8.5926}\\{7.6667}\\{4.7407}\end{array}} \right]\), \({{\bf{s}}_3} = \left[ {\begin{array}{*{20}{c}}{4.7407}\\{9.2469}\\{12.0602}\\{12.2610}\\{9.4222}\end{array}} \right]\), \({{\bf{t}}_3} = \left[ {\begin{array}{*{20}{c}}{3.5988}\\{6.0556}\\{6.9012}\\{6.0556}\\{3.5988}\end{array}} \right]\), \({{\bf{s}}_4} = \left[ {\begin{array}{*{20}{c}}{3.5988}\\{7.2551}\\{9.6219}\\{9.7210}\\{7.3104}\end{array}} \right]\), \({{\bf{t}}_4} = \left[ {\begin{array}{*{20}{c}}{2.7922}\\{4.7778}\\{5.4856}\\{4.7778}\\{2.7922}\end{array}} \right]\).

Thus, the temperature distributions are \({{\bf{t}}_1} = \left[ {\begin{array}{*{20}{c}}{6.5556}\\{9.6667}\\{10.4444}\\{9.6667}\\{6.5556}\end{array}} \right]\), \({{\bf{t}}_2} = \left[ {\begin{array}{*{20}{c}}{4.7407}\\{7.6667}\\{8.5926}\\{7.6667}\\{4.7407}\end{array}} \right]\), \({{\bf{t}}_3} = \left[ {\begin{array}{*{20}{c}}{3.5988}\\{6.0556}\\{6.9012}\\{6.0556}\\{3.5988}\end{array}} \right]\), \({{\bf{t}}_4} = \left[ {\begin{array}{*{20}{c}}{2.7922}\\{4.7778}\\{5.4856}\\{4.7778}\\{2.7922}\end{array}} \right]\).

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with 91Ó°ÊÓ!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

Let Abe an invertible \(n \times n\) matrix, and let B be an \(n \times p\) matrix. Show that the equation \(AX = B\) has a unique solution \({A^{ - 1}}B\).

When a deep space probe launched, corrections may be necessary to place the probe on a precisely calculated trajectory. Radio elementary provides a stream of vectors, \({{\bf{x}}_{\bf{1}}},....,{{\bf{x}}_k}\), giving information at different times about how the probe’s position compares with its planned trajectory. Let \({X_k}\) be the matrix \(\left[ {{x_{\bf{1}}}.....{x_k}} \right]\). The matrix \({G_k} = {X_k}X_k^T\) is computed as the radar data are analyzed. When \({x_{k + {\bf{1}}}}\) arrives, a new \({G_{k + {\bf{1}}}}\) must be computed. Since the data vector arrive at high speed, the computational burden could be serve. But partitioned matrix multiplication helps tremendously. Compute the column-row expansions of \({G_k}\) and \({G_{k + {\bf{1}}}}\) and describe what must be computed in order to update \({G_k}\) to \({G_{k + {\bf{1}}}}\).

Suppose block matrix \(A\) on the left side of (7) is invertible and \({A_{{\bf{11}}}}\) is invertible. Show that the Schur component \(S\) of \({A_{{\bf{11}}}}\) is invertible. [Hint: The outside factors on the right side of (7) are always invertible. Verify this.] When \(A\) and \({A_{{\bf{11}}}}\) are invertible, (7) leads to a formula for \({A^{ - {\bf{1}}}}\), using \({S^{ - {\bf{1}}}}\) \(A_{{\bf{11}}}^{ - {\bf{1}}}\), and the other entries in \(A\).

Let T be a linear transformation that maps \({\mathbb{R}^n}\) onto \({\mathbb{R}^n}\). Is \({T^{ - 1}}\) also one-to-one?

In Exercise 9 mark each statement True or False. Justify each answer.

9. a. In order for a matrix B to be the inverse of A, both equations \(AB = I\) and \(BA = I\) must be true.

b. If A and B are \(n \times n\) and invertible, then \({A^{ - {\bf{1}}}}{B^{ - {\bf{1}}}}\) is the inverse of \(AB\).

c. If \(A = \left( {\begin{aligned}{*{20}{c}}a&b\\c&d\end{aligned}} \right)\) and \(ab - cd \ne {\bf{0}}\), then A is invertible.

d. If A is an invertible \(n \times n\) matrix, then the equation \(Ax = b\) is consistent for each b in \({\mathbb{R}^{\bf{n}}}\).

e. Each elementary matrix is invertible.

See all solutions

Recommended explanations on Math Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.