【正文】 參考文獻[1]查光明．熊賢祚．擴頻通信lM]．西安：西安電子科技大學出版社．1990[2]田日才．擴頻通信[M]．北京：清華大學出版社．2007．[3]孫屹，李妍. MATLAB 通信仿真開發(fā)手冊[M]. 北京:國防工業(yè)工業(yè)出版社，2006,5.[4][J].. 23 [5]樊昌信，張甫翊，[M].北京:高等教育出版社，2001.[6],擴頻通信導論，電子工業(yè)出本社，2006附錄一 程序代碼產生gold序列的程序：r=6。N=2^r1。%移位寄存器的長度和序列的長度s1(1:6)=[1 0 0 0 0 1]。 %initial value 1s2(1:6)=[1 0 0 0 0 0]。 %initial value 1f1=[1 0 0 0 0 1 1]。 %特征多項式f＝x^6+x+1。f2=[1 1 0 0 1 1 1]。%特征多項式f＝x^6+x^5+x^2+x+1。for n=r+1:N。%進行循環(huán)s1(n)=mod(sum(s1(nr:n1).*s1(1:r)),2)。%產生m序列end。for n=r+1:N。%進行循環(huán)s2(n)=mod(sum(s2(nr:n1).*s2(1:r)),2)。%產生m序列end。for n=r+1:N。%進行循環(huán)s=mod(s1+s2,2)。%進行摩爾加gold_sequence_1=s。figure(1)。%產生gold序列stem(gold_sequence_1)。end。Gold序列自相關的分析程序：clear all%先要清除clcr=6。N=2^r1。 s1(1:6)=[1 0 0 0 0 0]。 %initial value 1s2(1:6)=[1 0 0 0 0 0]。 %initial value 1f1=[1 0 0 0 0 1 1]。 %特征多項式f＝x^6+x+1。f2=[1 1 0 0 1 1 1]。 %特征多項式f＝x^6+x^5+x^2+x+1。for n=r+1:N %creat pn 1s1(n)=mod(sum(s1(nr:n1).*f1(1:r)),2)。%產生m序列end%figure(0)。%plot(s1,n)。%grid on。for n=r+1:N %creat pn 2s2(n)=mod(sum(s2(nr:n1).*f2(1:r)),2)。%產生m序列ends=mod(s1+s2,2)。 %creat pn goldgold_sequence_1 = s。figure(1)。stem(gold_sequence_1)。grid on。gold_sequence_1=2*gold_sequence_11。 %變?yōu)殡p極性序列for j=0:N1 s3(j+1)=sum(gold_sequence_1.*[gold_sequence_1(1+j:N),gold_sequence_1(1:j)])/N。%自相關函數分析endj=N+1:N1。%進行循環(huán)rho=[fliplr(s3(2:N)),s3]。figure(2)plot(j,rho)。axis([30 30 1 ])。title(39。第一個gold序列的自相關函數39。)%y1=xcorr(gold_sequence_1)。%t=1:1:120。%figure(2)。%plot(y1)。%axis([0,120,1,40])。%gridGold序列互相關的程序：clear all%先要清除r=6。N=2^r1。 %移位寄存器的長度和序列的長度s1(1:6)=[1 0 0 0 0 0]。 %initial value 1s2(1:6)=[1 0 0 0 0 0]。 %initial value 1f1=[1 0 0 0 0 1 1]。 %特征多項式f＝x^6+x+1。f2=[1 1 0 0 1 1 1]。 %特征多項式f＝x^6+x^5+x^2+x+1。for n=r+1:N %creat pn 1s1(n)=mod(sum(s1(nr:n1).*f1(1:r)),2)。%產生m序列end%figure(0)。%plot(s1,n)。%grid on。for n=r+1:N %creat pn 2s2(n)=mod(sum(s2(nr:n1).*f2(1:r)),2)。%產生m序列ends=mod(s1+s2,2)。 %creat pn goldgold_sequence_1 = s。figure(1)。stem(gold_sequence_1)。grid on。r=6。M=2^r1。 s1(1:6)=[1 0 0 0 1 0]。 %initial value 1s2(1:6)=[1 0 0 0 1 0]。 %initial value 1f1=[1 0 0 0 0 1 1]。 %特征多項式f＝x^6+x+1。f2=[1 1 0 0 1 1 1]。 %特征多項式f＝x^6+x^5+x^2+x+1。for n=r+1:M %creat pn 3s1(n)=mod(sum(s1(nr:n1).*f1(1:r)),2)。%產生m序列end%figure(0)。%plot(s1,n)。%grid on。for n=r+1:M %creat pn 4s2(n)=mod(sum(s2(nr:n1).*f2(1:r)),2)。%產生m序列ends=mod(s1+s2,2)。 %creat pn goldgold_sequence_2 = s。figure(2)。stem(gold_sequence_2)。grid on。gold_sequence_1=2*gold_sequence_11。%變?yōu)殡p極性序列gold_sequence_2=2*gold_sequence_21。%變?yōu)殡p極性序列for j=0:M1 s3(j+1)=sum(gold_sequence_1.*[gold_sequence_2(1+j:M),gold_sequence_2(1:j)])/M。%互相關函數分析endj=M+1:M1。rho=[fliplr(s3(2:M)),s3]。%圖形figure(3)plot(j,rho)。axis([30 30 1 ])。title(39。第一個gold序列的互相關函數39。)%y1=xcorr(gold_sequence_1)。%t=1:1:120。%figure(3)。%plot(y1)。%axis([0,120,1,40])。%grid30