%產生二進制編碼數據
global dt t f df N T


   close all

   k=input('取樣點數＝2^k, k取10左右');

   if isempty(k), k=13; end

   z=input('每個信號取樣點數＝2^z,z<k');

   if isempty(z), z=5; end

   f0=input('f0=取1(kz)左右');

   if isempty(f0), f0=1; end

   N=2^k;  %一共的取樣點數

   L=2^z;  %每個碼元的取樣點數（偶數）
   M=N/L;  %信號數目 即碼元數目

   Rb=2;

   Ts=0.05                         %碼元寬度是0.05us=50ns

   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));            %產生隨機碼
   
%對隨機碼b每符號抽樣32次產生s
   s=zeros(1,N);

   for ii=1:L;

   s(ii+[0:M-1]*L)=b;

   end
   figure(1)                       %輸出信號波形畫圖

   plot(t,s,'b')

   grid

   axis([-3,3,0,1.3]);

   xlabel('t (ms)');ylabel('s(t) (V)');title('單極性RZ的時域圖')

   end
   %以下程序段產生