MATLAB Help - Translational Orbit Dynamics for a Low Earth Satellite using ode45

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

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

  • @giorgioharmony4824
    @giorgioharmony4824 3 года назад +6

    Just the video I needed to find to finish my final job, you’re a genius, thank you very much

  • @beebee_0136
    @beebee_0136 4 года назад +3

    @Monte.Carlos
    Hello Prof,
    I have been following your series on Low Orbit satellite dynamics. And it's been great!
    Regards.

  • @TheEdmaster87
    @TheEdmaster87 2 года назад

    Hi, love your stuff. You are a very good teacher. Can you tell me where in the code I need to change to get elliptical orbit?

    • @CarlosMontalvo251
      @CarlosMontalvo251  2 года назад

      You'd have to change the initial velocity. If you increase the velocity from the circular orbit velocity it will become elliptic. Does that make sense?

    • @TheEdmaster87
      @TheEdmaster87 2 года назад

      @@CarlosMontalvo251 of course yes thank you!

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

    great job! ... i tried to do that in Simulink but something's going wrong and i don't know how to fixe it, can you help if i send to you the simulink blocks..?... thank you!!

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

    Hello Carlos,
    thank you very much for this walkthrough
    how would you change this script in order to work with elliptical orbits?
    thank you in advance

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

      You just need to increase the initial velocity. You'll notice at 15:52 I define vcircular which is the velocity for a circular orbit. If you increase that velocity the orbit will be elliptic. In order to compute velocity as a function of eccentricity you will need to refer to the Fundamentals of Astrodynamics www.amazon.com/Fundamentals-Astrodynamics-Dover-Aeronautical-Engineering/dp/0486600610

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

      @@CarlosMontalvo251 thank you for the quick reply

  • @MrEngineer_
    @MrEngineer_ 4 месяца назад

    Hi sir ! The algorithm failed when I tried to simulate it for a year, also with increasing numbers of orbits the satellite sprirals down

    • @CarlosMontalvo251
      @CarlosMontalvo251  4 месяца назад

      Did you lower the timestep? Sometimes the simulation won't converge if the timestep is too large

    • @MrEngineer_
      @MrEngineer_ 4 месяца назад

      @@CarlosMontalvo251 But ode45 uses adaptive timestep right ? Do I have the control over timestep ?

    • @MrEngineer_
      @MrEngineer_ 4 месяца назад

      Actually the error said that the timestep got so reduced, it went below permissible values and caused the algorithm to stop.

    • @CarlosMontalvo251
      @CarlosMontalvo251  4 месяца назад

      @@MrEngineer_ This simulation uses RK4 so you have to set the timestep in the main.m

    • @CarlosMontalvo251
      @CarlosMontalvo251  4 месяца назад

      My apologies. The latest version of the software uses ode45. The simulation with ode45 probably threw an error when you went too close to the center of the planet and the gravity model became too nonlinear.

  • @joshuathompson6275
    @joshuathompson6275 5 месяцев назад

    Hi. Great video! Where in the repo are these codes? Thanks!

    • @CarlosMontalvo251
      @CarlosMontalvo251  5 месяцев назад

      github.com/cmontalvo251/aerospace/tree/main/adcs_seminar_series

    • @joshuathompson6275
      @joshuathompson6275 5 месяцев назад

      @@CarlosMontalvo251 Thank you! Do you have any code that involves spacecraft detumblization?

    • @CarlosMontalvo251
      @CarlosMontalvo251  4 месяца назад

      @@joshuathompson6275 The code above detumbles with magnetorquers and reaction wheels.

  • @aeronautics2740
    @aeronautics2740 2 года назад

    I have 2 questions about this and hope if you can help me with.
    1- I'm using Python and I faced an error in the integration tool: I did like this: stateout = sci.odeint(dstatedt,stateinitial,tspan) where dstatedt is a function that I called from the satellite file because I don't know how to call satellite file to the main. so I did like this : from satellite import dstatedt
    2- Also I'm receiving an error (TypeError: 'numpy.ndarray' object is not callable) when defining the initialstate like this: stateinitial = np.asarray([x0,y0,z0,xdot0,ydot0,zdot0])
    thanks a lot.

    • @CarlosMontalvo251
      @CarlosMontalvo251  Год назад +1

      Converting this to python is a pretty large endeavor. I would try watching one of my more simpler videos on python integration for help. You can check out this new playlist I have created. ruclips.net/video/ptgCdILud48/видео.html

    • @aeronautics2740
      @aeronautics2740 Год назад

      @@CarlosMontalvo251 thank you!!

  • @de_michael1222
    @de_michael1222 11 месяцев назад

    hey, which MatLab add-ons will I need for these simulations?

    • @CarlosMontalvo251
      @CarlosMontalvo251  11 месяцев назад

      You actually don't need any extra add ons. Just the standard MATLAB version will work.

    • @de_michael1222
      @de_michael1222 11 месяцев назад

      Thank you for the quick reply and helpful content Sir. Does this apply to the other ADCS in MatLab tutorials as well?

    • @CarlosMontalvo251
      @CarlosMontalvo251  11 месяцев назад

      ​@@de_michael1222 There is only 1 ADCS software on Github. There are multiple videos because every major edit that I've made has been done live on youtube. So if you downloaded the latest iteration of the software you wouldn't need any other packages.

    • @de_michael1222
      @de_michael1222 11 месяцев назад

      @@CarlosMontalvo251 alright, thank you very much!

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

    FAST

  • @diehardAMD
    @diehardAMD 2 года назад

    Get a damn noise cancellation headset. Your audio sucks. I have my audio all the way up and can't hear you.

    • @CarlosMontalvo251
      @CarlosMontalvo251  2 года назад

      I did. This is an old video with a bad mic. Perhaps I need to redo the video.

  • @tarekmahmoudkhalifa8067
    @tarekmahmoudkhalifa8067 Год назад

    Please help me, i need to plote 3d shape from run of Bvp4c code?
    Is there a way to send you the code to contact you

    • @CarlosMontalvo251
      @CarlosMontalvo251  Год назад +1

      I'm not familiar with bvp4c. Is bvp boundary value problem? What is 4C?

    • @CarlosMontalvo251
      @CarlosMontalvo251  Год назад +1

      @@tarekmahmoudkhalifa8067 I typically only help my watchers via youtube videos. If you describe the problem I can try and make a video about it.