linec.rar_微带线阻抗计算仿真_零极点数值仿真

function [ R , P ] = rootsr( X ,eps ) % 求零极点函数 [ R , P ] = rootsr( -2:0.001:2 ,0.000001 ); % X是离散自变量值，Y是离散因变量值，eps是允许的因变量误差 % R是零点值，P是极点值 Y = f(X); Z = double(1./Y); n = length(X); m1 = 0; m2 = 0; X(n+1) = X(n); Y(n+1) = Y(n); Z(n+1) = Z(n); for k = 2 : n %----------------------------------------------------------------- if Y(k-1) * Y(k) <= 0 && ( abs(Y(k-1)) <0.1 ) && ( abs(Y(k)) <0.1 ) m1 = m1 + 1; R(m1) = (X(k-1) + X(k)) / 2; end s = (Y(k) - Y(k-1)) * (Y(k+1) - Y(k)); if ( abs(Y(k)) < eps ) && ( s <= 0 ) m1 = m1 + 1; R(m1) = X(k); end %----------------------------------------------------------------- if ( Z(k-1) * Z(k) <= 0 ) && ( abs(Z(k-1)) <0.1 ) && ( abs(Z(k)) <0.1 ) m2 = m2 + 1; P(m2) = (X(k-1) + X(k)) / 2; end t = (Z(k) - Z(k-1)) * (Z(k+1) - Z(k)); if ( abs(Z(k)) < eps ) && ( t <= 0 ) m2 = m2 + 1; P(m2) = X(k); end %----------------------------------------------------------------- end if m1 == 0 R = 0;disp('Cant find roots.'); end if m2 == 0 P = 0;disp('Cant find poles.'); end end