function ysum=wave_sum(t,ampl,phase,period)
% Вычисление рЯда набором синусоид с заданными периодами
%
%            ysum=wave_sum(t,ampl,phase,period);
%
%     t - вектор абсцисс
%  ampl - вектор амплитуд аппроксимирующих гармоник
% phase - вектор фаз (в тех же единицах, что и t)
%period - вектор периодов аппроксимирующих синусоид.
%         Его длина определЯет кол-во гармоник, единица
%         измерениЯ - та же, что и длЯ t.
%
%  ysum - вектор вычисленных значений
%
%
% Сумма вычислЯетсЯ по формуле SUM{ampl(k)*SIN[f(k)*(t+ph(k))]}

%
% Хомутов С.Ю.,   12.04.96
%
freq=2*pi./period;      % частоты гармоник

for i=1:length(t);
   ysum(i)=sum(ampl.*sin(freq.*(t(i)+phase)));
end
ysum=ysum';
