代码: b = [1/3, 1/3, 1/3]; a = [1, -0.95, 0.9025]; % x(n) y(n) coefficient [R, p, C] = residuez(b,a) Mp = (abs(p))', Ap = (angle(p))'/pi % polar form Y = [-2, -3]; X =[1, 1]; xic = filtic(b,a,Y,X); %% ---------------------------------------------- %% S