物理層仿真 實驗

1、物理層仿真實驗系別：通信工程系專業(yè)：通信工程系x級學(xué)號：姓名： 呂XX實驗時間：2014年6月30日撰寫日期：2014年7月3日實驗一：隨機信號的產(chǎn)生1.1實驗?zāi)康恼莆崭鞣N偽隨機序列的產(chǎn)生方法。1.2實驗原理1.2.1 Wichmann-Hill算法產(chǎn)生均勻分布隨機變量該算法是通過將3個周期相近的隨機數(shù)發(fā)生器產(chǎn)生的數(shù)據(jù)序列進行相加，進而得到更大的周期的數(shù)據(jù)序列。定義三個隨機數(shù)發(fā)生器：以上三式中均需要設(shè)定一初始值，這三個初始值一般稱為種子。產(chǎn)生的三個序列的周期分別是：30269、30307、30323。將這三個序列組合相加即可得到一個周期更大的均勻分布隨機序列：1.2.2 逆變換發(fā)產(chǎn)生Rayle

2、igh分布隨機變量逆變換法的基本思想如圖1.1所示，條件是產(chǎn)生的隨機變量的分布函數(shù)具有閉合表達式。圖1.1 將一個不相關(guān)均勻分布的隨機序列U映射到一個具有概率分布函數(shù)Fx(x)的不相關(guān)序列隨機序列XRayleigh分布的分布函數(shù)：根據(jù)上面的逆變換方法有：因此，根據(jù)上式即可將均勻分布的隨機變量映射為Rayleigh分布的隨機變量。1.2.3 根據(jù)Rayleigh分布隨機變量產(chǎn)生Gaussian分布隨機變量基于Rayleigh隨機變量，可以方便的產(chǎn)生Gaussian分布隨機變量。關(guān)系如下：或者其中u1和u2分別是兩個（01）之間均勻分布的隨機變量，產(chǎn)生的X和Y均為高斯隨機變量。1.3 實驗內(nèi)容根據(jù)

3、實驗原理部分，完成以下實驗任務(wù)。任務(wù)一：采用Willmann-Hill算法產(chǎn)生10000個均勻分布的隨機變量，根據(jù)兩組種子做出隨機序列的直方圖（命令hist），兩組種子自行設(shè)定。任務(wù)2：使用逆變換法產(chǎn)生其它分布的隨機變量（Rayleigh）設(shè)定sigma2為0.5、1和2，分別做出Rayleigh隨機變量序列的直方圖，觀察圖形。任務(wù)3：通過Rayleigh分布隨機變量產(chǎn)生Gaussian分布隨機變量，設(shè)定sigma2為1產(chǎn)生標(biāo)準(zhǔn)正態(tài)分布的隨機變量，計算其方差并做直方圖。改變sigma2值觀察直方圖的變化。1.4 實驗結(jié)果與分析任務(wù)一：均勻分布兩組直方圖如下由上圖可看出，產(chǎn)生的隨機序列滿足（01

