《移動(dòng)通信原理與應(yīng)用》實(shí)驗(yàn)報(bào)告

1、重慶交通大學(xué)信息科學(xué)與工程學(xué)院綜合性設(shè)計(jì)性實(shí)驗(yàn)報(bào)告專 業(yè)： 通信工程專業(yè)12級(jí) 學(xué) 號(hào)： 631206040218 姓 名： 柴闖闖 實(shí)驗(yàn)所屬課程： 移動(dòng)通信原理與應(yīng)用 實(shí)驗(yàn)室(中心)： 信息技術(shù)軟件實(shí)驗(yàn)室 指 導(dǎo) 教 師 ： 譚晉 2014年11月教師評(píng)閱意見：簽名： 年 月 日實(shí)驗(yàn)成績(jī)：一、題目擴(kuò)頻通信系統(tǒng)仿真實(shí)驗(yàn)二、仿真要求 傳輸?shù)臄?shù)據(jù)隨機(jī)產(chǎn)生，要求采用頻帶傳輸（DPSK調(diào)制）； 擴(kuò)頻碼要求采用周期為63（或127）的m序列；  仿真從基站發(fā)送數(shù)據(jù)到三個(gè)不同的用戶，各不同用戶分別進(jìn)行數(shù)據(jù)接收；  設(shè)計(jì)三種不同的功率延遲分布，從基站到達(dá)三個(gè)不同的用戶

2、分別經(jīng)過多徑衰落（路徑數(shù)分別為2，3，4）；  三個(gè)用戶接收端分別解出各自的數(shù)據(jù)并與發(fā)送前的數(shù)據(jù)進(jìn)行差錯(cuò)比較。三、仿真方案詳細(xì)設(shè)計(jì)(1)通信系統(tǒng)的總體框圖如下：由上圖可以看出，整個(gè)設(shè)計(jì)由發(fā)送端、信道和接收機(jī)三個(gè)部分組成。 發(fā)射機(jī)原理發(fā)送端首先產(chǎn)生三組用戶數(shù)據(jù)和三組不同的m序列，并用三組m序列分別對(duì)用戶信息進(jìn)行擴(kuò)頻。再將擴(kuò)頻信號(hào)與載波進(jìn)行DPSK調(diào)制，得到高頻的已調(diào)調(diào)信號(hào)并將其送入無線的多徑信道。無線信道信道模擬成無線的多徑多用戶信道，在這個(gè)信道中有三個(gè)用戶進(jìn)行數(shù)據(jù)傳輸，每個(gè)用戶的數(shù)據(jù)分別通過三徑傳輸?shù)竭_(dá)接收端。三徑會(huì)有不同的延時(shí)，衰減。最終，還要將三徑用戶數(shù)據(jù)增加高斯白噪聲

3、。接收機(jī)原理接收端會(huì)接收到有燥的三徑信息的疊加。首先，要對(duì)接收到的三徑信息進(jìn)行解擴(kuò)，分離出三組用戶信息；其次，在將解擴(kuò)后的信息進(jìn)行帶通濾波去除帶外噪聲；最后，分別對(duì)三組用戶信息進(jìn)行解調(diào)得到原始數(shù)據(jù)，在對(duì)接收到的數(shù)據(jù)進(jìn)行誤碼率統(tǒng)計(jì)，得出系統(tǒng)的性能指標(biāo)。(2)功能模塊的詳細(xì)設(shè)計(jì)擴(kuò)頻碼（m序列）的產(chǎn)生擴(kuò)頻碼為偽隨機(jī)碼,可以m序列、Golden序列。本設(shè)計(jì)采用自相關(guān)特性好，互相關(guān)特性較差的m序列，為了節(jié)省運(yùn)算量，我選取了周期為63擴(kuò)頻序列，經(jīng)過計(jì)算易知要產(chǎn)生周期為63的m序列需要長(zhǎng)度為6的反饋系數(shù)，經(jīng)過查找資料得出三組反饋系數(shù)（八進(jìn)制）45、67、75，其對(duì)應(yīng)的二進(jìn)制為1000011、1100111、

