FIR Filter Design and Software Implementation - Phil's Lab #17
ฝัง
- เผยแพร่เมื่อ 25 มิ.ย. 2024
- FIR (Finite Impulse Response) filter theory, design, and software implementation. Real-time software implementation on a custom STM32-based PCB.
Overview of digital filtering, use-cases, mathematical theory, convolution, circular buffers, practical considerations, how to design FIR filters, filter design tools, implementation in real-time using C on an embedded system.
Visit jlcpcb.com/RHS for $2 for five 2-layer PCBs and $5 for five 4-layer PCBs.
Little Brain Hardware: github.com/pms67/LittleBrain-...
FIR Filter Design Tools: github.com/pms67/HadesFCS/tre...
and t-filter.engineerjs.com/
Circular Buffers: www.allaboutcircuits.com/tech...
If you're enjoying my content, please consider becoming a patron:
/ phils94
[TIMESTAMPS]
00:00 Preview of FIR Filter Implementation
00:34 JLCPCB Ordering (Custom STM32 PCB)
01:18 Introduction and Overview
01:57 Digital Filter Overview
02:40 Input Signal Representation (Discrete Time)
03:21 Filter Frequency Response
04:42 Output Signal
05:11 FIR Filter Theory (Impulse Response, Convolution)
08:20 Window-Sinc FIR Filter Design Procedure
11:17 Choices When Designing FIR Filters
12:19 Filter Design Tool
14:06 Implementation: Convolution and Circular Buffers
16:58 Header File
17:55 Filter Init and Update Functions
23:03 Filtering Accelerometer Data (DMA, Callbacks, USB, etc.)
25:13 Designing a Practical FIR Filter (Low-Pass)
28:37 Serial Oscilloscope Tool
28:50 Real-Time Accelerometer Filtering Example
ID: QIBvbJtYjWuHiTG0uCoK - วิทยาศาสตร์และเทคโนโลยี
![ขายสมบัติ Ferrari 488 ลูกรักก!! l [Nickynachat]](http://i.ytimg.com/vi/pwF1Il1Uevw/mqdefault.jpg)
New hardware and PCB design course on mixed-signal embedded systems just released!
⏵Course content: www.phils-lab.net/courses
⏵Course sign-up: phils-lab-shop.fedevel.education
Do you have any plans to look at implementing IIR and FIR filters using the STM32G series FMAC?
Thank you for your work! Also interested in the IIR Filter.
FWIW I googled the FMAC and found this informational video: st-onlinetraining.s3.amazonaws.com/STM32G4-Peripheral-Filter_Math_Accelerator_%28FMAC%29/index.html
Another great video Phil, thanks for clearing up how to actually obtain the filter coefficients, so far I have just trusted the calculators .. It makes complete sense that it would be the inverse FFT of the filters frequency response. It would have been nice to see the effects of scalloping prior to using the window function, just to really clarify what it does to a signal :)
-Keep 'em coming, love me some DSP.
Your vid was totally intelligibly explained. IMO the depth of the theory was on point 👍🏻. Also, i am looking forward to hear how IIR differs from FIR
I would also love some more videos on control theory! Starting at PID, which you have, and continuing to LQR and Kalman.
Took my undergrad DSP final 4 days ago, this video would have been incredible to have for my earlier days in that class.
To all students taking DSP, they are lucky to have your incredible videos, Phil
Thank you so much for your very kind comment, George! Hope all went well with your DSP final :)
This class is increadible, so clear. Wow just wow.
Thank you Paul, you've made this complex topic so much more approachable.
This is exactly the DSP video I needed to bridge the gap in my knowledge between theory and implementation. Keep these videos coming! I also like the fact that you share resources you found helpful.
Another amazing video. Great job! Your stm32 coverage is better than any I've seen on youtube, and you have such a wide breath of topics from initial board layout to fab, testing and implementing drivers and filters!
Really really precious tutorial!
The git repository is also a treasure as well!
Great job, Phi!
Just, Thank you for all your time and hard work to educate people. So refreshing to see complex concepts, explained thoroughly and implemented in a place where anyone can see the results. Merry Christmas mate.
Thanks for this quick FIR refresher and for mentioning t-filter. I found it to be extremely easy to understand and utilize for my STM32 projects.
Phil, your video libray of tutorials is absolutely amazing! Thank you for putting together these relatively short, practical, tutorials on extremely valuable information! 🙏
Thank you very much! I'm vey glad to hear your comments on these videos.
Cheers Phil! You make these topics very digestible and encouraging!
Thank you very much, very glad to hear that! :)
Wow, this will come in so handy on the automotive application I'm working on, which will be polling something like 10 signals from various analogue sensors. Thanks for yet another excellent video!
I just discovered this video now. Smooth as butter. It is like a good action movie, well constructed, fluent, interesting and also very informative. Thanks for the great effort.
Thank you very much for your kind comment!
Thank you so much for the content you create. These videos are immensely helpful and I appreciate the time and effort you put into this.
These videos are much clearer and more concise than many university courses and textbooks. Well done Phil!
I just thought about exactly the same!!!!!! Clear problem description, a bit of theory and an easy to understand implementation which makes the solution crystal clear, very very well done
this is the best video I have watched to explain FIR filters so far. thank you
Thanks for this - you've put a lot of effort into your video! Appreciated.
I'm stunned! Thanks for the class! Absolutly flawless! Keep up the great work!
Thank you very much!
A superb job in demonstrating from concept and theory to product.
I just wana thank you Phill to sharing your knowledge with us!
Appreciated that.
Thank you very much, Tomáš!
Excellent. Looking forward to your new course.
I'm in the last 6 weeks of MSc Digital Systems Engineering at the University of York, they have provided _LITERALLY_ zero DSP content... So thanks to Phil & his lab for actually providing this video.
Thank you once again for the amazing job you put in your videos. First digital control theory, then firmware configuration + hardware design and now this. It's like you are reading my mind or smth. Really amazing stuff. I'm all for a Q&A video!
Awesome, I'm glad you like them, Patrick! Let's see if I can read your mind with the next video :D
Thank you very much for these two superb an clear videos. Please keep teaching us about digital filters with this “hands-on” aproach.
Cheers!!
Thank you so much, Alberto! Many more videos to come :)
This is amazing. I am still very new and naive to dsp, but it’s amazing to see your content. Thank you so much for your willingness to share. It will not be forgotten ❤
Thank you very much :)
I’ve been looking for such good explanations for a long time!! Thanks!
You have done in 30 minutes what my brain has struggled with for hours - trying to bridge the gap between the practical FIR filtering code that I've examined with the theoretical/maths in Lyons book. This is one of the best tutorials I've seen on any subject - perfectly paced, clear and just the right level. Thanks!
Thank you very much, Robert!
@@PhilsLab One question after playing with your FIR Filter designer: what is the purpose of "Shift Samples"? Once you set the filter length (# of taps) isn't it implied that you will always shift by 1? Or is this a way to do decimation and filtering at the same time?
12 hours of lectures in DSP < this video. Thank you for doing this, it has been incredibly helpful!
Very glad to hear that, thank you for watching!
This is extremely useful as I'm currently working on an application based on imu data and require signals with very little noise. Thanks a ton !!
That's great to hear, thank you. Hope this can help out with your project!
This channel is a goldmine!
Amazing video and as an Eclipse contributor it’s nice to see other people using Eclipse for C/C++ as well :)
pure gold! I just learned at Special maths about convolution and Dirac impulse, Heavyside and so on but never though how to make the connection between theory and practice! As I said, pure gold! THANK YOU!!!!
That's great to hear, thank you very much, Lucian!
Thanks for another awesome video! 😁 I'm love watching your channel while working on my projects 😇
Awesome, thank you very much - very glad to hear that! :)
Great work as usual Phil! Keep up it mate 👍
Thank you very much!
Your channel is a treasure. Thank you
Thank you very much!
Thank you so much! this was so helpful, specially the firt part of the video, so I could better understand the concepts. This video was uploaded only 6 months ago, I'm glad I'm taking this course now
Thank you very much for watching!
Fantastic! Cant wait for the Kalman one!
I repeat what others said: WOW. GREAT STUFF!!
Thank you very much, Sanel!
Excellent explaination and excellent keboard sound.
really clear much more clear than my teacher used to explain this!
Thank you so much for those very clear explanation ! Your voice is perfect for non English speaker !
Thank you very much!
Nice content Phil, it helped me out a lot
As always love your work ❤️❤️
It was amazing to see this video. Keep up the good work. Please make more of such videos. :)
Great video man! FIR is really powerful, even my total bonehead implementation of 32 values moving average manages to increase DHT11 termometer resolution by more than order of magnitude.
Excellent work man. Thank you for sharing your knowledge in such a good way.
Thank you very much!
amazing video phil
another good educational video, good work phil really enjoyed this.
Thank you, very glad to hear that!
This is really gorgeous. Thank you so much. God bless you.
This is what I needed desperately. Skill upgraded :)
Awesome, glad to hear that, Frank!
Awesome, needed this ❤️
Thank you, glad to hear that!
Thank you very much, Phil!
Thank you for watching, Rodrigo!
Thanks for sharing this awesome thing.
This is awesome! Really hope for some more DSP videos!
Thank you, next video is on Z-transforms!
You are Awesome!! I have been looking at all of your videos, Just finished the ADXL355 build and took 20 pages of notes. I am trying so hard to understand all of this but it is hard for me because I am a self taught and just want to learn. You do a perfect job at teaching what you are doing.
Thank you! Really glad to hear that - many more videos to come :)
Very well explained, thank you. Waiting the next videos...
Thank you very much!
Thank you, useful content again. Additionally, I suggest SerialPlot Real Time Plotting Software to plot your data.
Excelente video. Muchas gracias Phil's
Thank you for watching, Wilson!
Clearer than sun.
Thanks you.
Awesome, thanks!
This is amazing !
Great video, thanks!
really great video sir!!
Welp, I better pick myself up some stationery then! Superb video :)
Really really great video, thank you for this. Great job
Thank you very much!
That was great Boss, thanks a lot!
In a previous attempt to understand FIR, I used Excel to model a Low Pass Filter. Rather than impulse, I used a step (square wave) and altered the response until it matched a sine wave.
Having the proper math would probably be better, but it worked as a one-off.
WOW PAUL! I am so happy with this one! You are such a good boy!
great video and great job !!!
Great video! Thanks a lot!
Thanks!
Kalman filter in discrete time is pretty interesting. May be you can describe in details Madgwick filter?
Also very interested in that! Cheers
Also interested. Cheers
And I'm gonna rewrite them using fixed point arithmetics to make life happier for low end MCUs ! go on mate !
Even though the main topic of these videos seems to be µC programming, they still give the most understandable explanations about digital filters i ever watched. Could you perhaps a video on Fast Fourier Transform FFT and its Inverse (IFFT) implementation in a similar way of your FIR and IIR videos ? Greetings
Really great video, watched the one after it about IIR as well. This is what TH-cam has been missing. Put more math. If people run away then they are your wrong audience. I loved it, I tried to show some how this was done and had difficulty with tools to bring a proper demonstration. This video will be my goto from now on.
Thank you so much! Yeah, I'm always unsure of how much maths to put in with it still being accessible to most people. I'll try and sprinkle in more and more and see when people have had enough :)
This is tremendous!
Thanks, Kadir!
very helpful... many thx
Amazing channel 👏👏👌
Amazing video Phil! The only thing that was getting to me was that your voice was competing with your clicky keyboard, which made it difficult to follow the software as you spoke. But that is by no means an indication of the quality content you share!!... So please keep on sharing your knowledge with us.
very good video. Summarizes a lot of hours of explanation by my "signal processing" professor :D Thanks!
Thanks, Joey!
Nicely presented. The C code can be made a lot more performant by storing two copies of the impulse response and dropping the branch in the inner loop.
Thanks for the video!
Very nice!! Thank you!
Thank you for sharing.
Where you format the accelerometer output before calling CDC_Transmit_FS(...) sprintf() returns the number of characters written to logBuf[] so you can save yourself the call to strlen() which is yet another pass over logBuf. I'll keep my comments about using sprintf() in ISR's for another day :-) ... Ok, well, this is relatively low rate stuff so formatting and transmitting from another ringbuffer in the main loop fed by the ISR would be another option to keep the ISR's lean and mean. Great stuff.
Thank you so so much for this!
Thank you for watching, Damola!
Thank you. Very instructive. I would've perform the convolution using the Discrete FFT.
Programming tip: If you keep buffer length a power of 2, i.e. 2^4=16, you can reset the index faster using: fir->bufIndex = ++fir->bufIndex&(FIR_FILTER_LENGTH-1)
One comment: if filter length is L=2^N, then the group delay of the output comes in-between the samples: (2^N - 1) / 2. For example, for N=3, the group delay is 3.5.
This can be important, if filter data need to be combined with the input data.
Great video, perhaps would be nice to see a fixed point Implementation, or maybe compare the two on uC with or without a FPU.
Wow I wish I saw this video a year ago! I hope you continue this series. It would be neat if you did a video using DSP to do hand motion detection on a microcontroller.
Thank you! I'll be making a couple more DSP videos. Next ones on IIR filtering and Extended Kalman Filters.
Cool. I learned a lot.
Awesome, glad to hear that!
Keep it up Subscribed.!
Your keyboard is very passionate.
I wish there was a bit of electronics demonstrations. I don't know actually where all you are talking about is in electronics world.
1. What hardware setup you are using?
is it the microcontroller dev board or it's another piece of hardware?
2. What level of electronics this happens? audio, wifi, bluetooth signals...etc. ?
What I liked:
1. The code explanation is very good especially following up other drivers libraries you've done in other videos
2. Theoretical information about is good even I didn't get a clear idea of where FIR should be implemented?
Two last days I've being reading articles about PCM DSD, and now you have published this video. And this is not a random vid but rather one of my subscriptions. I wonder what does it mean ))) And how do you have time for all these cool things for us?
Haha maybe I can read your mind :D It's getting a bit tough with time but I'm trying my best to get two vids out a month :)
Thanks a lot for the super interesting videos dude! Would you Say that trasforming a continuos time filter from s to Z domain, night be an equivalent workflow?
Nice video bro, tnx 👏👏👏👏
Thank you, Nicola!
You are amazing!
Thank you!
Very cool training course, so well done, I love it. Don’t forget to setup a PayPal URL to support your activities for me and others.
Here is some minor programming feedback.
- for the sprintf use always the snprintf function e.g.: snprintf(logBuf, sizeof(logBuf), …)
This is a good programming practise to avoid buffer overflows.
- For circular buffers the wrapping can be made easy via modulo operations.
e.g.: fir->buf[fir->firbufIndex++ % FIR_FILTER_LENGTH] = inp;
This way you never need to check if your are at the end of a buffer or at the beginning)
Again really great stuff you are doing here.
PS: I would write this as a C++ class where the constructor gets the length of the filter as a parameter and allocates the float array, however for the audience here keeping it simpler with basic C is better because everyone must understand the individual steps.
Regards from Hannover, Helmut
Amazing tutorial and very clean code. When i handle interrupts, what I do is set a flag when the interrupt is identified, then immediately clear the interrupt. Then that then triggers execution of the appropriate callback in the super loop. I do this to keep the interrupt vector free and clear. What do you think about this approach? Of course if I was using FreeRTOS, I would have some kind of thread safe semaphore or even a queue.
But damn man your explanation of FIR filter was spot on and probably the best I have ever seen.
Edit: also, subscribed! Keep up the amazing work!
Very good.
Thank you, John!
Thanks a lot, I will try implement your solution, I have problem in ntc sensor!