SolutionofTypicalmck-SystemODEviaMatlab.docVIP

Solution of Typical mck-System ODE via Matlab Last changed on 11 June 2007 A typical, unforced mk-system has the ODE: This can be algebraically manipulated by solving for the highest order derivative: Now we introduce two new (state) variables, x1 = x and x2 = dx1/dt = dx/dt. So the second-order equation can be rewritten as two first-order equations: This is the state-space representation of the second order ODE. Note that the system has two states, and the second-order ODE has been converted into a system of first-order ODEs. In this form the system can be input into a Matlab ODE solver to get the motion, i.e. x = x(t). Here’s how. The Matlab solution requires two steps: One writes a Matlab function in an m-file One invokes or calls this function from within Matlab. The function that contains the first-order ODEs should have a meaningful name, like mckSystem, for example. Then the m-file would be named mckSystem.m. This file would contain: function dx = mckSystem(t,x) m = 3.0; c = 6.93; k = 100.0; dx(1,1) = x(2); dx(2,1) = –(c/m)*x(2) - (k/m)*x(1); This describes the system and includes specific values of m, c, and k. To invoke this function from Matlab, one has to first define two things: The time that the resulting simulation should span The initial conditions for x1 and x2, i.e. x(0) and v(0) Once these are defined, the user can invoke the function from the Matlab command window to produce the simulation of the motion. tspan = [0,10]; x0 = [1.5,0]; [t,x] = ode45( ‘mckSystem’, tspan, x0 ); After these commands are issued, the results are stored in the vectors t and x. t is an array of time instances. x is an array of locations (x(*,1)) and velocities (x(*,2)) corresponding to each time instance. So, for example, if t(28) = 0.73 and x(28,1) = -0.043, that means that at t = 0.73 sec, the mass is located at -0.043 units. The units are inches, millimeters, or whatever uinits you decide to use in the problem. Though th