%單極性不歸零碼的時域圖及功率譜密度，“0”碼“1”碼等概出現，只有在零頻率處存在離散譜。
function demo_danjinrz
clear
close all
k=16; %取樣點數＝2^k
z=8;  %每個信號取樣點數＝2^z,z<k
f0=1; %f0=取1(kz)左右
N=2^k;
L=2^z;
M=N/L;
Ts=1/f0; %碼元寬度是1us
dt=Ts/L; %時域采樣間隔
df=1/(N*dt); %MHz
T=N*dt; %截短時間
Bs=N*df/2; %系統帶寬
f=[-Bs+df/2:df:Bs]; %頻域橫坐標
t=[-T/2+dt/2:dt:T/2]; %時域橫坐標
EP=zeros(size(f))+eps;
for ii1=1:30
    b=round(rand(1,M)); %產生隨機碼
    s=zeros(1,N);
    for ii=1:L;%for ii=1:L/2;碼型變為占空比1/2的NRZ.
        s(ii+[0:M-1]*L)=b;
    end
    S=fftshift(fft(s))*dt; %S是s的傅氏變換
    P=S.*conj(S)/T; %功率譜
    EP=(EP*(ii1-1)+P+eps)/ii1;
end
subplot(2,1,2); %輸出的功率畫圖
plot(f,10*log10(EP),'m');
grid
axis([-15,15,-100,20]);
xlabel('f (KHz)');ylabel('P (V/KHz)') ;title('功率譜圖')
subplot(2,1,1); %輸出信號波形畫圖
plot(t,s,'b')
grid
axis([-10,10,-0.2,1.2]);
xlabel('t (ms)');ylabel('s(t) (V)');title('單極性NRZ的時域圖')