Curve Fitting in Python (2022)
Vložit
- čas přidán 6. 07. 2024
- Check out my course on UDEMY: learn the skills you need for coding in STEM:
www.udemy.com/course/python-s...
In this video I show how to use the curvefit function in the scipy.optimize library. I also look at practical examples from physics.
Tutorial Playlist:
• The Full Python Tutorial
Code:
github.com/lukepolson/youtube...
Discord:
/ discord - Věda a technologie
That's the scientific channel that we need. Good sound (which is the most important), good picture, Chad-looking science guy (in movies smart people are often nerds with big glasses) and a good explanation - this is really nice! Thank you for video!
Sometimes one may both know coding and physics but it becomes really exhausting to combine them together; you may not able to know where to start... You are directly attacking to that tricky intersection perfectly and that is exactly what is needed for lots of science students/academics!
I just personally wanted to thank you for such great videos as a PhD student in physics who had worked with paper and pen for years but needed numerical computation eventually. Your videos are awesome and very informing! I hope you keep going.
If you ever try your hand at 3D Graphics Programming and write the underlying core components of the 3D Rendering Engine where you implement a Physics Engine with basic kinematics, collision detection systems, A.I. Path Finding Algorithms, etc. then you get into doing Animations, Particle Effects, Rigid Bodies, Rag Dolls, Liquid or Water Animations, etc... by using advanced shaders for many various lighting, shading, shadow, and other special effects... That's one of the best ways to merge the two fields.
Now, this kind of Physics programming is different from Data Science, Equation Solver type Physics system, but it's a great place to start and most of the topics and transitions would already be covered. A basic 3D Rendering Engine or 3D Game Engine uses mostly linear algebra and trigonometry, but there is still a bit of calculus here and there such as when basic integrals, performing interpolations of various kinds, using runge kutta as well as FFTs all come in handy. Especially FFTs when working with Audio Processing and 3D Sound Environments is another great benefit.
You end up incorporating many mathematical and physics, (maybe even a little bit of chemistry, biology, and environmental systems depending on the context of your Engine and what you are going to build with it) concepts into your Engine that stem across many software engineering and program application development disciplines by using almost every kind of container and associated algorithm there are out there. Other key features would be file management and file parsers, you can end up writing your own scripting or run time interpreted language, if doing multiplayer over networking then multithreading and multiprocessing principles, compression and decompression algorithms, encryption and decryption algorithms come in handy. Then you have many other topics such as batch processing, Storage or Assets Memory Managers, Many Various Shaders with an accompanying Shader Manager, and so much more. This doesn't even include things such as GUIs, HUDs and Font Management...
Now, once you have the Engine Built with all of its core components (and I don't mean using a 3D engine such as Unity or Unreal...) Then the fun begins when you start building your game that is seperate from then engine but requires it as a core underlying library to function properly. Now imagine building a game like Factorio or Minecraft... Both of those games are Turing Complete... or even a 2D with 3D layer slices such as in Oxygen Not Included or a semi-factory but defense tower game like Mindustry which are both Turing Complete... or something similar to those... Then I think you can build just about any kind of software application that you could imagine except maybe an actual robust compiler-debugger-linker suite, or an actual Operating System... Oh they could be done... but they're a challenge. I found it was actually a bit easier to learn how to build a Hardware Emulator in C++ with a focus on the 6502 for the NES... Now that was a fun project as well! So almost anything is possible to build but it does depend on the limitations of your hardware...
This is just a summary of about a 20 year journey of independent studying purely self taught hobbyist and some of the many topics I've learned over the years! There's many other things that I didn't mention as this is starting to become quite long... It's been a great journey. Even watching Ben Eaters 8-Bit Breadboard CPU series was a delight! So much so, that I actually implemented it in Logisim.
Wow! You explain very clearly! I was fighting with my code since I did not understand how to work with curve_fit and you fell out of the sky! Thank you so much!
Great video. Especially focusing and repeating multiple times that our goal is the parameters and not just having a good fit! Parameters = Context = Reason why we're doing the fit
The kind of detailed explanation i needed! The examples nailed it! Thanks a lot :)
Your initial guess was good enough 😂, great tutorial!
Thanks a lot! This was a really thorough explanation which gave me everything I need to plot my data in a nice way :) Keep it up!
Two weeks ago I had a physics assignment involving this method, which was a absolute disaster. This is why I should watch your videos as soon as they show up. Anyway, Thanks for the amazing job!
This is exactly the tutorial I needed given I had to perform many counting experiments . I wish I had found this channel sooner before going insane over tutorials that made zero sense to me.
Super helpful and so nicely explained! Thank you so much for making this video! The multiple examples helped me a lot and in fact covered what I was as looking for!
Absolutely love this intro to SciPy's 'curve_fit' function! Suggest you to make a video on spline fitting, which is included into SciPy's 'interpolate' module. With splines it is possible to find best fit to a really complex data, including those presented in this tutorial.
dude your videos rock, i love the no-nonsense type of attitude when tackling things, cheers! (;
This is excellent content. Clearly explained a problem I've been having trouble grasping for months. Thanks!
Your videos are sweet! I found you earlier today and have already learned a ton. Thanks!
Thank you so much for making such videos. Specifically in physics. It is of so much help. Please keep doing such videos.
Excellent video. Very clear unpacking of the logic. Thank you! I've subbed.
Dude! this video is so useful! Great explanation and great material!
Thank you for making this video I particularly liked that you gave more than one example.
superb video with respect to physics and maths! Best application based explanation. Period. Please I request you to post such concepts with such physics based applications visualisation. I just feel blessed. Thanks a lot. tons
Amazing! Thank you. This has been perfect help with plotting an exponential curve to my data.
Nice! Note that (in 19:40) you need to set absolute_sigma=True to get the absolute value errors on the parameters. A good example is using the known formulas to calculate the error on the parameters (A,B) of a strait line fit and compare with the results returning with absolute_sigma=False e absolute_sigma=True.
It's a good comment
Your Videos are really helping me. Thanks a lot.
Brilliant explanation of the covariance!
You give quality out there. Respect man.
Thanks for explaining this so well!! Super helpful!!
Thank you now I understood how the errors behave!
This is great! Glad to find your channel
great videos! continue your greatwork bro! supporting from philippines here.
Thanks a lot for this nice explanation! I wonder what you use to fix one or some of the fitted parameters using the Scipy curve fitting function?
i found your videos about 10 minutes ago and mate, what an amazing content, you remind me of keith galli, he is awesome too
excellent explanation, very detailed
This is so very helpful...thank you
Great video, thank you. Hope to see you @ UBC some day :)
I learnt a lot from thus vedio, thanks bro👍
Thank you so much !! You've been a really big help.
This is super clear!
Hi, first off EXCELLENT VIDEO THANK YOU SO MUCH
Just for those who need this video, there are a couple of clarifications (sorry if you made them and I didnt notice), but its popt, pcov = scipy.optimize.curve_fit
not just curve_fit () . Unless you made a shortcut beforehand
Other than that , thank you again so much help !
Great video, just a little correction, as the physicist in me is happy to show himself every now and then: It's actually Lennard-Jones potential, sometimes also called Van der Waals potentials (dipole-dipole potential) big brother. I don't think there's a thing called Leonard-Weibeck potential reguarding atomic repulsion, but if there is and it's related to the LJ potential, let me know!
Oh I might have misspoke! Thank you for this correction 😂
pretty nice explanation, thanks! would it give some sense to use R or Scilab to achieve this, would it be more efficient?
Thank you very much for such beautiful tutorials.
Can you show fitting hysteresis models
good video. How would I estimate the data based on those parameters?
Thank you for the content,
I would like to see a video on hysteresis curvefitting with multiple conditions
Very helpful content. The error(as obtained from experiment) is sometimes not symmetric about the mean y position. Is it possible to encode that into curve_fit function?
Hello, Thanks a lot for yours tutorials! They help me a lot. I have a question. How to use curve_fit (or how to cope with the problem) if I have a system of ODEs. How to use it as a model for curve_Fit function?
Like 25 years ago i wrote code to fit multiple gauss function to some spectrum comming from Cherenkov radiation emitted from cristals. Those days everythibg had to be done in C/C++. I use the "Numerical recipies in C" book and lib and the Marquard Levenberg algorithm. Took 3 weeks to progam that (including some GUI and other stuff). Today with Scipy maybe one or two days... Some stuff actually get's better.
Great video! Thank you!
Hello Mr. P Solver,
First of all a great video on optimization. I have a question regarding the optimization and identification of values in a given curve. Since I am completely new to this, I am not much aware of reverse engineering. I have an experimental Time-Temperature curve of a heat transfer system and based on it I need to identify the convective heat transfer coefficient. Is it possible to make it? if yes could you please give me an example. Thanks in advance.
With Regards,
Preshit
Have you done any videos on FFTs or any Runge Kutta methods?
Sir, can u make video on Heisenberg spin chain and correlation between spins.
These videos are really great
could you please show how to add bound to parameter?
Please consider doing some stuff on Bayesian inference using Markov Chains, maybe using the PyMC library?
Can you tell something about fitting multiple gaussians in a signal?
Is there a way to do R^2?
Was expecting a good description of what's going on with the curve-fitting functions, albeit with boring data. But actually I found exactly what I was looking for. Trying to get the photopeaks and their resolutions in a similar spectroscopy. Also the walking through cell-by-cell helps for digestion. Used to begrudingly use python, but after a semester I much more prefer it over finicky excel sheets.
Very helpful video...👍👍👍
perfect please do more in physcis and more intuitive
Hey could you do a tutorial for a CFD simulation in python?
Great video!
In the given case that you want to work multiple samples, and perform simultaneous fits. Is that possible? It is something I use often in XPS.
Same question here. I have tried clustering without luck.
For the Gaussian fit again. If I add a +c factor at the end in the gauss_f definition, will the plot work if I expand popt with c_opt so that the measurement points are taken into account which do not look "parabolic"?
You definitely making a good thing, killing a Matlab!
Really great video!
awesome video as always. Learned a lot - especially the cov-Matrix was super helpful
One note: It is the Lennard-Jones potential, not the Leonard Weibeck potential
I had a question, what if the function had constants in its definition that were an array that take certain values at specific points of xdata?
Very good ,very Useful very : )
Thanks for the great tutorial. Since in real life we have only the data, how can we find the best non-linear model function that fits to our data?
Great video ! I just finished a lab report last Saturday, this would been useful jajaja. I would like you to make a video about animations in python tool
Same. I wish I had found this tutorial a couple of weeks ago for my lab reports. But better late than never.
The heatmap of the covariance is brilliant, I'm definitely going to experiment with that. You mentioned the importance of the diagonal, but what do the other values in the covariance matrix mean?
That’s how one parameter affects another if at all.
Please, make a numba tutorial :)
Your videos are awsome
Say I have a histogram plot with bars. I would like to do a skew norm fit to it. How do I proceed with this? can anyone help? like what initial data should I have
Can we solve a multi-variable function through curve-fitting , for example - y = f(x1,x2,x3,x4)?
Good job!
Hello Mr. P solver, is there a way to solve for a complex function- meaning that the function that I have defrined returns complex(both real and imaginary parts) simultaneously, I tried doing it this with Scipy.optimize.curve_fit() but does not help If you have any idea pls share, thanks :)
My values in y axis are all less than 1. So taking the sq root considering Poison statistics of photons gives very high error. What should I do ? Should i take some kind of scaling factor?
Simply WOW!
how to get that sample data, if i need
OH GOD!!, thank you very much!!!!!
Perfect
What is that algorith? I mean that is working in the background.
I am getting the following error. I tried multiple methods like using numpy.asarray, np.array...... Is there any other solution?
setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (5,) + inhomogeneous part.
thank you..
great vid
Amazing
where is the csv file ? @13:47
Thanks 👏👏
Hmmm i was expecting something more clicking this video. When you have object returned by the fit method it has method called plot. It shows data with fitted model as well as plots the residuals.
Great video as always!
May I ask you to recommend some good books to dive further into this stuff? Maybe the ones you learned from, if any? It would be great!
Keep doing these great videos, sir!
ye im also wondering how did he reach all that knowledge
@@jay89boy Same here, he might and actually should do a short video about giving some advices on how to dig deeper into the topic, especially book recommendations. He is a great guy and I think he will read this and respond in one way or another.
👍💯
you left out a 2 in the gaussian formula
but when i plug in 2 it gives wrong answer. i tried it on my dataset. without it it gives right answer. why?
This isn't at all related to the video content, but I have the same exact scar on my forehead
Your style annoys me but the content is good
Great content. Just some feedback, I feel it would be less distracting if you didn’t show yourself talking while showing content on screen.