%??單極性非歸零碼及其功率譜
?
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）;?