p_of_rz.m

单极性归零码的波形图及功率谱 一般信道,无噪声

global dt df t f N
close all
%初始化
N=2^16;%总取样点数
L=64;%每码元的采样点数
M=N/L;%码元数
Rb=2;%码率
Rt=0.5;%占空比
Ts=0.5;%码元宽度
dt=Ts/L;%系统时域采样间隔
df=1/(N*dt);%系统频域采样间隔
T=N*dt;%截断时间
Bs=N*df/2;

t=[-T/2+dt/2:dt:T/2];
f=[-Bs+df/2:df:Bs];

figure(1)
  set(1,'position',[10,50,300,200]);
figure(2)
  set(2,'position',[400,50,300,200]);
  
EP=zeros(size(f));
for ii=1:50
    signal_rz=round(rand(1,M));%产生信息码元
    rz=zeros(1,N);%全零序列以备附值
    for jj=1:(Ts*Rt)/dt, rz(jj+[0:M-1]*L)=signal_rz;end%产生单极性RZ码
    RZ=t2f(rz);%傅立叶变换
    P=RZ.*conj(RZ)/T;%求功率谱
    EP=(EP*(ii-1)+P)/ii;
    
    figure(1)
    aa=30+10*log10(EP+eps);
    plot(f,aa,'b');
    grid
    axis([-20,+20,-50,+50]);
    xlabel('f (MHz)');
    ylabel('Ps(f)  (dbm/MHz)');%作出功率谱密度随频率变化图
    
    figure(2)
    plot(t,rz,'g')
    grid
    axis([-2,+2,-1.5,+1.5]);
    xlabel('t (us)');
    ylabel('s(t) (V)');%作出RZ信号时域波形图
end