基于matlab的擴頻通信仿真

1、<p>　　基于MATLAB的擴頻通信仿真</p><p><b>　　仿真原理</b></p><p>　　擴展頻譜通信具有很強的抗干擾性，其多址能力、保密、抗多徑等功能也倍受人們的關(guān)注，被廣泛的應用于軍事通信和民用通信中。擴頻技術(shù)，將信號擴展到很寬的頻譜上，在接收端對擴頻信號進行相關(guān)處理即帶寬壓縮，恢復成窄帶信號。對干擾信號而言，由于與擴頻信號不相關(guān)，則

2、被擴展到一個很寬的頻帶上，使之進入信號通頻帶內(nèi)的干擾功率大大的降低，具有很強的對抗能力。</p><p>　　MATLAB是由美國mathworks公司發(fā)布的主要面對科學計算、可視化以及交互式程序設計的高科技計算環(huán)境。它將數(shù)值分析、矩陣計算、科學數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中，為科學研究、工程設計以及必須進行有效數(shù)值計算的眾多科學領(lǐng)域提供了一種全面的解決方案，

3、并在很大程度上擺脫了傳統(tǒng)非交互式程序設計語言（如C、Fortran）的編輯模式。使用MATLAB進行通信仿真有兩種方式，分別是使用m文件編寫程序和用SIMULINK進行可視化建模。</p><p>　　長期以來,人們總是想法使信號所占頻譜盡量的窄,以充分利用十分寶貴的頻譜資源。為什么要用這樣寬頻帶的信號來傳送信息呢? 簡單的回答就是主要為了通信的安全可靠。擴頻通信的基本特點是傳輸信號所占用的頻帶寬度(W )遠大于

4、原始信息本身實際所需的最小</p><p>　　(有效)帶寬(DF) ,其比值稱為處理增益Gp:</p><p>　　Gp =W /△F (1)</p><p>　　眾所周知,任何信息的有效傳輸都需要一定的頻率寬度,如話音為1. 7kHz～3. 1kHz,電視圖像則寬到數(shù)兆赫。為了充分利用有限的頻率資源,增加通路數(shù)目,人們廣泛選擇不同調(diào)制方式,采用寬頻信道(同

5、軸電纜、微波和光纖等) ,和壓縮頻帶等措施,同時力求使傳輸?shù)拿浇橹袀鬏數(shù)男盘栒加帽M量窄的帶寬。因現(xiàn)今使用的電話、廣播系統(tǒng)中,無論是采用調(diào)幅、調(diào)頻或脈沖編碼調(diào)制制式, Gp值一般都在十多倍范圍內(nèi),統(tǒng)稱為“窄帶通信”。而擴頻通信的Gp值,高達數(shù)百、上千,稱為“寬帶通信”。</p><p>　　擴頻通信的可行性,是從信息論和抗干擾理論的基本公式中引伸而來的。</p><p>　　信息論中關(guān)于信息

6、容量的香農(nóng)( Shannon)公式為:</p><p>　　C =Wlog2 (1 + P /N ) (2)</p><p>　　其中:C——信道容量(用傳輸速率度量) ,W ——信號頻帶寬度, ——信號功率, N ——白噪聲功率。式( 2)說明,在給定的傳輸速率C不變的條件下, 頻帶寬度W和信噪比P /N 是可以互換的。即可通過增加頻帶寬度的方法,在較低的信噪比P /N (S /N

7、)情況下, 傳輸信息。擴展頻譜換取信噪比要求的降低,是擴頻通信的重要特點,并由此為擴頻通信的應用奠定了基礎。</p><p>　　擴頻通信與一般的通信系統(tǒng)相比，主要是在發(fā)射端增加了擴頻調(diào)制，而在接收端增加了擴頻解調(diào)的過程，擴頻通信按照其工作方式不同主要分為直接序列擴頻系統(tǒng)、跳頻擴頻系統(tǒng)、跳時擴頻系統(tǒng)、現(xiàn)行調(diào)頻系統(tǒng)和混合調(diào)頻系統(tǒng)?，F(xiàn)以直接序列擴頻系統(tǒng)為例說明擴頻通信的實現(xiàn)方法。下圖1為直接序列擴頻系統(tǒng)的原理框圖。&

