1 2 3 4 5 6 7 8 9 10 11 12 13
1; %EM field definition Ep=0.01; Epc=conj(Ep); Ed=1; Edc=conj(Ed); wd=w13; wp=w12; modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; %modulation_freq=[0, wp, -wp]; %E_field =[0, Ep, Epc];