【正文】 plot(5,1,2)plot(q,F)。 n6=1。 n6=n6+1。 E(n5)=B(unum,n4).*n3(n6)。n6=1。) n2=128/N。)ylabel(39。 %生成用戶數(shù)據(jù)的波形axis([1 1024 ])title(39。endq=1:1:128*8。 end k=k+1。 u=u+1。 for u=(k1)*128+1:k*128。i=1。 end i=i+1。 data_2(i)=1。 n=fix(da)。for i=1:8。data=dec2base(n,2)。)。 otherwise disp(39。 case 64 A=H64。 case 16 A=H16。 case 4 A=H4。%128階哈達(dá)碼NH128=1*H128。 H128=[H64 H64。H32 NH32]。%32階哈達(dá)碼NH32=1*H32。 H32=[H16 H16。H8 NH8]。%8階哈達(dá)碼NH8=1*H8。 H8=[H4 H4。H2 NH2]。%2階哈達(dá)碼NH2=H2*(1)。 (2) function A=walsh(x) H2=[1 1。39。202。184。194。206。for i=1:32 if send(i)~=receive(i) count=count+1。receiveD。receiveB。nD_2。nB_2。247。253。198。237。238。 end i=i+1。error39。 receiveD=receive_module(N,y,m)。 case 4 m=i+1。 figure(8)。 receiveB=receive_module(N,y,m)。 case 2 m=i+1。 figure(6)。) for i=1:4。197。197。213。211。axis([1 1024 ])title(39。 subplot(2,1,2)。measured39。)axis([1 1024 5 5]) snr=。197。197。228。162。title(39。subplot(2,1,1)。Channel=ChannelA+ChannelB+ChannelC+ChannelD。 end i=i+1。error39。 [nD_2,ChannelD]=send_module(nD,N,m)。 m=i+1。186。208。187。203。189。246。187。 nD=fix(unifrnd(0,255))。nD=39。190。D202。187。211。 [nC_2,ChannelC]=send_module(nC,N,m)。 m=i+1。186。208。187。203。189。246。187。 nC=fix(unifrnd(0,255))。nC=39。190。C202。187。211。 [nB_2,ChannelB]=send_module(nB,N,m)。 m=i+1。186。208。187。203。189。246。187。 nB=fix(unifrnd(0,255))。nB=39。190。B202。187。211。 [nA_2,ChannelA]=send_module(nA,N,m)。 m=i+1。186。208。187。203。189。246。187。 nA=fix(unifrnd(0,255))。nA=39。190。A202。187。211。N=16。N=39。202。235。***********clearclc%N=input(39。186。4. 仿真程序（1）%***********214。在數(shù)字通信系統(tǒng)中，積分器轉(zhuǎn)換為累加器，我們?cè)贑DMA通信原理中找到了理論依據(jù)，通過(guò)求規(guī)格化內(nèi)積，求一個(gè)用戶碼片內(nèi)所有擴(kuò)頻碼累加以后求平均值。（2）Walsh解擴(kuò)我們知道，在現(xiàn)代通信系統(tǒng)中，大多采用數(shù)字通信。%y為接收信號(hào)，wc為載波頻率，t為時(shí)間，x為解調(diào)信號(hào)，即y=x/2+[cos(2*wc*t)*xsin(2*wc*t)*sqrt(1x^2)]/2。這里，我們利用MATLAB強(qiáng)大的運(yùn)算功能，采用符號(hào)函數(shù)進(jìn)行求反運(yùn)算。為了方便計(jì)算，我們舍去了噪聲部分，然后將接收到的信號(hào)直接代入計(jì)算。 于是，利用MATLAB強(qiáng)大的運(yùn)算能力，用數(shù)學(xué)的方法反推出。有遞減的趨勢(shì)，當(dāng)信噪比接近于0時(shí)，系統(tǒng)誤碼率接近于0；當(dāng)信噪比接近負(fù)無(wú)窮時(shí)，系統(tǒng)誤碼率接近于100％。measured39。信噪比計(jì)算是數(shù)學(xué)表達(dá)式為： (6)其中為信碼發(fā)射功率，為噪聲功率， 為信源碼，為信宿碼，計(jì)算結(jié)果單位為dB。在實(shí)際測(cè)量中，往往需要用到帶限高斯白噪聲。誤碼率的數(shù)學(xué)表達(dá)式為： (5)其中S是信碼個(gè)數(shù)，e是誤碼個(gè)數(shù)，E就是誤碼率。擴(kuò)頻序列通過(guò)終端機(jī)和信道后，輸出仍為擴(kuò)頻序列。通常認(rèn)為二進(jìn)制信號(hào)中0和1是以等概率隨機(jī)出現(xiàn)的，所以測(cè)量誤碼率時(shí)最理想的信源應(yīng)是隨機(jī)信號(hào)發(fā)生器。 3仿真結(jié)果分析 誤碼率在數(shù)字通信中，誤碼率是一項(xiàng)主要的性能指標(biāo)。誤碼個(gè)數(shù)39。%統(tǒng)計(jì)誤碼個(gè)數(shù)for i=1:32 if send(i)~=receive(i) count=count+1。receiveC。nD_2]%發(fā)送端用戶原始數(shù)據(jù)receive=[receiveA。nB_2。end 差錯(cuò)計(jì)數(shù)器將4個(gè)用戶的發(fā)送函數(shù)和接收函數(shù)的返回值逐位進(jìn)行比較。)。 otherwise disp(39。 figure(9)。 receiveC=receive_module(N,y,m)。 case 3 m=i+1。 figure(7)。 receiveA=receive_module(N,y,m)。 switch i case 1 m=i+1。end如圖9所示，是接收端信號(hào)解調(diào)過(guò)程中各個(gè)子模塊的波形圖。)。 otherwise disp(39。 figure(9)。 receiveC=receive_module(N,y,m)。 case 3 m=i+1。 figure(7)。 receiveA=receive_module(N,y,m)。 switch i case 1 m=i+1。返回值receiveA表示接收對(duì)用戶A發(fā)送的信號(hào)進(jìn)行解調(diào)后還原的信號(hào)，作為最后差錯(cuò)計(jì)數(shù)器還原比較信號(hào)。接收信號(hào)39。plot(q,y)。) 。%信道的信噪比，單位dBy = awgn(Channel,snr,39。發(fā)射信號(hào)39。plot(q,Channel)。q=1:1:128*8。信道傳輸?shù)脑创a如下：figure(5)。measured39。將用戶A、B、C、D已調(diào)信號(hào)相加后，送入傳輸信道中，用高斯加性白噪聲對(duì)信道中的信號(hào)進(jìn)行干擾。 end i=i+1。error39。 [nD_2,ChannelD]=send_module(nD,N,m)。%產(chǎn)生一個(gè)字節(jié)隨機(jī)信號(hào) m=i+1。)。 case 4 %nD=input(39。 figure(3)。 nC=fix(unifrnd(0,255))。用戶C數(shù)據(jù)nC=39。 [nB_2,ChannelB]=send_module(nB,N,m)。%產(chǎn)生一個(gè)字節(jié)隨機(jī)信號(hào) m=i+1。)。 case 2 %nB=input(39。 figure(1)。 nA=fix(unifrnd(0,255))。用