Доплеровский сдвиг частоты

function [sigBPSK,t]=BPSK(tFirst,nSymbFirst,phaseInit,Fd,N,PRNseq,periodPRNseq,M,amplSig,F0) % функция формирования сигнала с расширением спектра 
%M - количество символов ПСП 
%amplSig - амплитуда сигнала 
%PRNseq - ПСП 
%TSymb - длина символа ПСП 
%Fseq - частота сигнала 
%phaseInit - начальная фаза 
%N - количество отсчетов 
%nSymbFirst - номер символа, с которого начинается формирование сигнала 
% tfirst - начал время смещения внутри символа 
%sigBPSK - вектор отсчетов 
%phase - фаза последнего отсчета 
%tSymb - временное смещение текущего отсчета внутри текущего символа ПСП 
%nSymb - номер текущего символа ПСП 
%t - Вектор моментов взятия отсчетов сигнала. 
Td=1/Fd; 
t=[0:Td:(N-1)*Td]; 
tSymb=tFirst; 
phase=phaseInit; 
nSymb=nSymbFirst; 
TSymb=periodPRNseq/M; 
sigBPSK = zeros (1,N); 
for i=1:N 
if tSymb>TSymb % проверка пересечения границы символа 
tSymb=tSymb-TSymb; % вычисление смещения внутри символа ПСП 
nSymb=nSymb+1; % номер следующего символа ПСП 
if PRNseq(nSymb)==1 %проверка необходимости сдвига фазы 
phase=phase+pi; % сдвиг фазы 
end 
end 
tSymb=tSymb+Td; 
sigBPSK(i)=amplSig*cos(2*pi*F0*t(i)+phase); 
end