cplxcomp.m

来自「介绍了一些数字信号处理方面的matlab程序」· M 代码 · 共 17 行

function I = cplxcomp(p1,p2)
%  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';