單極性歸零(RZ)波形及其Matlab仿真

1、精選優(yōu)質(zhì)文檔-傾情為你奉上單極性歸零（RZ）波形及其仿真1. 仿真原理及思路功率信號的平均功率：。該信號的雙邊功率譜密度為：，其中是截短后的傅氏變換，是的能量譜，是在截短時間內(nèi)的功率譜。對于仿真系統(tǒng)，若x是時域取樣值矢量，X是對應(yīng)的傅氏變換，那么x的功率譜便為矢量。對于采用歸零（RZ）及不歸零（NRZ）矩形脈沖波形的數(shù)字信號，可以用簡單的方法信號矢量s。設(shè)a是碼元矢量，N是總?cè)狱c數(shù)，M是總碼元數(shù)，L是每個碼元內(nèi)的點數(shù)，那么NRZ信號可這樣獲得：s=zeros(1,N);for ii=1:L, s(ii+0:M-1*L)=a; end對于，若Rt是要求的占空比，dt是仿真系統(tǒng)的時域采樣間隔，則

2、RZ信號的產(chǎn)生方法是：s=zeros(1,N);for ii=1:Rt/dt, s(ii+0:M-1*L)=a; end首先，利用rand函數(shù)產(chǎn)生一個尺寸為1*M的矩陣，其元素按均勻分布隨機取值于區(qū)間0,1，并用round函數(shù)對其四舍五入，得到一個隨機產(chǎn)生的0,1序列；其次，利用for循環(huán)產(chǎn)生碼元長度為L點，碼元為0,1，且占空比為50%的單極性歸零碼，并畫出其波形圖；最后，計算該RZ的功率譜密度，并畫出其波形。2. 程序流程圖產(chǎn)生0,1分布的隨機序列產(chǎn)生碼元為0,1的RZ計算功率譜，并畫出波形畫出波形RZ波形及其功率譜密度仿真流程圖3. 仿真程序及運行結(jié)果仿真程序：%實驗二：畫出單極性歸零碼

3、及其功率譜close allclear allglobal dt t df NN=214;%采樣點數(shù)L=64;%每碼元的采樣點數(shù)M=N/L;%碼元數(shù)Rb=2;%碼速率為2Mb/sTs=1/Rb;%碼元間隔dt=Ts/L;%時域采樣間隔Rt=0.5; %占空比df=1/(N*dt);%頻域采樣間隔T=N*dt; %截短時間Bs=N*df/2;%系統(tǒng)帶寬t=linspace(-T/2,T/2,N);%時域橫坐標(biāo)f=linspace(-Bs,Bs,N);%頻域橫坐標(biāo)EP=zeros(1,N);for jj=1:100 a=round(rand(1,M); %產(chǎn)生M個取值0，1等概的隨機碼 s=zer

4、os(1,N); %產(chǎn)生一個N個元素的零序列 for ii=1:Rt*Ts/dt s(ii+0:M-1*L)=a;%產(chǎn)生單極性歸零碼 end Q=t2f(s); %付氏變換 P=Q.*conj(Q)/T; %P為單極性歸零碼的功率 EP=(EP*(jj-1)+P)/jj; %累計平均end aa=30+10*log10(EP+eps);%加eps以避免除以零figure(1)set(1,'Position',10,350,600,200)%設(shè)定窗口位置及大小 figure(2)set(2,'Position',10,50,600,200)%設(shè)定窗口位置及大小fi

5、gure(1)plot(f,aa,'g')xlabel('f(MHZ)')ylabel('Ps(f)')axis(-15,+15, 1.1*min(aa), 1.1*max(aa)grid onfigure(2)plot(t,s,'r')xlabel('t(ms)')ylabel('s(t)(V)')axis(-10,+10,1.1*min(s),1.1*max(s)grid on4. 程序結(jié)果：輸入取樣點數(shù)=2k,k=14，得到如下波形：5. 實驗結(jié)果分析由上圖結(jié)果可以看單極性歸零碼序列的功率譜密度不僅含有離散的直流分量及連續(xù)譜,而且還包含離散的時鐘分量和奇次諧波分量，功率譜主瓣寬度為碼元速率2Rb Mb/s。根據(jù)數(shù)字PAM信號功率譜密度公式得到：功率譜分為兩個部分，第一部分是連續(xù)譜，形狀取決于G