单极性归零(RZ)波形及其Matlab仿真.docVIP

单极性归零（RZ）波形及其 1. 仿真原理及思路 功率信号的平均功率：。该信号的双边功率谱密度为：，其中是截短后的傅氏变换，是的能量谱，是在截短时间内的功率谱。对于仿真系统，若x是时域取样值矢量，X是对应的傅氏变换，那么x的功率谱便为矢量。 对于采用归零（RZ）及不归零（NRZ）矩形脉冲波形的数字信号，可以用简单的方法信号矢量s。设a是码元矢量，N是总取样点数，M是总码元数，L是每个码元内的点数，那么NRZ信号可这样获得： s=zeros(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的功率谱密度，并画出其波形。 2. 程序流程图 3. 仿真程序及运行结果 仿真程序： %实验二：画出单极性归零码及其功率谱 close all clear all global dt t df N N=2^14; %采样点数 L=64; %每码元的采样点数 M=N/L; %码元数 Rb=2; %码速率为2Mb/s Ts=1/Rb; %码元间隔 dt=Ts/L; %时域采样间隔 Rt=0.5; %占空比 df=1/(N*dt); %频域采样间隔 T=N*dt; %截短时间 Bs=N*df/2; %系统带宽 t=linspace(-T/2,T/2,N); %时域横坐标 f=linspace(-Bs,Bs,N); %频域横坐标 EP=zeros(1,N); for jj=1:100 a=round(rand(1,M)); %产生M个取值0，1等概的随机码 s=zeros(1,N); %产生一个N个元素的零序列 for ii=1:Rt*Ts/dt s(ii+[0:M-1]*L)=a;%产生单极性归零码 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])%设定窗口位置及大小 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 on figure(2) plot(t,s,r) xlabel(t(ms)) ylabel(s(t)(V)) axis([-10,+10,1.1*min(s),1.1*max(s)]) grid on 4. 程序结果： 输入取样点数=2^k,k=14，得到如下波形： 5. 实验结果分析 由上图结果可以看单极性归零码序列的功率谱密度不仅含有离散的直流分量及连续谱,而且还包含离散的时钟分量和奇次谐波分量，功率谱主瓣宽度为码元速率2Rb Mb/s。根据数字PAM信号功率谱密度公式得到：功率谱分为两个部分第一部分是连续谱形状取决于GT(f);第二部分是离散线谱相邻线谱频率间隔为1/Ts若序列的均值ma为零,则第二部分为零即离散线谱消失单极性码的均值不为零故都存在直流分量。 生 产生RZ 计算功率谱，并画出波形 画出波形 RZ波形及其功率谱密度仿真流程图