I Generated Guitar Audio in python using NUMBA
Вставка
- Опубліковано 28 вер 2024
- Here we use the python package NUMBA to solve the FULL wave equation and create both animations and audio outputs.
Code located in the link below. Go to "Python Metaphysics tutorials" and then "Vid 5"
github.com/luk...
You have no clue how enjoyable these videos are. Finally someone who can connect my mechanics classes to my computational ones. Could you solve a problem with continuum mechanics next?
Currently working on a 2D schrodinger equation one but send me a problem and I'll consider it!!!
28:55 that sounded WAY more realistic than what i was expecting.
thank you for such an instructive and creative video!
Great effort! And just a suggestion: you solved the diff equation and got the time domain information of the string. Since FFT kills the time domain information of the simulated sound, I think it would be much more realistic that taking every smallest partion of the string as a seperate signal source (movement of the smallest string element in time) and combining all of the signals together. Or maybe instead of using FFT, using wavelet transformations would be more realistic.
this has quickly become my favorite channel, i love the videos. you do an amazing job!
Thank you, that means a lot! Look forward to lots more exciting topics in the future as well. And definitely let me know if you have any suggested topics
thanks for this! I'm new to python and this is a great start.
One thing: to make it audible, instead of the FFT, you can imagine using a virtual guitar pickup. Then your timeseries is just a chain of same positions in different rows. If you average out more positions, you can also get a bigger pickup.
That way you can easily model the pickup for your virtual guitar...
Complex topic explained nicely. Quite elegant in python with numba. I did something similar with the heat flow equation in fortran about 45 years ago. I set up the finite difference equations in a tridiagonal matrix then used Crank Nicholson to solve iteratively. Upper and lower boundary condtions could change over time as could physical properties of the media. I'd program the IBM punch cards, run them through the card reader in batch at night and pick up the printout the next morning. Each run cost about $25. One day a semi truck and crane turned up at the computer science building and they lifted a huge box through the roof and added 5MB of ram to the mainframe. We were in hog heaven, imagine 200x200 arrays! Can't wait to see your other videos. Thanks.
Wonderful memories...i also have similar experience in 90's when i created some indian music using foxpro 2.6....cherishing .....
You're killing it with these videos! Now I'm curious how the body of the instrument would change this, if you had FEM models for different "guitar" shapes... I understand that's a much more difficult simulation though :)
Damn, I clicked on this video just because, and then I watched through the end because it was really fascinating. I really enjoyed it, thank you!
When I watch your videos, I get happiness, excitment all over my body. Really interesting!
This is the second time I watch this tutorial. Still very amazing to me. Thx.
A new star is rising.
Woah !!! Just tippin' my pinkie into Python. Man !! It sure is an ocean. Your content is amazing. Thanks.
Subs from 🇮🇳 🇮🇳 🇮🇳 India
This is so cool and well explained ! Thank you !
This is really cool! Just one thing though: you set L=0.7 m (=700 mm) and Nx=101, so shouldn't delta x be 7 mm, not 0.7 mm?
Yes, good catch!I I will update the code.
Loved the video!
this video has accomplished 2 things for me:
1st: It's awesome you programmed a WAV.
2nd: Watching this made me realize I'm a moron when it comes to programming in Python.
Still interesting to watch though.
Amazing work my friend! Next I wanna see is the same thing, but using the wave function without the small angle approximation. You probably need a super computer for that, but should be interesting! :)
We don't need a supercomputer to solve the nonlinear wave equation, the same approach he used (finite difference) should be enough. What must be careful it's about the stability of the solver.
Excellent video again! I have one question though, why did you sum the hms list to get the waveform? The hms list just gives the amplitude associated with each harmonic right, so why should summing their amplitudes have any significance at all?
If the solution to the wave equation is f(x,t), then the integral of sin(n pi x/L)f(x,t) tells us "how much" of the wave is at the harmonic n. This can be expressed as g(t) (we integrate over x), and g(t) will have an associated frequency that gives rise to a certain waveform.
Let sum_n be summing over all harmonics and integral_L be integrating. Note that sum_n(integral_L(f(x,t)sin(n pi x/L))) = integral_L(sum_n(f(x,t)sin(n pi x /L))) = integral_L(f(x,t)) by Fourier theory, so if we wanted to include ALL frequencies in our final waveform, we would just sum f(x,t) over x (equivalent to integrating) at each time to get the waveform. This is also proportional to taking the average of f(x,t) over x for each time.
@@MrPSolver It should be sum(integral_L*sin(c*n*pi/L * t))
It seems like you forget adding sin(c*n*pi/L * t) the term. This should be the wave that it's added to the air, not just sum(integral_L)
@@adammusk4530 Maybe I'm wrong but I was thinking about the same. Glad somebody pointed it out but I couldn't understand his explanation. In the end I tried both, with the sin term there is a high-frequency envelope, so the two sounds similar. Then I recorded an actual harmonic on guitar for comparison, sadly if we zoom in, the waveform seems completely random, while our generated ones are very regular. Which means even if for some reason there shouldn't be a sin term, I won't get it, at least from experiment.
wow man! You are sooooo clean! Well done sir!
Super interesting!!
higher notes are unstable. any idea why?
Thanks for this! Might try a 2D membrane and see if I can nail down a snare sound. Get some Navier-Stokes sims of wind in a tube with holes and we can start a band 🤣
I'm so down with trying this out LETS GO
why did you choose NUMBA over cython?
was it faster or just what you were used to using?
24:55 would probably been faster if you didn’t use loops. something like
return (sin_arr*sol).sum(axis=0)
maybe sin_arr needs to be manually broadcasted though, forgot the exact syntax
numba is faster than vectorization
Super cool!
No way... clap clap clap youuuu youuu... yes youuu... awesome !
how can I use pencil codes in python?
Nice
If you care to measure the actual lengths of the strings on your classical guitar, you will find these to be 0.65m, rather than 0.7m as is the case with a Cello.
This is not much of an error at the end of a long day understanding the given algorithm, but an error nonetheless.
You seem to be fudging a lot to make a sound similar to an A3 note, however this A3 note measured to be 224Hz. 😢
wow😶👍
oooh
Subscribed bro....
i cant find the source on git, sorry i am dumb
I really enjoy your videos. This kinda video is what I'm looking for. Python + real application science. Thank you very much for this work.
Can you also make a video to simulate tsunami that generated by an earthquake of certain location? For example earthquake in Japan. If we can find the ETA of that particular Tsunami in some area close to that earthquake, that would be awesome
These videos are so good man, please keep it up! This is what makes physics such a beautiful subject :)
S T R I N G S
I am getting sol = compute_d(sol, Nt, Nx, dt, dx, l, gamma) as IndexError: Getitem out of range. any ideas.
Where did your equation for the stress-strain term come from? May I have the source?
When you approximate the fourth derivative, shouldn't you use the following scheme: [y(x+2h)- 4y(x+h) + 6y(x) - 4y(x-h)+ y(x-2h)]/h^4 instead [y(x+2h)- 4y(x+h) + 4y(x) - 4y(x-h)+ y(x-2h)]/h^4 ? Can you clarify?
Must be a typo because it is correct in the code
Love it bro
Where did you get that wave equation? I can't seem to find it anywhere online.
The damping is just added based on the fact that there will be air resistance (see air resistance formulas). The stress strain force (1D) is common among materials that bend (i.e metal strings)
@@MrPSolver okay I see equations where the fourth order term is multiplied by EI (the elastic modulus and second moment of area). Is that equivalent to the characteristic length?
A bit scary that you don't blink, but it's a great content! Thumbs up!
Damn numba seems very simple. Maybe I should update my simuls that way
I would like to see your python solution to sound attenuation. In the building industry, they use baffled insulated in ducting to reduce noise from the plant room. But, there seems to be NO equation that governs the air gap between baffles. From this video I'm think that the sound wave could be modelled with different baffle gaps.
Is there any way I can find further information about these things you call "baffles"?
❤️❤️❤️❤️