4、1101101。并將二進(jìn)制與移位寄存器級(jí)數(shù)對(duì)應(yīng)，以1000011為例，設(shè)初始化各寄存器單元內(nèi)容為1，其具體的寄存器結(jié)構(gòu)圖如下所示：產(chǎn)生m序列的matlab程序如下： %m序列的產(chǎn)生,s為初始狀態(tài)，抽頭系數(shù)為（1000011）2function PN=mseq(s)N=length(s); c=; D=s; for i=1:2N-1 c1=rem(D(N)+D(N-1),2); c=c,D(N);D=c1,D(1:N-1);end c=c*2-1;%變?yōu)?，-1的序列 PN=c;end擴(kuò)頻 擴(kuò)頻的主要思想是每一位數(shù)據(jù)位都擴(kuò)展成長(zhǎng)度為m序列長(zhǎng)的信息，其具體做法是將數(shù)據(jù)

5、信息中的1用m序列代替，而對(duì)于-1用-m序列代替，這樣對(duì)每一個(gè)數(shù)據(jù)位都進(jìn)行擴(kuò)展就實(shí)現(xiàn)了對(duì)原始數(shù)據(jù)的擴(kuò)頻。其結(jié)構(gòu)框圖如下:擴(kuò)頻仿真代碼如下：%擴(kuò)頻function expandsignal=expand(signal,PN)expand=;for i=1:length(signal) if signal(i)=1 expand=expand,PN; else expand=expand,-1*PN; endendexpandsignal=expand;%擴(kuò)頻后的序列end多徑信道的仿真對(duì)多徑信道的仿真主要通過對(duì)不同用戶產(chǎn)生的數(shù)據(jù)經(jīng)三徑進(jìn)行傳輸，每一徑信號(hào)又有三個(gè)用戶數(shù)據(jù)疊加而成，并且每一徑信號(hào)具

