信道編碼與擴頻通信

1、信道編碼與擴頻通信直接擴頻技術(shù)是二十一世紀(jì)通信的主要潮流之一。特別是在移動通信領(lǐng)域，碼分多址已經(jīng)成為公認(rèn)的二十一世紀(jì)蜂窩式移動通信基本特征。故而通過仿真手段討論擴頻通信與信道編碼的關(guān)系，具有相當(dāng)重要的現(xiàn)實意義。從某種意義上說，直接擴頻技術(shù)也可以看作一種信道編碼的分組編碼形式，編譯碼都是通過碼序列分組和同一個PN碼的乘運算實現(xiàn)。應(yīng)用這一思想，可以進(jìn)行二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真。二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真模型均勻隨機數(shù)發(fā)生器差 錯 計 數(shù) 器 1比較累加判決器PN碼1（gold碼）同步序列發(fā)生器二進(jìn)制數(shù)據(jù)源 輸出數(shù)據(jù)1高斯隨機數(shù)發(fā)

2、生器 ±E r n A sin（wn）正弦干擾序列發(fā)生器 ±E r累加判決器PN碼2（gold碼）同步序列發(fā)生器二進(jìn)制數(shù)據(jù)源 輸出數(shù)據(jù)2比較均勻隨機數(shù)發(fā)生器差 錯 計 數(shù) 器 2雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真流程圖開始外部輸入：信噪比（dB）初始化變量：信噪比、雙用戶發(fā)送接收信號變量數(shù)組、噪聲變量數(shù)組、仿真序列長度、循環(huán)步長、數(shù)據(jù)塊長、臨時變量清零： 循環(huán)計數(shù)器和差錯計數(shù)器、數(shù)據(jù)塊長計數(shù)器循環(huán)計數(shù)器置1生成發(fā)送用戶1數(shù)據(jù)塊模塊生成發(fā)送用戶2數(shù)據(jù)塊模塊發(fā)送用戶2數(shù)據(jù)塊擴頻模塊發(fā)送用戶1數(shù)據(jù)塊擴頻模塊循環(huán)計數(shù)器加1生成疊加噪聲與人為正弦干擾信號模塊解擴判決接收用戶

3、2數(shù)據(jù)塊模塊解擴判決接收用戶1數(shù)據(jù)塊模塊接收用戶2數(shù)據(jù)塊差錯計數(shù)模塊接收用戶1數(shù)據(jù)塊差錯計數(shù)模塊循環(huán)計數(shù)器次數(shù)已滿？ N Y顯示：接收用戶1差錯計數(shù)器數(shù)值、接收用戶2差錯計數(shù)器數(shù)值，計算：接收用戶1誤碼率、接收用戶2誤碼率結(jié)束返回值：接收用戶1誤碼率、接收用戶2誤碼率為了簡化二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真模型，使我們把主要精力真正集中于編碼仿真本身，本次仿真用戶數(shù)為2，PN碼采用長度為30的gold碼。如上圖所示，雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真模型中還加入了正弦干擾序列發(fā)生器，以用來測試擴頻技術(shù)對人為正弦干擾的抑制作用。在雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真流程圖

4、中，正弦干擾被視做噪聲的一部分置于疊加噪聲模塊。MATLAB源程序function p1,p2=ss_Pe94_2105(snr_in_dB,A,w)%ss_Pe94_2105.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真%snr_in_dB 信噪比%p 誤碼率%A 正弦干擾幅度%w 正弦干擾角頻率%PN碼采用長度Lc為30的gold碼,Lc=length(gold)=30snr=10(snr_in_dB/10);sgma=1;Eb=2*sgma2*snr;N=105; %仿真序列長度105，運行時間約3分鐘p1=0;p2=0;gold1=0 1 1 1 0 0 0 0 1 0 0 0 0