4、）均勻分布。任務(wù)二：Rayleigh分布直方圖由這三張圖可看出，圖像大致滿足瑞利分布，且方差不同，瑞利分布的橫坐標(biāo)不同。方差越大，圖像在橫坐標(biāo)軸上的跨度越大任務(wù)三：通過瑞利分布隨機變量產(chǎn)生高斯隨機分布變量，方差為1計算其方差得：varx =1.0033；vary =1.0079與所設(shè)定的sigma2=1的值相接近。實驗二：QPSK系統(tǒng)仿真2.1 實驗?zāi)康恼莆蘸唵握{(diào)制方法的基帶仿真實現(xiàn)，以及AWGN信道和Rayleigh平坦衰落信道的建模，并完成在這兩種信道下的誤碼率仿真。2.2 實驗原理QPSK發(fā)射信號的等效基帶格式和星座圖Sk = A(cosk+jsink)其中A為信號幅度，k為調(diào)制相位，取

5、值范圍為(/4, 3/4,5/4,7/4),比特與相位之間的映射關(guān)系如圖2.1所示。圖1.2 QPSK、8PSK、16QAM信號星座圖2.3 實驗內(nèi)容任務(wù)1：根據(jù)參考教材第一章的BPSK示例，采用實驗一產(chǎn)生隨機數(shù)的方法重新改寫程序，并仿真的BER結(jié)果，將結(jié)果與示例仿真結(jié)果進行比較。（提示：主要修改的內(nèi)容有比特的產(chǎn)生、AWGN噪聲的產(chǎn)生，即只要涉及需要隨機數(shù)的地方均改用實驗一中的方法得到）。任務(wù)2：參考BPSK程序，搭建QPSK、8PSK、16QAM的基帶仿真程序，仿真在AWGN信道和Rayleigh平坦衰落信道下的誤碼率性能，比較它們的誤比特率和帶寬效率，解釋它們誤碼率性能差別的原因。 2.4

6、 實驗結(jié)果與分析任務(wù)一：用實驗一的方法產(chǎn)生的隨機數(shù)和用matlab自帶函數(shù)產(chǎn)生的隨機數(shù)經(jīng)過BPSK通信系統(tǒng)仿真的BER結(jié)果的比較。如圖，對這兩種方法產(chǎn)生的隨機數(shù)進行BPSK仿真，其誤碼率幾乎一致。任務(wù)二：BPSK、QPSK、8PSK、16QAM在AWGN信道下的仿真：以上四圖分別是BPSK、QPSK、8PSK、16QAM在偽隨機數(shù)下AWGN信道中的仿真曲線，可見其結(jié)果與理論曲線幾乎重疊，這是在滿足大數(shù)定理下的仿真結(jié)果。BPSK、QPSK、8PSK、16QAM在Rayleigh信道下的仿真曲線：可見各個調(diào)制在AWGN和Rayleigh信道下的仿真結(jié)果幾乎一致。接下來為節(jié)省仿真時間，在比較他們之間

7、的性能差異將直接利用理論值曲線進行比較。BPSK、QPSK、8PSK、16QAM仿真在AWGN信道下的誤碼率性能比較。如上圖所示，在相同信噪比下，BPSK和QPSK的誤碼率性能最好，8PSK 次之，QAM最差，8PSK和QAM的誤碼率性能大致在同一數(shù)量級。但是從帶寬效率看，QAM的帶寬效率最高為4,8PSK的帶寬效率為3，QPSK為2。說明在相同信噪比下，帶寬效率和誤碼率是一對矛盾互補的，它們分別對應(yīng)著通信系統(tǒng)的有效性和可靠性。QPSK、8PSK、6QAM在Rayleigh衰落信道下的誤碼率性能，并比較它們的比特率。如上圖可看出，在Rayleigh平坦衰落信道下BPSK、QPSK和8PSK的性

8、能較在AWGN信道下有較大的減少，要得到相同的誤碼率，Rayleigh信道下要使用比AWGN信道下更大的信噪比。在相同信噪比下，QPSK的誤碼率會比8PSK小，但小的不是很多，而8PSK的帶寬效率比QPSK高，所以在Rayleigh平坦衰落信道下，采用8PSK調(diào)制比QPSK更好。實驗三：直接序列擴頻抗ISI研究3.1 實驗?zāi)康恼莆罩苯有蛄袛U頻技術(shù)(DSSS)原理及基帶仿真實現(xiàn)，理解擴頻碼設(shè)計對直接序列擴頻抗多徑干擾的影響。3.2 實驗原理定義：擴展頻譜技術(shù)一般是指用比信號帶寬寬得多的頻帶寬度來傳輸信息的技術(shù)。直接序列擴頻就是用比信息速率高很多倍的偽隨機噪聲碼（PN）與信號相乘來達到擴展信號的帶

9、寬。上圖是直接序列擴頻的系統(tǒng)框圖，從圖可以看到，與傳統(tǒng)窄帶系統(tǒng)相比主要引入了擴頻和解擴過程(通過擴頻碼實現(xiàn))，廣義來說可以把擴頻/解擴也堪稱一種調(diào)制/解調(diào)。通過擴頻碼的作用窄帶信號變成了寬帶信號。因此，擴頻系統(tǒng)具有良好的抗窄帶干擾和抗多徑(ISI)干擾的能力。下圖展示了DSSS系統(tǒng)抗窄帶干擾的原理，其中S(f)是調(diào)制信號的頻譜，經(jīng)過擴頻進入接收端的信號頻譜為S(f)*Sc(f)，由于擴頻碼具有寬頻特性，因此，信號帶寬被展寬，如圖1.3(b)所示，I(f)表示窄帶干擾的頻譜，其幅度超過有用信號，但經(jīng)過解擴后，如圖1.3(c)所示，由于擴頻碼的相關(guān)性，調(diào)制信號的功率得到恢復(fù)，而干擾信號的頻譜被擴散

