cplxcomp.m

Digital Signal Processing Using MATLAB V.4随书代码

function I = cplxcomp(p1,p2)

%  I = cplxcomp(p1,p2)

% Compares two complex pairs which contain the same scalar elements

%  but (possibly) at differrent indices.  This routine should be

%  used after CPLXPAIR routine for rearranging pole vector and its

%  corresponding residue vector.

%      p2 = cplxpair(p1)

%

I=[];

for j=1:1:length(p2)

    for i=1:1:length(p1)

	if (abs(p1(i)-p2(j)) < 0.0001)

	   I=[I,i];

        end

    end

end

I=I';