- LAB 4.3.1: (Undamped, forced harmonic oscillator) First consider solutions to ...
- LAB 4.3.2: (Damped, forced oscillator) Using the same value of k1 as in Part 1...
- LAB 4.3.3: (Small-amplitude periodic solutions for the forced mass-spring syst...
- LAB 4.3.4: (Large-amplitude periodic solutions for the forced mass-spring syst...
Solutions for Chapter LAB 4.3: The Tacoma Narrows Bridge
Full solutions for Differential Equations 00 | 4th Edition
Tv = Av + Vo = linear transformation plus shift.
Augmented matrix [A b].
Ax = b is solvable when b is in the column space of A; then [A b] has the same rank as A. Elimination on [A b] keeps equations correct.
Upper triangular systems are solved in reverse order Xn to Xl.
cond(A) = c(A) = IIAIlIIA-III = amaxlamin. In Ax = b, the relative change Ilox III Ilx II is less than cond(A) times the relative change Ilob III lib II· Condition numbers measure the sensitivity of the output to change in the input.
Cross product u xv in R3:
Vector perpendicular to u and v, length Ilullllvlll sin el = area of parallelogram, u x v = "determinant" of [i j k; UI U2 U3; VI V2 V3].
Dot product = Inner product x T y = XI Y 1 + ... + Xn Yn.
Complex dot product is x T Y . Perpendicular vectors have x T y = O. (AB)ij = (row i of A)T(column j of B).
Gram-Schmidt orthogonalization A = QR.
Independent columns in A, orthonormal columns in Q. Each column q j of Q is a combination of the first j columns of A (and conversely, so R is upper triangular). Convention: diag(R) > o.
A symmetric matrix with eigenvalues of both signs (+ and - ).
Left inverse A+.
If A has full column rank n, then A+ = (AT A)-I AT has A+ A = In.
Left nullspace N (AT).
Nullspace of AT = "left nullspace" of A because y T A = OT.
The pivot row j is multiplied by eij and subtracted from row i to eliminate the i, j entry: eij = (entry to eliminate) / (jth pivot).
Orthogonal matrix Q.
Square matrix with orthonormal columns, so QT = Q-l. Preserves length and angles, IIQxll = IIxll and (QX)T(Qy) = xTy. AlllAI = 1, with orthogonal eigenvectors. Examples: Rotation, reflection, permutation.
Positive definite matrix A.
Symmetric matrix with positive eigenvalues and positive pivots. Definition: x T Ax > 0 unless x = O. Then A = LDLT with diag(D» O.
Projection matrix P onto subspace S.
Projection p = P b is the closest point to b in S, error e = b - Pb is perpendicularto S. p 2 = P = pT, eigenvalues are 1 or 0, eigenvectors are in S or S...L. If columns of A = basis for S then P = A (AT A) -1 AT.
Reflection matrix (Householder) Q = I -2uuT.
Unit vector u is reflected to Qu = -u. All x intheplanemirroruTx = o have Qx = x. Notice QT = Q-1 = Q.
Row space C (AT) = all combinations of rows of A.
Column vectors by convention.
If x gives the movements of the nodes, K x gives the internal forces. K = ATe A where C has spring constants from Hooke's Law and Ax = stretching.
Sum V + W of subs paces.
Space of all (v in V) + (w in W). Direct sum: V n W = to}.
Symmetric matrix A.
The transpose is AT = A, and aU = a ji. A-I is also symmetric.
v + w = (VI + WI, ... , Vn + Wn ) = diagonal of parallelogram.