function sum_wave=wave_val(t,ampl,period,phase)
%Сумма синусоид с заданными периодами, амплитудами и фазами
%
%Вызов:		sum_wave=wave_val(t,ampl,period,phase);
%
%        t - вектор, дл? которого вычисл?етс? сумма синусоид;
%     ampl - вектор амплитуд синусоид размерности n;
%   period - вектор периодов синусоид размерности n (в тех же единицах,
%	     что и t);
%    phase - вектор фаз синусоид размерности n (в тех же единицах,
%	     что и t);
%
% sum_wave - выходной вектор размерности t (в единицах ampl).
%
%См. также: WAVE_FIT

if nargin<4
   error ('!!! Мало аргументов !!!')
end
% ------------------------- Array initilizing
sum_wave=zeros(size(t));
tmp=zeros(size(t));
pi2=2*pi;
for k=1:length(ampl),
    tmp=ampl(k)*sin((pi2./period(k)).*(t+phase(k)));
    sum_wave=sum_wave+tmp;
end
