用马赫-曾德尔调制器（MZM）进行信号调制，实现了OOK，BPSK，DPSK，QPSK，DQPSK，两种16QAM的调制方式

phase=pi./2; %调制器2 e_out2是RZ波形的输出 ff=Rb/2; s=sin(2*pi*ff*t./1000+phase); %注意t的单位,以及正弦信号的初始相位和振荡频率 %%%%%%%%%%%%% Vbias2=0; Vin2=V_pi.*s; %%%%%%%%%%%%% e_out2=e_out1.*cos((pi*(Vin2-Vbias2)/(2*V_pi)))*exp(1i*pi*Vbias2/(2*V_pi)); % e_out2是RZ的输出 RZ-33 Vbias=0; Vin=V_pi.*sin(2*pi*ff*t./1000+phase); RZ-50 Vbias=V_pi/2; Vin=V_pi/2.*sin(2*pi*ff*t./1000+phase); RZ-66 Vbias=V_pi; Vin=V_pi.*sin(2*pi*ff*t./1000+phase); %%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%% Vbias1=V_pi/2; Vin1=V_pi/2.*nrz2； %%%%%%%%%%%%%%%%%%% e_out1=A_modulator*cos(pi*(Vin1-Vbias1)/(2*V_pi))*exp(1i*pi*Vbias1/(2*V_pi)); NRZ-OOK Vbias=V_pi/2; Vin=V_pi/2.*nrz2; NRZ-BPSK Vbias=V_pi; Vin=V_pi.*nrz2; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%差分计算 for ii=1:M-1 inp(ii+1)=bitxor(inp(ii+1),inp(ii)); %数据位异或计算，返回值不是0就是1