noisecodeacqure.asv

来自「本人做的GPS并行码捕获的算法.用的是1024点FFT,原始数据可以假如不同功率」· ASV 代码 · 共 53 行

%-------------------------------------------------------------------------
% File:NoiseCodeAcqure.m
% Author:Zhu Jiezuo
% Data: 2006-11-02
% ver:  v1.0
% Description: 研究在噪声环境下的码并行的相关特性
%------------------------------------------------------------------------
clear;
clc;
load AllSate32CA.mat;
SateNum1CA=[AllSate32CA(1,:)*2-1 0];   %选择特定卫星的PN码
Shift10ChipNum1CA=[SateNum1CA(512:1024),SateNum1CA(1:511)];
NoiseDBW(1,:) = wgn(1,1024,3);      %2  DBW的噪声
NoiseDBW(2,:) = wgn(1,1024,6);      %6  DBW的噪声
NoiseDBW(3,:)= wgn(1,1024,10);      %10  DBW的噪声
NoiseDBW(4,:)= wgn(1,1024,14);      %15 DBW的噪声
NoiseDBW(5,:)= wgn(1,1024,18);      %18 DBW的噪声
NoiseDBW(6,:)= wgn(1,1024,19);      %19  DBW的噪声
NoiseDBW(7,:)= wgn(1,1024,12);      %12  DBW的噪声
NoiseDBW(8,:)= wgn(1,1024,15);      %15 DBW的噪声
NoiseDBW(9,:)= wgn(1,1024,18);      %18 DBW的噪声
NoiseDBW(10,:)= wgn(1,1024,19);      %19  DBW的噪声

FftSateNum1CA=fft(SateNum1CA);
ConjFftSateNum1CA=conj(FftSateNum1CA);
for  i=1:6
    NoiseTem=NoiseDBW(i,:);
    NoiseShift10ChipNum1CA(i,:)=Shift10ChipNum1CA+NoiseTem;
    FftShift10ChipNum1CA(i,:)=fft(NoiseShift10ChipNum1CA(i,:));
    
    PrData(i,:)=ConjFftSateNum1CA.*FftShift10ChipNum1CA(i,:);
    IffData(i,:)=ifft(PrData(i,:));
end  
figure(1);
subplot(2,1,1);
plot(abs(IffData(1,:)));
subplot(2,1,2);
plot(abs(IffData(2,:)));
figure(2);
subplot(2,1,1);
plot(abs(IffData(3,:)));
subplot(2,1,2);
plot(abs(IffData(4,:)));
figure(3);
subplot(2,1,1);
plot(abs(IffData(5,:)));
subplot(2,1,2);
plot(abs(IffData(5,:)));
%------------------end of the file---------------------------------------