北斗B1C B2a伪码生成的MATLAB代码

北斗B1C、B2a导频路和数据路伪码生成代码 generateB1CDcode.m B1C数据路主码 generateB1CPcode.m B1C导频路主码 generateB1COcode.m B1C导频路子码 generateB2aDcode.m B2a数据路主码 generateB2aPcode.m B2a导频路主码 generateB2aOcode.m B2a导频路子码

function code = generateB1CDcode（prn）
% B1C Data code
% Phase difference （w）
w=[2678 4802 958  859  3843 2232 124  4352 1816 1126 1860 4800 2267 424  4192 ...
   4333 2656 4148 243  1330 1593 1470 882  3202 5095 2546 1733 4795 4577 1627 ... 
   3638 2553 3646 1087 1843 216  2245 726  1966 670  4130 53   4830 182  2181 ...
   2006 1080 2288 2027 271  915  497  139  3693 2054 4342 3342 2592 1007 310 ...
   4203 455  4318];
% Truncation point （p） 
p=[699  694  7318 2127 715  6682 7850 5495 1162 7682 6792 9973 6596 2092 19 ...   
  10151 6297 5766 2359 7136 1706 2128 6827 693  9729 1620 6805 534  712  1929 ... 
   5355 6139 6339 1470 6867 7851 1162 7659 1156 2672 6043 2862 180  2663 6940 ... 
   1645 1582 951  6878 7701 1823 2391 2606 822  6403 239  442  6769 2560 2502 ... 
   5072 7268 341];

L=zeros（110243）;
for x=1:5121
    L（1+mod（x*x10243））=1;
end

code=zeros（110230）;
for n=0:10229
    k=mod（n+p（prn）-110243）;
    k_w=mod（k+w（prn）10243）;
    code（n+1）=mod（L（k+1）+L（k_w+1）2）;
end

% Convert 1 to -1 0 to 1
for i =1:length（code）
    if code（i） == 1
        code（i） = -1;
    else
        code（i） = 1;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1174  2018-01-28 17:20  generateB1COcode.m

     文件       1179  2018-01-28 17:11  generateB1CPcode.m

     文件       2513  2018-01-28 17:09  generateB2aDcode.m

     文件       1028  2018-01-28 17:32  generateB2aOcode.m

     文件       2512  2018-01-28 17:16  generateB2aPcode.m

     文件       1165  2018-01-28 14:08  generateB1CDcode.m

----------- ---------  ---------- -----  ----

                 9571                    6