%-------------------------------------------------------------------------
% 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---------------------------------------