Reference no: EM132761062
Question: A coupled system of ODEs representing a 2 mass 3 spring system has the following equations of motion and initial conditions:
x1" = -3xl + 2 x2 + 9 cos(2t), x2"=4x1 - 5x2,x1(0) = 7, x2(0) = 1, x1' (0) = 0, x2'(0) = 0.
a) The solution of this IVP has integer coefficients. Verify this using Maple and write down its solutions for the two unknowns. [Use function notation for all variables in Maple: x1(t), x1"(t). Ask bob if something goes wrong.]
b) Identify the exact and numerical values of the two natural frequencies ω1 < ω2 which appear in these solutions together with the driving frequency ω3 = 2. What are the exact and numerical values of the corresponding 3 periods T1, T2, T3? Give 3 decminal place accuracy. Because the frequency ratios are not all rational numbers, there is no common period.
c) Rewrite this system of DEs and its initial conditions in explicit matrix form x"→ = A x→ + F→ for the vector variable x = (x1→, x2→), identifying the coefficient matrix A and the driving vector F→.
d) Use Maple to write down its choice of eigenvalues and eigenvectors of A, ordered so that |λ1| < |λ2|.
e) By hand showing all (matrix method) steps, find the smallest integer component eigenvectors b1→, b2→ of the coefficient matrix produced by the solution algorithm after rescaling of the standard results by positive multiples if necessary, ordered so that the corresponding eigenvalues satisfy |λ1| < |λ2| i.e., ω1 < ω2. Write down this matrix B =<b1→| b2→> and use Maple to evaluate its inverse and the diagonalized matrix AB = B-1 AB. [Make sure your results agree with Maple's eigenvectors modulo rescaling and/or permutation.]
f) What are the slopes m1, m2 of the lines through the origin containing the two eigenvectors (remember, use integer component eigenvectors)? On the grid provided, draw in those two lines, labeling them by their corresponding coordinates y1, y2 at the ends in the positive direction (include arrow heads) determined by the eigenvectors and then indicate by thicker arrows both eigenvectors b1→, b2→, labeled by their symbols. Indicate the unit tickmarks marked off along each new axis. Recall x→ = By→, y→ = B-1 x→, where y→ = <y1, y2>. Also label the x1, x2 axes.
g) Evaluate y→(0) = B-1x→(0), B-1 F→(0) to find the new components of these two vectors.
h) On the grid provided, draw in the vector x→(0) and label this vector. On your graph, draw in exactly the parallelograms parallel to the new coordinate axes which project this vector along those axes and lightly shade it in pencil (pen?). Are the part g) new components numerically consistent with your plot parallelogram? Explain.
i) Find by hand the general solution of the corresponding decoupled system of DEs y→" = ABy→+ B-1F→. First write these equations out in explicit matrix form, then obtain the two equivalent scalar DEs which are its components. Then solve them to find their general solutions using the method of undetermined coefficients. State your general solutions in scalar form and box them: y1(t) = ...., y2(t), identifying the homogeneous and particular parts of each solution: y1 = y1h + y1p, y2 = y2h + y2p.
j) Then express the general solution for x→ = By→ in explicit matrix form (without multiplying matrix factors) and impose the initial conditions using matrix methods to solve the linear systems. Write out and box the final scalar solutions: x1(t) = .....x2(t). Do they agree with Maple's solution from part a)? If not, look for your error. Did you input the equations correctly?
k) Express the (correct) solution as a sum of the two eigenvector modes and the response mode in the form:
x→ = y1h b1→ + y2h b2→ + cos(2t) a3→ thus identifying the particular solution xp→ (last term), the response vector coefficient a3→ and the homogeneous solution xh→ (first two terms), as well as the final expressions for the two decoupled variables y1h and y2h. Is the response term a tandem or accordian mode? Identify the vectors a1→ = y1h(0) b1→ + y2h(0) b2→. Include these vectors in your plot (together with a3→) and use them to create the sides of the bounding box enclosing the solution with four endpoints ±a1→ ±a2→.
l) Plot the homogeneous solution for large t to "see" the parallelogram box containing it. Is it consistent with your plot on the grid?
Attachment:- eigenvalues and eigenvectors.rar