samplecacode.m

来自「利用自己仿真产生的卫星数据」· M 代码 · 共 43 行

function SampSequence=sampleCAcode(svnum,sampleFrequence,sampleTime)

%sample a given C/A sequence with a high sample frequence

%svnum                      is the ID of SV's PRN
%sampleFrequence      is the sample frequence
%sampleTime             is the millisecond number

%Created 4/15/05, Hao Yu, [e-mail] 
%Copy right olinkstar corp. Bejing China

%*********C/A of sat # sampled at sampleFrequence for sampleTime ms

% %use for debug
% clear;clc;
% svnum=12;   sampleFrequence=16367667;   
% sampleTime=2;

UnitTime=1.0;                                                     % time's unit is 1 second;
CACyc=UnitTime/(1023*1000);                          %c/a code period
SampCyc=UnitTime/(sampleFrequence);          %sampled period
totalTime=sampleTime;                                  %total time 

n=totalTime;

catemp=CAcode(svnum);
caseque=catemp;
for counter=2:n,
    caseque=[caseque catemp];
end

SampPoint=fix(sampleFrequence*n/1000);          %the sampled points'number in n ms
SampSequence=zeros(1,SampPoint);              %define the sample sequence

%sample the C/A sequence by a higher frequence
CASNo=1;
for SSNo=1:1:SampPoint,
    if SSNo*SampCyc>CASNo*CACyc,
        CASNo=CASNo+1;
    end    
    SampSequence(SSNo)=caseque(CASNo);
end