commGSP/commGSP_filter_poly.m
234a68af8299master

commGSP_filter_poly.m
View Options

			function y = commGSP_filter_poly(A, signal, p, type)
			% Filters a graph signal with a polynomial filter
			% INPUT:
			% A - (symmetric) adjacency matrix
			% signal - column-vector signal to filter set of signals organized in a matrix
			% p - vector of polynomial coefficients of the filter
			% (first entry corresponds to the highest order term, i.e.
			% to K-1 if K is the length of p)
			% type - string defining what shift operator to use when filtering
			% 'laplacian' - Laplacian matrix
			% 'modularity' - modularity matrix
			% OUTPUT
			% y - filtered signal vector

			X = commGSP_build_matrix(A, type);

			ord=length(p)-1;
			y = p(end).*signal;
			tmp=signal;
			for iter=ord:-1:1
			tmp = X*tmp;
			y = y + tmp.*p(iter);
			end

			end