function I = cplxcomp(p1,p2)
%《數(shù)字信號處理教程——MATLAB釋義與實現(xiàn)》子程序
% 向量重新排序后的序號計算
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
%  I = cplxcomp(p1,p2)
% -------------------------------------------------------------------
%  計算復數(shù)極點p1變?yōu)閜2后留數(shù)的新序號
%  本程序必須用在CPLXPAIR 程序之后以便重新頻率極點向量
%  及其相應(yīng)的留數(shù)向量:
%      p2 = cplxpair(p1)
%
I=[];                               % 設(shè)一個空的矩陣
for j=1:1:length(p2)                % 逐項檢查改變排序后的向量p2
    for i=1:1:length(p1)            % 把該項與p1中各項比較
	if (abs(p1(i)-p2(j)) < 0.0001)  % 看與哪一項相等
	   I=[I,i];                     % 把此項在p1中的序號放入I
        end
    end
end
I=I';                       % 最后的I表示了p2中各元素在p1中的位置