8、lt;/p><p>　　圖1 直接序列擴頻系統(tǒng)原理圖</p><p>　　由直擴序列擴頻系統(tǒng)原理圖可以看出，在發(fā)射端，信源輸出的信號與偽隨機碼產(chǎn)生器產(chǎn)生的偽隨機碼進行模2加，產(chǎn)生一速率與偽隨機碼速率相同的擴頻序列，然后再用擴頻序列去調(diào)制載波，這樣得到已擴頻調(diào)制的射頻信號。在接收端，接受到的擴頻信號經(jīng)高放和混頻后，用與發(fā)射端同步的偽隨機序列對擴頻調(diào)制信號進行相關(guān)解擴，將信號的頻帶恢復為信息序列的

9、頻帶，然后進行解調(diào)，恢復出所傳輸?shù)男畔ⅰ?lt;/p><p>　　2．SIMULINK模型建立</p><p>　　基于MATLAB/Simulink所建立的擴頻通信系統(tǒng)的仿真模型，能夠反映擴頻通信系統(tǒng)的動態(tài)工作過程，可進行波形觀察、品剖分析和性能分析等，同時能根據(jù)研究和設計的需要擴展仿真模型，實現(xiàn)以擴頻通信為基礎的現(xiàn)代通信的模擬仿真，為系統(tǒng)的研究和設計提供強有力的平臺。圖2為基于MATLAB

10、/Simulink的擴頻通信系統(tǒng)的仿真模型。</p><p>　　圖2 基于Simulink 的擴頻通信系統(tǒng)的仿真模型</p><p>　　2.1 隨機整數(shù)發(fā)生器（Random Integer genarator）</p><p>　　隨機整數(shù)發(fā)生器(Randon Integer generator)作為仿真系統(tǒng)的信源，隨機整數(shù)發(fā)生器產(chǎn)生二進制隨機信號，采用時間、初

11、始狀態(tài)可自由設置，從而滿足擴頻通信系統(tǒng)所需信源的要求。</p><p>　　2.2 PN序列發(fā)生器（PN Sequence Generator）</p><p>　　用PN序列發(fā)生器（PN Sequence Generator）產(chǎn)生偽隨機碼，對基帶信號進行擴頻。擴頻過程通過信息碼與PN碼進行雙極性變換后相乘加以實現(xiàn)。解擴過程與擴頻過程相同，即將接收的信號用PN碼進行第二次擴頻處理。<

12、/p><p>　　2.3 最小相移鍵控調(diào)制（M-PSK Modulator Passband）</p><p>　　對基帶信號進行數(shù)字調(diào)制，將頻譜從基帶搬移至頻帶，以實現(xiàn)遠距離傳輸?shù)男Ч?lt;/p><p>　　2.4 加性高斯白噪聲信道（AWGN Channel）</p><p>　　信號通過加性高斯白噪聲信道，以達到對信號的傳輸后進行再進行解擴

13、解調(diào)的目的。改信道信噪比參數(shù)設置為負數(shù)，可驗證擴頻通信系統(tǒng)具有良好的抗噪聲能力。</p><p>　　2.5 最小相移鍵控解調(diào) （M-PSK Demodulator Passband）</p><p>　　對解擴后的信號進行解調(diào)，恢復基帶信號。</p><p>　　2.6 誤碼率分析儀（Error Rate Calculation）</p><p

14、>　　在誤碼率計算中，接收到的信號，由于經(jīng)過擴頻解擴、調(diào)制解調(diào)、相關(guān)統(tǒng)計等處理，會存在一個延遲，在誤碼儀模塊的對話框中要設置一個合適的延遲。</p><p><b>　　2.7 其他設置</b></p><p>　　2.7.1 RELAY 的設置</p><p>　　2.7.1 product的設置</p><p>

15、;　　2.7.3 display 的設置</p><p>　　2.7.4 頻譜分析儀（spectrum scope）的設置</p><p>　　系統(tǒng)中有四個頻譜分析儀器，分別命名為：baseband、mod、SS和DSS。其中，baseband能觀察基帶信號的頻譜；mod能觀察到基帶信號調(diào)制到頻帶后的頻譜；SS能觀察擴頻后的頻譜；DSS能觀察解擴解調(diào)后的接受端輸出信號的頻譜。這四個頻譜分析

