function [am,ph,tap]=am_ph(x,Ifr,Step,F,tau);
%Вычисление изменений со временем амплитуды и фазы гармоники с заданной
%частотой (скольз?ща? оценка)
%
%Вызов:      [am,ph,[t]]=am_ph(x,Ifr,Step,F,tau); 
%
%    x - исходный вектор;
%  Ifr - длина фрагмента, по которой оцениваетс? амплитуда и фаза;
% Step - шаг между последовательными оценками;
%    F - частота оцениваемой гармоники (д.б. согласована с единицами
%        измерени? tau);
%  tau - дискретность измерений (д.б. согласована с единицами измерени? F);
%
%   am - вектор амплитуд;
%   ph - вектор фаз (в единицах tau);
%    t - вектор моментов оценок (в единицах tau);
%
%ВНИМАНИЕ !!!  1. Врем? не контролируетс?. Пропуск записей будет
%                приводить к смене фазы.
%              2. Первое измерение соответствует моменту t=0.

% Хомутов С.Ю.,  12.11.97

x=x(:);			% делаем столбцы, если были строки
t=(0:length(x)-1)'*tau;	% фиктивное врем? (в единицах tau)
T=1/F;
%-------------------------------------------------------
am=[]; ph=[]; tap=[]; k=1;
i1=1; i2=Ifr;
while i2<length(x),
   [am(k),ph(k)]=wave_fit(t(i1:i2),x(i1:i2),T);
   tap(k)=(Ifr/2+(k-1)*Step)*tau;
   k=k+1;
   i1=(k-1)*Step+1;
   i2=i1+Ifr-1;
end
am=am(:); ph=ph(:); tap=tap(:);