5、 1 1 0 0 1 0 0 1 0 1 1 1 1 0 0 0 0 0;gold2=0 1 0 1 1 0 1 1 1 0 0 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 1 0 0 0 1;Lc=length(gold1); %Lc=30;pn_seq1=2*gold1-1;pn_seq2=2*gold2-1;E_chip=Eb/Lc;temp=0;data1=0;data2=0;repeated_data1=zeros(1,Lc);repeated_data2=zeros(1,Lc);trans_sig1=zeros(1,Lc);trans_sig2=zeros(1,Lc

6、);noise=zeros(1,Lc);n=1:Lc;interference=zeros(1,Lc);rec_sig=zeros(1,Lc);temp1=zeros(1,Lc);temp2=zeros(1,Lc);decision_variable1=0;decision_variable2=0;decision1=0;decision2=0;num_of_err_1=0;num_of_err_2=0;time=0;for i=1:N temp=rand; if (temp<0.5) data1=-1; else data1=1; end for j=1:Lc repeated_dat

7、a1(j)=data1; end trans_sig1=sqrt(E_chip)*repeated_data1.*pn_seq1; temp=rand; if (temp<0.5) data2=-1; else data2=1; end for j=1:Lc repeated_data2(j)=data2; end trans_sig2=sqrt(E_chip)*repeated_data2.*pn_seq2; noise=sgma*randn(1,Lc); time=time+1, n=(time-1)*Lc+1:time*Lc; interference=A*sin(w*n); re

8、c_sig=trans_sig1+trans_sig2+noise+interference; temp1=rec_sig.*pn_seq1; decision_variable1=sum(temp1); if (decision_variable1<0) decision1=-1; else decision1=1; end if (decision1=data1) num_of_err_1=num_of_err_1+1; end temp2=rec_sig.*pn_seq2; decision_variable2=sum(temp2); if (decision_variable2&

9、lt;0) decision2=-1; else decision2=1; end if (decision2=data2) num_of_err_2=num_of_err_2+1; endendnum_of_err_1,num_of_err_2,p1=num_of_err_1/N;p2=num_of_err_2/N;%ss_Pe94_2105.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真echo on%cd94_2_030_105.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真繪圖%PN碼采用長度Lc為30的gold碼,Lc=length(gold)=30%仿真序列長度105，運行時間

