代写Math 466 Spring 2024 Assignment 9代写Matlab语言

Math 466

Spring 2024

Assignment 9

1.   Consider the linear input output system:

where ∈ Rn,  ∈ Rp , and  ∈ Rm. It is desired to input signals of the form.   = −K +  where the

feedback term −K serves to stabilize the system so that resonance can be averted.  The term  represents an external disturbance or forcing function. Unfortunately, however, we are unable to measure the entire    state, , and can only observe  as given by the output equation  = C .    To remedy this, we design an

observer system the state of which approximates  .  We then input a signal of the form.

and attempt to design the feedback gain Kin such away that we obtain the desired result.  The observer system takes the form.

where the matrix L is an appropriately chosen gain.

a.   What are the dimensions of the design matrices Kand L?

b.   Let   =  − , and use (1)-(3), (5),(6) to show that  satisfies the homogeneous linear system given by

c.   Give a design condition on the gain L that will guarantee that lim (t) = 0 (and therefore that t→∞

t(l)(t) − (t)) = 0.

d.   Now substitute the input given by (4) with (t) = 0(⃗⃗) into (1) and (5) and use (1)-(8) to find a 2n

dimensional first order homogeneous linear system describing the evolution of the 2n

dimensional  vector function given by

e.   Give a design condition on the gain K, that when used in conjunction with the design criteria

given in part (c) will guarantee that lim (t) = 0 and therefore the design objective has been met. t→∞

(Hint: For a block upper triangular matrix of the form. T  =  [ 0(u)    v(w)] , it follows that

det T = detudetv .   If you do not see how this hint is helpful, then you may have made an error in parts ((a) –(d)).

f.    In the case of the undamped simple harmonic oscillator with both mass and stiffness set equal to one (i.e. m = k = 1) and the ability to measure only displacement and not velocity, provide

explicit conditions on the elements in the gains L and kthat will ensure that the conditions you set out in parts (c) and (e) are satisfied.

g.   Use Simulink to demonstrate that your design in part (f) works by choosing values for the gains L and kthat satisfy the conditions you set out in part (f) for the cases (i) v(t) = 0, and (ii) v(t) =    sint.

2. A United States Navy Contractor is designing the next generation of sono-buoy.  A sono-buoy

is a device used to track enemy submarines.  It is dropped into the ocean from low flying aircraft  and sends out low frequency sound waves which travel through the water and bounce off of large metallic objects such as a submarine.  The new generation of sono-buoy is cylindrical in shape

with circular bases and floats on the surface of the ocean with the circular bases at the top and

bottom.  See the figure on the next page. The radius of the circular bases is rand the height of the cylinder ish.  The average mass density of the sono-buoy isp and the mass density of seawater   is pw.

a.   What is the mass m of the sono-buoy in terms of p,T, and ℎ.

b.   The buoyancy (which is a force!) of an object is equal in magnitude to the weight of  water being displaced.  In terms of l, at equilibrium with gravity, what is the buoyancy?

c.   Based on part (b), what must lbe equal to?

d.   In order to be able to accurately process the data it receives from the sono-buoy, the Navy requires a mathematical model for the vertical motion of the device as it bobs up and

down on the surface of the ocean. Let u(t) denote any external forces acting on the

device at time t ≥ 0, and use Newton’s second law to derive a differential equation for the quantity measured in reference to the coordinate system shown in the figure.  Write your model as a first order linear input/output system.

e.   Using either Laplace transforms or the matrix exponential and the variation of constant formula, provide an expression for the output, y, in terms of the input u.

f.    If incoming ocean waves aremodeled in the form of a cosine function of the form,

u(t) = R cos幼t

where 幼 is never larger than some frequency 幼max,  in designing the sono-buoy, what values of the sono-buoy dimensions,h, r and p will the Navy engineers want to avoid and why?

 

g.   If in their design, the Navy decides to include a damping device that can be modeled via a term that is proportional to the vertical velocity of the sono-buoy, with a positive constant of proportionality b, where the value of b depends on the design of the damper, how

should the damper be designed (i.e. how should the gain b be chosen) so as to achieve critical damping?

h.   Choose reasonable values for p,T, and ℎ, lookup pw, and simulate your model using

Simulink for different values of R, b and 幼.  In your simulations, be certain to include    examples where you choose values for 幼 that cause the system response to exhibit beats and resonance.  Be certain to pay careful attention to units and make conversions where  necessary.

3.    Why are marching soldiers told to “break step” when they cross a bridge?  In this problem we will see how a

marching column of soldiers can de-stabilize a flexible bridge.   Model the bridge as an undamped simple harmonic oscillator with mass m and stiffness k, assume that the bridge is initially at rest when the soldiers begin marching across it, and model the marching soldiers as a sequence of pulses or impulses.   That is, set

where δ(t) is the unit impulse (Dirac delta function) at time t  = 0 (and therefore its translate, δ(t − c), is the unit impulse at time t = c) and τ > 0 is a constant that represents the time it takes the soldiers to take a single step.  Show that for an appropriate choice of τ (that of course will depend on the values of mandk) a resonance-like situation will result and the amplitude of the oscillations of the bridge will grow without bound.  (Hint: Use the impulse response function and the convolution form. of the solution you found as part of Problem 3 in Assignment 4, the integration properties of the Dirac delta function, and the fact that the system being linear with zero initial conditions, implies that the response of the system to a sum of inputs, will be the sum of the responses of the system to each of the individualsummands that make up the input. Or equivalently that integration is a linear operation.  It may also help to realize that u can be written as

where uc (t) is the unit step function at t  = c; uc (t) = 1, fort ≥ c  anduc (t) = 0, fort < c .)

4.    Use Simulink to simulate your results in Problem 1 for the case where a) τ is such that resonance occurs, b) τ is such that resonance does not occur, and c) the soldiers have “broken step” and are casually walking across the  bridge  (this  last one  will  require your  modeling skills  and  possibly  Matlab’s  random  number generator! )