单极性归零波形及其Matlab仿真

1、单极性归零(RZ)波形及其仿真1.仿真原理及思路P. = I P(f)df = lim f(t) dt功率信号f(t)的平均功率：1T*2TJ_T V)。该信号的双边功率谱密度P(f)为：P(f) = liniJ52L其中Sy(/)是$(/)截短后的傅氏变换, |Sr(/)|2是5应)的能量谱，已是片在截短时间内的功率谱。对于仿真系 统，若x是时域取样值矢量，X是对应的傅氏变换，那么x的功率谱便为矢量 P=(X.*conj(X)/T o对于采用归零(RZ)及不归零(NRZ)矩形脉冲波形的数字信号，可以用简 单的方法信号矢量s。设a是码元矢量，N是总取样点数，M是总码元数，L是 每个码元内的点数

2、，那么NRZ信号可这样获得：s=zeros(l,N);for ii=l:L, s(ii+O:M-l*L)=a; end对于，若Rt是要求的占空比，dt是仿真系统的时域采样间隔，则RZ信号的 产生方法是：s=zeros(l,N);for ii=l:Rt/dt/ s(ii+O:M-l*L)=a; end首先，利用rand函数产生一个尺寸为1*M的矩阵，其元素按均匀分布随机 取值于区间0,1,并用round函数对其四舍五入，得到一个随机产生的0,1序列; 其次，利用for循环产生码元长度为L点，码元为0?且占空比为50%的单极性 归零码，并画出其波形图；最后，计算该RZ的功率谱密度，并画出其波形。2

3、.程序流程图产生0,1分布的随机序列 卜RZ波形及其功率谱密度仿真流程图3仿真程序及运行结果 仿真程序： %实验二：画岀单极性归零码及其功率谱 close allclear allglobal dt t df NN=2A14;%釆样点数L=64;%每码元的釆样点数M=N/L;%码元数Rb=2;%码速率为2Mb/sTs=l/Rb;%码元间隔dt=Ts/L;%时域采样间隔/Rt=;%占空比df=l/(N*dt);%频域采样间隔 T=N*dt;%截短时间Bs=N*df/2; %系统带宽 t=linspace(-T/2,T/2,N);% 时域横坐标f=linspace(-Bs,Bs/N); %频域横坐

4、标 EP=zeros(l,N);forjj=l:100a=round(rand(lzM); %产生M个取值0, 1等概的随机码 s=zeros(l,N);%产生一个N个元素的零序列for ii=l:Rt*Ts/dts(ii+0:M-l*L)=a;%产生单极性归零码endQ=t2f(s);%付氏变换P=Q.*conj(Q)/T;%P为单极性归零码的功率EP=(EP*(jj-l)+P)/jj; %累计平均endaa=30+10*logl0(EP+eps);%加 eps 以避免除以零 figure(l)set(l,Position,10,350,600,200)% 设定窗口位置及大小figure(2

5、)setfPositionJlO,50,600,200)% 设定窗口位置及大小 figure(l)plotftaa/g)xlabel(f(MHZ)ylabel(Ps(f)axis(-15,+15z *min(aa), *max(aa)grid onfigure(2)plot(t,s,F)(xlabel(t(ms)ylabefs(t)(V)axis(-10,+10,*min(s),*max(s)grid on4.程序结果：输入取样点数=2Ak,k=14,得到如下波形:山上图结果可以看单极性归零码序列的功率谱密度不仅含有离散的直流分 量及连续谱，而且还包含离散的时钟分量和奇次谐波分量，功率谱主瓣宽度为码 元速率2Rb Mb/so根据数字PAM信号功率谱密度公式&=单叩卅+竽 G岸）得到：功率谱分为两个部分，第一 部分是连续谱，形状