10、約1小時A1=3;A2=7;A3=12;A4=0;%A1 A2 A3 A4 正弦干擾幅度w0=1;%w0 正弦干擾角頻率SNRindB=0:2:30;for i=1:length(SNRindB) smld_err_prb11(i),smld_err_prb12(i)=ss_Pe94_2105(SNRindB(i),A1,w0); smld_err_prb21(i),smld_err_prb22(i)=ss_Pe94_2105(SNRindB(i),A2,w0); smld_err_prb31(i),smld_err_prb32(i)=ss_Pe94_2105(SNRindB(i),A3,w0

11、); smld_err_prb41(i),smld_err_prb42(i)=ss_Pe94_2105(SNRindB(i),A4,w0);end;%繪圖函數(shù)semilogy(SNRindB,smld_err_prb11,'mx-');title('雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真')axis(0,30,10(-5),1)hold onsemilogy(SNRindB,smld_err_prb12,'mo-');semilogy(SNRindB,smld_err_prb21,'gx-');semilogy(SNRind

12、B,smld_err_prb22,'go-');semilogy(SNRindB,smld_err_prb31,'kx-');semilogy(SNRindB,smld_err_prb32,'ko-');semilogy(SNRindB,smld_err_prb41,'rx-');semilogy(SNRindB,smld_err_prb42,'ro-');%繪圖函數(shù)SNRindB2=0:0.1:30;theo_err_prb=zeros(1,length(SNRindB2);for i=1:length(SNRi

13、ndB2) SNR=exp(SNRindB2(i)*log(10)/10); theo_err_prb(i)=Qfunct(sqrt(2*SNR); %Qfunct y=(1/2)*erfc(x/sqrt(2); 理論誤碼率公式endsemilogy(SNRindB2,theo_err_prb,'b:'); %繪圖函數(shù)%cd94_2_030_105.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真繪圖運行結(jié)果下圖（cd94_2_030_105的運行結(jié)果）給出了不同信噪比條件（0dB-30dB）和不同幅度（A1=3、A2=7、A3=12、A4=0）正弦干擾下，發(fā)送100000比特

14、的雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真結(jié)果以及理論值曲線。從仿真結(jié)果來看，無正弦干擾下的蒙特卡羅仿真結(jié)果與理論值曲線吻合程度好的令人驚嘆，仿真中兩用戶在同一信道真正實現(xiàn)了和平共處，因此用蒙特卡羅仿真模型來模擬實際雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)，是可以被接受的；同時可以看出，擴頻技術(shù)對人為的正弦干擾具有抑制能力，并且對于同樣的人為正弦干擾，不同gold碼的抑制干擾能力相對來說也有高下之分，在本次仿真中，用戶1使用的gold碼抑制正弦干擾能力更強。級聯(lián)碼對雙用戶二進(jìn)制雙極性擴頻系統(tǒng)誤碼性能的改善無正弦干擾下的仿真結(jié)果和擴頻理論都已經(jīng)說明，從單個用戶角度來看，雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)和

15、原來的單用戶二進(jìn)制雙極性通信系統(tǒng)是相同的。于是我們可以推論，在雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)中對上文各種信道編碼性能測試結(jié)果，也將和前面的測試結(jié)果相同，即按照誤碼率性能指標(biāo)，信道編碼的最好方案是級聯(lián)碼。另外，考慮到實際應(yīng)用信道編碼時通常采用級聯(lián)碼，本次雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)中的信道編碼仿真采用了級聯(lián)碼。在雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真模型框圖中加入級聯(lián)碼的編譯碼器，就構(gòu)成了測試級聯(lián)碼對雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)誤碼性能改善的蒙特卡羅仿真模型，系統(tǒng)框圖如下：雙用戶二進(jìn)制雙極性擴頻級聯(lián)碼通信系統(tǒng)的蒙特卡羅仿真模型均勻隨機數(shù)發(fā)生器差 錯 計 數(shù) 器 1比較級聯(lián)碼譯碼器級聯(lián)碼

16、編碼器累加判決器PN碼1（gold碼）同步序列發(fā)生器二進(jìn)制數(shù)據(jù)源 數(shù)據(jù)1高斯隨機數(shù)發(fā)生器 ±E r n A sin（wn）正弦干擾序列發(fā)生器 ±E r級聯(lián)碼譯碼器級聯(lián)碼編碼器累加判決器PN碼2（gold碼）同步序列發(fā)生器二進(jìn)制數(shù)據(jù)源 數(shù)據(jù)2比較均勻隨機數(shù)發(fā)生器差 錯 計 數(shù) 器 2雙用戶二進(jìn)制雙極性擴頻級聯(lián)碼通信系統(tǒng)的蒙特卡羅仿真流程圖開始外部輸入：信噪比（dB）初始化變量：信噪比、雙用戶發(fā)送接收信號變量數(shù)組、噪聲變量數(shù)組、仿真序列長度、循環(huán)步長、數(shù)據(jù)塊長、臨時變量清零： 循環(huán)計數(shù)器和差錯計數(shù)器、數(shù)據(jù)塊長計數(shù)器循環(huán)計數(shù)器置1生成發(fā)送用戶2數(shù)據(jù)塊模塊生成發(fā)送用戶1數(shù)據(jù)塊模塊發(fā)

17、送用戶2數(shù)據(jù)塊編碼模塊發(fā)送用戶1數(shù)據(jù)塊編碼模塊發(fā)送用戶2數(shù)據(jù)塊擴頻模塊發(fā)送用戶1數(shù)據(jù)塊擴頻模塊生成疊加噪聲與人為正弦干擾信號模塊循環(huán)計數(shù)器加1解擴判決接收用戶2數(shù)據(jù)塊模塊解擴判決接收用戶1數(shù)據(jù)塊模塊接收用戶2數(shù)據(jù)塊譯碼模塊接收用戶1數(shù)據(jù)塊譯碼模塊接收用戶2數(shù)據(jù)塊差錯計數(shù)模塊接收用戶1數(shù)據(jù)塊差錯計數(shù)模塊 N 循環(huán)計數(shù)器次數(shù)已滿？ Y顯示接收用戶1、2差錯計數(shù)器數(shù)值，分別計算接收用戶1、2誤碼率結(jié)束返回值：接收用戶1誤碼率、接收用戶2誤碼率MATLAB源程序function p1,p2=ss_Pe94_2(snr_in_dB,A,w)%ss_Pe94_2.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡

