轨道角动量仿真

clear;
%----------------constants-------------------
c=1;
%-------------things to change---------------
l_start=0; %start OAM
lam=1.5; %Wavelength
A=0.5; %Amplitude
N=100; %
t_stop=10; %Amount of timesteps
z_range=2*lam; %Length of z values
%--------------------------------------------
k=2*pi/lam;
w=k*c;
k_vector=[1 0 0]*k;
%---------vector x--------
z_value=linspace（0z_rangeN）; %z on x-axis in plot
L_z=length（z_value）;
r=A*ones（1N）;
%----- --------- -------
t_steg=lam/（2*2*pi*c）; %Not to big timesteps
t_vector=0:t_steg:t_stop;
%------------------------
z_antal=z_range/lam;
steg=zeros（1z_antal）; %s? att den plana v?gen upprepas
%d? den passerar utanf?r bild.
t=1;
n=1;
l=l_start;
while t<=t_stop
    while l<=l_start+3
        if l==0 %Plane wave
            clf;
            subplot（22n）;
            hold on;
            axis（[0 z_range -2 2 -2 2]）;
            view（5130）;
            axis fill;
            axis square;
            grid on;
%