1、单极性归零（ RZ）波形及其仿真1. 仿真原理及思路功率信号 的平均功率： 。该信号的双f(t) TfP(f)dlimftd212边功率谱密度 为：，其中 是 截短后的傅氏变换，PTF()li2SfT()st是 的能量谱， 是 在截短时间内的功率谱。对于仿真SfT()2stT()SfstT()系统，若 x 是时域取样值矢量，X 是对应的傅氏变换，那么 x 的功率谱便为矢量 。P=.*conj/对于采用归零（RZ）及不归零（NRZ）矩形脉冲波形的数字信号，可以用简单的方法信号矢量 s。设 a 是码元矢量，N 是总取样点数，M 是总码元数，L是每个码元内的点数，那么 NRZ 信号可这样获得：s=z
2、eros(1,N);for ii=1:L, s(ii+0:M-1*L)=a; end对于，若 Rt 是要求的占空比，dt 是仿真系统的时域采样间隔，则 RZ 信号的产生方法是：s=zeros(1,N);for ii=1:Rt/dt, s(ii+0:M-1*L)=a; end首先，利用 rand 函数产生一个尺寸为 1*M 的矩阵，其元素按均匀分布随机取值于区间0,1 ，并用 round 函数对其四舍五入，得到一个随机产生的 0,1 序列；其次，利用 for 循环产生码元长度为 L 点，码元为 0,1，且占空比为 50%的单极性归零码，并画出其波形图；最后，计算该 RZ 的功率谱密度，并画出其波
3、形。2. 程序流程图产生 0,1 分布的随机序列 产生码元为 0,1 的RZ 计算功率谱，并画出波形画出波形RZ 波形及其功率谱密度仿真流程图3. 仿真程序及运行结果仿真程序：%实验二：画出单极性归零码及其功率谱close allclear allglobal dt t df NN=214; %采样点数L=64; %每码元的采样点数M=N/L; %码元数Rb=2; %码速率为 2Mb/sTs=1/Rb; %码元间隔dt=Ts/L; %时域采样间隔Rt=0.5; %占空比df=1/(N*dt); %频域采样间隔T=N*dt; %截短时间Bs=N*df/2; %系统带宽t=linspace(-T/
4、2,T/2,N); %时域横坐标f=linspace(-Bs,Bs,N); %频域横坐标EP=zeros(1,N);for jj=1:100a=round(rand(1,M); %产生 M 个取值 0，1 等概的随机码s=zeros(1,N); %产生一个 N 个元素的零序列for ii=1:Rt*Ts/dts(ii+0:M-1*L)=a;%产生单极性归零码endQ=t2f(s); %付氏变换P=Q.*conj(Q)/T; %P 为单极性归零码的功率 EP=(EP*(jj-1)+P)/jj; %累计平均end aa=30+10*log10(EP+eps);%加 eps 以避免除以零figure
5、(1)set(1,Position,10,350,600,200)%设定窗口位置及大小 figure(2)set(2,Position,10,50,600,200)%设定窗口位置及大小figure(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. 程序结果：输入取样点数=2k,k=14，得到如下波形：5. 实验结果分析由上图结果可以看单极性归零码序列的功率谱密度不仅含有离散的直流分量及连续谱,而且还包含离散的时钟分量和奇次谐波分量，功率谱主瓣宽度为码元速率 2Rb Mb/s。根据数字 PAM 信号功率谱密度公式得到：功率谱分为两个部分，第222()()()aaSTTmssssmPfGf f一部分是连续谱，形状取决于 GT(f);第二部分是离散线谱，相邻线谱频率间隔为 1/Ts。若序列的均值 ma 为零,则第二部分为零，即离散线谱消失，单极性码的均值不为零，故都存在直流分量。