16、儀的參數(shù)設置大致一致，其前方的零階保持器的采樣時間略有不同。</p><p><b>　　M文件編寫仿真程序</b></p><p>　　用m文件編寫的MATLBA程序來實現(xiàn)觀察系統(tǒng)在不同信噪比下的誤碼率變化情況。此仿真要實現(xiàn)擴頻碼的產(chǎn)生、發(fā)射機部分、接收機部分、信道仿真部分以及誤碼分析部分，此外還有進行循環(huán)仿真以達到繪出信噪比與誤碼率見的函數(shù)圖。</p>

17、<p>　　3.1 擴頻碼的產(chǎn)生</p><p>　　選取m序列作為仿真系統(tǒng)的擴頻碼，產(chǎn)生函數(shù)如下：</p><p>　　function [mout] = mseq(stg, taps, inidata, n) </p><p><b>　　%變量含義說明 </b></p><p>　　% stg

18、 m序列階數(shù)</p><p>　　% taps 線性移位寄存器的系數(shù)</p><p>　　% inidata 序列的初始化</p><p>　　% n 輸出序列的數(shù)目 </p><p>　　% mout 輸出的m序列</p><p>　　if nargin < 4 &l

19、t;/p><p><b>　　n = 1; </b></p><p><b>　　end </b></p><p>　　mout = zeros(n,2^stg-1); </p><p>　　fpos = zeros(stg,1); </p><p>　　fpos(taps) =

20、 1; </p><p>　　for ii=1:2^stg-1 </p><p>　　mout(1,ii) = inidata(stg); % 輸出數(shù)據(jù)的存儲 </p><p>　　num = mod(inidata*fpos,2); % 反饋數(shù)據(jù)的計算</p><p>　　inidata(2:stg) = inid

21、ata(1:stg-1); % 線形移位寄存器的一次移位</p><p>　　inidata(1) = num; % 返回反饋值 </p><p><b>　　end </b></p><p><b>　　if n > 1 </b></p><p>　　for

22、 ii=2:n </p><p>　　mout(ii,:) = shift(mout(ii-1,:),1,0); </p><p><b>　　end </b></p><p><b>　　end</b></p><p><b>　　3.2 發(fā)射機部分</b></p>

23、;<p>　　發(fā)射機部分包括擴頻和調(diào)制兩部分。</p><p>　　3.2.1 QPSK調(diào)制部分</p><p>　　此處采用QPSK調(diào)制，相比SIMULINK的M-PSK調(diào)制的實現(xiàn)要求要低一些。</p><p>　　function [iout,qout]=qpskmod(paradata,para,nd,ml) </p><p

24、><b>　　%各變量含義說明 </b></p><p>　　% paradata 輸入數(shù)據(jù)</p><p>　　% iout 輸出的實部數(shù)據(jù)</p><p>　　% qout 輸出的虛部數(shù)據(jù)</p><p>　　% para 并行信道數(shù)</p>

25、<p>　　% nd 輸入數(shù)據(jù)個數(shù)</p><p>　　% ml 調(diào)制階數(shù) </p><p>　　m2=ml./2; </p><p>　　paradata2=paradata.*2-1; </p><p>　　count2=0; </p><p>　　for j

26、j=1:nd </p><p>　　isi = zeros(para,1); </p><p>　　isq = zeros(para,1); </p><p>　　for ii = 1 : m2 </p><p>　　isi = isi + 2.^( m2 - ii ) .* paradata2((1:para),ii+count2); &

27、lt;/p><p>　　isq = isq + 2.^( m2 - ii ) .* paradata2((1:para),m2+ii+count2); </p><p><b>　　end </b></p><p>　　iout((1:para),jj)=isi; </p><p>　　qout((1:para),jj)=i

28、sq; </p><p>　　count2=count2+ml; </p><p><b>　　end</b></p><p>　　3.2.2 擴頻部分</p><p>　　function [iout, qout] = spread(idata, qdata, code1) </p><p>

29、<b>　　%變量含義說明 </b></p><p>　　% idata 輸入序列實部 </p><p>　　% qdata 輸入序列虛部 </p><p>　　% iout 輸出序列實部 </p><p>　　% qout 輸出序列虛部</p><p>　　%

30、 code1 擴頻碼序列</p><p>　　switch nargin </p><p>　　case { 0 , 1 } </p><p>　　error('lack of input argument'); </p><p><b>　　case 2 </b></p><

31、p>　　code1 = qdata; </p><p>　　qdata = idata; </p><p><b>　　end </b></p><p>　　[hn,vn] = size(idata); </p><p>　　[hc,vc] = size(code1); </p><p>　

