convolution_sum.cpp

#include "AnalyticDrawing.h"
#include <raylib.h>
#include <vector>


class ConvolutionSum : public AnalyticDraw
{
    private:
        float altSinir;
        float ustSinir;
        float artisMiktari;
        std::vector<Vector2> points;

    public:
        ConvolutionSum(float altSinir = 0.0f, float ustSinir = 2.0f, float artisMiktari = 0.001f)
        {
            this -> altSinir = altSinir;
            this -> ustSinir = ustSinir;
            this -> artisMiktari = artisMiktari;
        }
        
        template<typename Func1, typename Func2, typename... Args>
        void CalcConvSum(Func1 firstSignal, Func2 secondSignal)
        {
              float convValue;
              points.clear();     
              for(float i = 0; i < ustSinir; i += artisMiktari )
              {

                 auto convFunc = [=](auto time, auto a, auto f) {return firstSignal(time, a, f) * secondSignal(i-time,f);}; 
                 convValue = Integrate(altSinir, i, 1000, convFunc, 5.0f, 10.0f);                       
                 points.push_back({10 * (i * 100) + (800), -((convValue) * 100) + (450)});  
                
             }
              
        }
        void DrawConvolutionGraph()
        {
            DrawSplineBezierCubic(points.data(), points.size(), 2, GREEN);
        }
        template<typename Func, typename... Args>
        float Integrate(float a, float b, int steps, Func func, Args... args)
        {
            float h = (b - a) / steps; 
            float area = 0.5f * (func(a, args...) + func(b, args...));

            for(int i = 1; i < steps; ++i)
            {
                area += func(a + i * h, args...);
               
            }
            
            return area * h;
        }
};