【正文】 e=(s_bpska+rez).*cos(2*pi*fs*ts)。 plot(ts,re)。 axis([,])。 xlabel(39。t39。)。 title(39。擴頻后接收信號乘以恢復載波 39。)。 xlabel(39。t/s39。)。 subplot(2,1,2) yre=fft(re,N)。 magre=abs(yre)。 plot(fb,magre(1:N/2)*2/N)。 title(39。碼分復用擴頻后乘以恢復載波后信號頻譜 39。)。 axis([0,5000,0,])。 xlabel(39。f/Hz39。)。 44 信號經凱薩爾窗低通濾波器 figure(7) subplot(2,1,1) magrel=bs.*magre。 plot(fb,magrel(1:N/2)*2/N)。 title(39。碼分復用信號經過凱薩爾窗函數低通濾波 39。)。 axis([0,200,0,])。 xlabel(39。f/Hz39。)。 subplot(2,1,2) tm=(1:N)/N*4。 yrel=real(ifft(bs.*yre,400000))。 plot(tm,yrel)。 xlabel(39。t/s39。)。 title(39。經過凱薩爾窗函數濾波后時域波形 39。)。 可以看出，信號經過濾波后的時域波形前后部分畸變較大，說明用改善的walsh 碼進行碼分復用時，不同用戶間仍然后一定的碼間干擾。 45 對濾波后信號分別用 m1 和 m2 進行解擴 figure(8) subplot(2,1,1) j1=rectpulse(M1,1000)。 j2=rectpulse(M2,1000)。 yrej1=j1.*yrel。 plot(ts,yrej1)。 title(39。時域波形與解擴碼 m1 相乘波形 39。)。 xlabel(39。t/s39。)。 subplot(2,1,2) yrej2=j2.*yrel。 plot(ts,yrej2)。 xlabel(39。t/s39。)。 title(39。時域波形與解擴碼 m2 相乘波形 39。)。 對比碼分復用系統(tǒng)和無碼分復用系統(tǒng)，可以看出，解擴后波形出現了 0 項，每個有效樣值的持續(xù)時間比無碼分 復用的系統(tǒng)減少了一半。而且解擴信號時域波形前后部分的干擾比解擴前還要嚴重，說明碼分復用引入了碼間干擾。 46 對兩路信號分別采樣，判決 figure(9) for i=1:1:50 S1(i)=0。S2(i)=0。 end for i=1:1:800 ij=i*500 ss1(i)=yrej1(ij)。ss2(i)=yrej2(ij)。 end for i=1:1:800 if ss1(i) ss1(i)=1。 elseif ss1(i) ss1(i)=1。 else ss1(i)=0。 end end for i=1:1:800 if ss2(i) ss2(i)=1。 elseif ss2(i) ss2(i)=1。 else ss2(i)=0。 end end n1=1。n2=1。 for ii=1:1:100 i=ii*87。 k1=ss1(i)+ss1(i+1)+ss1(i+2)+ss1(i+3)+ss1(i+4)+ss1(i+5)+ss1(i+6)+ss1(i+7)。 k2=ss2(i)+ss2(i+1)+ss2(i+2)+ss2(i+3)+ss2(i+4)+ss2(i+5)+ss2(i+6)+ss2(i+7)。 if k1=5 S1(n1)=0。n1=n1+1。 elseif k1=5 S1(n1)=1。n1=n1+1。 end if k2=5 S2(n2)=0。n2=n2+1。 elseif k2=5 S2(n2)=1。n2=n2+1。 end end subplot(2,1,1) 47 stem(S1)。 title(39。序列 1 判決后的最終信號 39。)。 subplot(2,1,2) stem(S2)。 title(39。序列 2 判決后的最終信號 39。)。 比較原信息序列，可以看出最終判決信號與原信號相同，但是有時會出現誤碼。所以，為了防止多用戶帶來的碼間干擾， CDMA 除了進行碼分復用以外還應采用其他措施抗干擾。至此， walsh 正交碼的碼分復用仿真完成。 48 八、 產生隨 機序列 Gold 碼和正交 Gold 碼 Gold 碼是由 m 序列派生出的一種偽隨機碼，它具有類似于 m 序列具有的偽隨機性質，但其同長度不同序列的數目比 m 序列的多得多。 Gold 碼發(fā)生器框圖： Gold 碼是由 m 序列的優(yōu)選對移位模二加構成， m1和 m2 為同長度的兩個不同m 序列并為優(yōu)選對。 Gold 碼的自相關特性： Gold 證明過 Gold 碼序列的自相關函數的所有非最高峰的取值為三值，如下式所示，其中 p=2n1，為 Gold 碼序列的周期。 ?????????????ptptpR21 當 n 為奇數時， 212?? nt +1，當 n 為偶數且不是 4 的整數倍時， 12 21 ?? ?nt 。 實驗中，我采用了長度 n=7 的一對優(yōu)選對構成 Gold 碼。采用 matlab 里面的gfprimfd(7,39。all39。)命令得到全部的 18 個本原多項式表達式，如下圖： 之后選取的本原多項式的八進制數表示為 m1（ 211）， m2（ 217）。 m1 m2 Gold 碼 49 產生 Gold 碼并仿真其自相關函數 clear all。 clc。 X1=1。X2=1。X3=1。X4=1。X5=1。X6=1。X7=1。 m=2^71。 for i=1:1:m Y7=X7。Y6=X6。Y5=X5。Y4=X4。Y3=X3。Y2=X2。Y1=X1。 X7=Y6。X6=Y5。X5=Y4。X4=Y3。X3=Y2。X2=Y1。 X1=Y3*Y7。 L1(i)=Y1。 end for i=1:1:m Y7=X7。Y6=X6。Y5=X5。Y4=X4。Y3=X3。Y2=X2。Y1=X1。 X7=Y6。X6=Y5。X5=Y4。X4=Y3。X3=Y2。X2=Y1。 X1=Y1*Y2*Y3*Y7。 L2(i)=Y1。 end for i=1:1:m L(i)=L1(i)*L2(i)。 L(i)=1L(i)。 L(i)=*L(i)。 end figure(1) subplot(2,1,1) i=1:127。 stairs(i,L)。 axis([0,130,])。 title(39。n=7Gold 碼 39。)。 xlabel(39。k39。)。 subplot(2,1,2) a=L。b=a。Nb=127。N=2*m。 for k=1:N c=xor(a,b)。 D=sum(c)。 A=NbD。 R(k)=(AD)/(A+D)。%計算相關系數 b=[b(Nb),b(1:(Nb1))]。 end k=0:(N1)。 plot(k,R,39。r.39。)。 axis([64,191,])。 title(39。n=7Gold 碼自相關函數 39。)。 50 xlabel(39。k39。)。 可以根據?????????????ptptpR21計算出，除了在 k=127*n 即序列周期的整數倍時自相關函數為 1 以外，其余自相關數值應為三值函數，取值分別為 , 和。而從圖上看出， Gold 自相關函數滿足以上條件，具有尖銳的自相關特性，說明仿真正確。 51 產生正交 Gold 碼并仿真其互相關函數 Gold 碼的長度等于對應的 m 序列的長度，是奇數，上述仿真中為 127，因此其互相關不為 0，屬于準正交碼。可以證明，若在同一優(yōu)選對產生的 Gold 碼末尾加一個 0，則構成的偶位 Gold 碼相互正交，其股相關函數 Rij（ 0） =0，稱為正交 Gold 碼或偶位 Gold 碼。 Rij≠0， ?? ?? ,0 。 figure(2) for i=1:1:m Lz1(i)=12*L(i)。 end Lz1(m+1)=1。 X1=1。X2=1。X3=1。X4=1。X5=1。X6=1。X7=1。 for i=1:1:m Y7=X7。Y6=X6。Y5=X5。Y4=X4。Y3=X3。Y2=X2。Y1=X1。 X7=Y6。X6=Y5。X5=Y4。X4=Y3。X3=Y2。X2=Y1。 X1=Y1*Y2*Y3*Y7。 L3(i)=Y1。 end for i=1:1:m L4(i)=L2(i)*L3(i)。 Lz2(i)=L4(i)*L1(i)。 end Lz2(m+1)=1。 subplot(2,1,1) i=1:128。 stairs(i,Lz1)。 axis([0,130,])。 title(39。n=7 正交 Gold 碼 39。)。 xlabel(39。k39。)。 subplot(2,1,2) k=128。xk=fft(Lz1,2*k)。yk=fft(Lz2,2*k)。rm=real(ifft(conj(xk).*yk))。 rm=ifftshift(rm)。 m=k:k1。 stem(m,rm,39。r.39。)。 grid on。 axis([10,10,20,20])。 title(39。n=7 正交 Gold 碼互相關函數 39。)。 xlabel(39。k39。)。 52 如圖，正交 Gold 碼在碼型上只是比 Gold 碼多了一位（單極性為 0，雙極性為 1），使得其成為了真正的隨機序列，滿 足了隨機序列的 3 個基本特性。從其互相關函數的截選段可以看出，互相關函數只在 0 時取 0，正確。至此， Gold碼仿真完成。 九、 實驗心得體會 本次實驗我進行了直接序列擴頻系統(tǒng)的仿真工作，不僅完成了基本的 DSSS仿真，還在其基礎上增加了窄帶強干擾，碼分多址技術應用， Gold 碼和正交 Gold碼仿真，并第一次使用 simulink 工具仿真了誤碼率。通過實驗，我更詳細地了解了直接序列擴頻系統(tǒng)的工作原理。由于我沒有選修移動通信，通信原理書上對于直擴系統(tǒng)講得也不是很詳細，很多地方我只好自己探索，自己查資料，慢慢編程。尤其是在后來 的 Gold 碼方面，書上幾乎一帶而過，我只好查閱了許多文獻來確定優(yōu)選對的尋找方法。而且 simulink 仿真中，由于沒有任何經驗，我遇到過許多建模問題和錯誤，有時不得不全部重來。總體來看，本次課設既更深入學習了直擴系統(tǒng)， Walsh 函數碼分多址， Gold 碼，窄帶干擾等知識點，并通過 matlab 將知識點化為圖像，更加直觀地掌握了所學內容，還讓我更加熟練地使用了 matlab和 simulink 工具，收獲頗豐。