function [ cavity_pole ] = estimate_cavitypole( ff_cavpole, abs_cavpole, ang_cavpole, varargin)
%local function for FSS noise budget
[b,a] = invfreqs(abs_cavpole.*exp(1i*ang_cavpole*pi/180),ff_cavpole*2*pi,0,1);
cavity_pole = a(2)/a(1)/2/pi;
if isempty(varargin)
    plotflag = false;
else
    plotflag = varargin{1};
end
if plotflag
h = freqs(b,a,ff_cavpole*2*pi);
mnbode(ff_cavpole,abs_cavpole,ang_cavpole,ff_cavpole,abs(h),angle(h)*180/pi,'Legend',{'measured','fitted'});
['cavity pole : ' num2str(cavity_pole)]
end
end