Solving 2nd Order Equations with Maple: Analyzing Anemometer & Economics Models, Assignments of Differential Equations

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

Pre 2010

Uploaded on 08/19/2009

koofers-user-pg7
koofers-user-pg7 🇺🇸

10 documents

1 / 3

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Math 308 - Differential Equations Fall 2003
Homework Assignment 2
Due Wednesday, October 29.
1. Using Maple for Second Order Equations. Work through the examples shown in the handout
“Solving Second Order Differential Equations”, and then use Maple to solve the following problems.
(You should review the earlier Maple handouts, too.)
(a) Consider the differential equation
y00 + 3y0+y=t2etsin(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 yhto 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 yhand 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
y00 +y0+ 8y=te2tcos(t), y(0) = 1, y0(0) = 1,
and plot y(t) for 1< t < 8.
2. A Deflection Anemometer. Consider a spring-mass system as shown in Figure 3.8.10 on page 200
of the text (but, for the moment, we use τfor time and yfor the displacement). The mass is m, the
spring constant is k, the damping coefficient is γ, and y(τ) is the displacement of the mass at time
τ. Assume that wind blows horizontally, parallel to the direction of motion of the mass, and that the
force exerted on the mass by the wind is proportional to the wind speed v(τ). Then Newton’s Second
Law gives
my00(τ) + γ y0(τ) + ky(τ) = cv(τ),(1)
where cis the proportionality constant that relates the wind speed to the force.
(a) Show that by changing to the new function udefined by
y(τ) = c
ku(ω0τ),where ω0=rk
m,
and then changing the independent variable from τto a new variable t=ω0τ, we can rewrite the
differential equation (1) as
u00(t) + bu0(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(τ) = c
ku(ω0τ), then
dy
=d
³c
ku(ω0τ)´=0
ku0(ω0τ) = c
mk u0(t).
1
pf3

Partial preview of the text

Download Solving 2nd Order Equations with Maple: Analyzing Anemometer & Economics Models and more Assignments Differential Equations in PDF only on Docsity!

Math 308 - Differential Equations Fall 2003

Homework Assignment 2 Due Wednesday, October 29.

  1. Using Maple for Second Order Equations. Work through the examples shown in the handout “Solving Second Order Differential Equations”, and then use Maple to solve the following problems. (You should review the earlier Maple handouts, too.)

(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.

  1. A Deflection Anemometer. Consider a spring-mass system as shown in Figure 3.8.10 on page 200 of the text (but, for the moment, we use τ for time and y for the displacement). The mass is m, the spring constant is k, the damping coefficient is γ, and y(τ ) is the displacement of the mass at time τ. Assume that wind blows horizontally, parallel to the direction of motion of the mass, and that the force exerted on the mass by the wind is proportional to the wind speed v(τ ). Then Newton’s Second Law gives my′′(τ ) + γy′(τ ) + ky(τ ) = cv(τ ), (1) where c is the proportionality constant that relates the wind speed to the force.

(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.