18、羅仿真函數(shù)%snr_in_dB 信噪比%p 誤碼率%A 正弦干擾幅度%w 正弦干擾角頻率%PN碼采用長度Lc為30的gold碼,Lc=length(gold)=30snr=10(snr_in_dB/10);sgma=1;Eb=2*sgma2*snr;N=16; %Nloop=104; %loopNs=N*loop1*loop2; %仿真序列長度loop*N=16*104，運行時間約5分鐘p1=0;p2=0;gold1=0 1 1 1 0 0 0 0 1 0 0 0 0 1 1 0 0 1 0 0 1 0 1 1 1 1 0 0 0 0 0;gold2=0 1 0 1 1 0 1 1 1 0 0

19、 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 1 0 0 0 1;Lc=length(gold1); %Lc=30;pn_seq1=2*gold1-1;pn_seq2=2*gold2-1;E_chip=Eb/Lc;temp=0;data1=0;data2=0;repeated_data1=zeros(1,Lc);repeated_data2=zeros(1,Lc);trans_sig1=zeros(1,Lc);trans_sig2=zeros(1,Lc);noise=zeros(1,Lc);n=1:Lc;interference=zeros(1,Lc);rec_sig=zero

20、s(1,Lc);temp1=zeros(1,Lc);temp2=zeros(1,Lc);decision_variable1=0;decision_variable2=0;decision1=0;decision2=0;num_of_err_1=0;num_of_err_2=0;time=0;for h=1:loop for i=1:N temp=rand; if (temp<0.5) data1=-1; else data1=1; end for j=1:Lc repeated_data1(j)=data1; end trans_sig1=sqrt(E_chip)*repeated_d

21、ata1.*pn_seq1; temp=rand; if (temp<0.5) data2=-1; else data2=1; end for j=1:Lc repeated_data2(j)=data2; end trans_sig2=sqrt(E_chip)*repeated_data2.*pn_seq2; noise=sgma*randn(1,Lc); time=time+1, n=(time-1)*Lc+1:time*Lc; interference=A*sin(w*n); rec_sig=trans_sig1+trans_sig2+noise+interference; tem

22、p1=rec_sig.*pn_seq1; decision_variable1=sum(temp1); if (decision_variable1<0) decision1=-1; else decision1=1; end if (decision1=data1) num_of_err_1=num_of_err_1+1; end temp2=rec_sig.*pn_seq2; decision_variable2=sum(temp2); if (decision_variable2<0) decision2=-1; else decision2=1; end if (decis

23、ion2=data2) num_of_err_2=num_of_err_2+1; end endendnum_of_err_1,num_of_err_2,p1=num_of_err_1/Ns;p2=num_of_err_2/Ns;%ss_Pe94_2.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真函數(shù)function p1,p2=ss_Pe94_2_cnv_jz_74(snr_in_dB,A,w)%ss_Pe94_2_cnv_jz_74.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真函數(shù)%snr_in_dB 信噪比%p 誤碼率%A 正弦干擾幅度%w 正弦干擾角頻率%PN碼采用長度Lc為30

