Convolution is a mathematical operation that blends two functions relative to the overlap of one function as it is shifted over another. Although MATLAB contains a pre-built convolution function, it is possible to calculate the discrete convolution integral yourself. The discrete convolution of two functions f and g is defined as the sum over the range 0 to j of f(j) * g(k-j).
Define two vectors, f and g, containing the two functions you want to convolve. The lengths of f and g do not have to be equal. The length of result of the convolution, k, will be one less than the sum of the length of f and g:
m = length(f); n = length(g); k = m + n - 1;
Define the range j over which the convolution will occur. The value of j is the range where subscripts of the two functions to be convolved, f(j) and g(k+1-n), are legal. The value of 1 added to k is to account for the fact that MATLAB begins indexing vectors at 1 rather than 0:
j = max(1,k+1-n):min(k,m)
Preallocate space for the result of the convolution:
my_result = zeros(k);
Write a for loop to iterate through values of k:
for result_index = 1:k
Calculate the convolution for all values of j:
my_result(k) = sum(f(j) .* g(k-j+1));
Close the for loop with the "end" command.