clear?all;
Ts???????=?1;????????????????%碼元周期
N_sample?=?1000;??????????????%每個碼元抽樣點數
fs=N_sample;

dt???????=?Ts/N_sample;??????%抽樣的時間間隔
N????????=?2000;???????????????%碼元數
t????????=?0:dt:（N*N_sample-1）*dt;

%!!!!!!!!!!!!!!!!!!!!!!!!!
df???????=?fs/length（t）;
f????????=?-fs/2:df:fs/2-df;?%定義頻率矢量（頻譜圖的橫坐標）

%?基本碼元?g（t）?=?1
gt0??????=?zeros（1N_sample）;
gt1??????=?ones（1N_sample）;??????????????????????????%單極性不歸零
gt2??????=?[ones（1N_sample/2）zeros（1N_sample/2）];??%單極性歸零，占空比為0.5
gt3??????=?gt0-gt1;???????????????????????????????????%雙極性不歸零
gt4??????=?[zeros（1N_sample/2）-ones（1N_sample/2）zeros（1N_sample/2）];%雙極性歸零


%生成隨機碼元
random_code?=?randi（[01]1N）;
%單極性波形的最后結果
yt1??????=?zeros（1N*N_sample-1）;???????%單極性不歸零NRZ波形
yt2??????=?zeros（1N*N_sample-1）;???????%單極性不歸零NRZ波形
yt3??????=?zeros（1N*N_sample-1）;???????%雙極性不歸零NRZ波形
yt4??????=?zeros（1N*N_sample-1）;???????%雙極性歸零RZ波形

for?i?=?1:N
????if?random_code（i）?==?1?????????????????????????????%得到波形圖
????????yt1（（i-1）*N_sample+1:i*N_sample）?=?gt1;????????%單極性不歸零NRZ波形
????????yt2（（i-1）*N_sample+1:i*N_sample）?=?gt2;????????%單極性歸零RZ波形
????????yt3（（i-1）*N_sample+1:i*N_sample）?=?gt1;????????%雙極性不歸零NRZ波形
????????yt4（（i-1）*N_sample+1:i*N_sample）?=?gt2;????????%雙極性歸零RZ波形
????else
????????yt1（（i-1）*N_sample+1:i*N_sample）?=?gt0;????????
????????yt2（（i-1）*N_sample+1:i*N_sample）?=?gt0;
????????yt3（（i-1）*N_sample+1:i*N_sample）?=?gt3;
????????yt4（（i-1）*N_sample+1:i*N_sample）?=?gt4;
????end
end
%*********************************************************************
%求單極性不歸零時域波形于波功率譜密度

figure（1）
plot（tyt1‘LineWidth‘1.5）;
axis（[0?20?-0.1?1.1]）;
title（‘二進制單極性不歸零NRZ波形‘）;
xlabel（‘時間t/s‘）;?????????
grid;

figure（2）
fmt=fft（yt1）;?%?對時域信號進行FFT變換，計算其頻譜
fmt=fftshift（fmt）;
fmt=abs（fmt）;
fmt_dB=fmt.^2/Ts;
maxF=max（fmt_dB）;
fmt_dB=fmt_dB/maxF;
fmt_dB=?10*log10（fmt_dB?+eps）;
plot（ffmt_dB）;grid?on;
axis（[-6?6?-80?0]）;
xlabel（‘頻率（Hz）‘）;ylabel（‘功率譜幅度值（dB）‘）;
title（‘二進制單極性不歸零NRZ波形功率譜密度（db）‘）;

%*********************************************************************
%*********************************************************************
%求單極性歸零時域波形與波功率譜密度
figure（3）
plot（tyt2‘LineWidth‘1.5）;
title（‘二進制單極性歸零RZ波形‘）;
xlabel（‘時間t/s‘）;
axis（[020-0.091.1]）;??????????????%畫出單極性歸零RZ波形功率譜密度
grid;

figure（4）
fmt=fft（yt2）;?%?對時域信號進行FFT變換，計算其頻譜
fmt=fftshift（fmt）;
fmt=abs（fmt）;
fmt_dB=fmt.^2/Ts;
maxF=max（fmt_dB）;
fmt_dB=fmt_dB/maxF;
fmt_dB=?10*log10（fmt_dB?+eps）;
plot（ffmt_dB）;grid?on;
axis（[-6?6?-80?0]）;
xlabel（‘頻率（Hz）‘）;ylabel（‘功率譜幅度值（dB）‘）;
title（‘二進制單極性歸零RZ波形功率譜密度（db）‘）;


%*********************************************************************
%*********************************************************************
%求雙極性不歸零時域波形與波功率譜密度
figure（5）
plot（tyt3‘LineWidth‘1.5）;
title（‘二進制雙極性不歸零波形‘）;
xlabel（‘時間t/s‘）;
axis（[0?20?-1.1?1.1]）?
grid;

figure（6）
fmt=fft（yt3）;?%?對時域信號進行FFT變換，計算其頻譜
fmt=fftshift（fmt）;
fmt=abs（fmt）;
fmt_dB=fmt.^2/T