drz.m

包括RZ

%试验一 RZ的功率谱
%本程序中时间单位是微秒
%频率单位为MHz
%码元速率单位是Mb/s
global dt df t f N
close all
N=2^13;   %采样点数
L=64;    %每码元的采样点数
M=N/L     %码元数
Rb=2;
Ts=0.5		%码元宽度是0.5us
dt=Ts/L       
df=1/(N*dt)     %MHz
T=N*dt          %截短时间
Bs=N*df/2       %系统带宽
Rt=0.5;
t=[-T/2+dt/2:dt:T/2]; %时域横坐标
f=[-Bs+df/2:df:Bs];   %频域横坐标
figure(1)
set(1,'Position',[10,50,500,290])                %设定窗口位置及大小 
figure(2)
 set(2,'Position',[520,50,500,290])                %设定窗口位置及大小 
EP=zeros(size(f));
for jj=1:50
s=zeros(size(f));
a=sign(randn(1,M));
for ii=1:Rt*Ts/dt,s(ii+[0:M-1]*L)=a; end
  RZ=t2f(s);       %RZ信号的傅氏变换
P=RZ.*conj(RZ)/T;        %RZ信号的功率谱
EP=(EP*(jj-1)+P)/jj;       %累计平均
 figure(1)
 aa=30+10*log10(EP+eps);%加eps以避免除以零
 plot(f,aa,'g');
 grid
 axis([-10,+10,-70,70])
 xlabel('f  (MHz)')
 ylabel('Ps(f)   (dBm/MHz)')
 figure(2)
 plot(t,s,'g')
 grid
 axis([-2,+2,-1.5,1.5])
 xlabel('t (us)')
 ylabel('s(t) (V)')
end