单极性不归零码及频谱

%  单极性非归零码及其功率谱
 
global dt df t f N   %全局变量
close all
k=input（‘取样样点=2^kk>=[14]‘）
if k==[]k=14;end
N=2^k;				 %采样点数
L=64;			%每码元的采样点数
M=N/L;			%码元数
Rb=2;			%码速率是2Mb/s
Ts=1/Rb;		%码元间隔
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];	 %频域横坐标

EP=zeros（size（f））;
for jj=1:100
a=round（rand（1M））;	 %码元矢量
s=zeros（1N）;    %产生一个N个元素的零序列
for ii=1:L
   s（ii+[0:M-1]*L）=a;
end					%产生NRZ码
S=T2F（s）;