(a) Modify the function ex_with_2eqs to solve the IVP (4) for 0 t 45 using the...
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
(a) Modify the function ex_with_2eqs to solve the IVP (4) for 0 t 45 using the MATLAB routine ode45. Call the new function LAB04ex1. Let [t,Y] (note the upper case Y) be the output of ode45 and y and v the unknown functions. Use the following commands to define the ODE: function dYdt= f(t,Y) y=Y (1); v=Y (2); dYdt = [v; -6* sin(t)-6*v-2*y]; Plot y(t) and v(t) in the same window (do not use subplot), and the phase plot showing v vs y in a separate window. Add a legend to the first plot. (Note: to display v(t) = y(t), use 'v(t)=y'' (t)'). Add a grid. Use the command ylim ([-3.9,3.9]) to adjust the y-limits for both plots. Adjust the x-limits in the phase plot so as to reproduce the pictures in Figure 3 2 1 0 -1 www -2 -3 0 10 20 t y(t) v(t) 30 40 A-a 3 2 1 0 -1 -2 -3 -2 0 Y 2 Figure Time series y = y(t) and v=v(t) = y'(t) (left), and phase plot v = y' vs. y for (4). (b) By reading the matrix Y and the vector t, find (approximately) the last three values of t in the interval 0 t 45 at which y reaches a local maximum. Note that, because the M-file LAB04ex1.m is a function file, all the variables are local and thus not available in the Command Window. To read the matrix Y and the vector t, you need to modify the M-file by adding the line [t, Y(:,1), Y(:,2)]. Do not include the whole output in your lab write-up. Include only the values necessary to answer the question, i.e. just the rows of [t, y, v] with local y-maxima and the adjacent rows. To quickly locate the desired rows, recall that the local maxima of a differentiable function appear where its derivative changes sign from positive to negative. (Note: Due to numerical approximations and the fact that the numerical solution is not necessarily computed at the exact t-values where the maxima occur, you should not expect v(= y) to be exactly 0 at local maxima, but only close to 0). (c) What seems to be the long term behavior of y? (d) Modify the initial conditions to y(0) = 1.1, v(0) = -1.2 and run the file LAB04ex1.m with the modified initial conditions. Based on the new graphs, determine whether the long term behavior of the solution changes. Explain. Include the pictures with the modified initial conditions to support your answer. (a) Modify the function ex_with_2eqs to solve the IVP (4) for 0 t 45 using the MATLAB routine ode45. Call the new function LAB04ex1. Let [t,Y] (note the upper case Y) be the output of ode45 and y and v the unknown functions. Use the following commands to define the ODE: function dYdt= f(t,Y) y=Y (1); v=Y (2); dYdt = [v; -6* sin(t)-6*v-2*y]; Plot y(t) and v(t) in the same window (do not use subplot), and the phase plot showing v vs y in a separate window. Add a legend to the first plot. (Note: to display v(t) = y(t), use 'v(t)=y'' (t)'). Add a grid. Use the command ylim ([-3.9,3.9]) to adjust the y-limits for both plots. Adjust the x-limits in the phase plot so as to reproduce the pictures in Figure 3 2 1 0 -1 www -2 -3 0 10 20 t y(t) v(t) 30 40 A-a 3 2 1 0 -1 -2 -3 -2 0 Y 2 Figure Time series y = y(t) and v=v(t) = y'(t) (left), and phase plot v = y' vs. y for (4). (b) By reading the matrix Y and the vector t, find (approximately) the last three values of t in the interval 0 t 45 at which y reaches a local maximum. Note that, because the M-file LAB04ex1.m is a function file, all the variables are local and thus not available in the Command Window. To read the matrix Y and the vector t, you need to modify the M-file by adding the line [t, Y(:,1), Y(:,2)]. Do not include the whole output in your lab write-up. Include only the values necessary to answer the question, i.e. just the rows of [t, y, v] with local y-maxima and the adjacent rows. To quickly locate the desired rows, recall that the local maxima of a differentiable function appear where its derivative changes sign from positive to negative. (Note: Due to numerical approximations and the fact that the numerical solution is not necessarily computed at the exact t-values where the maxima occur, you should not expect v(= y) to be exactly 0 at local maxima, but only close to 0). (c) What seems to be the long term behavior of y? (d) Modify the initial conditions to y(0) = 1.1, v(0) = -1.2 and run the file LAB04ex1.m with the modified initial conditions. Based on the new graphs, determine whether the long term behavior of the solution changes. Explain. Include the pictures with the modified initial conditions to support your answer.
Expert Answer:
Answer rating: 100% (QA)
Solutions function dYdt ftY y Y1 v Y2 dYdt v cost4v3y end function ty LAB04ex1 ... View the full answer
Related Book For
Posted Date:
Students also viewed these databases questions
-
The following additional information is available for the Dr. Ivan and Irene Incisor family from Chapters 1-5. Ivan's grandfather died and left a portfolio of municipal bonds. In 2012, they pay Ivan...
-
You are required to write a Python program that will manage character (heroes and villain) information. Character (hero and villain) information will be stored in a text file that will be read in...
-
August 8, 2019, was a grim day at the San Francisco headquarters of Uber, the global leader in ride-hailing companies. That day it announced a staggering loss of $5.2 billion in just the previous 3...
-
When your clients need tax-sheltered life insurance coverage, you always recommend a particular insurance agent. This agent appreciates your referrals and frequently refers clients who need...
-
Prepare a report for the managing director of your company explaining how costs may be classified by their behaviour, with particular reference to the effects both on total and on unit costs. Your...
-
What stakeholder affected by their behavior was Gilead weighting very lightly when it decided what to tell the FDA about the medicines it wished to have permission to sell?
-
Wells Technical Institute (WTI), a school owned by Tristana Wells, provides training to individuals who pay tuition directly to the school. WTI also offers training to groups in off-site locations....
-
What other options exist if you do not want to use the relational model? Researchon the Internet and find one other type of database model, describe that model, and post the link to the source that...
-
Use anonymous FTP to connect to ftp.rfc-editor.org (directory innotes), and retrieve the RFC index. Also retrieve the protocol specifications for TCP, IP, and UDP.
-
Pete Maverick was invited to his friend Bradley Bradshaw's housewarming dinner party. The dinner party took place in the wooden cabin located in the garden of Bradley's new home. When Pete arrived at...
-
Explain how an effective performance management system along with compensation and benefits can attract, develop and retain talented employees. Discuss labor and employment laws and regulations that...
-
You have been asked to send an email to Stacey Mintern. The email contained some background notes attached about an important meeting she has scheduled for tomorrow with the other departments. You...
-
A firm that is currently locating facilities in a large number of other countries in order to capitalize on lower production and distribution costs is at what level of global participation?...
-
Describe the legal, social, and psychological employment contracts in organisations. Citing examples, explain why it is important for HRM to manage each of these three contracts at the workplace....
-
What special compensation and benefits issues arise in international HRM ?
-
You work as a director of leadership and learning for a company that makes prosthetics. The company has been in business since 1999 and employs over 350 people in its Tampa, Florida, headquarters and...
-
Refer to the data for problem 13-36 regarding Long Beach Pharmaceutical Company. Required: Compute each division's residual income for the year under each of the following assumptions about the...
-
Evaluate the integral. 2.x (1/2) xe fo (1 + 2x) Jo dx
-
Differentiate. f(x) = 1 xe x /x + e x
-
Calculate y. y = e 1/x /x 2
-
For coordinates \(\left(x^{1}, x^{2} ight)\) and metric \(g=\operatorname{diag}\left(g_{11}, g_{22} ight)\), the Gaussian curvature is For a sphere with coordinates defined in the following figure,...
-
Consider the holonomic basis defined in Box 26.1 . Using that the tangent vector for a curve can be written \(t=t^{\mu} e_{\mu}=\left(d x^{\mu} / d \lambda ight) e_{\mu}\), show that Thus, \(g_{\mu...
-
The Lie bracket of vector fields \(A\) and \(B\) is defined as their commutator, \([A, B]=\) \(A B-B A\). The Lie bracket of two basis vectors vanishes for a coordinate basis but not for a...
Study smarter with the SolutionInn App