noisecodeacqure.m

本人做的GPS并行码捕获的算法.用的是1024点FFT,原始数据可以假如不同功率的噪声.

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

FftSateNum1CA=fft(SateNum1CA);
ConjFftSateNum1CA=conj(FftSateNum1CA);
for  i=1:10
    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(6,:)));
figure(4);
subplot(2,1,1);
plot(abs(IffData(7,:)));
subplot(2,1,2);
plot(abs(IffData(8,:)));
figure(5);
subplot(2,1,1);
plot(abs(IffData(9,:)));
subplot(2,1,2);
plot(abs(IffData(10,:)));

%------------------end of the file---------------------------------------