Coding Challenge #93: Double Pendulum
Вставка
- Опубліковано 27 тра 2024
- In this coding challenge, I create a double pendulum simulation in Processing. Code: thecodingtrain.com/challenges...
🕹️ p5.js Web Editor Sketch: editor.p5js.org/codingtrain/s...
🎥 Previous video: • Coding Challenge #92: ...
🎥 Next video: • Coding Challenge #94: ...
🎥 All videos: • Coding Challenges
References:
🔴 Physics Lab Double Pendulum: www.myphysicslab.com/pendulum...
Videos:
🚂 Pendulum Simulation: • 3.4: Pendulum Simulati...
🔴 Coding Train Live 120: • Live Stream #120: MNIS...
Related Coding Challenges:
🚂 #86 Cube Wave by Bees and Bombs: • Coding Challenge #86: ...
🚂 #159 Simple Pendulum Simulation: • Coding Challenge #159:...
Timestamps:
0:00:00 Double Pendulum Simulation
0:03:45 Creating The Double Pendulum
0:08:35 Adding The Tracing of The Path
0:12:51 Implementing The Double Pendulum Formulas
0:24:45 Visualizing It
0:28:21 Adding Damping
0:29:16 Conclusions And Suggestions
Editing by Mathieu Blanchette
Animations by Jason Heglund
Music from Epidemic Sound
🚂 Website: thecodingtrain.com/
👾 Share Your Creation! thecodingtrain.com/guides/pas...
🚩 Suggest Topics: github.com/CodingTrain/Sugges...
💡 GitHub: github.com/CodingTrain
💬 Discord: thecodingtrain.com/discord
💖 Membership: ua-cam.com/users/thecodingtrainjoin
🛒 Store: standard.tv/codingtrain
🖋️ Twitter: / thecodingtrain
📸 Instagram: / the.coding.train
🎥 Coding Challenges: • Coding Challenges
🎥 Intro to Programming: • Start learning here!
🔗 p5.js: p5js.org
🔗 p5.js Web Editor: editor.p5js.org/
🔗 Processing: processing.org
📄 Code of Conduct: github.com/CodingTrain/Code-o...
This description was auto-generated. If you see a problem, please open an issue: github.com/CodingTrain/thecod...
#doublependulum #pendulumharmonicmotion #processing
![[柴犬ASMR]曼玉Manyu&小白Bai 毛发护理Spa asmr](http://i.ytimg.com/vi/0TsXQ7z2Dh4/mqdefault.jpg)
13:15
Note: the following proof is trivial and has been left as an exercise to the reader
lmfao love those videos
Physics books be like
Lagrange bra.
ua-cam.com/video/tc2ah-KnDXw/v-deo.html
in case you are interested ;)
@ bill bassichis
Fun fact in case nobody else has mentioned it: updating the positions using the "current velocity" vs the previous one changes the numerical method from Euler's method to the Euler-Cromer method, which is quite a lot more accurate. The regular Euler method is prone to inadvertently increasing the energy of the system, while the Euler-Cromer method conserves energy quite decently. The two methods seem almost identical, but the behind the scenes math is quite interesting and shows that they are actually of different orders.
"Okay, I'm gonna see if I can just copy and paste this enormous differential equation right into my code."
i have no idea wats going on. yet im enjoying this video.
Story of my life online
Same
Same
Like 2 livers connected and if we swing it from the free end, it creates a beautiful art. I guess that's what he did. Not sure.
😂😂😂 Me too
“G is the universal gravitational constant, I assume. That’s just gonna be one.”
If only it became true at the moment he wrote it ! Some physics genius would have re-died. ^^
Missing Semicolon I died. Small g is acceleration of gravity on earth = 9.81 m/s^2. Big G is the universal gravitational constant, 6.67 x 10^-11. TRIGGERED!!!
Most of the time you don't use actual real values in code because in order to have them work properly ud have to do a lot more work.
ħ = c = G = 1
This was golden hahahah thank you so much for this
To those who want to make this simulation more real do this.
Currently Dan is setting g = 1. Although its far from the actual value the simulation looks very real.
This is because one second in the simulation is one frame. And by default processing or p5.js is set to 60 fps.
This means when one second passes in real world, 60 seconds would have passed in the simulation world.
Hence the correct value of g would be 9.8/ ( 60 * 60 ) .
There's still one more thing to adjust. The masses are actually in kilograms hence you should set the value of bob masses to around 0.05 which is 50 grams. Also, the lengths are in meters so make that they aren't very large.
Thanks for these comments!
This Guy has made a magic. But when we give the real numbers within no time, the software is getting into errors. he has very intelligently demonstrated the capabilities given the limitation of a 30 min video. The world appreciated him for this.
Nope
I changed and now it's really really slow
OMG, I've watched a 30 min video without even noticing it. Thanks for existing, Dan!
João Vitor Bonadiman
Same...
When I noticed the video would end soon from the way he talked I looked up the runtime because I hoped there was more to come with messing with the values. Then I noticed I've watched 30 min of coding feeling like 5...
I really REALLY laughed when the pendulum just started warping into another dimension at 12:30
*"See, this is kind of the idea"*
you just made my day hahahaha
Hahahaha
He just made the flower of life
Was my thought too :D
24:10 you are changing the integration method from simple explicit Euler to symplectic Euler which is numerically much more stable.
Top comment!! This explains your astute observation at 23:33 "It's got quite a bounce to it." The simple/forward Euler method tends to artificially add energy to this system, while the symplectic/semi-implicit Euler method is better at conserving energy. There are many other variant methods, including the so-called midpoint method (angle is updated by the average of present and future velocity steps).
W h a t
I found some math gods here...
I think a cool feature to add would be a fade of the trace line, so as a segment of line is drawn, say 5 seconds later it begins to fade, so you only ever have the last 5 or so seconds of trace at any given time
Something like this in draw:
fill(255,10);
rect(0,0,width,height);
I never knew i could be so engaged and laugh at someone doing something like this. Great video! :)
17:49
"Oh I forgot a semi colon"
Me: (Laughs in Python)
Imagine running your programs at 60 fps
Me: laughs in c++
@@killereks imagine not having Dynamic Memory Allocation
Me: also laughs in python
@@thewild2334 imagine spending too much memory on a slow list.
Me: laughs in C
@@killereks Imagine running your program in the browser.
Me: laughs in JavaScript
why did you guys stop??
You’re one of those special people in the world that make it a better place. You can always tell, and it just radiates from you.
I'm so happy you've done this; I started making one of these in P5 using a class, and I just figured you'd use basic Trig to calculate the angles. Had literally no idea the lre was a huge algorithm derived for it! Gonna give this a looking soon!
We need more teachers like you, I'm starting to enjoy coding thanks to you!
this is the first of your videos i have seen, i am so glad people actally go through these sort of things for others to see
keep up the good work :)
"I'm not going to derive those formulas" *cries in physics major
@Max Beer more easy to use Hamiltonian
@@termitori Hamiltonian is advanced Lagrange
@@mr.mirror1213 but it makes it so easy. the hamiltonian destroys everything you thought you knew about classical mechanics and it's fucking beautiful
Exactly man, after using lagrangian it bacame easior
Idk why this was recommenced to me, I never even coded, but you're so full of energy and entertaining, thumbs up dude
First video I watch from this genius. You’re funny as hell, your channel’s gonna get me hooked. Thanks.
you are a blessing !!!! a true blessing. I cant thank you enough for what you are doing.
The first time i see, why i had to learn all the Math in school xD
All those physics courses really made me understand this, I'd be lost without taking those courses tbh.
It's so satisfying to see the joy on his face as if he became a father😄great man to be convinced to get into coding
THIS IS A PERFECT EXAMPLE OF HOW TO DO AND HANDLE COMPLEX PROJECTS AND COMPLEX ALGORITHM. YOU'RE GREAT!!
This is the best UA-cam channel I've ever seen
Pretty cool! A few things: the dampening should be a function of velocity, proportional using some constant. It would be harder/maybe fun to calculate the forces along the bars and have them extend elastically under the load.
thanks for this feedback!
There is a small error here in the equations of motion. Because the canvas draws the origin in the top left, the derivation used here is slightly wrong. This is accounted for in drawing the pendulum. However the error propagates into the angular acceleration equation.
There are two ways to fix this:
1) Fix the equations of motion by doing a coordinate transformation of the solution
g -> -g, and cos(angle1) and cos(angle2) also go to negatives. This is actually a bit harder than it looks because there are a lot of trig identities you'd need to pull out to di this.
2) Fix the coordinate system
When I implemented this myself using plain old JavaScript the Canvas Context object allows you to scale it by a negative number. Then fix how you draw the pendulums so y = -cos(angle) * length. I believe a way to accomplish this in Processing would be the rotateX() function.
Thank you for this feedback!
hey, I've encountered the exact same problem you pointed out, could you help me with it in more detail ? I'm kinda lost. For g--> - g, easy to understand but what do I do for the trig identities ? Thanks
You do great videos, congrats!!!
Make a pendulum wave, multiple pendulums with differents arm length, that effect is pretty awesome !!
awesome :D i always get so excited when i see a new video on this channel. I even managed to attend some live streams.
I love how weird you are :D and how you can explain these pretty complicated concepts in a way that is so easy to understand.
I've tried the same with the matter.js library but failed miserably :p nicely done!
Quick wallpaper I did using this: i.imgur.com/BSb37so.png
Cédric Reber nice wallpaper
Too Good
You are the one Level 99 developer that I've spent my entire life looking for. If there was a 10th Dan for programmers, it'd be yours. Thanks for existing.
And by the way, if you or anyone else want to solve the problem at the end ( I mean that the pendulum is getting higher and higher and even faster) you need to use another way to solve numerically the differential equations.
Here you used the Euler method (explicit) ( v += a and x+=v).
There is a loose of accuracy because of the integration step.
In order to fix it (or to make it better) you can use the RK4-method (Runge Kutta 4) which will make you gain in accuracy without reducing the integration step of time (which is expensive in terms of calculation time).
I know it's a bit more complicated but it's worth it ( I personally have tested it with pendulums and with gravity simulation with planets and orbits which can degenerate).
man , watching this smart and fun guy is another level of happiness :D !
that is VERY hard to create a double pendulum, just because it's so random and hard! I didn't watched the video yet, because as always; I'm first trying the challenge myself and then trying your code. This one was IMPOSSIBLE!! good job Dan!!!
You should try to create a simulation with vectors as well, it would be interesting.
(Calculate the motions without formulas, just using vectors and addition)
Wow, Dan! I didn't know you are so good at object oriented programming!
22:50 The most beautiful sound you can hear...
"I'm right about the angular acceleration, because it's going in a perfect circle"
- Dan 2018
i've literally never watched a coding video before (but i've at least taken calc and physics haha) but somehow this made sense, so props to you as a teacher!
Glad to hear! You can use either of these playlists to start from the beginning!
ua-cam.com/users/shiffmanplaylists?view=50&shelf_id=14&sort=dd
ua-cam.com/users/shiffmanplaylists?view=50&shelf_id=2&sort=dd
This was my second year BSc Physics coding project (circa 2016). It was hardcore, but definitly do-able :D
Dude you are a gifted individual.
15:00 you hurt the physicist in me so darn hard XD
René Carannante Actually, that's exactly what a physicist would have done. We set constants to 1 all the time. Like in relativity, c is always 1
First of all, c is simply the speed of light, it is NOT 1. In special relativity we treat it as 1 because that's the the ratio speed of light/speed of light = 1 , and other speeds are written in terms of c so it's easier to do lorentz transformations.
Physicists don't set constants. Ever. You can't just set a constant to whatever you want.
Pati, clearly you aren't a physicist. There is nothing fundamental about the meter or the second that makes them absolute in any way. All dimensional constants are dependent entirely on the measuring system you use them with. Who sets the measuring system? Physicists. You CAN set the speed of light to whatever you want, just by using whatever length you want. It's not just a notational trick, the speed of light IS one lightyear per year, or one lightsecond per second, etc.
The only FUNDAMENTAL constants are the dimensionless ones, like the fine structure constant or the strong coupling constant, dimensional constants like the speed of light and gravitational constants are all secondary.
Pati I think you must understand that anything can be referenced in relation to everything, regarding certain vectors of location and assuming small scale locality and infinite spatial dimensions of said locality. Therefore, symbology and meaning in this world is manifested solely based on our reality and thought processes from the point of which we stand.
There’s a reason why everyone else besides the US uses the metric system. Because it’s easier to work with. Why would you want to make things harder for yourself?
@@pati7132 I agree
This was awesome. Thanks!
Absolutely love you videos. I also found out that if you try to compute the total energy of the system, it is NOT conserved, I wonder what is wrong, everything seems to work fine except that
Go Dan! GO!
I just found my new favorite channel.
that was awesome man!!
That probably the craziest formula I've ever seen
You are simply awesome!!!
I go to know so many things related to coding bcuz of you thanks a lot for making such amazing tutorials. Now I can kinda convert Java into Javascript(as thats the only language I know). They are so fun to follow along!!!
Now me as a physics student want to make 3D version where all of the formulas are calculated by me! Funny thing is that last semester I had classes called Classical Mechanics and one of the tasks on final exam was to calculate this formulas in double pendulum ;>
Sounds hella fun. Would love to see a 3d version!
Morszczux Lagrangians save lives
Mr. Fluffypantz True!
Amazing... Chaos moviment.. and the butterfly effect. You are great!
Cara, muito pai dégua teu trabalho!
Professor de Física Ele aceitou na programação mas errou na física
You are surely the best coding tutor available
Thanks for introducing me to Processing....
You can check your code by making both angles kept equal throughout and starting out with a small angle such that it recreates a rigid rotor case.
Great Video, thanks for sharing.
For those who does not understand why at 24:30, switching the velocity before position fixed the problem, here is the technical explanation:
it seems like its the same, but if you work out the math, position first is called Euler's method whose error is not bounded. As time goes by, the error will be big enough(if you plot the energy of the system, it continuously goes up) , it seems like it getting more energetic. If you put the velocity before the position, the error is bounded. He accidentally used a method called Verlet method whose error is bounded in this case. therefore, as time goes by, the energy of the system fluctuate within in a range, therefore, it looks real.
29:03 actually you don't need that dampening since you have in the g value, since it's 1 it won't really change the result of the formula but if you tweak it to be 9,807 m/s² (earth gravity) it'd change, I like moon, mars and sun values to play with.
So much fun watching this
Hey! I'd love to see you implement the '4th order runge kutta' method for that. What you did (v+=a; x+=v) is called 'Forward Euler Integration', 'Runge Kutta Integration' is basically the same, but with multiple intermediate steps. That achieves an accuracy that goes down two more taylor iterations, removing the 'energy gain' over time.
Yes, it would be nice to cover this sometime!
I did that here for multiple pendula: github.com/kevinhp/multiplePendula
See here for a working implementation:
www.openprocessing.org/sketch/507468
Cheers, Daniel
The reason that updating the velocities and then the positions works better is because it's a symplectic forward Euler integration scheme which more closely captures the symplectic geometric structure of the (Hamiltonian) double pendulum system, and so you get less numerical integration error. You can also write down the general equations of motion for a generalized pendulum fairly easily in the case where the lengths and masses of the pendulums are the same, and I'm sure it's just a tedious generalization to consider distinct masses and lengths.
Thanks for this information!
The double pendulum is a system of four first order non-linear differential equations. A better method of approximating the solution would be an ODE solving algorithm like Runge-Kutta. Simply adding acceleration to velocity and velocity to position will only be valid for a very short time interval after release due compounding error.
It was interesting, thank you)
The pendulum going crazy is my favoirite part
Rather than calculate the angles, I calculated the force in the connecting rods and used a pair of standard movers, with the child pendulum applying its rod force to the parent pendulum. Seems to work fairly well. I also call canvas.background() with a very transparent color so the trace fades over time.
Next challenge
Gravitational Waves
Purifier Phoenix The Mecca hell yeah
wow, that's pretty awesome.
g is acceleration due to gravity
G is the universal gravitational constant
This code helped me make a spirograph generator.
One tiny thing I would change is making the width and height of the masses proportional to the square root of the mass, that way the area is proportional to the mass.
7:36 You could delete that second fill function call just before the ellipse function call you are in this moment, sine the you are just resetting it to the same value. Indeed, the program also works if you move the fill function call to setup()
For those interested, formula comes from approximating Runge Kutta numerical solution. You might as well program it so you do not have the need to write the "analytical" solution.
Awesome !
25 mins in and I'm subbing. You're awesome. If you need help with understanding these equations I would be happy to help out. I'm an currently a master engineering student in the Netherlands. And I know python language coding
"Somebody crazy calculated it"
Yes, intro phys students are absolutely insane
You are really awesome man.
Dan, i would like to challenge you to recreate this but in 3D. Still a double pendulum but can rotate along a z axis aswell. Think this would make a cool project, since it seems like a very difficult thing to recreate in the real world.
Two strings Two masses and a conection to the celling. but a simulation would be cool either way.
Hi , I have a doubt. Without getting the function called , how is there constant incrementing and decrementing in the first third of the video?
12:43 hmm that pendulum looks like a famous logo... that's a new brand I believe?
95/5000
it will not be the same pattern, the inverted pendulum is a non-linear and chaotic system ... great work
It could be a good beginner project, specifically for people learning to deal with neural networks, to train a neural network to predict a double pendulum.
What is the font that you use to get those 'i' and 'l'?
I like it a lot.
The way to calculate the double pendulum is actually not incredibly complex. There's really 2 main things you have to know. 1.) you have to know the physics in order to get a formula for the kinetic and potential energy for each of the masses. 2.) you have to know the Euler-Lagrange equation. you use the formulas for the kinetic and potential energy to get two differential equations, and really the complicated part is the algebra for actually solving for the angular acceleration for theta1 and theta2 because the equations is a little long and both accelerations show up in both equations. The original differential equation is (relatively) simpler and not too complex to see where each part of the equation comes from.
I love your work so much
"that seemed to have fixed it" -every coder working on bugs
Loved it
One question, what video editor do you use to overlap you so smoothly to your screen? Thx
You're using the 1st order Euler method to solve the differential equations, it is very rough and you get the increasing total energy to your system with time. Check please the sum of kinetic and potential energy of the pendulum. To fix the problem you should damp the velocity in such a way that the energy is kept constant.
What font is that which is used in the text editor in this video, it looks really nice?
I love this guy
In the limit as _r2_ , _m2_ → 0 and _r1_ , _m1_ get large the behaviour should in theory also start to approach that of a single pendulum
Update the velocity before the position, it is a more stable integration method (semi-implicit euler integration)
I know absolutely nothing about coding and all of this is foreign to me but I love watching you lol
Mathematical art's great for introducing kids to math, physics and programming. This 2D representation of a 3D object captures interesting properties of a double pendulum, but I'd imagine that 2DoF robotic arms (Scara configurations?) probably treat it better. What I'd like to see is a colored map that shows where the pendulum weights spend most of their time and I'd like to see implementations of quantum algorithms to generate it. I'd guess it would resemble a blobby v shape. I bet this has a name and a figure, but I have yet to come across it.
Also, if the video outlined coding concepts it covers, it would be more useful for structured learning, for example in school.
P.S. While looking for previous Processing tutorials about double pendulums, this guy popped up and I think he looks like Metro/Music Man; no offense, just an observation.
scripps.ucsd.edu/profiles/hstaudigel
BTW, those 'crazy' expressions of accelerations can be derived using euler-langrange equation (a part of langrangian mechanics). Infact, It's very simple, not crazy at all :)
Just, wow.
Any idea why the translate function is not translating to the coordinate (300, 50)? It draws the tracing points at the origin (0, 0). I tried with the translate func at different places, no effect. I am at 12:40
Hi, this is a great video! Could you please tell me the environment you have for coding. Like language and libraries, etc. ?
This video uses Processing (which is built on top of the Java programming language). For more info, visit processing.org and also this video might help ua-cam.com/video/AmlAiKsiy0o/v-deo.html.
What environment do you code in? How can you have all the methods without having to import a library? Also is that java?