10、，達到了抗窄帶干擾的效果。類似的由于擴頻碼的自相關(guān)特性，DSSS系統(tǒng)還有有效的抗ISI干擾，其頻域原理如圖3.3所示，其中H(f)=+e-j2f 表示ISI信道的頻域傳遞函數(shù)。在直接序列擴頻系統(tǒng)中，其擴頻碼的設(shè)計是整個系統(tǒng)的關(guān)鍵，擴頻碼的選擇要求具有好的互相關(guān)性和自相關(guān)性，即要能有效的抗多用戶干擾和抗ISI干擾。為了抗多用戶干擾，需要盡量選擇互相正交的擴頻碼，為了抗ISI干擾，需要擴頻碼的自相關(guān)特性盡量滿足一定的特性。常用的擴頻碼序列有m序列、Golden序列等等。接收機可采用相干接收機、RAKE接收機。3.3 實驗內(nèi)容根據(jù)DSSS技術(shù)原理搭建仿真系統(tǒng)，并完成下面的任務(wù)。任務(wù)一：參考后面提供的

11、程序，采用m序列替代程序中的隨機序列作為擴頻碼，仿真其性能，完成與隨機序列擴頻性能之間的比較；任務(wù)二：在ISI信道下，研究m序列與golden序列抗ISI的性能。3.4 實驗結(jié)果與分析任務(wù)一：比較m序列與隨機序列的擴頻性能高斯信道下采用m序列擴頻和采用隨機序列擴頻，信噪比相同時，它們的誤碼率是大致也是相同的。說明在無碼間串?dāng)_的信道下m序列擴頻與一般隨機序列的性能沒有顯著差別。任務(wù)二：在ISI信道下，m序列與golden序列抗ISI性能比較。如上圖所示，在ISI信道下，采用m序列和golden序列擴頻都能夠減少通信系統(tǒng)的誤碼率。但是golden序列擴頻的性能明顯優(yōu)于m序列的性能。m序列需要在足夠

12、大信噪比下才能得到滿意的誤碼率。4實驗小結(jié)4.1實驗一實驗小結(jié)在通信系統(tǒng)中的各種信號以及噪聲、信道衰落都可以看作相應(yīng)的隨機過程，通信系統(tǒng)仿真就是要通過計算機產(chǎn)生各種隨機信號，并依據(jù)通信原理對這些隨機信號進行相應(yīng)的處理，獲取期望的結(jié)果。這次的實驗中，我掌握了Wichmann-Hill算法，并利用此算法產(chǎn)生的均勻分布的隨機變量變換得到了Rayleigh分布的隨機變量，以及Gussian分布的隨機變量，為后續(xù)實驗奠定了基礎(chǔ)。4.2實驗二實驗小結(jié)在本次實驗中，我利用上一個實驗一產(chǎn)生的偽隨機數(shù)，對BPSK、QPSK、8PSK和16QAM通信系統(tǒng)進行了仿真，并且得到了和理論曲線一致的仿真曲線。在實驗過程中

13、的第一個問題：我一開始采用10000個隨機數(shù)作為信號源，發(fā)現(xiàn)仿真的結(jié)果與理論相距甚遠。后來通過分析，明白在誤碼率達到10(-3)10(-5)時，仿真的結(jié)果在滿足大數(shù)定理的情況下即錯誤個數(shù)盡可能的大，實驗數(shù)據(jù)至少要10(7)數(shù)量級。因此后來的實驗中，我都采用另一種仿真模式，即設(shè)定錯誤限，在每次達到錯誤限之前一直產(chǎn)生偽隨機源和偽隨機噪聲，直到達到大數(shù)定理，才完成仿真。通過設(shè)定錯誤限來控制仿真的曲線穩(wěn)定性。最后得到幾乎與理論曲線一致的仿真曲線。在實驗過程中的第二個問題：在進行仿真的過程中，我發(fā)現(xiàn)16QAM的仿真曲線與理論曲線不一致，雖然仿真曲線和理論曲線在數(shù)量級上一致，但仿真曲線的誤碼率總是高于理論