32、　if hn > hc </p><p>　　error('lack of spread code sequences'); </p><p><b>　　end </b></p><p>　　iout = zeros(hn,vn*vc); </p><p>　　qout = zeros(hn,vn

33、*vc); </p><p>　　for ii=1:hn </p><p>　　iout(ii,:) = reshape(rot90(code1(ii,:),3)*idata(ii,:),1,vn*vc); </p><p>　　qout(ii,:) = reshape(rot90(code1(ii,:),3)*qdata(ii,:),1,vn*vc); <

34、/p><p><b>　　end</b></p><p>　　3.3 信道仿真部分</p><p>　　實現(xiàn)高斯白噪聲代碼函數(shù)</p><p>　　function [iout, qout] = comb2(idata, qdata, attn) </p><p><b>　　%各變量含義說

35、明</b></p><p>　　% idata 輸入序列實部</p><p>　　% qdata 輸入序列虛部 </p><p>　　% iout 輸出序列實部 </p><p>　　% qout 輸出序列虛部</p><p>　　% attn 根據(jù)信噪比得到的信號衰減水

36、平 </p><p>　　v = length(idata); </p><p>　　h = length(attn); </p><p>　　iout = zeros(h,v); </p><p>　　qout = zeros(h,v); </p><p>　　for ii=1:h </p><p

37、>　　iout(ii,:) = idata + randn(1,v) * attn(ii); </p><p>　　qout(ii,:) = qdata + randn(1,v) * attn(ii); </p><p><b>　　end </b></p><p><b>　　3.4 接收機部分</b></

38、p><p>　　3.4.1 接收機的解擴部分函數(shù)</p><p>　　function [iout, qout] = despread(idata, qdata, code1) </p><p><b>　　%各變量含義說明</b></p><p>　　% idata 輸入序列實部 </p><p

39、>　　% qdata 輸入序列虛部 </p><p>　　% iout 輸出序列實部 </p><p>　　% qout 輸出序列虛部 </p><p>　　% code1 擴頻碼序列</p><p>　　switch nargin </p><p>　　case { 0 ,

40、 1 } </p><p>　　error('lack of input argument'); </p><p><b>　　case 2 </b></p><p>　　code1 = qdata; </p><p>　　qdata = idata; </p><p><b

41、>　　end </b></p><p>　　[hn,vn] = size(idata); </p><p>　　[hc,vc] = size(code1); </p><p>　　vn = fix(vn/vc); </p><p>　　iout = zeros(hc,vn); </p><

42、p>　　qout = zeros(hc,vn); </p><p>　　for ii=1:hc </p><p>　　iout(ii,:) = rot90(flipud(rot90(reshape(idata(ii,:),vc,vn)))*rot90(code1(ii,:),3)); </p><p>　　qout(ii,:) = rot90(flipu

43、d(rot90(reshape(qdata(ii,:),vc,vn)))*rot90(code1(ii,:),3)); </p><p><b>　　end</b></p><p>　　3.4.2 接收機解調(diào)部分的函數(shù)</p><p>　　function [demodata]=qpskdemod(idata,qdata,para,nd,ml)

44、 </p><p><b>　　%各變量含義說明</b></p><p>　　% idata 輸入數(shù)據(jù)的實部 </p><p>　　% qdata 輸入數(shù)據(jù)的虛部</p><p>　　% demodata 解調(diào)后的數(shù)據(jù)</p><p>　　% para 并行的信道數(shù)&l

45、t;/p><p>　　% nd 輸入數(shù)據(jù)個數(shù)</p><p>　　% ml 調(diào)制階數(shù)</p><p>　　demodata=zeros(para,ml*nd); </p><p>　　demodata((1:para),(1:ml:ml*nd-1))=idata((1:para),(1:nd))>=0; <

46、/p><p>　　demodata((1:para),(2:ml:ml*nd))=qdata((1:para),(1:nd))>=0; </p><p>　　3.5 誤碼計算和繪圖部分</p><p><b>　　%誤碼率分析 </b></p><p>　　noe2 = sum(sum(abs(data-demodat

47、a))); </p><p>　　nod2 = user * nd * ml; </p><p>　　noe = noe + noe2; </p><p>　　nod = nod + nod2; </p><p>　　% fprintf('%d\t%e\n',ii,noe2/nod2); </p>&

