北郵通信原理軟件實(shí)驗(yàn)報(bào)告2013

1、通信原理軟件實(shí)驗(yàn)報(bào)告學(xué)院：信息與通信工程學(xué)院班級(jí)：一、 通信原理Matlab仿真實(shí)驗(yàn)實(shí)驗(yàn)八一、實(shí)驗(yàn)內(nèi)容假設(shè)基帶信號(hào)為m(t)=sin(2000*pi*t)+2cos(1000*pi*t),載波頻率為20kHz，請(qǐng)仿真出AM、DSB-SC、SSB信號(hào)，觀察已調(diào)信號(hào)的波形和頻譜。二、 實(shí)驗(yàn)原理1、具有離散大載波的雙邊帶幅度調(diào)制信號(hào)AM該幅度調(diào)制是由DSB-SC AM信號(hào)加上離散的大載波分量得到，其表達(dá)式及時(shí)間波形圖為：應(yīng)當(dāng)注意的是，m(t)的絕對(duì)值必須小于等于1，否則會(huì)出現(xiàn)下圖的過(guò)調(diào)制：AM信號(hào)的頻譜特性如下圖所示： 由圖可以發(fā)現(xiàn)，AM信號(hào)的頻譜是雙邊帶抑制載波調(diào)幅信號(hào)的頻譜加上離散的大載波分量。

2、2、雙邊帶抑制載波調(diào)幅（DSBSC AM）信號(hào)的產(chǎn)生 雙邊帶抑制載波調(diào)幅信號(hào)s(t)是利用均值為0的模擬基帶信號(hào)m(t)和正弦載波 c(t)相乘得到，如圖所示： m(t)和正弦載波s(t)的信號(hào)波形如圖所示：若調(diào)制信號(hào)m(t)是確定的，其相應(yīng)的傅立葉頻譜為M(f)，載波信號(hào)c(t)的傅立葉頻譜是C(f),調(diào)制信號(hào)s(t)的傅立葉頻譜S(f)由M(f)和C(f)相卷積得到，因此經(jīng)過(guò)調(diào)制之后，基帶信號(hào)的頻譜被搬移到了載頻fc處，若模擬基帶信號(hào)帶寬為W，則調(diào)制信號(hào)帶寬為2W，并且頻譜中不含有離散的載頻分量，只是由于模擬基帶信號(hào)的頻譜成分中不含離散的直流分量。3、 單邊帶條幅SSB信號(hào) 雙邊帶抑制載波

3、調(diào)幅信號(hào)要求信道帶寬B=2W,其中W是模擬基帶信號(hào)帶寬。從信息論關(guān)點(diǎn)開(kāi)看，此雙邊帶是有剩余度的，因而只要利用雙邊帶中的任一邊帶來(lái)傳輸，仍能在接收機(jī)解調(diào)出原基帶信號(hào)，這樣可減少傳送已調(diào)信號(hào)的信道帶寬。 單邊帶條幅SSB AM信號(hào)的其表達(dá)式：或其頻譜圖為：三、 仿真設(shè)計(jì)1、流程圖：作圖傅立葉變換AM信號(hào)SSB信號(hào)加上已調(diào)DSB-SC信號(hào)與載頻的希爾伯特變換相乘DSB-SC信號(hào)加上載波信號(hào)希爾伯特變換與載波的余弦信號(hào)相乘與載波的余弦信號(hào)相乘基帶信號(hào)m(t)2、實(shí)驗(yàn)結(jié)果&分析討論實(shí)驗(yàn)仿真結(jié)果從上至下依次是AM信號(hào)、DSB信號(hào)、SSB信號(hào)。從仿真結(jié)果看，AM調(diào)制信號(hào)包絡(luò)清晰，可利用包絡(luò)檢波恢復(fù)原

