數字信号處理應用

%%
% 式(1.4)
N = 2;
alpha = 0.9;
k = 0:N;
b = zeros(1,N*D+1);
b(1+k*D) = alpha.^k;
x = filter(b,1,y);
           

%%
% 式(1.5)
alpha = 0.9;
b = [alpha, zeros(1,D-1), 1];
a = [1, zeros(1,D-1), alpha];
x = filter(b,a,y);