Псевдослучайный код спутников ГЛОНАСС

fc=4*fDat ; %несущая частота

fd=32*fDat; %частота дискретизации

fd_fDat=fd/fDat; %отношение частоты дискретизации к символьной скорости

t=(0:length(Dat)*fd_fDat-1)/fd; %дискретное время

t1=(0:length(codGL)*fd_fDat-1)/fd; %дискретное время

BPSK=cos(2*pi*fc*t + pi*Dat(floor(fDat*t)+1)); %модулированный сигнал

y=BPSK;

%преобразование частоты

x1=y.*cos(2*pi*fc*t);

[b1,a1]=butter(3,fc*2/fd); %фильтр нижних частот с максимально-плоской характеристикой

[b,a]=ellip(10,2,60,fc*2/(32*fDat)); %фильтр нижних частот с эллиптической характеристикой

%фильтрация высокочастотных составляющих

x1=filtfilt(b1,a1,x1);

x=filtfilt(b,a,x);

m=3.20;

%Графика

fig1=figure;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Сигналы спутников ГЛОНАСС

%Несущая

subplot(3,1,1),plot(t,cos(2*pi*fc*t))

xlabel('а','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

%Несущая, модулированная М-последовательностью

subplot(3,1,2),plot(t,BPSK)

xlabel('б','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

%М-последовательность

subplot(3,1,3),stem(t,Dat(floor(fDat*t)+1))

xlabel('в','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Фильтрация демодулированного сигнала ГЛОНАСС фильтром нижних частот с максимально плоской и эллиптической характеристиками

m=5;

fig2=figure;

subplot(3,1,1), plot(t,Dat(floor(fDat*t)+1))

xlabel('а','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

subplot(3,1,2), plot(t,x1)

xlabel('б','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

subplot(3,1,3), plot(t,x)

xlabel('в','FontSize',12,'FontName', 'TimesNewRoman')

grid on

xlim([3.07*10^(-4) m*10^(-4)])

ylim([-1.5 1.5])

%Фильтры

fig3=figure;

freqz(b1,a1,512,32*fDat),

xlabel('а','FontSize',12,'FontName', 'TimesNewRoman')

fig4=figure;

freqz(b,a,512,32*fDat),

xlabel('а','FontSize',12,'FontName', 'TimesNewRoman')

Контрольные вопросы.

1.      Запишите математическую процедуру «сложение по модулю 2».

2.      При каких условиях «сложение по модулю 2» можно заменить умножением?

3.      Есть ли разница между псевдослучайным сигналом и M-последовательностью?

4.      Для чего кодируются сигналы навигационных спутников?

5.      Какая связь между регистрами сдвига, формирующими коды сигналов спутников и образующими полиномами?

6.      Как идентифицируются спутники ГЛОНАСС?

7.      Как понимать термин «отношение сигнал/шум»?