Matlab Training Session Iv Simulating Dynamic Systems: Sampling The Solution Equation

Simulating Dynamic Systems
MATLAB TRAINING SESSION IV SIMULATING DYNAMIC SYSTEMS

MATLAB provides a variety of techniques for solving dynamics problems numerically. On-line help for each of the commands shown below is available by typing help command_name or doc command_name at the MATLAB prompt.
Sampling the Solution Equation:

Suppose you already have the solution to a dynamics problem and you want to obtain a plot of the response. The task would reduce to implementing the equation in MATLAB and sampling the function over some domain of interest. For a high delity reproduction of the function, the function must be sampled with a spacial or time increment that captures many samples for each frequency that comprises the Fourier series for the function. In simpler terms, we can tell if we have used a small enough increment if the resulting plot appears smooth. The example shows how the same function can look jagged or smooth in the plot depending on the sampling increment.
1 of 9
Suppose we nd:
y(t ) = 8 5e
2 t
( sin ( 2 t ) + 0.5 cos ( 2 t ) )
Linear Systems:
There are three general purpose methods for solving dynamics problems; transfer function method, state space integration method, and the MATLAB function lsim. Transfer Function Method: The transfer function method uses partial fraction expansion of the transfer function multiplied by the Laplace transform of the forcing function to obtain the solution. The solution of the differential equation is expressed in terms of a linear combonation of complex exponentials that are determined by the poles and the corresponding residues in the partial fraction expansion. Thus the form of the solution is: y(t ) =
r ( i ) exp ( p ( i ) t )
i=1
Suppose we have Then we can write
y + 2 y + 5 y = 5u(t ) + u ( t )
and
u ( t ) = input
2 of 9
s Y ( s ) sy ( 0 ) y ( 0 ) + 2 sY ( s ) 2 y ( 0 ) + 5 Y ( s ) = 5 U ( s ) + sU ( s ) + u ( 0 ) Let y(0) = 0 y ( 0 ) = 0 u(0) = 0 s+5 -U (s) Then Y ( s ) = ------------------------2 s + 2s + 5 Let the input u(t) be the step function and to solve for y(t) do With 1 U ( s ) = -s then s+5 Y ( s ) = --------------------------------2 s(s + 2s + 5)
Let u(t) be sin(2t) and to solve for y(t) do
2 U ( s ) = ------------2 s +4
3 of 9
Integration Method: The integration method utilizes the state space representation of the system and the MATLAB ode23 or the ode45 integrator. One advantage of this method is that it is much easier to incorporate the initial conditions in the solution without having to nd the new transfer function. Lets look at the last example given above. x ( t ) = 0 1 x ( t ) + 0 u ( t ) x = Ax + Bu 5 2 1 The state space representation is y = 5 1 x(t ) To implement the model create the M-le eqn1.m as shown y = Cx + Du
4 of 9
Then at the MATLAB prompt, enter:
Special Purpose MATLAB Commands: The MATLAB control systems toolbox provides several built in functions for solving for the solution to linear dynamic equations. All functions except either the transfer function or the state space representations. When the functions are invoked with left hand arguments, the variable are lled and a no plot is made. Each function has a variable number of possible valid argument congurations. Most function support both transfer function (num/den) and state space (A,B,C,D) forms. More information about each function can be found by typing help function_name. A list is provided below.
initial(A,B,C,D,x0) impulse(num,den,t)
impulse(A,B,C,D,IU,t) step(num,den,t) step(A,B,C,D,IU,t) lsim(num,den,U,t) lsim(A,B,C,D,U,t,x0)
Initial condition response of cont-time LTI systems Impulse response of continuous time LTI systems IU = input channel number; t = sample time vector Step response of continuous time LTI systems Simulation of c-time LTI systems to arbitrary inputs U = input time history (# of columns = # of inputs)
5 of 9
For example lets work with the series resonant circuit: di 1 - i d = V s Ri + L + --dt C
0 t
THEN
2
LC
di di dV + RC + i = C dt dt dt I (s) Cs ---------- = ------------------------------------------2 V (s) ( LCs + RCs + 1 ) C
Vs
( RCs + LCs + 1 ) I ( s ) = CsV ( s )
OR
x1 = i
x2 = x 1
x 1 x 2
0 1 x 1 = + 0 u 1 R ------- -- x 1 LC L 2
1 x1 y = 0 -L x2
With Matlab:
There is a corresponding function for discrete time dynamics for the functions mentioned above . The function would start with a d like dstep for example. Suppose we nd the z-transform to be Then 0.4673 z 0.3393 G ( z ) = ---------------------------------------------------2 z 1.5327 z + 0.6607
6 of 9
Nonlinear Systems: Nonlinear dynamics are solved by the integration method. The system equations are expressed in state form described in a M-le function and the ode integrators are used to solve for the state trajectories. 1 - [ x1 + x2 g ( x2 ) ] x 1 = -----------C 1 R1 Consider the Wien-Bridge Oscillator: 1 1 - [ x 1 + x 2 g ( x 2 ) ] ------------x x 2 = -----------C 2 R1 C 2 R2 2 g ( v ) = 3.234 v 2.195 v + 0.666 v Create the M-le wien.m as shown (Let R1 = R2 = C1 = C2 = 1.0)
3 5
Then with MATLAB
7 of 9
Consider the Cart and Pendulum Problem 2 sin = 0 (M + m)x + ml cos ml sin + mgl sin = 0 ml + mlx cos mlx Use a change of variables to y then we can write y1 x(t ) y x ( t ) = 2 (t ) y3 (t) y4 Write the M-le cart.m 1 0 0 0 0 ( M + m ) 0 ml cos y 3 0 0 1 0 ml
2 2
y 1 y 2 y 3 y 4 =
y2 mly 4 sin y 3 y4 ml ( y 2 y 4 g ) sin y 3

2
0 ml cos y 3 0
Then with MATLAB
8 of 9
9 of 9

Matlab Training Session Iv Simulating Dynamic Systems: Sampling The Solution Equation

Uploaded by

Copyright:

Available Formats

Matlab Training Session Iv Simulating Dynamic Systems: Sampling The Solution Equation

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Matlab Training Session Iv Simulating Dynamic Systems: Sampling The Solution Equation

Uploaded by

Copyright:

Available Formats

Simulating Dynamic Systems