48、lt;p><b>　　end </b></p><p><b>　　%數(shù)據(jù)文件 </b></p><p>　　ber = noe / nod; </p><p>　　fprintf('%d\t%d\t%d\t%e\n',ebn0,noe,nod,noe/nod); &

49、lt;/p><p>　　fid = fopen('BER.dat','a'); </p><p>　　fprintf(fid,'%d\t%e\t%f\t%f\t\n',ebn0,noe/nod,noe,nod); </p><p>　　fclose(fid); </p><p>

50、　　err_rate_final(ebn0+6)=ber;</p><p><b>　　end</b></p><p><b>　　%性能仿真圖 </b></p><p><b>　　figure</b></p><p>　　semilogy(SNR,err_rate_fina

51、l,'b-*'); </p><p>　　xlabel('信噪比/dB')</p><p>　　ylabel('誤碼率')</p><p>　　axis([-5,10,0,1])</p><p><b>　　grid on</b></p><p>　

52、　3.6 初始化設置以及仿真循環(huán)設置</p><p><b>　　Main函數(shù)為</b></p><p>　　clear all;</p><p><b>　　clc </b></p><p>　　sr = 256000.0; % 符號速率</p><p&g

53、t;　　ml = 2; % 調(diào)制階數(shù) </p><p>　　br = sr * ml; % 比特速率 </p><p>　　nd = 100; % 符號數(shù)</p><p>　　SNR=-5:1:10; % Eb/No </p><

54、p>　　%**************************擴頻碼初值設定************** </p><p>　　user = 1; % 用戶數(shù)</p><p>　　stage = 3; % 序列階數(shù) </p><p>　　ptap1 = [1 3]; % 第一個

55、線性移位寄存器的系數(shù)</p><p>　　disp('--------------start-------------------');</p><p>　　%******************擴頻碼的產(chǎn)生******************</p><p>　　for ebn0=-5:1:10 % m序列<

56、/p><p>　　code = mseq(stage,ptap1,regi1,user); </p><p>　　code = code * 2 - 1; </p><p>　　clen = length(code); </p><p>　　%******************仿真運算開始****************** </p>

57、;<p>　　nloop = 1000; % 仿真循環(huán)次數(shù)</p><p>　　noe = 0; </p><p>　　nod = 0; </p><p>　　for ii=1:nloop </p><p>　　%******************發(fā)射機*******

58、***********</p><p>　　data = rand(user,nd*ml) > 0.5; </p><p>　　[ich, qch] = qpskmod(data,user,nd,ml); % QPSK 調(diào)制</p><p>　　[ich1,qch1] = spread(ich,qch,code); %

59、 擴頻 </p><p>　　ich2 = ich1; </p><p>　　qch2 = qch1; </p><p>　　%******************接收機****************** </p><p>　　[ich3,qch

60、3] = comb2(ich2,qch2,attn); % 添加高斯白噪聲(AWGN) </p><p>　　sampl = irfn * IPOINT + 1; </p><p>　　ich4 = ich3(:,sampl:IPOINT:IPOINT*nd*clen+sampl-1); </p><p>　　qch4 = qch3(:,sam

61、pl:IPOINT:IPOINT*nd*clen+sampl-1); </p><p>　　[ich5 qch5] = despread(ich4,qch4,code); % 解擴</p><p>　　demodata = qpskdemod(ich5,qch5,user,nd,ml); % QPSK解調(diào)</p><p>　　%*

62、*****************誤碼率分析******************</p><p>　　noe2 = sum(sum(abs(data-demodata))); </p><p>　　nod2 = user * nd * ml; </p><p>　　noe = noe + noe2; </p><p>　　nod = no

63、d + nod2; </p><p>　　%******************數(shù)據(jù)文件****************** </p><p>　　ber = noe / nod; </p><p>　　fprintf('%d\t%d\t%d\t%e\n',ebn0,noe,nod,noe/nod); </p

64、><p>　　fid = fopen('BER.dat','a'); </p><p>　　fprintf(fid,'%d\t%e\t%f\t%f\t\n',ebn0,noe/nod,noe,nod); </p><p>　　fclose(fid); </p><p>　　err

65、_rate_final(ebn0+6)=ber;</p><p>　　%******************性能仿真圖******************</p><p><b>　　figure</b></p><p>　　semilogy(SNR,err_rate_final,'b-*');</p><p

