MATLAB ode45: How To Solve a System of Ordinary Differential Equations (ODE

No video

MATLAB ode45: How To Solve a System of Ordinary Differential Equations (ODE - with discrete data)

Рет қаралды 44,593

VDEngineering

Күн бұрын

Пікірлер: 39

@DEFIB7 10 ай бұрын

6 years after posting this and you saved my LIFE on my final lab for my mechatronics class, thank god man! you are the BEST

@Mtmonaghan 3 ай бұрын

Top class work, well done and thank you.

@liamgerritsen6040 2 жыл бұрын

you have no idea how much you've helped me, thank you so so much. Extremely well done video

@Bu.T.F.I 2 жыл бұрын

Well done! SSooooo useful, and I understood it! Very much appreciated.

@josipsacher 2 жыл бұрын

Thank you. You can also define constants as global variables (no need to pass it as an argument)

@funnysele 6 жыл бұрын

Thank you so much for this! It's very well explained :)

@frederickajustagosto4765 5 жыл бұрын

great video thanks for posting. covers many topics that are encountered.

@susramu 6 жыл бұрын

Many many many thanks you helped me alot... again thanks to you.

@mikal_1 5 жыл бұрын

thanks brother

@felipebm13 6 жыл бұрын

Hi again VDE, could you tell me if there is any way of setting a minimum "ground" value for y(2), such as 0? If there is no minimum value and the initial force values are less than the ones of gravity, it means that the rocket would go down instead of simply not taking off.

@madhuryahazarika3266 4 жыл бұрын

Please discuss a problem using keller box method with MATLAB code

@danielfrantz1815 2 жыл бұрын

Hi From MAE 484

@joseanchieta1706 4 жыл бұрын

Good night, I have a doubt. If the FT were dependent on y (1), in this case, the displacement, what should change and how to change?

@funnyjone9473 5 жыл бұрын

Thanks Brother

@freddycaban3306 5 жыл бұрын

Was this solved numerically? I'm asking becauseyour using quadratic drag in your equation which gives V^2=Vx^2+Vy^2 and will lead to a coupled system of non-linear differential equations.

@robin9740 4 жыл бұрын

ODE45 solves ODE's according to an explicit Runge-Kutta order 4. It is a numerical way of solving ODE's indeed.

@tarekmahmoudkhalifa8067 4 жыл бұрын

Can this code to solve five equation coupled nonlinear ??

@ajinkyadalvi4761 4 жыл бұрын

hello sir, My name is ajinkya dalvi, I have acceleration data in excel format and i want to convert that data into time displacement data and use to solve a ODE. so can you please guide me in that.

@ThomasHaberkorn 4 жыл бұрын

Great video, but how to solve a higher order system of ODEs?

@Tasty_triangles Жыл бұрын

same query please reply

@sudhansumtripathy 7 жыл бұрын

do we have any built in tool in matlab , simulink, simscape or any thing where I can build my own propeller and simulate it.

@VDEngineering 7 жыл бұрын

Yes, it is called SimMechanics Link

@chenzofunari2924 6 жыл бұрын

If thrust and density weren't constant but rather functions of altitude, how would I implement that?

@VDEngineering 6 жыл бұрын

see this: kzfaq.info/get/bejne/m7J2qqR9vKfGqGw.html You can also create additional excel columns with thrust and density as a function of height.

@chenzofunari2924 6 жыл бұрын

Thanks for the response. I am still having trouble relating the functions of heights back to the functions of time. For example, my problem started off by calculating thrust at every altitude by relating sea level and altitude thrust coefficients using values from the standard atmosphere model. I could then calculate mass flow which is a function of time. That makes my dv/dt = (T(h) - D(h,v))/m(t) - g*sin(gamma). My time vector and altitude vectors are different lengths and I am not sure how to normalize them to get altitude as a function of time.

@VDEngineering 6 жыл бұрын

The easiest thing would be to truncate the time vector (or extend it), Maybe use an add array operation

@chenzofunari2924 6 жыл бұрын

I was able to clean up the discrepancy of lengths. However, after think more about it, I don’t think you can use Euler’s method like this for my specific problem. In the for loop, you will have T(i) in the numerator and m(i) in the denominator for the dv/vt equation. With thrust based on altitude in my case, you can’t just simply have the same indices because T(1) will be the thrust at one height step used and m(1) will be the mass at the time step used. You cannot assume that after the one time step, the vehicle will be at the altitude assumed after one height step. This assumes a constant velocity of delta-h/delta-t with the changes being your steps chosen.

@felipebm13 7 жыл бұрын

Hey VDE, how can I define a time value as a parameter and find the y value at that time? Thanks for the video, very nice

@VDEngineering 7 жыл бұрын

Your y output is an array so use array ops

@felipebm13 7 жыл бұрын

What is the command I should type? I am new to Matlab and I want to find the velocity at the end of the burn so I can calculate the total height of my rocket

@felipebm13 7 жыл бұрын

To be more clear, I want to know the total altitude of the rocket, and since my excel file shows the time until the burnout only, I can't see past the altitude when the burnout happens. What should I do? I have tried to expand the timespan but it did not work... Is there any tool in which I can input a value of y'(0) and it gives me y(t) at that moment? Thanks again for the amazing videos

@VDEngineering 7 жыл бұрын

X = y(:,1) h = length(X) X(h) try this, if that doesn't work, just manipulate it

@felipebm13 7 жыл бұрын

Thank you! it worked.

@nemanipatil8370 3 жыл бұрын

Sir I am begginred please help me matlab program

@mayraalvarez5879 3 жыл бұрын

Hi! Is it possible to have your email to contact you for a doubt on this topic? It is a code that I am doing and it is running for too long. If you could help me I can join Patreon in exchange or alike. I hope you can help me! All the best

@VDEngineering 3 жыл бұрын

Msg on Instagram! sorry just saw this, I don't get comment notifications lol