Question

Problem 8 [5 points The following system of ODEs, formulated by Lorenz, represents are crude model of atmospheric circulation

0 0
Add a comment Improve this question Transcribed image text
Answer #1

In MATLAB:

function yp = programa (t,y)

y1=15;

y2=15;

y3=36

yp1=10*(y2-y1);

yp2=28*y1-y2-y1*y3;

yp3=y1*y2-(8/3)*y3;

yp=[yp1;yp2;yp3];

create new script with the following:

[t,y]= ode45(@programa, [0,100], [1,4]); #(change the initial values)

plot (t,y (:,1), t,y(:,2)); #(change for (y1,y3) and (y2,y3))

grid on

title (exercise);

xlabel(´t´);

xlabel(´y´);

Add a comment
Know the answer?
Add Answer to:
Problem 8 [5 points The following system of ODEs, formulated by Lorenz, represents are crude mode...
Your Answer:

Post as a guest

Your Name:

What's your source?

Earn Coins

Coins can be redeemed for fabulous gifts.

Not the answer you're looking for? Ask your own homework help question. Our experts will answer your question WITHIN MINUTES for Free.
Similar Homework Help Questions
  • Q5 (30 pts). A well-studied first order system of ODEs is given by which are called the Lorenz eq...

    Q5 (30 pts). A well-studied first order system of ODEs is given by which are called the Lorenz equations, and they were derived by the Meteorologist Edward Lorenz in the early 1960's as a simple model of weather (more precisely as a model of a pair of coupled convection cells in the atmosphere). In the model o, β, ρ > 0 are positive real parameters. Fix the classical parameter values to be σ = 10, β-8/3 and ρ 28. Fix...

  • Use MATLAB’s ode45 command to solve the following non linear 2nd order ODE: y'' = −y'...

    Use MATLAB’s ode45 command to solve the following non linear 2nd order ODE: y'' = −y' + sin(ty) where the derivatives are with respect to time. The initial conditions are y(0) = 1 and y ' (0) = 0. Include your MATLAB code and correctly labelled plot (for 0 ≤ t ≤ 30). Describe the behaviour of the solution. Under certain conditions the following system of ODEs models fluid turbulence over time: dx / dt = σ(y − x) dy...

  • Problem Set A Problem 6. (20%) A ordinary differential equation for a mass-damper-spring system is following....

    Problem Set A Problem 6. (20%) A ordinary differential equation for a mass-damper-spring system is following. The mass m 1, damping coetfic e initial position y(o) O, and the initial velocity i constant k 3 and force 10, all are in appropriate units. Th 1, spring zero, within the time range of O to 20 unit of time, use Matlab find the solution of function y(t)? Hint: you need to convert the 2nd order ODE into two 1st order ODEs....

  • Instructions: For your lab write-up follow the instructions of LAB 1 1. (a) Modify the function e...

    Instructions: For your lab write-up follow the instructions of LAB 1 1. (a) Modify the function exvith2eqs to solve the IVP (4) for o St S40 using the MATLAB routine ode45. Call the new function LAB04exl Let [t,Y] (note the upper case Y) be the output of ode45 and y and v the unknown functions. Use the following commands to define the ODE function dYdt f(t.Y) y-Y();Y(2) Plot y(t) and e(t) in the same window (do not use subplot), and...

  • t0 = 0; tf = 20; y0 = [10;60]; a = .8; b = .01; c...

    t0 = 0; tf = 20; y0 = [10;60]; a = .8; b = .01; c = .6; d = .1; [t,y] = ode45(@f,[t0,tf],y0,[],a,b,c,d); u1 = y(:,1); u2 = y(:,2); % y in output has 2 columns corresponding to u1 and u2 figure(1); subplot(2,1,1); plot(t,u1,'b-+'); ylabel('u1'); subplot(2,1,2); plot(t,u2,'ro-'); ylabel('u2'); figure(2) plot(u1,u2); axis square; xlabel('u_1'); ylabel('u_2'); % plot the phase plot %---------------------------------------------------------------------- function dydt = f(t,y,a,b,c,d) u1 = y(1); u2 = y(2); dydt = [ a*u1-b*u1*u2 ; -c*u2+d*u1*u2 ]; end Only...

ADVERTISEMENT
Free Homework Help App
Download From Google Play
Scan Your Homework
to Get Instant Free Answers
Need Online Homework Help?
Ask a Question
Get Answers For Free
Most questions answered within 3 hours.
ADVERTISEMENT
ADVERTISEMENT
ADVERTISEMENT