6、有不同的延遲和衰減。最后在信道的另一端三徑信號(hào)進(jìn)行疊加并對(duì)其加燥。多徑信道仿真代碼如下：function multiS=channels(modusignal,k) a=length(modusignal);signal1=0,modusignal;signal2=0,0,modusignal;signal3=0,0,0,modusignal;r1=0.01*abs(randn(1,a)+j*randn(1,a);r2=0.1*abs(randn(1,a+1)+j*randn(1,a+1);r3=0.2*abs(randn(1,a+2)+j*randn(1,a+2);r4=0.5*abs(ra

7、ndn(1,a+3)+j*randn(1,a+3);if k=2 r=modusignal.*r1,0,0,0+signal1.*r2,0,0;endif k=3 r=modusignal.*r1,0,0,0+signal1.*r2,0,0+signal2.*r3,0;end if k=4r=modusignal.*r1,0,0,0+signal1.*r2,0,0+signal2.*r3,0+signal3.*r4; end multiS=r;end解擴(kuò)在擴(kuò)頻碼與接收信號(hào)同步的情況下, 可對(duì)接收信號(hào)進(jìn)行解擴(kuò), 其解擴(kuò)的過程與擴(kuò)頻的過程一樣, 也是利用擴(kuò)頻碼與接收信

8、號(hào)進(jìn)行相乘即可。其原理圖如下：解擴(kuò)代碼如下：function deexpandsignal=deexpand(signal,k)expandsignal1=expand(round(rand(1,10),mseq(1,0,0,1,1,0);expandsignal2=expand(round(rand(1,10),mseq(1,1,0,1,1,0);expandsignal3=expand(round(rand(1,10),mseq(1,0,0,1,1,1);if k=2 for i=1:length(expandsignal1) deexpPN1(i-1)*100+1:i*100)=expa

9、ndsignal1(i); end for i=1:length(signal)-3 deexpsignal1(i)=signal(i).*deexpPN1(i); end deexpandsignal=deexpsignal1;endif k=3 for i=1:length(expandsignal1) deexpPN2(i-1)*100+1:i*100)=expandsignal2(i); end for i=1:length(signal)-3 deexpsignal2(i)=signal(i).*deexpPN2(i); end deexpandsignal=deexpsignal2

10、;endif k=4 for i=1:length(expandsignal1) deexpPN3(i-1)*100+1:i*100)=expandsignal3(i); end for i=1:length(signal)-3 deexpsignal3(i)=signal(i).*deexpPN3(i); end deexpandsignal=deexpsignal3;endEnd調(diào)制與解調(diào) 為了使低頻信號(hào)能夠在高頻中傳輸，并且增強(qiáng)系統(tǒng)的抗噪聲性能，我們必須采用一定的調(diào)制解調(diào)技術(shù)，這里我采用DPSK對(duì)擴(kuò)頻信號(hào)進(jìn)行調(diào)試，并在接收端對(duì)解擴(kuò)信號(hào)進(jìn)行DPSK解調(diào)以恢復(fù)出原有信號(hào)。DPSK是

11、利用前后相鄰碼元的載波相對(duì)相位變化傳遞數(shù)字信息。DPSK的主要思路利用前后相鄰碼元的載波相對(duì)相位變化進(jìn)行調(diào)制，以將低頻的信號(hào)調(diào)制到較高頻率，在接收端再將調(diào)制信號(hào)乘以載波做相干解調(diào)，便得到倍頻信號(hào)和原信號(hào)的疊加，再低通濾波器將倍頻信號(hào)濾掉就能得到原始信號(hào)。I、調(diào)制框圖如下：調(diào)制仿真代碼如下：%差分編碼function Y = difference(X)a=length(X);for i=1:a if X(i)=-1 X(i)=0; endendY=zeros(1,a);b=0.01:0.01:a;for(i=1:a) is(i-1)*100+1:i*100)=X(i);endif (X(1)=1

12、) Y(1)=1;else Y(1)=0;endfor i=2:a Y(i)=rem( Y(i-1)+X(i),2);endY;t=0.01:0.01:a;for(i=1:a) st(i-1)*100+1:i*100)=Y(i);End%DPSK調(diào)制function z=dpsk(X)fs=2000;w=2;%調(diào)制頻率n=length(X);t=1:100;nn=1:100;n1=0.01:0.01:n;for(i=1:n) is(i-1)*100+1:i*100)=X(i); st(i-1)*100+1:i*100)=0;endfor(nn=1:100) oss(nn)=sin(2*pi*w

13、*(t(nn)/100); osc(nn)=sin(2*pi*w*(t(nn)/100)+pi);endfor(i=1:n) if(X(i)=1) for(nn=1:100) st(100*(i-1)+nn)=oss(nn); end else for(nn=1:100) st(100*(i-1)+nn)=osc(nn); end endendz=st;II、解調(diào)框圖如下：解調(diào)仿真代碼如下： function y=dpskdemod(X) fs=200; w=2; n=length(X)/100; tt=0.01:0.01:n; dt=X.*sin(2*pi*w*tt); %低通濾波器 N,W

14、n=buttord(2*pi*25,2*pi*50,3,25,'s');%臨界頻率采用角頻率表示 b,a=butter(N,Wn,'s'); bz,az=impinvar(b,a,fs);%映射為數(shù)字的 dt=filter(bz,az,dt); %抽樣判決逆碼變換部分 dt1=0; dt2=0; for(i=1:n) for(j=1:100) dt1=dt1+dt(100*(i-1)+j); end dt2=dt1/100; if(dt2>0.05) dtt(100*(i-1)+1:100*i)=1; dt2=0; dt1=0; else dtt(100*

15、(i-1)+1:100*i)=0; dt2=0; dt1=0; end end for i=1:n y(i)=dtt(i-1)*100+50); end四、仿真結(jié)果及結(jié)論第一組用戶數(shù)據(jù)：第二組用戶數(shù)據(jù)：第三組用戶數(shù)據(jù)：信道中三徑數(shù)據(jù)的疊加：仿真說明及結(jié)論：為了觀察到較好的圖示效果，仿真時(shí)只產(chǎn)生了10個(gè)數(shù)據(jù)，并對(duì)一些波形圖進(jìn)行了部分區(qū)間的顯示。由上圖可知接收端數(shù)據(jù)與發(fā)送端數(shù)據(jù)基本一致，但仍存在一些誤碼如：用戶1接收端的數(shù)據(jù)在第一個(gè)比特和第二個(gè)比特之間有一個(gè)負(fù)脈沖，這是誤碼的存在引起的。由于時(shí)間原因，在系統(tǒng)中沒有對(duì)誤碼進(jìn)行統(tǒng)計(jì)，總體而言基本上實(shí)現(xiàn)了要求。五、總結(jié)與體會(huì)本次實(shí)驗(yàn)是集擴(kuò)頻通信、DPSK

16、調(diào)制與解調(diào)、多徑衰落于一體的綜合性實(shí)驗(yàn)。在實(shí)驗(yàn)過程中遇到很多困難，如濾波延時(shí)引起經(jīng)過濾波器之后的信號(hào)在濾波系數(shù)長(zhǎng)度的2倍加1段的幅值接近為0，多徑衰落引起信號(hào)的延時(shí)，DPSK在解調(diào)之后還需要與m序列相乘才能恢復(fù)出原始數(shù)據(jù)等等，后來經(jīng)過認(rèn)真思考和上網(wǎng)查閱相關(guān)資料等途徑終于解決了這些問題。通過本次實(shí)驗(yàn)設(shè)計(jì)，讓我對(duì)以前學(xué)習(xí)過的MATLAB軟件編程系統(tǒng)地復(fù)習(xí)了一遍，重新溫習(xí)了對(duì)于用rand來產(chǎn)生隨機(jī)序列的編寫方法及過程，還有其中的如DPSK調(diào)制函數(shù)、高斯白噪聲函數(shù)、濾波器的設(shè)計(jì)等等的用法。在系統(tǒng)的鞏固了MATLAB知識(shí)的基礎(chǔ)上，也讓我學(xué)習(xí)了擴(kuò)頻碼的產(chǎn)生原理及擴(kuò)頻碼的要求(自相關(guān)性高、互相關(guān)性低)、多徑

17、延遲信道的特性。這次實(shí)驗(yàn)有以下幾點(diǎn)收獲： （1）實(shí)驗(yàn)的核心部分是接收機(jī)部分，接收部分主要是對(duì)已調(diào)的信號(hào)進(jìn)行解擴(kuò)和解調(diào)，解擴(kuò)要注意的就是需要一個(gè)與接收端完全一樣的同步PN序列，把已調(diào)信號(hào)解擴(kuò)出來，解擴(kuò)信號(hào)再經(jīng)過DPSK解調(diào)，得到的信號(hào)波形應(yīng)與發(fā)端信號(hào)相同，說明完全恢復(fù)了所發(fā)送數(shù)據(jù)波形。 （2）擴(kuò)頻是用來傳輸信號(hào)的帶寬遠(yuǎn)遠(yuǎn)大于信息本身帶寬的一種傳輸方式，頻帶的擴(kuò)展由獨(dú)立于信息的擴(kuò)頻碼(PN碼)實(shí)現(xiàn)，與所傳數(shù)據(jù)無關(guān)，在接收端用同步接收實(shí)現(xiàn)解擴(kuò)和數(shù)據(jù)恢復(fù)。擴(kuò)頻通信的理論基礎(chǔ)為香農(nóng)定理：。由香農(nóng)定理可得出重要結(jié)論：對(duì)于給定的信息傳輸速率，可以用不同的帶寬和信噪比的組合來傳輸。擴(kuò)頻通信

18、系統(tǒng)就是利用香農(nóng)定理將信道帶寬擴(kuò)展許多倍以換取低信噪比，增強(qiáng)了系統(tǒng)的抗干擾能力。直接序列擴(kuò)頻(DS)，是直接利用具有高碼率的擴(kuò)頻碼序列在發(fā)送端擴(kuò)展信號(hào)的頻譜，而在接收端，用相同的擴(kuò)頻碼序列進(jìn)行解擴(kuò)，把展寬的擴(kuò)頻信號(hào)還原成原始的信息，是一種數(shù)字調(diào)制方法。（3）本次實(shí)驗(yàn)基本完成了實(shí)驗(yàn)的預(yù)期效果，達(dá)到了實(shí)驗(yàn)的要求，使我對(duì)M序列的生成、直接序列擴(kuò)頻、信道的多徑多用戶數(shù)據(jù)疊加、解擴(kuò)過程有了更加升入的理解和更熟練的應(yīng)用，對(duì)以后的學(xué)習(xí)和使用有很大的幫助。（4）本次實(shí)驗(yàn)讓我明白對(duì)于一個(gè)綜合性實(shí)驗(yàn)而言，首先要學(xué)會(huì)查找資料，把實(shí)驗(yàn)涉及到的基本原理搞懂，然后要把要寫程序的基本思路弄清，把整個(gè)程序分成幾個(gè)小模塊，再逐

19、個(gè)模塊去解決，最后才是寫程序和調(diào)試程序。（5）本次實(shí)驗(yàn)，讓我對(duì)無線通信原理理論和MATLAB軟件有了更深的了解，對(duì)其中很多函數(shù)的使用有了更深一步的掌握，加強(qiáng)了自己動(dòng)手能力和學(xué)習(xí)技能的提升。六、主要仿真代碼主要程序代碼如下：function PN=mseq(s)%m序列子函數(shù)function expandsignal=expand(signal,PN)%擴(kuò)頻子函數(shù)function Y = difference(X)%差分編碼子函數(shù)function z=dpsk(X)%DPSK調(diào)制子函數(shù)function multiS=channels(modusignal,k)%無線信道子函數(shù)function d

20、eexpandsignal=deexpand(signal,k)%解擴(kuò)子函數(shù)function y=dpskdemod(X)%解調(diào)子函數(shù)function y=dpskdemod(X)%解差分子函數(shù)function y=f(signal,PN)%原始信號(hào)恢復(fù)子函數(shù)%主函數(shù)user1=round(rand(1,10);%生成用戶1的0,1隨機(jī)數(shù)據(jù)user2=round(rand(1,10);%生成用戶2的0,1隨機(jī)數(shù)據(jù)user3=round(rand(1,10);%生成用戶3的0,1隨機(jī)數(shù)據(jù)x=1/100:1/100:10;for i=1:length(x)yuser1(i)=user1(ceil(

21、x(i);yuser2(i)=user2(ceil(x(i);yuser3(i)=user3(ceil(x(i);endt1=1,0,0,1,1,0;%用戶1的m序列的初始狀態(tài)t2=1,1,0,1,1,0;%用戶2的m序列的初始狀態(tài)t3=1,0,0,1,1,1;%用戶3的m序列的初始狀態(tài)PN1=mseq(t1);%生成用戶1的m序列PN2=mseq(t2);%生成用戶2的m序列PN3=mseq(t3);%生成用戶3的m序列expandsignal1=expand(user1,PN1);%生成用戶1的擴(kuò)頻序列expandsignal2=expand(user2,PN2);%生成用戶2的擴(kuò)頻序列e

22、xpandsignal3=expand(user3,PN3);%生成用戶3的擴(kuò)頻序列yy=difference(expandsignal1);y1=dpsk(difference(expandsignal1);%用戶1擴(kuò)頻調(diào)制后的信號(hào)y2=dpsk(difference(expandsignal2);%用戶2擴(kuò)頻調(diào)制后的信號(hào)y3=dpsk(difference(expandsignal3);%用戶3擴(kuò)頻調(diào)制后的信號(hào)ys1=channels(y1,2);%用戶1經(jīng)過多徑衰落后的信號(hào)ys2=channels(y2,3);%用戶2經(jīng)過多徑衰落后的信ys3=channels(y3,4);%用戶3經(jīng)過多

23、徑衰落后的信號(hào)yn=ys1+ys2+ys3;yt1=deexpand(ys1,2);%用戶1解擴(kuò)后的信號(hào)yt2=deexpand(ys2,3);%用戶2解擴(kuò)后的信號(hào)yt3=deexpand(ys3,4);%用戶3解擴(kuò)后的信號(hào)yd1=dpskdemod(yt1);%用戶1解調(diào)后的信號(hào)yd2=dpskdemod(yt2);%用戶2解調(diào)后的信號(hào)yd3=dpskdemod(yt3);%用戶3解調(diào)后的信號(hào)yi1=indifference(yd1);%用戶1解差分yi2=indifference(yd2);%用戶2解差分yi3=indifference(yd3);%用戶3解差分yr1=f(yi1,PN1)

24、;%用戶1數(shù)據(jù)恢復(fù)yr2=f(yi2,PN2);%用戶2數(shù)據(jù)恢復(fù)yr3=f(yi3,PN3);%用戶3數(shù)據(jù)恢復(fù)%用戶1部分%figuresubplot(3,2,1)plot(x,yuser1);title('用戶1的發(fā)送數(shù)據(jù)');axis(0 10 -0.2 1.2);gridsubplot(3,2,2)plot(1/63:1/63:10,yr1);axis(0,10,-0.1,1.1);title('接收到用戶1數(shù)據(jù)');gridsubplot(3,2,3)plot(expandsignal1);title('用戶1擴(kuò)頻后的時(shí)域信號(hào)');axi

25、s(0 630 -1.2 1.2);subplot(3,2,4)plot(0.01:0.01:length(yy),yt1);title('用戶1解擴(kuò)后的時(shí)域信號(hào)');axis(0 20 -0.3 0.3);subplot(3,2,5)plot(0.01:0.01:length(yy),y1);title('用戶1DPSK調(diào)制后的時(shí)域信號(hào)');axis(0 20 -1.2 1.2);subplot(3,2,6)plot(yi1);title('用戶1DPSK解調(diào)后的時(shí)域信號(hào)');axis(0 630 0.2 1.2);figuresubplot(

26、2,2,1)plot(0:length(expandsignal1)-1)*100/length(expandsignal1)-100/2,abs(fftshift(fft(expandsignal1);title('用戶1擴(kuò)頻后的頻譜圖');subplot(2,2,2)plot(0:length(yt1)-1)*100/length(yt1)-100/2,abs(fftshift(fft(yt1);title('用戶1解擴(kuò)后的頻譜圖');axis(-15 15 0 400);subplot(2,2,3)plot(0:length(y1)-1)*100/leng

27、th(y1)-100/2,abs(fftshift(fft(y1);title('用戶1調(diào)制后的頻譜圖');axis(-20 20 0 3000);subplot(2,2,4)plot(0:length(yi1)-1)*100/length(yi1)-100/2,abs(fftshift(fft(yi1);title('用戶1解調(diào)后的頻譜圖');axis(-20 20 0 400);%用戶2部分%figuresubplot(3,2,1)plot(x,yuser2);title('用戶2的發(fā)送數(shù)據(jù)');axis(0 10 -0.2 1.2);gri

28、dsubplot(3,2,2)plot(1/63:1/63:10,yr2);axis(0,10,-0.1,1.1);title('接收到用戶2數(shù)據(jù)');gridsubplot(3,2,3)plot(expandsignal2);title('用戶2擴(kuò)頻后的時(shí)域信號(hào)');axis(0 630 -1.2 1.2);subplot(3,2,4)plot(0.01:0.01:length(yy),yt2);title('用戶2解擴(kuò)后的時(shí)域信號(hào)');axis(0 20 -0.8 0.8);subplot(3,2,5)plot(0.01:0.01:lengt

29、h(yy),y2);title('用戶2DPSK調(diào)制后的時(shí)域信號(hào)');axis(0 20 -1.2 1.2);subplot(3,2,6)plot(yi2);title('用戶2DPSK解調(diào)后的時(shí)域信號(hào)');axis(0 630 0.2 1.2);figuresubplot(2,2,1)plot(0:length(expandsignal2)-1)*100/length(expandsignal2)-100/2,abs(fftshift(fft(expandsignal2);title('用戶2擴(kuò)頻后的頻譜圖');subplot(2,2,2)pl

30、ot(0:length(yt2)-1)*100/length(yt2)-100/2,abs(fftshift(fft(yt2);title('用戶2解擴(kuò)后的頻譜圖');axis(-15 15 0 1500);subplot(2,2,3)plot(0:length(y2)-1)*100/length(y2)-100/2,abs(fftshift(fft(y2);title('用戶2調(diào)制后的頻譜圖');axis(-20 20 0 3000);subplot(2,2,4)plot(0:length(yi2)-1)*100/length(yi2)-100/2,abs(f

31、ftshift(fft(yi2);title('用戶2解調(diào)后的頻譜圖');axis(-20 20 0 400);%用戶3部分%figuresubplot(3,2,1)plot(x,yuser3);title('用戶3的發(fā)送數(shù)據(jù)');axis(0 10 -0.2 1.2);gridsubplot(3,2,2)plot(1/63:1/63:10,yr3);axis(0,10,-0.1,1.1);title('接收到用戶3數(shù)據(jù)');gridsubplot(3,2,3)plot(expandsignal3);title('用戶3擴(kuò)頻后的時(shí)域信號(hào)&

32、#39;);axis(0 630 -1.2 1.2);subplot(3,2,4)plot(0.01:0.01:length(yy),yt3);title('用戶3解擴(kuò)后的時(shí)域信號(hào)');axis(0 20 -1.6 1.6);subplot(3,2,5)plot(0.01:0.01:length(yy),y3);title('用戶3DPSK調(diào)制后的時(shí)域信號(hào)');axis(0 20 -1.2 1.2);subplot(3,2,6)plot(yi3);title('用戶3DPSK解調(diào)后的時(shí)域信號(hào)');axis(0 630 0.2 1.2);figuresubplot(2,2,1)plot(0:length(expandsignal3)-1)*100/length(expandsignal3)-100/2,abs(fftshift(fft(expandsignal3);title('用戶3擴(kuò)頻后的頻譜圖');subplot(2,2,2)plot(0:length(yt3)-1)*100/length(yt3)-100/2,abs(fftshift(fft(yt3);title('用戶3解擴(kuò)后的頻譜圖');axis(-15 15 0 5000);subplo