cplxcomp.m

来自「digital signal processing-MATLAB」· M 代码 · 共 22 行

function I = cplxcomp(p1,p2)
%《数字信号处理教程——MATLAB释义与实现》子程序
% 向量重新排序后的序号计算
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
%  I = cplxcomp(p1,p2)
% -------------------------------------------------------------------
%  计算复数极点p1变为p2后留数的新序号
%  本程序必须用在CPLXPAIR 程序之后以便重新频率极点向量
%  及其相应的留数向量:
%      p2 = cplxpair(p1)
%
I=[];                               % 设一个空的矩阵
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中的位置