66、>　　xlabel('信噪比/dB')</p><p>　　ylabel('誤碼率')</p><p>　　axis([-5,10,0,1])</p><p><b>　　grid on</b></p><p>　　disp('--------------end-------

67、------------'); </p><p><b>　　仿真結(jié)果與分析</b></p><p>　　Simulink仿真結(jié)果</p><p><b>　　基帶信號的頻譜如下</b></p><p><b>　　經(jīng)調(diào)制后頻譜如下</b></p><

68、;p><b>　　擴頻后頻譜如下</b></p><p>　　接收端恢復后的輸出信號頻譜如下</p><p><b>　　誤碼計算儀的結(jié)果</b></p><p>　　Simulink仿真結(jié)果分析</p><p>　　可以看出，基帶信號的頻率分量主要是在-100Hz到100Hz，這與信號源所設

69、置的頻率相符合?；鶐盘栠M過M-PSK調(diào)制后，正頻域的頻率分量集中在2.9KHz-3.1KHz處，可看出載波頻率為3KHz，這個結(jié)果與M-PSK調(diào)制器所設置的載波頻率一致。頻帶信號經(jīng)過與PN碼相乘擴頻后，觀察其信號的頻譜可以看出，正頻域的頻率分量主要集中在0Hz到20KHz這主瓣中，帶寬擴展了約20倍。此時信號將經(jīng)過模擬信道的AWGN信道模塊，來到接收端。在接收端處，先經(jīng)過解擴（即與擴頻PN碼同步的PN碼相乘），得到解擴后的頻帶信號，再

70、經(jīng)過解調(diào)模塊，恢復成基帶信號。觀察恢復后的信號的頻譜，與信源處的信號基本一致。</p><p>　　再看系統(tǒng)的誤碼率計算分析。 輸出的數(shù)據(jù)是一個n 行（與輸入數(shù)據(jù)數(shù)目相等）3列的矩陣。 第1列是差錯率， 第2列是差錯碼的數(shù)量， 第3列是碼元總數(shù)即前述的n。 可看出，差錯率為0.002，差錯碼的數(shù)量為2，碼元總數(shù)為1000。由此可見，系統(tǒng)具有較好的抗噪聲性能。企抗噪聲原理如圖3所示。</p><

71、p>　　圖3 擴頻系統(tǒng)抗噪聲原理示意圖</p><p>　　4.2.1 m文件的仿真結(jié)果</p><p>　　仿真結(jié)果如圖4和圖5所示。其中，</p><p>　　圖4是信號比和誤碼率的值；圖5 是繪制</p><p><b>　　的函數(shù)關(guān)系圖像。</b></p><p>　　圖4 仿

72、真結(jié)果的值</p><p>　　圖5 信噪比和誤碼率的關(guān)系函數(shù)圖象</p><p>　　4.2.2 m文件的仿真結(jié)果分析</p><p>　　通過觀察可知，信噪比越大，系統(tǒng)的誤碼率越小。這符合通信原理的規(guī)律。從圖中還可得知，在負的信噪比條件下，擴頻通信系統(tǒng)還能有比較小的誤碼率?？芍獢U頻系統(tǒng)的抗干擾性能很好，能在負信噪比條件下工作。這是因為擴頻系統(tǒng)特有的頻譜擴展特點

73、和特有的解擴技術(shù)，能使有用信號淹沒在偽噪聲之中，通過相關(guān)解擴將之還原。</p><p><b>　　實驗總結(jié)</b></p><p>　　此次仿真通過MATLAB的SIMULINK仿真實現(xiàn)了擴頻系統(tǒng)從PN碼產(chǎn)生到擴頻和調(diào)制后發(fā)送到通過加性高斯白噪聲然后接收后解擴和解調(diào)還原的過程。觀察了幾個階段的頻譜，對比了各頻譜之間的變化，分析了擴頻系統(tǒng)抗干擾能力及其抗干擾的原因。&

74、lt;/p><p>　　除此之外，還通過m文件的編程，第一步觀察了不同信噪比條件下擴頻通信系統(tǒng)的誤碼率的大小。分析了擴頻通信系統(tǒng)抗干擾能力的優(yōu)越性。</p><p>　　使用SIMULINK仿真和M文件仿真有各自的好處。SIMULINK是可視化的，有模塊庫的支持，使用十分方便；M文件仿真有很好的計算能力。他們都是通信仿真里面的兩把利劍，有些大的系統(tǒng)可能會用到兩者的混合仿真。所以，熟悉這兩種仿真