fir1.m

来自「包含有各种各样的数字信号处理经典算法源代码,很有用的.」· M 代码 · 共 16 行

% fir.m - sample processing algorithm for FIR filter.
%
% [y, w] = fir(M, h, w, x)
%
% h = order-M filter (row vector)
% w = filter state (row vector)
% x = scalar input
% y = scalar output
% based on fir2.c

function [y, w] = fir(M, h, w, x)

w(1) = x;                          % read input
y = h * w';                        % compute output
w = delay(M, w);                   % update delay