4、信號(hào)，接收設(shè)備較為簡(jiǎn)單。其頻譜含有離散大載波，從理論分析可知，此載波占用了較多發(fā)送功率，使得發(fā)送設(shè)備功耗較大。DSB-SC信號(hào)波形和頻譜，其時(shí)域波形有相位翻轉(zhuǎn)，頻譜不含離散大載波。SSB信號(hào)比DSB信號(hào)節(jié)省一半帶寬，適合于語(yǔ)聲信號(hào)的調(diào)制，因?yàn)槠錄](méi)有直流分量，也沒(méi)有很低頻的成分。3、結(jié)果分析：根據(jù)通原理論課的知識(shí)可知，信號(hào)的AM調(diào)制比較容易實(shí)現(xiàn)，但其功率譜中有相當(dāng)大一部分是載頻信號(hào)，效率非常低；DSB-SC調(diào)制解決了AM信號(hào)效率低下的問(wèn)題，但仍然存在的問(wèn)題是調(diào)制信號(hào)的帶寬為基帶信號(hào)的兩倍，頻譜利用率較低；SSB調(diào)制方式在頻譜利用上又做出了改進(jìn)，為原先的一半，但其可靠性降低了，總之，可靠性與有效性

5、是難以?xún)扇涿赖?，為一?duì)矛盾體。四、 程序代碼fs = 800;%kHzT = 200;N = T*fs;t = linspace(-T/2,T/2,N);f = linspace(-fs/2,fs/2,N);fm1 = 1;%kHzfm2 = 0.5;%kHzfc = 20;%kHzmt = sin(2*pi*fm1*t)+2*cos(2*pi*fm2*t);%AM調(diào)制信號(hào)波形和頻譜A = 2;a = 0.3;AM = (1 + a*mt).*cos(2*pi*fc*t);figure(1);subplot(2,1,1);plot(t,AM);xlabel('時(shí)間t(ms)'

6、);ylabel('AM');axis(-3,3,-3,3);title('AM波形');grid on;FAM = t2f(AM,fs);subplot(2,1,2);plot(f,abs(FAM);xlabel('頻率f(kHz)');ylabel('FAM');axis(-25,25,0,150);title('AM幅頻特性');grid on;%DSB-SC信號(hào)波形和頻譜DSB = mt.*cos(2*pi*fc*t);figure(2);subplot(2,1,1);plot(t,DSB);xlabel(

7、'時(shí)間t(ms)');ylabel('DSB');axis(-3,3,-3,3);title('DSB波形');grid on;FDSB = t2f(DSB,fs);subplot(2,1,2);plot(f,abs(FDSB);xlabel('頻率f(kHz)');ylabel('FDSB');axis(-25,25,0,150);title('DSB幅頻特性');grid on;%SSB信號(hào)波形和頻譜M = t2f(mt,fs);MH = -j*sign(f).*M;mh = real(f2t(

8、MH,fs);SSB = mt.*cos(2*pi*fc*t)- mh.*sin(2*pi*fc*t);FSSB = t2f(SSB,fs);figure(3);subplot(2,1,1);plot(t,SSB);xlabel('時(shí)間t(ms)');ylabel('SSB');axis(-3,3,-5,5);title('SSB波形');grid on;subplot(2,1,2);plot(f,abs(FSSB);xlabel('頻率f(kHz)');ylabel('FSSB');axis(-25,25,0,1

9、50);title('SSB幅頻特性');grid on;實(shí)驗(yàn)九一、 實(shí)驗(yàn)內(nèi)容假設(shè)基帶信號(hào)為m(t)=sin(2000t)+2cos(1000t)+4sin(500t+/3),載波頻率為40kHz，仿真產(chǎn)生FM信號(hào)，觀察波形與頻譜。FM的頻偏常數(shù)為5kHz/V。二、 實(shí)驗(yàn)原理1、調(diào)頻信號(hào)表達(dá)式為：2、 調(diào)頻信號(hào)的頻譜特征： 單頻調(diào)制信號(hào)包括無(wú)窮多頻率分量，在實(shí)際應(yīng)用中，我們運(yùn)用等效帶寬這個(gè)概念。等效帶寬的定義為：包含98%或者99%的已調(diào)信號(hào)總功率的帶寬，根據(jù)卡松公式有：三、 仿真設(shè)計(jì)1、流程圖作圖（時(shí)域圖&頻域圖）傅立葉變換將累計(jì)求和的結(jié)果作為瞬時(shí)相位偏移，得到已調(diào)信

10、號(hào)s(t)的表達(dá)式，即FM信號(hào)利用scilab中的cumsum函數(shù)對(duì)基帶信號(hào)m(t)進(jìn)行求累積和，相當(dāng)于對(duì)其進(jìn)行積分基帶信號(hào)m(t)2、實(shí)驗(yàn)結(jié)果&分析討論3、結(jié)果分析由于max(m(t)=6.4833，信號(hào)帶寬fm=1KHz，Kf=5KHz/V，由卡松公式 ,計(jì)算得帶寬為66.8320KHz，與頻域仿真結(jié)果基本相符。（事實(shí)上，由于調(diào)制信號(hào)不是單頻信號(hào)，卡松公式不滿(mǎn)足）四、 程序代碼fs = 800;%kHzT = 200;N = T*fs;dt = 1/fs;t = linspace(-T/2,T/2,N);f = linspace(-fs/2,fs/2,N);fm1 = 1;%kHz

11、fm2 = 0.5;%kHzfm3 = 0.25;%kHzfc = 40;%kHzk = 5;%kHzmt = sin(2*pi*fm1*t)+2*cos(2*pi*fm2*t)+4*sin(2*pi*fm3*t+pi/3);phi = 2*pi*k*cumsum(mt)*dt;FM = cos(2*pi*fc*t+phi);subplot(2,1,1);plot(t,FM);xlabel('時(shí)間t(ms)');ylabel('FM');axis(-2,2,-3,3);title('FM波形');grid on;FFM = t2f(FM,fs);

12、subplot(2,1,2);plot(f,abs(FFM);xlabel('頻率f(kHz)');ylabel('FFM');axis(-80,80,0,20);title('FM幅頻特性');grid on;實(shí)驗(yàn)十一一、 實(shí)驗(yàn)內(nèi)容通過(guò)仿真測(cè)量占空比為50%、75%以及100%的單、雙極性歸零碼波形以及其功率譜，分析不同占空比對(duì)仿真結(jié)果的影響。二、 實(shí)驗(yàn)原理1、時(shí)域特性1.1二進(jìn)制單極性歸零碼：1.2二進(jìn)制雙極性歸零碼：2、頻域特性（功率譜）：數(shù)字基帶信號(hào)s(t)的功率譜密度為：數(shù)字基帶信號(hào)s(t)的功率譜密度與隨機(jī)序列的功率譜特性以及發(fā)送濾波

13、器的頻率特性有關(guān)。在實(shí)隨機(jī)序列的各符號(hào)互不相關(guān)時(shí)，s(t)的功率譜為2.1單極性歸零碼：2.4雙極性歸零碼：三、 仿真設(shè)計(jì)1、流程圖：作圖（波形&功率譜）利用公式P=abs(S1).2/T，求樣本信號(hào)的功率譜密度，然后對(duì)各個(gè)樣本的功率譜密度求數(shù)學(xué)期望，得到功率譜密度設(shè)置占空比（分別為25%、50%、75%、100%）利用randn(1,M)以及邏輯判斷，產(chǎn)生單（雙）極性數(shù)據(jù)2、實(shí)驗(yàn)結(jié)果&分析討論2.1單極性從上至下依次是占空比為50%、75%、100%。從仿真結(jié)果可以看出，單極性歸零碼的頻譜主瓣寬度隨占空比增加而減小，且含有沖激2.2雙極性從上至下依次是占空比50%、75%、1

14、00%。從仿真結(jié)果可以看出，隨占空比增加，頻譜主瓣寬度減小，且不含沖激。3、結(jié)果分析：3.1單極性 由于單極性碼含有直流分量，所以表現(xiàn)在頻域內(nèi)為在直流處奇級(jí)次諧波處有一個(gè)沖激，而其功率譜主瓣寬度隨著占空比的不同而不同，對(duì)于50%、75%、100%的占空比的帶寬分別對(duì)應(yīng)為碼元速率B的4倍、2倍、1.33倍和1倍。3.2雙極性 由于雙極性碼不含直流分量，所以沒(méi)有單極性碼所具有的直流處和奇次諧波處的沖激，但在帶寬上與單極性碼表現(xiàn)一致。3.3對(duì)比 因?yàn)楹芏鄠鬏斁€(xiàn)路不能傳送直流分量，而且在誤碼特性上雙極性馬要比單極性碼好很多，故此，一般都優(yōu)先選用雙極性碼。四、 程序代碼%通過(guò)仿真測(cè)量占空比為50%、75

15、%以及100%的單雙極性歸零碼波形及其功率譜，分析不同占空比對(duì)仿真結(jié)果的影響clear all;xdel(winsid(); /關(guān)閉所有圖形窗口exec t2f.sci; exec f2t.sci;L=32; /每個(gè)碼元間隔內(nèi)的采樣點(diǎn)數(shù)N=213; /總采樣點(diǎn)數(shù)M=N/L; /總碼元數(shù)Rb=2; /碼元速率Ts=1/Rb; /比特間隔 fs=L/Ts; /采樣速率T=N/fs; /截短時(shí)間 Bs=fs/2; /系統(tǒng)帶寬 t=-T/2+0:N-1/fs; /時(shí)域采樣點(diǎn)f=-Bs+0:N-1/T; /頻域采樣點(diǎn)L0=input('請(qǐng)輸入占空比(01)：')EP=zeros(1,N)

16、;ch=input('請(qǐng)選擇要觀察的碼型：1-單極性；其他-雙極性')for loop=1:1000/1000次樣本函數(shù)取平均 if ch=1 a=(rand(1,M)>0.5)+0; /生成單極性序列 else a=sign(rand(1,M)>0.5)-0.5); /生成雙極性序列 end tmp=zeros(L,M);/一個(gè)碼元的歸零部分取零 L1=L*L0;/占空比，求出一個(gè)碼元不歸零部分的取樣點(diǎn)數(shù) tmp(1:L1,:)=ones(L1,1)*a;/將一個(gè)碼元不歸零部分的取樣點(diǎn)值置為1 s=tmp(:)' S=t2f(s,fs); /傅里葉變化 P

17、=abs(S).2/T; /樣本的功率譜密度 EP=EP*(1-1/loop)+P/loop; /隨機(jī)過(guò)程的功率譜是各個(gè)樣本的功率譜的數(shù)學(xué)期望end xset("window",1)plot(t,s)set(gca(),"grid",1,1) title('時(shí)域圖')xlabel('t')ylabel('S(t)')mtlb_axis(-3,3,-1.5,1.5);xset("window",2)plot(f,abs(EP+%eps)set(gca(),"grid",

18、1,1) title('功率譜圖形')xlabel('f')ylabel('功率') mtlb_axis(-35,35,-5,max(EP+%eps); xset("window",3) plot(f,10*log10(EP+%eps)set(gca(),"grid",1,1) title('功率譜圖形(dB)')xlabel('f')ylabel('功率') 實(shí)驗(yàn)十二一、實(shí)驗(yàn)內(nèi)容仿真測(cè)量滾降系數(shù)為a=0.25的根升余弦滾降系統(tǒng)的發(fā)送功率譜密度及眼圖。二、實(shí)驗(yàn)

19、原理1. 理想限帶情況下的最佳基帶傳輸：a的存在，使得傳輸頻帶增加，降低了頻譜利用率。A愈小，波形振蕩起伏愈大，傳輸頻帶擴(kuò)展??；反之，a愈大，波形振蕩起伏愈小，頻帶擴(kuò)展增大。實(shí)際中一般選取a不小于0.2。2、眼圖實(shí)際通信系統(tǒng)中，數(shù)字信號(hào)經(jīng)過(guò)非理想的傳輸系統(tǒng)產(chǎn)生畸變，總是在不同程度上存在碼間干擾的，系統(tǒng)性能很難進(jìn)行定量的分析，常常甚至得不到近似結(jié)果。而眼圖可以直觀地估價(jià)系統(tǒng)碼間干擾和噪聲的影響，是常用的測(cè)試手段。眼圖分析中常用結(jié)論：(1)最佳取樣時(shí)刻應(yīng)選擇在眼睛張開(kāi)最大的時(shí)刻；(2)眼睛閉合的速率，即眼圖斜邊的斜率，表示系統(tǒng)對(duì)定時(shí)誤差靈敏的程度，斜邊愈陡，對(duì)定位誤差愈敏感；(3)在取樣時(shí)刻上，陰

20、影區(qū)的垂直寬度表示最大信號(hào)失真量；(4)在取樣時(shí)刻上，上下兩陰影區(qū)的間隔垂直距離之半是最小噪聲容限，噪聲瞬時(shí)值超過(guò)它就有可能發(fā)生錯(cuò)誤判決；(5)陰影區(qū)與橫軸相交的區(qū)間表示零點(diǎn)位置變動(dòng)范圍，它對(duì)于從信號(hào)平均零點(diǎn)位置提取定時(shí)信息的解調(diào)器有重要影響。三、仿真設(shè)計(jì)1、流程圖產(chǎn)生M個(gè)0、2等概隨機(jī)碼定義升余弦滾降系統(tǒng)系統(tǒng)函數(shù)序列經(jīng)過(guò)升余弦滾降濾波器作出功率譜和眼圖傅里葉反變換得時(shí)域波形2、實(shí)驗(yàn)結(jié)果&分析討論發(fā)送功率譜密度：實(shí)驗(yàn)結(jié)果：仿真結(jié)果與計(jì)算公式=2(1+0.25)/2=1.25基本相等。眼圖：五、 程序代碼clear allglobal dt t df N close all N=214;

21、 %采樣點(diǎn)數(shù) L=32; %每碼元的采樣點(diǎn)數(shù) M=N/L %碼元數(shù) Rb=2; %碼速率是2Mb/s Ts=1/Rb; %碼元間隔 dt=Ts/L; %時(shí)域采樣間隔 df=1/(N*dt) %頻域采樣間隔 T=N*dt %時(shí)域截短時(shí)間 Bs=N*df/2 %系統(tǒng)帶寬 Na=4; %設(shè)示波器掃描寬度為4個(gè)碼元alpha= 0.25; t=-T/2+dt/2:dt:T/2; %時(shí)域橫坐標(biāo) f=-Bs+df/2:df:Bs; %頻域橫坐標(biāo) g1=sin(pi*t/Ts)./(pi*t/Ts); g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).2); g=g1.*g2

22、; %升余弦脈沖波形 G=t2f(g,df); %升余弦脈沖波形的傅氏正變換 figure(1) %set(1,'Position',100,50,300,200) %設(shè)定窗口位置及大小 %set(2,'Position',350,50,300,200); %設(shè)定窗口位置及大小 figure(2) hold on %圖形疊加，即重復(fù)掃描 grid xlabel('t in us') ylabel('s(t) in V') EP=zeros(size(f)+eps; %加eps是為了躲零 for ii=1:2 % 設(shè)累計(jì)平均的循環(huán)

23、a=sign(randn(1,M); imp=zeros(1,N); %產(chǎn)生沖激序列 imp(L/2:L:N)=a/dt;S=t2f(imp,df).*G; %升余弦信號(hào)的傅氏變換 s=f2t(t2f(imp,df).*G,df); %升余弦信號(hào)的時(shí)域波形 s=real(s); P=S.*conj(S)/T; %升余弦信號(hào)的功率譜 EP=(EP*(ii-1)+P+eps)/ii; %累計(jì)平均 figure(1) plot(f,30+10*log10(EP),'g'); grid axis(-3,+3,-200,200) xlabel('f (MHz)') yla

24、bel('Ps(f) (dBm/MHz)')figure(2) tt=0:dt:Na*L*dt; %設(shè)定畫(huà)眼圖的坐標(biāo)的位置 for jj=1:Na*L:N-Na*L %Na*L是每次畫(huà)眼圖的間隔步長(zhǎng) ，N-Na*L是畫(huà)眼圖的 %最后一點(diǎn)的起始點(diǎn)。 plot(tt,s(jj:jj+Na*L); %畫(huà)升余弦信號(hào)的眼圖 end end % 對(duì)應(yīng)累計(jì)平均的循環(huán)實(shí)現(xiàn)數(shù)字基帶傳輸系統(tǒng)實(shí)驗(yàn)?zāi)康模豪镁幊虒?shí)現(xiàn)數(shù)字基帶傳輸系統(tǒng)1）畫(huà)出發(fā)送端輸入碼序列波形和功率譜、發(fā)送濾波器輸出波形和功率譜2）畫(huà)出接收端采樣判決后碼序列波形和功率譜、接收濾波器輸出波形和功率譜3）畫(huà)出接收濾波器輸出信號(hào)眼圖（在升余

25、弦滾將系數(shù)分別為0、0.5、1的情況下）4）分別畫(huà)出升余弦滾將系數(shù)為0、0.5、1，采樣判決點(diǎn)在眼圖最大處的系統(tǒng)的實(shí)際誤碼曲線(xiàn)（Pes/n0曲線(xiàn)），并在同坐標(biāo)系中畫(huà)出理論誤碼曲線(xiàn)5）改變采樣點(diǎn)重復(fù)1）4）。程序代碼：主程序：clear allexec t2f.sci;exec f2t.sci;exec eyes.sci;k=13;N=2k;/采樣點(diǎn)數(shù)L=16;/每碼元的采樣點(diǎn)數(shù)M=N/L;/碼元數(shù)Rb=2; /碼速率是2Mb/sTs=1/Rb; /碼元間隔dt=Ts/L; /時(shí)域采樣間隔fs=1/dt;/采樣速率df=1/(N*dt); /頻域采樣間隔T=N*dt; /截短時(shí)間Bs=N*df/

26、2; /系統(tǒng)帶寬f=-Bs+df/2:df:Bs; /頻域橫坐標(biāo)t=-T/2+dt/2:dt:T/2; /時(shí)域橫坐標(biāo)alpha=1; /升余弦滾降系數(shù)Hcos=zeros(1,N);ii=find(abs(f)>(1-alpha)/(2*Ts)&abs(f)<=(1+alpha)/(2*Ts);Hcos(ii)=Ts/2*(1+cos(%pi*Ts/alpha*(abs(f(ii)-(1-alpha)/(2*Ts);ii=find(abs(f)<=(1-alpha)/(2*Ts);Hcos(ii)=Ts;GT=sqrt(Hcos);GR=GT; /最佳系統(tǒng)的發(fā)送接收濾

27、波器的傅式變換EPC=zeros(1,N);EPB=zeros(1,N);EP1=zeros(1,N);EP2=zeros(1,N);RECT=ones(L,1);for l1=1:20; Eb_N0(l1)=(l1-1); /信噪比 eb_n0(l1)=10(Eb_N0(l1)/10); Eb=1; n0=Eb/eb_n0(l1); /信道的噪聲譜密度 sita=n0*Bs; /信道中噪聲功率 n_err=0; /誤碼計(jì)數(shù) for loop=1:100 a=sign(rand(1,M,"normal"); tmp1=RECT*a; seq_send=tmp1(:)'

28、; /發(fā)送端產(chǎn)生序列！ S1=t2f(seq_send,fs); P1=abs(S1).2/T;/樣本信號(hào)的功率譜密度 EP1=EP1*(1-1/loop)+P1/loop;/隨機(jī)過(guò)程的功率譜是各個(gè)樣本的功率譜的數(shù)學(xué)期望 EP11=20+10*log10(EP1+%eps); sp=zeros(1,N);/產(chǎn)生沖激序列 sp(1:L:N)=a/dt; SP=t2f(sp,fs); SH=SP.*GT; sh=real(f2t(SH,fs); /通過(guò)發(fā)送濾波器后的信號(hào)！ PB=abs(SH).2/T; /過(guò)升余弦后能量譜密度 EPB=EPB*(1-1/loop)+PB/loop;/求平均 EPB

29、B=20+10*log10(EPB+%eps); nw=sqrt(sita)*rand(1,N,"normal"); /信道噪聲 SY=SH.*GR; sy=real(f2t(SY,fs); r=sy+nw; /接收信號(hào) R=t2f(r,fs); PC=abs(R).2/T; EPC=EPC*(1-1/loop)+PC/loop;/收端 EPCC=20+10*log10(EPC+%eps); sd=r(L/4:L:N);/*取樣！ sdp=sign(sd);/判決 tmp2=RECT*sdp; seq_reci=tmp2(:)' /接收端信號(hào)！ S2=t2f(seq

30、_reci,fs); P2=abs(S2).2/T;/樣本信號(hào)的功率譜密度 EP2=EP2*(1-1/loop)+P2/loop; EP22=20+10*log10(EP2+%eps); n_err=n_err+length(find(sdp=a);/錯(cuò)誤累計(jì)endPe(l1)=n_err/(M*loop); xset("window",5) plot(Eb_N0,log10(Pe+%eps),'g'); /PeEb/N0曲線(xiàn)畫(huà)圖 xlabel('Eb/N0');ylabel('Pe');title("PeEb/N0

31、曲線(xiàn)"); eb_n0=10.(Eb_N0/10); set(gca(),"auto_clear","off") k = log10(0.5*erfc(sqrt(eb_n0); plot(Eb_N0,k) mtlb_axis(0,15,-3.5,0); xlabel('Eb/N0') ylabel('Pe') legend('實(shí)際的','理論的'); set(gca(),"auto_clear","on")endxset("wind

32、ow",1)subplot(2,2,1)plot(t,seq_send)title("最初產(chǎn)生序列輸出波形")mtlb_axis(0,L,-1.5,1.5)xgridsubplot(2,2,2)plot(t,sh)title("通過(guò)發(fā)送濾波器后的信號(hào)波形")xlabel("t(ms)")ylabel("s2(t)(V)")mtlb_axis(0,L,-2.5,2.5)xgridsubplot(2,2,3)plot(t,sy)title("采樣前的信號(hào)波形")xlabel("t

33、(ms)")ylabel("y(t)(V)")mtlb_axis(0,L,-2,2)xgridsubplot(2,2,4)plot(t,seq_reci)title("判決后輸出序列")mtlb_axis(0,L,-1.5,1.5)xgridxset("window",2)subplot(2,2,1)plot(f,EP11)xlabel("f/kHz")ylabel("|FMf|/(V/Hz)")mtlb_axis(-10,10,-20,max(EP11)title("發(fā)送序列頻率譜")xgridsubplot(2,2,2)plot(f,EPBB)title("發(fā)送濾波器輸出功率譜")xlabel("f(kHz)")ylabel(&