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

Поделиться
HTML-код
  • Опубликовано: 3 окт 2024

Комментарии • 49

  • @JousefM
    @JousefM  4 года назад +8

    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
      @sriti_hikari 4 года назад +2

      Good.

    • @zenianggara4958
      @zenianggara4958 4 года назад +1

      thanks

    • @zenianggara4958
      @zenianggara4958 4 года назад

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

    • @JousefM
      @JousefM  4 года назад

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

    • @zenianggara4958
      @zenianggara4958 4 года назад

      @@JousefM no. name function is ode45.

  • @zofiaorlowski5451
    @zofiaorlowski5451 3 года назад

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

    • @JousefM
      @JousefM  3 года назад

      Thank you so much!!

  • @brianconway9085
    @brianconway9085 2 года назад +1

    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.

  • @MachineLearningwithPhil
    @MachineLearningwithPhil 4 года назад +4

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

    • @JousefM
      @JousefM  4 года назад

      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 🙂😉

  • @ThijsNL1
    @ThijsNL1 3 года назад +2

    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?

  • @techlead3743
    @techlead3743 4 года назад +2

    Nice!

    • @luizfelipecarvalho1469
      @luizfelipecarvalho1469 4 года назад

      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...]

  • @vargalaszlo1989
    @vargalaszlo1989 3 года назад

    Weekly Regelungstechnik practice when the entire class was shaking with fear

    • @JousefM
      @JousefM  3 года назад

      I feel you 😄

  • @simranjeetsinghdhillon6571
    @simranjeetsinghdhillon6571 3 года назад +1

    How can I plot for acceleration vs time

  • @andreitanasache5019
    @andreitanasache5019 4 года назад

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

  • @shanmukharao1683
    @shanmukharao1683 4 года назад

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

    • @JousefM
      @JousefM  3 года назад

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

  • @OatCmor
    @OatCmor 3 года назад

    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
      @JousefM  3 года назад +1

      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!

  • @pabloalmagro5181
    @pabloalmagro5181 4 года назад

    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
      @JousefM  4 года назад +1

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

    • @luizfelipecarvalho1469
      @luizfelipecarvalho1469 4 года назад

      @@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?

  • @nikamirfirdaus2516
    @nikamirfirdaus2516 4 года назад

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

    • @JousefM
      @JousefM  4 года назад

      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!

  • @shanmukharao1683
    @shanmukharao1683 4 года назад

    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
      @JousefM  4 года назад

      Maybe in the future! Stay tuned :)

    • @shanmukharao1683
      @shanmukharao1683 4 года назад

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

  • @shanmukharao1683
    @shanmukharao1683 3 года назад

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

    • @JousefM
      @JousefM  3 года назад

      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
      @shanmukharao1683 3 года назад

      @@JousefM thanku sir, i will try this

  • @kid-vf4lu
    @kid-vf4lu 3 года назад

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

    • @JousefM
      @JousefM  3 года назад

      Yes, this is Runge-Kutta!

  • @philippebergna5783
    @philippebergna5783 4 года назад

    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
      @JousefM  4 года назад

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

    • @veereshguruguntla5204
      @veereshguruguntla5204 4 года назад

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

  • @malikzeeshanawan3542
    @malikzeeshanawan3542 4 года назад

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

    • @JousefM
      @JousefM  4 года назад +1

      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
      @malikzeeshanawan3542 4 года назад

      @@JousefM thank you

  • @MS-ld3jn
    @MS-ld3jn 4 года назад

    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
      @JousefM  4 года назад

      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
      @MS-ld3jn 4 года назад

      @@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!

  • @olakappassov6788
    @olakappassov6788 3 года назад

    hi. where did u take initial conditions?

    • @JousefM
      @JousefM  3 года назад

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