14、曲線，并與其幾乎平行。這與其他的仿真中仿真曲線幾乎與理論曲線重合的現(xiàn)象不一致，讓我很費解。后來分析，覺得這種現(xiàn)象產(chǎn)生的原因是，16QAM的星座圖不是最好的星座圖，在譯碼過程中，產(chǎn)生錯誤時，譯碼結(jié)果沒有最小限度的減小錯誤的比特，也就是星座圖上相鄰兩個碼點距離過大。通過修改星座圖，我發(fā)現(xiàn)仿真曲線與理論曲線又一致了，下圖是我修改后的星座圖。所有相鄰的碼字距離皆為1。這樣能最小化錯誤的比特數(shù)目。4.3實驗三實驗小結(jié)這個實驗中，我又重新上網(wǎng)搜索了m序列和golden序列的概念，并通過網(wǎng)上的一些描述原理，產(chǎn)生了兩組m序列和golden序列作為擴頻碼，完成了實驗的仿真。通過修改實驗書上的代碼，我對一些mat

15、lab矩陣實驗平臺上工具又有了進一步的了解。比如ISI信道的建模，可以直接通過卷積實現(xiàn)，而且實驗書上的例程是利用錯誤限來控制仿真數(shù)據(jù)盡量的滿足大數(shù)定理。這也給了我很多的啟發(fā)。雖然實驗中的一些原理不是很清晰，但是讓我有了總體的把握，并加深了我對通信系統(tǒng)的理解。附錄：代碼實驗一：lab1.mclc; clear all;%lab1的實驗一%種子設(shè)定x(1)=22;y(1)=21;z(1)=23;u(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%隨機數(shù)發(fā)生器for i=2:10000 x(i)=mod(171*x(i-1),30269); y(i)=mod

16、(170*y(i-1),30307); z(i)=mod(172*z(i-1),30323); u(i)=mod(x(i)/30269+y(i)/30307+z(i)/30323),1);end%種子設(shè)定x(1)=9;y(1)=8;z(1)=10;v(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%隨機數(shù)發(fā)生器for i=2:10000 x(i)=mod(171*x(i-1),30269); y(i)=mod(170*y(i-1),30307); z(i)=mod(172*z(i-1),30323); v(i)=mod(x(i)/30269+y(i)/

17、30307+z(i)/30323),1);end%直方圖figure(1);subplot(1,2,1);hist(u);subplot(1,2,2);hist(v);%實驗二%th為0.5時。th=0.5;r=sqrt(-2*th2*log(u);figure(2)subplot(2,2,1);hist(r);%th為1時。th=1;r=sqrt(-2*th2*log(u);subplot(2,2,2);hist(r);%th為2時。th=2;r=sqrt(-2*th2*log(u);subplot(2,2,3);hist(r);%實驗三%通過不同的種子產(chǎn)生不同的01偽隨機均勻分布的隨機變量

18、x(1)=22;y(1)=21;z(1)=23;u1(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%隨機數(shù)發(fā)生器for i=2:10000 x(i)=mod(171*x(i-1),30269); y(i)=mod(170*y(i-1),30307); z(i)=mod(172*z(i-1),30323); u1(i)=mod(x(i)/30269+y(i)/30307+z(i)/30323),1);endx(1)=11;y(1)=12;z(1)=10;u2(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%隨機

19、數(shù)發(fā)生器for i=2:10000 x(i)=mod(171*x(i-1),30269); y(i)=mod(170*y(i-1),30307); z(i)=mod(172*z(i-1),30323); u2(i)=mod(x(i)/30269+y(i)/30307+z(i)/30323),1);end%產(chǎn)生高斯分布的隨機變量,th=1th=1;X=sqrt(-2*th2*log(u1).*cos(2*pi*u2);%直方圖figure(3)subplot(2,2,1);hist(X);%產(chǎn)生高斯分布的隨機變量,th=0.5th=0.5;X=sqrt(-2*th2*log(u1).*cos(2

20、*pi*u2);subplot(2,2,2);hist(X);%產(chǎn)生高斯分布的隨機變量,th=2th=2;X=sqrt(-2*th2*log(u1).*cos(2*pi*u2);subplot(2,2,3);hist(X);實驗二代碼BPSK系統(tǒng)AWGN信道自帶隨機數(shù)lab2_BPSK_N.mclc;clear all;close all;%實驗一%BPSK仿真-系統(tǒng)自帶隨機數(shù)%EbN0dB=1:10;%設(shè)定信噪比errLimit=50;N=100;for ii=1:length(EbN0dB) errCount=0; totalN=0; sigma2(ii)=1/(2*(10(EbN0dB(

21、ii)/10); while errCounterrLimit bits = randi(1,1,N); %調(diào)制 syms = 1-2*bits; receiver = syms +sqrt(sigma2(ii)*randn(1,N); %解調(diào) for k = 1:N if receiver(k)0 demodata(k)=1; else demodata(k)=0; end end temp = sum(abs(bits-demodata); errCount = errCount + temp; totalN = totalN + N; end ber(ii)=errCount/total

22、N;endfigure(1);semilogy(EbN0dB,ber,r-o);hold on;grid on;%理論性能曲線berthe=berawgn(EbN0dB,psk,2,nodiff);semilogy(EbN0dB,berthe,b-*);axis(0 10 0.000001 1);xlabel(SNR in dB);ylabel(BER);legend(Simulation,Theoretical);BPSK系統(tǒng)AWGN信道偽隨機數(shù)lab2_BPSK_PN.mclc;clear all;close all;%實驗一%BPSK仿真-偽隨機數(shù)%產(chǎn)生高斯信道噪聲X種子x_g1seed

23、=22;y_g1seed=21;z_g1seed=23;x_g2seed=11;y_g2seed=12;z_g2seed=10;%產(chǎn)生偽隨機信源種子x_seed=5;y_seed=7;z_seed=6;EbN0dB=1:10;%設(shè)定信噪比errLimit=50;N=100;for ii=1:length(EbN0dB) errCount=0; totalN=0; sigma2(ii)=1/(2*(10(EbN0dB(ii)/10); while errCounterrLimit %產(chǎn)生偽隨機信源并轉(zhuǎn)換為0和1的形式 x_s(1)=x_seed;y_s(1)=y_seed;z_s(1)=z_se

24、ed;u1_s(1)=mod(x_seed/30269+y_seed/30307+z_seed/30323),1); for i=2:100 x_s(i)=mod(171*x_s(i-1),30269); y_s(i)=mod(170*y_s(i-1),30307); z_s(i)=mod(172*z_s(i-1),30323); u1_s(i)=mod(x_s(i)/30269+y_s(i)/30307+z_s(i)/30323),1); end x_seed=x_s(100);y_seed=y_s(100);z_seed=z_s(100); u1bits=round(u1_s); bits

25、 = u1bits; %調(diào)制 syms = 1-2*bits; %產(chǎn)生偽隨機高斯噪聲 x_g1(1)=x_g1seed;y_g1(1)=y_g1seed;z_g1(1)=z_g1seed;u1_g1(1)=mod(x_g1(1)/30269+y_g1(1)/30307+z_g1(1)/30323),1); for i=2:100 x_g1(i)=mod(171*x_g1(i-1),30269); y_g1(i)=mod(170*y_g1(i-1),30307); z_g1(i)=mod(172*z_g1(i-1),30323); u1_g1(i)=mod(x_g1(i)/30269+y_g1(

26、i)/30307+z_g1(i)/30323),1); end x_g2(1)=x_g2seed;y_g2(1)=y_g2seed;z_g2(1)=z_g2seed;u2_g2(1)=mod(x_g2(1)/30269+y_g2(1)/30307+z_g2(1)/30323),1); for i=2:100 x_g2(i)=mod(171*x_g2(i-1),30269); y_g2(i)=mod(170*y_g2(i-1),30307); z_g2(i)=mod(172*z_g2(i-1),30323); u2_g2(i)=mod(x_g2(i)/30269+y_g2(i)/30307+z_

27、g2(i)/30323),1); end x_g1seed=x_g1(100);y_g1seed=y_g1(100);z_g1seed=z_g1(100); x_g2seed=x_g2(100);y_g2seed=y_g2(100);z_g2seed=z_g2(100); th=1; X=sqrt(-2*th2*log(u1_g1).*cos(2*pi*u2_g2); receiver = syms +sqrt(sigma2(ii)*X; %解調(diào) for k = 1:N if receiver(k)0 demodata(k)=1; else demodata(k)=0; end end tem

28、p = sum(abs(bits-demodata); errCount = errCount + temp; totalN = totalN + N; end ber(ii)=errCount/totalN;endsemilogy(EbN0dB,ber,r-o);hold on;grid on;%理論性能曲線berthe=berawgn(EbN0dB,psk,2,nodiff);semilogy(EbN0dB,berthe,b-*);title(BPSK在AWGN信道下的仿真曲線與理論曲線);axis(0 10 0.000001 1);xlabel(SNR in dB);ylabel(BER

29、);legend(Simulation,Theoretical);QPSK系統(tǒng)AWGN信道偽隨機數(shù)lab2_QPSK_AWGN_PN.mclc;clear all;close all;%實驗一%QPSK仿真-偽隨機數(shù)%產(chǎn)生高斯信道噪聲X種子x_g1seed=22;y_g1seed=21;z_g1seed=23;x_g2seed=11;y_g2seed=12;z_g2seed=10;%產(chǎn)生偽隨機信源種子x_seed=5;y_seed=7;z_seed=6;EbN0dB=1:10;%設(shè)定信噪比errLimit=100;N=100;for ii=1:length(EbN0dB) errCount=0

30、; totalN=0; sigma2(ii)=1/(2*(10(EbN0dB(ii)/10); while errCount=0&ang=pi/2&ang=pi&ang=3*pi/2&ang4*pi/2) demodata(i:i+1)=1 0; end end temp = sum(abs(bits-demodata); errCount = errCount + temp; totalN = totalN + N; end ber(ii)=errCount/totalN;endsemilogy(EbN0dB,ber,r-o);hold on;grid on;%理論性能曲線berthe=be

31、rawgn(EbN0dB,psk,4,nodiff);semilogy(EbN0dB,berthe,b-*);title(QPSK在AWGN信道下的仿真曲線與理論曲線);axis(0 10 0.000001 1);xlabel(SNR in dB);ylabel(BER);legend(Simulation,Theoretical);8PSK系統(tǒng)AWGN信道偽隨機數(shù)lab2_8PSK_AWGN_PN.mclc;clear all;close all;%實驗一%QPSK仿真-偽隨機數(shù)%產(chǎn)生高斯信道噪聲X種子x_g1seed=22;y_g1seed=21;z_g1seed=23;x_g2seed=

32、555;y_g2seed=456;z_g2seed=364;%產(chǎn)生偽隨機信源種子x_seed=5;y_seed=7;z_seed=6;EbN0dB=1:11;%設(shè)定信噪比errLimit=200;N=102;for ii=1:length(EbN0dB) errCount=0; totalN=0; sigma2(ii)=1/(3*(10(EbN0dB(ii)/10); while errCount=0&ang=15*pi/8&ang=13*pi/8&ang=11*pi/8&ang=9*pi/8&ang=7*pi/8&ang=5*pi/8&ang=3*pi/8&ang=pi/8&ang3*pi/

33、8) demo=0 0 1; end demodata(i:i+2)=demo; end temp = sum(abs(bits-demodata); errCount = errCount + temp; totalN = totalN + N; end ber(ii)=errCount/totalN;endsemilogy(EbN0dB,ber,r-o);hold on;grid on;%理論性能曲線berthe=berawgn(EbN0dB,psk,8,nondiff);semilogy(EbN0dB,berthe,b-*);title(8PSK在AWGN信道下的仿真曲線與理論曲線);a

34、xis(0 11 0.000001 1);xlabel(SNR in dB);ylabel(BER);legend(Simulation,Theoretical);16QAM系統(tǒng)AWGN信道偽隨機數(shù)lab2_16QAM_AWGN_PN.mclc;clear all;close all;%實驗一%QPSK仿真-偽隨機數(shù)%產(chǎn)生高斯信道噪聲X種子x_g1seed=29;y_g1seed=21;z_g1seed=23;x_g2seed=555;y_g2seed=456;z_g2seed=364;%產(chǎn)生偽隨機信源種子x_seed=10;y_seed=7;z_seed=6;EbN0dB=1:10;%設(shè)定信

35、噪比errLimit=500;N=100;for ii=1:length(EbN0dB) errCount=0; totalN=0; sigma2(ii)=1/(4*(10(EbN0dB(ii)/10)*10; while errCounterrLimit %產(chǎn)生偽隨機信源并轉(zhuǎn)換為0和1的形式 x_s(1)=x_seed;y_s(1)=y_seed;z_s(1)=z_seed;u1_s(1)=mod(x_seed/30269+y_seed/30307+z_seed/30323),1); for i=2:100 x_s(i)=mod(171*x_s(i-1),30269); y_s(i)=mod(170*y_s(i-1),30307); z_s(i)=mod(172*z_s(i-1),30323); u1_s(i)=mod(x_s(i)/30269+y_s(i)/30307+z_s(i)/30323),1); end x_seed=x_s(100);y_seed=y_s(100);z_seed=z_s(100); u1bits=round(u1_s); bits = u1bits; %編碼調(diào)制 for i=1:4:97 d=bits(i:i+3); if d = 1 1 1 1 xd = -3-3*j; elseif d = 1 1 1 0 xd= -3 - j; e