CA碼生成原理及matlab程序?qū)崿F(xiàn)

1、作業(yè)：用Matlab寫C/A碼生成器程序，并畫生成碼的方波圖。C/A碼生成原理C/A 碼是用 m 序列優(yōu)選對組合形成的 Gold 碼。Gold碼是由兩個長度相同而互相關(guān)極大值為最小的 m 序列逐位模2 相加所得到的碼序列。它是由兩個10 級反饋移位寄存器組合產(chǎn)生的，其產(chǎn)生原理如圖1 所示。圖1 C/A碼生成原理G1發(fā)生器的抽頭號為3和10，G2發(fā)生器的抽頭號為2、3、6、8、9、10；G1發(fā)生器的第10位輸出的數(shù)字即為G1碼，而G2碼是由G2的兩個抽頭的輸出結(jié)果進行模2相加得到。 衛(wèi)星的PRN碼與延時的量是相關(guān)聯(lián)的,對C/A碼來說,每顆衛(wèi)星都有特別的延時,如第1顆GPS衛(wèi)星的G2 抽為2、6,

2、第2顆為3、7,第3 顆為4、8,第4 顆為5、9 等，如圖2所示。通過G2 相位選擇可以產(chǎn)生結(jié)構(gòu)不同的偽隨機碼,從而可以實現(xiàn)不同衛(wèi)星之間的碼分多址技術(shù)與衛(wèi)星識別。圖2 prn序號與G2抽頭、時延對應關(guān)系基于MATLAB的GPS信號實現(xiàn)編寫成“codegen”程序，輸入ca_used=codegen(svnum)，其中svnum為衛(wèi)星號，ca_used為得到的C/A碼序列。程序具體實現(xiàn)流程如下：在程序中定義一個數(shù)組，使得衛(wèi)星號與G2的碼片延時一一對應。gs2=5;6;7;8;17;18;139;140;141;251;252;254;255;256;257;258;469;470;471;47

3、2;473;474;509;512;513;514;515;516;859;860;861;862;定義兩個 11 023 的數(shù)組 g1、g2 用來存放生成的Gold 碼。定義一個全 1 的 10 位數(shù)組，作為移位寄存器，相當于G1、G2 生成模塊的初值均置為全“1”。按原理式來生成兩個 Gold 碼序列。設(shè)定兩個 Gold 碼產(chǎn)生器的反饋抽頭：save1=reg(3)*reg(10); % g1 碼的反饋值；save2=reg(2)*reg(3)*reg(6)*reg(8)*reg(9)*reg(10);%g2 碼的反饋值。在移位寄存器數(shù)組reg 中逐位移動，將末端值reg(10)送入g1（g2），將反饋抽頭的模2 值輸出 save1、save2 送入寄存器首位reg(1)。g1 碼這樣重復運算1 023 次后即可生成。g2 碼在此基礎(chǔ)上要引入與衛(wèi)星號相對應的碼片延時g2s(s)，s 為輸入的衛(wèi)星號。以碼片延時時間作為分界點，將其后半段數(shù)據(jù)前移，前半段數(shù)據(jù)后移，這樣就得到 g2 碼。將產(chǎn)生的 g1 和g2 碼逐點相乘，反向，即得 C/A 碼。程序?qū)崿F(xiàn)結(jié)果分別輸入codegen(22)和code