MATLAB's ode45 Solver - Single Degree-of-Freedom Oscillator

📝 Blog Post: www.jousefmurad.com/engineeri...
📈 APEX Consulting: theapexconsulting.com
🌎 Website: jousefmurad.com
In this video, we will have a look at MATLAB's ode45 solver and how to solve a 2nd-order differential equation for a single degree of freedom oscillator.
ONLINE PRESENCE
================
🌍 My website - jousefmurad.com/
💌 My weekly science newsletter - jousef.substack.com/
📸 Instagram - / jousefmrd
🐦 Twitter - / jousefm2
SUPPORT MY WORK
=================
👕 Science Merch: engineered-mind.creator-sprin...
👉 Become my Patreon: / theengiineer
🧠 Subscribe for more free videos: bit.ly/2RLmMxq
RECOMMENDATIONS
==================
🎥 1 Month Free of Skillshare Premium: bit.ly/3E3bXwC
📕 Books on AI, Data Science & Programming: bit.ly/2I1lUD1 - use "podengineered20" to get a 40% discount!
🎵 Tomsize & Simeon - Jump (Fetre Retwerk)
Contact:
If you need help or have any questions or want to collaborate feel free to reach out to me via email: [e-mail]
#MATLAB
#ODE45
#DifferentialEquation

[@JousefM]

Here you can find all the code: github.com/jousefm/ModelingAndSimulation - let me know if this video helped you and what application you are solving with ode45! 🎓 My Science Courses - courses.jousefmurad.com/ ✍️ Latest blog posts: jousefmurad.com 📥 My Newsletter - jousef.substack.com/

[@sriti_hikari]

Good.

[@zenianggara4958]

thanks

[@zenianggara4958]

sorry, function "[T,Y] = ode45(@StateSpaceForm, Time, IC, SolverOptions, m, d, c, f);" is not working. why ?

[@JousefM]

Is your function even called StateSpaceForm? Also name the file for the function the same!

[@zenianggara4958]

@@JousefM no. name function is ode45.

[@zofiaorlowski5451]

This video and the blog post instructions were so helpful and insightful. Thank you for making them!

[@JousefM]

Thank you so much!!

[@brianconway9085]

Great stuff, please do a quarter car and half car model. it would be amazing to find out using an if statement can you for loop and find the correct damping for a critically damped system.

[@ThijsNL1]

Thanks for the video, first time I managed to create such a thing in Matlab. But how can I create a sine wave as force input, so that f=sin(t) for example?

[@MachineLearningwithPhil]

Channel's blowing up! Congrats man. You posting these to reddit?

[@JousefM]

Thanks mate! Yes I do but the potential of being banned is quite high, got banned in 2 subreddits. If I were you I would give it a try in the ML/DL/RL subreddits. I will give a shoutout to your channel once I have the time to create a beginner‘s guide to machine learning/RL! Or you create one beforehand yourself 🙂😉

[@andreitanasache5019]

Niceee, how can i extrapolate the last value of ode45? like as y(end,1) ??? thanks

[@techlead3743]

Nice!

[@luizfelipecarvalho1469]

Nice work! I Wonder if it is possible to use external forces changing in time, but randomly like in time series? e.g. F(t) = [10 3 -10...]

[@pabloalmagro5181]

Nice work!!! In the case of a harmonic displacement for the same mass-spring-damper SDOF model where the goal would be to find the generated force due to the harmonic displacement. Do you have some recommendations to carry out such a case?

[@JousefM]

Does that help? www.mathworks.com/help/symbolic/examples/physics-damped-harmonic-oscillator.html :)

[@luizfelipecarvalho1469]

@@JousefM Nice work!! I Wonder if it is possible to have an force varying in time (time series, not described for a fubction). Is it possible? How would it be in code?

[@simranjeetsinghdhillon6571]

How can I plot for acceleration vs time

[@nikamirfirdaus2516]

Hi, thanks for this video! Can you do a system with 2 masses? I would really appreciate it!

[@JousefM]

Thanks for your comment! Sure, I can :) currently working on something else but feel free to subscribe and activate the notification bell to stay up-to-date with the latest uploads Nik!

[@MS-ld3jn]

Hey man thanks for the video! Also, would you know how to plot the force transmissibility of the 1 DOF system on matlab? Thanks.

[@JousefM]

You're welcome! This should help: vortarus.com/creating-a-vibration-transmissibility-plot/ - the formulas are straightforward, simply put them into MATLAB and let it do the job, :)

[@MS-ld3jn]

@@JousefM Thanks so much but I tried plotting damped force transmissibility but it didn't work :-{. But I think i found a solution from a textbook I found. Thanks anyways!

[@philippebergna5783]

I have a similar question to this one, however, I need to find the force (force * acceleration), so I need to find the acceleration of the mass, how would I be able to do it?

[@JousefM]

You can plot and calculate the acceleration as shown in the video. If you only have the velocity take the derivative

[@veereshguruguntla5204]

hi i also want to draw acceleration vs time and R.M.S acceleration and frequency plots ..can you help me out

[@shanmukharao1683]

thanku very much for your videos , if mass matrix is having 2x2 size, how can apply this code?? can u pls explain this

[@JousefM]

Just give it a try with 2x2 and see what error it throws. Always use the workspace to see what dimensions you’re working with and if these even make sense. „Debugging“ using the workspace is very helpful sometimes.

[@shanmukharao1683]

@@JousefM thanku sir, i will try this

[@OatCmor]

Hi. Thanks for your video. However, I visited your website "engineered mind". It says Page is under construction - great stuff coming soon!. When will it be back online?

[@JousefM]

Hey and thanks for asking :) will slowly build it in the coming weeks and make it more „bloggy“ to educate more and make it more user friendly. Stay tuned!

[@shanmukharao1683]

Hai sir, video is really good.... ca u pls explain how to solve for 2-DOF systems with damped forced vibrations.... thanks in advance

[@JousefM]

Maybe in the future! Stay tuned :)

[@shanmukharao1683]

@@JousefM Thanks for your response sir, if u do within a few days really is helpful for me... thanks in advance

[@olakappassov6788]

hi. where did u take initial conditions?

[@JousefM]

Initial Conditions are indicated by „IC“ in the Video - Hope That helps! :)

[@malikzeeshanawan3542]

why did you not include the weight (mg) in your diff. eq. ? Kindly tell

[@JousefM]

Hey Malik! The mass was included as the small parameter m but not gravity as you might have seen. But you could include it of course - the equation would look slightly different :)

[@malikzeeshanawan3542]

@@JousefM thank you

[@kid-vf4lu]

Is ODE45 a runge kutta method? I know that it is adaptive, so it's not the standard runge kutta 4th

[@JousefM]

Yes, this is Runge-Kutta!

[@vargalaszlo1989]

Weekly Regelungstechnik practice when the entire class was shaking with fear

[@JousefM]

I feel you 😄

[@shanmukharao1683]

kindly do 2-dof...soon sir...it will be helpfull

[@JousefM]

Working on it at the moment, should hopefully be done soon.