《单极性归零波形及其Matlab仿真.docx》由会员分享，可在线阅读，更多相关《单极性归零波形及其Matlab仿真.docx（3页珍藏版）》请在三一办公上搜索。
1、单极性归零(RZ)波形及其仿真1. 仿真原理及思路P =。P(f)df = lim +T f (t)2dt功率信号f(t)的平均功率：f _8j/T _t 。该信号的双边功率谱密度P(f)为：P(f) = limlT(I)2其中s (f)是s(t)截短后的傅氏变换， T T ， tJS (f )|2是s (t)的能量谱，冬整是s (t)在截短时间内的功率谱。对于仿真系 tTTT统，若X是时域取样值矢量，X是对应的傅氏变换，那么x的功率谱便为矢量P=(X.*conj(X)/T 。对于采用归零(RZ)及不归零(NRZ)矩形脉冲波形的数字信号，可以用 简单的方法信号矢量s。设a是码元矢量，N是总取样
2、点数，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%的单极 性归零码，并画出其波形图；
3、最后，计算该RZ的功率谱密度，并画出其波形。2. 程序流程图RZ波形及其功率谱密度仿真流程图3. 仿真程序及运行结果仿真程序：%实验二：画出单极性归零码及其功率谱close allclear allglobal dt t df NN=2A14;%采样点数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/2,T/2,N); % 时域横坐标f=linspace(-Bs,Bs
4、,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; % 累计平均endaa=30+10*log10(EP+eps);%加 eps 以避免除以零figure(1)set(1,Position,10,350,600,200)% 设定窗口位置及
5、大小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，得到如下波形:由上图结果可以看单极性归零码序列的功率谱密度不仅含有离散的直流分 量及连续谱，而且还包含离散的时钟分量和奇次谐波分量，功率谱主瓣宽度为码 元速率2Rb Mb/s。根据数字PAM信号功率谱密度公式P (f)=姓G (f )|2 + m党G (m)25 (f -竺)得到：功率谱分为两个部分，第一 s T TT 2T TTss m =3ss部分是连续谱，形状取决于6书);第二部分是离散线谱，相邻线谱频率间隔为 1/Ts。若序列的均值ma为零，则第二部分为零，即离散线谱消失，单极性码的均 值不为零，故都存在直流分量。
