

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Instructions for solving second order differential equations using maple and analyzing an anemometer model and an economics model. It includes homework problems for the math 308 - differential equations course in the fall of 2003. The anemometer model involves a spring-mass system and the effects of wind on the mass, while the economics model deals with the relationship between the price of a commodity and its supply.
Typology: Assignments
1 / 3
This page cannot be seen from the preview
Don't miss anything!


Math 308 - Differential Equations Fall 2003
Homework Assignment 2 Due Wednesday, October 29.
(a) Consider the differential equation
y′′^ + 3y′^ + y = t^2 e−t^ sin(t).
i. (By hand, not with Maple.) Write down the correct form of the guess that you would use to find a particular solution to this problem by using the method of undetermined coefficients. Do not solve for the coefficients! Just write down the final guess–with all its undetermined coefficients–that is sure to work. Remember that this means you must find yh to determine if your guess will work. ii. Use the dsolve command in Maple to find the general solution. Look at the solution carefully, and indicate which terms in the solution belong to yh and which belong to Yp. (You can print the Maple session and make notes on the print-out that you hand in.) (b) Use Maple to solve the initial value problem
y′′^ + y′^ + 8y = te−^2 t^ cos(t), y(0) = 1, y′(0) = − 1 ,
and plot y(t) for − 1 < t < 8.
(a) Show that by changing to the new function u defined by
y(τ ) =
c k u (ω 0 τ ) , where ω 0 =
k m
and then changing the independent variable from τ to a new variable t = ω 0 τ , we can rewrite the differential equation (1) as u′′(t) + bu′(t) + u(t) = w(t) (2) where b =
γ √ mk
, and w(t) = v(t/ω 0 ) = v(τ ).
(By choosing our coordinates in this way, we see that the problem is greatly simplified; there is really just one parameter that is important, not four.) We will use (2) for the rest of this problem. Hint. If y(τ ) = ck u (ω 0 τ ), then
dy dτ
d dτ
( (^) c k
u (ω 0 τ )
cω 0 k
u′(ω 0 τ ) =
c √ mk
u′(t).
(b) Find uh(t), the solution to the homogeneous equation associated with (2). To give your answer, you will have to find bc, the value of b for which the homogeneous system is critically damped. Then there are three cases to consider: b < bc, b = bc and b > bc. (c) Suppose that for t < 0, there is no wind, and the mass is at rest. At t = 0, the wind speed w(t) suddenly jumps to 1 and remains there. Solve the corresponding initial value problem for u(t) for three cases: b = bc/3, b = bc and b = 3bc (where bc is the value found in the previous part). Plot all three solutions on one set of axes for 0 < t < 20. Label each curve with the appropriate value of b. What is the steady-state behavior of u(t) in each case? (d) It should be clear from your plots in the previous part that if an anemometer displays u(t) as an approximation to the wind speed w(t), there will be a large error at first, but the error will approach zero asymptotically. Let’s suppose that an error of 5% is “acceptable”. The goal of this part of the problem is to design the anemometer so that the duration of the transient error is as short as possible. Since there is only one parameter to consider, design means find an appropriate value of b. Define T to be the time t at which u(t) is within 5% of the correct steady-state value for all t > T. The plots of the previous part should make clear that T depends on the parameter b. The following plot shows an example of T in an underdamped case. Note that u(t) crosses the “acceptable error” region several times, but it is only after t = T ≈ 10 .8 that the error remains within 5%. (The horizontal lines are u = 0.95 and u = 1.05.)
T 0
0.
0.
0.
0.
1
1.
1.
u
2 4 6 8 10 12 14 16 18 20 t Find the value of b (to three significant digits^1 ) that gives the smallest value of T for the initial value problem above. You can do this analytically, or by trial and error. Plot the solution for 0 < t < 20.
(^1) Remember that numbers such as 7.89 and 1. 23 × 10 − (^2) have three significant digits, but 0.003 has only one significant digit.