%MZM 仿真

%使用推挽式MZM，上下两路使用幅度，相位可控的交流信号

clear all;

nu=20;                       %所取的能带阶数

m=1;

N=100;                       %交流信号幅值范围取100个点

Beitam=linspace(1,3,N);   %Beita为交流信号幅值对应的β值范围1-3

Xm=linspace(0,2*pi,N);    %X为交流信号相位差范围0°到360°

P=[];                        %P为二阶矩阵，第一行为取0阶到20阶边带的功率，第二行为对应的阶数

save=struct('A',[],'B',[],'Beita',[],'X1step',[],'X2step',[],'X3step',[],'X4step',[],'X1',[],'X2',[],'X3',[],'X4',[],'X',[]);

save1=struct('A',[],'B',[],'Beita',[],'X1step',[],'X2step',[],'X3step',[],'X4step',[],'X1',[],'X2',[],'X3',[],'X4',[],'X',[]);

%A,B为射频电压对应β范围,Beita为交流信号幅值对应的β值范围,X1到X4为第一大到第四大的功率，X1step