用matlab实现利用双MZI结构产生NRZ-OOK信号

clear
%% 定义参数
Define_Parameters

global Vpi BitNum SampBit N T F
global  Pcw VmoduN VbiasN

%% 输入信号
Code_In=randint（1BitNum）;
a=Code_In*2-1;

%% NRZ
%Modulation
VinN = VmoduN*a;
E_Out_NRZ=ones（1N）;
for k=1:BitNum
    for j=1:SampBit
        E_Out_NRZ（（k-1）*SampBit+j）=sqrt（Pcw）*cos（pi/2/Vpi*（VinN（k）-VbiasN））...
            *exp（1i*pi*VbiasN/2/Vpi）;
    end
end
plot（T*10^9abs（E_Out_NRZ））                                                %时域图
axis（[0.5 1.5 -0.5 2.5]）
title（‘NRZ-OOK 时序图‘）
xlabel（‘t（ns）‘）
ylabel（‘E‘）
E_F_NRZ=fftshift（fft（E_Out_NRZ））;
figure
plot（F10*log10（abs（E_F_NRZ）.^2））                                          %频域图
% axis（[-200 200 5 -250]）
title（‘NRZ-OOK 频谱图‘）
xlabel（‘f（GHz）‘）
ylabel