cplxcomp.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 13 行

function I=cplxcomp(p1,p2)
%比较两个包含同样标量元素，但可能有不同下标的复数对，返回其数组下标
%这个函数必须用在cplxpair函数之后，以便重新排序频率极点向量及其响应的留数向量
%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';