cplxcomp.m

数字信号处理及MATLAB实现 具体的例子与码源 很能帮助初学者学习

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';