【正文】 subplot(2,1,1)rez=real(ifft(dz,N))。re=(s_bpska+rez).*cos(2*pi*fs*ts)。plot(ts,re)。axis([,])。xlabel(39。t39。)。title(39。擴(kuò)頻后接收信號(hào)乘以恢復(fù)載波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。碼分復(fù)用擴(kuò)頻后乘以恢復(fù)載波后信號(hào)頻譜39。)。axis([0,5000,0,])。xlabel(39。f/Hz39。)。 45 信號(hào)經(jīng)凱薩爾窗低通濾波器figure(7)subplot(2,1,1)magrel=bs.*magre。plot(fb,magrel(1:N/2)*2/N)。title(39。碼分復(fù)用信號(hào)經(jīng)過凱薩爾窗函數(shù)低通濾波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。經(jīng)過凱薩爾窗函數(shù)濾波后時(shí)域波形39。)。 可以看出，信號(hào)經(jīng)過濾波后的時(shí)域波形前后部分畸變較大，說明用改善的walsh 碼進(jìn)行碼分復(fù)用時(shí)，不同用戶間仍然后一定的碼間干擾。 46 對(duì)濾波后信號(hào)分別用 m1 和 m2 進(jìn)行解擴(kuò)figure(8)subplot(2,1,1)j1=rectpulse(M1,1000)。j2=rectpulse(M2,1000)。yrej1=j1.*yrel。plot(ts,yrej1)。title(39。時(shí)域波形與解擴(kuò)碼 m1 相乘波形39。)。xlabel(39。t/s39。)。subplot(2,1,2)yrej2=j2.*yrel。plot(ts,yrej2)。xlabel(39。t/s39。)。title(39。時(shí)域波形與解擴(kuò)碼 m2 相乘波形39。)。 對(duì)比碼分復(fù)用系統(tǒng)和無碼分復(fù)用系統(tǒng)，可以看出，解擴(kuò)后波形出現(xiàn)了 0 項(xiàng)，每個(gè)有效樣值的持續(xù)時(shí)間比無碼分復(fù)用的系統(tǒng)減少了一半。而且解擴(kuò)信號(hào)時(shí)域波形前后部分的干擾比解擴(kuò)前還要嚴(yán)重，說明碼分復(fù)用引入了碼間干擾。 47 對(duì)兩路信號(hào)分別采樣，判決figure(9)for i=1:1:50 S1(i)=0。S2(i)=0。endfor i=1:1:800 ij=i*500 ss1(i)=yrej1(ij)。ss2(i)=yrej2(ij)。endfor i=1:1:800 if ss1(i) ss1(i)=1。 elseif ss1(i) ss1(i)=1。 else ss1(i)=0。 endendfor i=1:1:800 if ss2(i) ss2(i)=1。 elseif ss2(i) ss2(i)=1。 else ss2(i)=0。 endendn1=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 endsubplot(2,1,1) 48 stem(S1)。title(39。序列 1 判決后的最終信號(hào)39。)。subplot(2,1,2)stem(S2)。title(39。序列 2 判決后的最終信號(hào)39。)。 比較原信息序列，可以看出最終判決信號(hào)與原信號(hào)相同，但是有時(shí)會(huì)出現(xiàn)誤碼。所以，為了防止多用戶帶來的碼間干擾，CDMA 除了進(jìn)行碼分復(fù)用以外還應(yīng)采用其他措施抗干擾。至此，walsh 正交碼的碼分復(fù)用仿真完成。 49 產(chǎn)生隨機(jī)序列 Gold 碼和正交 Gold 碼Gold 碼是由 m 序列派生出的一種偽隨機(jī)碼，它具有類似于 m 序列具有的偽隨機(jī)性質(zhì)，但其同長(zhǎng)度不同序列的數(shù)目比 m 序列的多得多。Gold 碼發(fā)生器框圖：m1m2Gold 碼Gold 碼是由 m 序列的優(yōu)選對(duì)移位模二加構(gòu)成，m 1 和 m2 為同長(zhǎng)度的兩個(gè)不同 m 序列并為優(yōu)選對(duì)。Gold 碼的自相關(guān)特性：Gold 證明過 Gold 碼序列的自相關(guān)函數(shù)的所有非最高峰的取值為三值，如下式所示，其中 p=2n1，為 Gold 碼序列的周期。 ??????ptR21當(dāng) n 為奇數(shù)時(shí)， +1，當(dāng) n 為偶數(shù)且不是 4 的整數(shù)倍時(shí)， 。21??nt 12??nt實(shí)驗(yàn)中，我采用了長(zhǎng)度 n=7 的一對(duì)優(yōu)選對(duì)構(gòu)成 Gold 碼。采用 matlab 里面的gfprimfd(7,39。all39。)命令得到全部的 18 個(gè)本原多項(xiàng)式表達(dá)式，如下圖：之后選取的本原多項(xiàng)式的八進(jìn)制數(shù)表示為 m1（211） ，m 2（217） 。 50 產(chǎn)生 Gold 碼并仿真其自相關(guān)函數(shù)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。endfor 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。endfor i=1:1:m L(i)=L1(i)*L2(i)。 L(i)=1L(i)。 L(i)=*L(i)。endfigure(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)。%計(jì)算相關(guān)系數(shù) b=[b(Nb),b(1:(Nb1))]。endk=0:(N1)。plot(k,R,39。r.39。)。axis([64,191,])。title(39。n=7Gold 碼自相關(guān)函數(shù)39。)。 51 xlabel(39。k39。)?？梢愿鶕?jù) 計(jì)算出，除了在 k=127*n 即序列周期的整數(shù)倍時(shí)自相關(guān)??????ptR21函數(shù)為 1 以外，其余自相關(guān)數(shù)值應(yīng)為三值函數(shù)，, 和。而從圖上看出，Gold 自相關(guān)函數(shù)滿足以上條件，具有尖銳的自相關(guān)特性，說明仿真正確。 52 產(chǎn)生正交 Gold 碼并仿真其互相關(guān)函數(shù)Gold 碼的長(zhǎng)度等于對(duì)應(yīng)的 m 序列的長(zhǎng)度，是奇數(shù)，上述仿真中為 127，因此其互相關(guān)不為 0，屬于準(zhǔn)正交碼?？梢宰C明，若在同一優(yōu)選對(duì)產(chǎn)生的 Gold 碼末尾加一個(gè) 0，則構(gòu)成的偶位 Gold 碼相互正交，其股相關(guān)函數(shù) Rij（0）=0，稱為正交 Gold 碼或偶位 Gold 碼。Rij≠0， 。???,0figure(2)for i=1:1:m Lz1(i)=12*L(i)。endLz1(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。endfor i=1:1:m L4(i)=L2(i)*L3(i)。 Lz2(i)=L4(i)*L1(i)。endLz2(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 碼互相關(guān)函數(shù)39。)。xlabel(39。k39。)。 53 如圖，正交 Gold 碼在碼型上只是比 Gold 碼多了一位（單極性為 0，雙極性為 1） ，使得其成為了真正的隨機(jī)序列，滿足了隨機(jī)序列的 3 個(gè)基本特性。從其互相關(guān)函數(shù)的截選段可以看出，互相關(guān)函數(shù)只在 0 時(shí)取 0，正確。至此，Gold碼仿真完成。實(shí)驗(yàn)心得體會(huì)本次實(shí)驗(yàn)我進(jìn)行了直接序列擴(kuò)頻系統(tǒng)的仿真工作，不僅完成了基本的 DSSS仿真，還在其基礎(chǔ)上增加了窄帶強(qiáng)干擾，碼分多址技術(shù)應(yīng)用，Gold 碼和正交Gold 碼仿真，并第一次使用 simulink 工具仿真了誤碼率。通過實(shí)驗(yàn)，我更詳細(xì)地了解了直接序列擴(kuò)頻系統(tǒng)的工作原理。由于我沒有選修移動(dòng)通信，通信原理書上對(duì)于直擴(kuò)系統(tǒng)講得也不是很詳細(xì)，很多地方我只好自己探索，自己查資料，慢慢編程。尤其是在后來的 Gold 碼方面，書上幾乎一帶而過，我只好查閱了許多文獻(xiàn)來確定優(yōu)選對(duì)的尋找方法。而且 simulink 仿真中，由于沒有任何經(jīng)驗(yàn)，我遇到過許多建模問題和錯(cuò)誤，有時(shí)不得不全部重來。總體來看，本次課設(shè)既更深入學(xué)習(xí)了直擴(kuò)系統(tǒng)，Walsh 函數(shù)碼分多址，Gold 碼，窄帶干擾等知識(shí)點(diǎn)，并通過 matlab 將知識(shí)點(diǎn)化為圖像，更加直觀地掌握了所學(xué)內(nèi)容，還讓我更加熟練地使用了 matlab 和 simulink 工具，收獲頗豐。