1、1、 m 序列產(chǎn)生及自相關和互相關函數(shù)曲線 序列產(chǎn)生及自相關和互相關函數(shù)曲線function PN=makem(x) %m 序列產(chǎn)生函數(shù)ss1=num2str(x); ss2=dec2bin(base2dec(ss1,8)); %先把八進制轉換為十進制，再把十進制轉換為二進制G=2^(length(ss2)-1)-1; %最大周期 sd=[];for j=1:(length(ss2)-2)sd=[sd 0]; endsd

2、=[sd 1]; %寄存器初始狀態(tài) 0 0 0...0 1PN=[]; for j=1:GPN=[PN sd(length(sd))]; %m 序列輸出的第一位onenum=[]; for jj=1:length(ss2)if str2num(ss2(jj))==1onenum=[onenum jj-1]; %存儲二進制反饋系數(shù)里面“1”的位置-1，即進行異或的位置endendtemp=s

3、d(onenum(2));for jj=3:length(onenum) %根據(jù)“1”的位置進行異或運算temp=xor(temp,sd(onenum(jj)));endfor jj=length(ss2)-1:-1:2 %移位（序列后一位值等于前一位值）sd(jj)=sd(jj-1);endsd(1)=temp; %序列第一位等于反饋出來的值 en

4、dfunction mandzi(ss) %m 序列曲線及自相關函數(shù)曲線繪圖函數(shù)ss1=num2str(ss); ss2=dec2bin(base2dec(ss1,8)); %轉換為二進制G=2^(length(ss2)-1)-1; %最大周期PN=makem(ss); %調用函數(shù)計算 m 序列pp=(-2).*PN+1; %0→1 1→-12、Rake 接收機仿真 接收機仿真clea

5、r all; clc Numusers=1;Nc=16; %擴頻因子ISI_Length=1; %每徑延時為 ISI_Lengh/2EbN0db=[0:1:30]; %信噪比，單位 dBTlen=8000;

6、 %數(shù)據(jù)長度Bit_Error_Number1=0; %誤比特率初始值 Bit_Error_Number2=0;Bit_Error_Number3=0;power_unitary_factor1=sqrt(6/9); %每徑功率因子power_unitary_factor2=sqrt(2/9); power_unitary_factor3=sqr

7、t(1/9);s_initial=randsrc(1,Tlen); %數(shù)據(jù)源wal2=[1 1;1 -1]; %產(chǎn)生 walsh 矩陣wal4=[wal2 wal2;wal2 wal2*(-1)];wal8=[wal4 wal4;wal4 wal4*(-1)];wal16=[wal8 wal8;wal8 wal8*(-1)];s_spread=z

8、eros(Numusers,Tlen*Nc); %擴頻ray1=zeros(Numusers,2*Tlen*Nc);ray2=zeros(Numusers,2*Tlen*Nc); ray3=zeros(Numusers,2*Tlen*Nc); for i=1:Numusersx0=s_initial(i,:).'*wal16(8,:);x1=x0.';s_spread(i,:)=(x1(:)).&#

9、39;; end%將每個擴頻后的輸出重復為兩次，有利于下面的延遲（延遲半個碼元） ray1(1:2:2*Tlen*Nc-1)=s_spread(1:Tlen*Nc);ray1(2:2:2*Tlen*Nc)=ray1(1:2:2*Tlen*Nc-1);%產(chǎn)生第二徑和第三徑信號 ray2(ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-ISI_Length); ray2(2*ISI_Length+1:2*T