24、的gold碼,Lc=length(gold)=30snr=10(snr_in_dB/10);sgma=1;Eb=2*sgma2*snr;N=16; %Nloop=104; %loopNs=N*loop; %仿真序列長度Ns=loop*N=16*104，運行時間約15分鐘dsource1=zeros(1,N);dsource2=zeros(1,N);output1_h74=hamming74(dsource1);output2_h74=hamming74(dsource2);G=1 0 1 1 ; 1 1 1 1 ;k=1;output_cnv_encd1=cnv_encd(G,k,output

25、1_h74);output_cnv_encd2=cnv_encd(G,k,output2_h74);channel_output1=zeros(1,length(output_cnv_encd2);channel_output2=zeros(1,length(output_cnv_encd2);decoder_output1,survivor_state1,cumulated_metric1=viterbi(G,k,channel_output1);decoder_output2,survivor_state2,cumulated_metric2=viterbi(G,k,channel_out

26、put2);h4output1,h7output1=h47(decoder_output1);h4output2,h7output2=h47(decoder_output2);p1=0;p2=0;numoferr1=0;numoferr2=0;num_of_err_1=0;num_of_err_2=0;cnv_jz_74_time=0;gold1=0 1 1 1 0 0 0 0 1 0 0 0 0 1 1 0 0 1 0 0 1 0 1 1 1 1 0 0 0 0 0;gold2=0 1 0 1 1 0 1 1 1 0 0 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 1 0 0

27、 0 1;Lc=length(gold1); %Lc=30;pn_seq1=2*gold1-1;pn_seq2=2*gold2-1;E_chip=Eb/Lc;temp=0;data1=0;data2=0;repeated_data1=zeros(1,Lc);repeated_data2=zeros(1,Lc);trans_sig1=zeros(1,Lc);trans_sig2=zeros(1,Lc);noise=zeros(1,Lc);n=1:Lc;interference=zeros(1,Lc);rec_sig=zeros(1,Lc);temp1=zeros(1,Lc);temp2=zero

28、s(1,Lc);decision_variable1=0;decision_variable2=0;decision1=0;decision2=0;for h=1:loop for i=1:N temp=rand; if (temp<0.5) dsource1(i)=0; else dsource1(i)=1; end temp=rand; if (temp<0.5) dsource2(i)=0; else dsource2(i)=1; end end %length(dsource) 28 output1_h74=hamming74(dsource1); output2_h74=

29、hamming74(dsource2); output1_jz28=jiaozhi(output1_h74,28); output2_jz28=jiaozhi(output2_h74,28); output_cnv_encd1=cnv_encd(G,k,output1_jz28); output_cnv_encd2=cnv_encd(G,k,output2_jz28); %length(output_cnv_encd),pause for m=1:length(output_cnv_encd1)%=length(output_cnv_encd2) data1=2*output_cnv_encd

30、1(m)-1; data2=2*output_cnv_encd2(m)-1; for j=1:Lc repeated_data1(j)=data1; repeated_data2(j)=data2; end trans_sig1=sqrt(E_chip)*repeated_data1.*pn_seq1; trans_sig2=sqrt(E_chip)*repeated_data2.*pn_seq2; noise=sgma*randn(1,Lc); cnv_jz_74_time=cnv_jz_74_time+1; n=(cnv_jz_74_time-1)*Lc+1:cnv_jz_74_time*

31、Lc; interference=A*sin(w*n); rec_sig=trans_sig1+trans_sig2+noise+interference;% temp1=rec_sig.*pn_seq1; decision_variable1=sum(temp1); if (decision_variable1<0) decision1=-1; else decision1=1; end if (decision1=data1) num_of_err_1=num_of_err_1+1; end% temp2=rec_sig.*pn_seq2; decision_variable2=su

32、m(temp2); if (decision_variable2<0) decision2=-1; else decision2=1; end if (decision2=data2) num_of_err_2=num_of_err_2+1; end% channel_output1(m)=(decision1+1)/2; channel_output2(m)=(decision2+1)/2; % end decoder_output1,survivor_state1,cumulated_metric1=viterbi(G,k,channel_output1); decoder_outp

33、ut2,survivor_state2,cumulated_metric2=viterbi(G,k,channel_output2); decoder_output1_jjz28=jiejiaozhi(decoder_output1,28); decoder_output2_jjz28=jiejiaozhi(decoder_output2,28); h4output1,h7output1=h47(decoder_output1_jjz28); h4output2,h7output2=h47(decoder_output2_jjz28); for i=1:N%=length(dsource1)=

34、length(dsource2) if(h4output1(i)=dsource1(i) numoferr1=numoferr1+1; end if(h4output2(i)=dsource2(i) numoferr2=numoferr2+1; end endendnum_of_err_1,numoferr1,num_of_err_1-numoferr1,num_of_err_2,numoferr2,num_of_err_2-numoferr2,p1=numoferr1/Ns;p2=numoferr2/Ns;%ss_Pe94_2_cnv_jz_74.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅

35、仿真函數(shù)echo on%cd94_2_0216_jz.m 雙用戶二進(jìn)制雙極性擴頻通信系統(tǒng)的蒙特卡羅仿真繪圖%仿真序列長度Ns=loop*N=16*104，運行時間約12小時 A1=3;A2=7;A3=12;A4=0;%A1 A2 A3 A4 正弦干擾幅度w0=1;%w0 正弦干擾角頻率SNRindB=0:2:16;smld_err_prb110=zeros(1,length(SNRindB);smld_err_prb120=zeros(1,length(SNRindB);smld_err_prb210=zeros(1,length(SNRindB);smld_err_prb220=zeros(

36、1,length(SNRindB);smld_err_prb310=zeros(1,length(SNRindB);smld_err_prb320=zeros(1,length(SNRindB);smld_err_prb410=zeros(1,length(SNRindB);smld_err_prb420=zeros(1,length(SNRindB);smld_err_jz_prb11=zeros(1,length(SNRindB);smld_err_jz_prb12=zeros(1,length(SNRindB);smld_err_jz_prb21=zeros(1,length(SNRin

37、dB);smld_err_jz_prb22=zeros(1,length(SNRindB);smld_err_jz_prb31=zeros(1,length(SNRindB);smld_err_jz_prb32=zeros(1,length(SNRindB);smld_err_jz_prb41=zeros(1,length(SNRindB);smld_err_jz_prb42=zeros(1,length(SNRindB);for i=1:length(SNRindB) i,SNRindB(i), smld_err_prb110(i),smld_err_prb120(i)=ss_Pe94_2(

38、SNRindB(i),A1,w0); smld_err_prb210(i),smld_err_prb220(i)=ss_Pe94_2(SNRindB(i),A2,w0); smld_err_prb310(i),smld_err_prb320(i)=ss_Pe94_2(SNRindB(i),A3,w0); smld_err_prb410(i),smld_err_prb420(i)=ss_Pe94_2(SNRindB(i),A4,w0);end;for i=1:length(SNRindB) i,SNRindB(i), smld_err_jz_prb11(i),smld_err_jz_prb12(

39、i)=ss_Pe94_2_cnv_jz_74(SNRindB(i),A1,w0); smld_err_jz_prb21(i),smld_err_jz_prb22(i)=ss_Pe94_2_cnv_jz_74(SNRindB(i),A2,w0); smld_err_jz_prb31(i),smld_err_jz_prb32(i)=ss_Pe94_2_cnv_jz_74(SNRindB(i),A3,w0); smld_err_jz_prb41(i),smld_err_jz_prb42(i)=ss_Pe94_2_cnv_jz_74(SNRindB(i),A4,w0);end;%繪圖函數(shù)figure;

40、semilogy(SNRindB,smld_err_jz_prb42,'r>:');axis(0,16,1e-5,1);xlabel('Eb/N0 in dB');ylabel('Pe');title('雙用戶級聯(lián)編碼及未編碼擴頻系統(tǒng)仿真誤比特率曲線比較');hold onsemilogy(SNRindB,smld_err_jz_prb11,'m<:');semilogy(SNRindB,smld_err_jz_prb12,'m>:');semilogy(SNRindB,smld_err_jz_prb21,'g<:');semilogy(SNRindB,smld_err_jz_prb22,'g>:');semilogy(SNRindB,smld_err_jz_prb31,'k<:');semilogy(SNRindB,smld_err_jz_prb32,'