Teilchenbewegungen in el./magn. Feldern (Visualisierung)

D.2 eqMotion MATLAB
case ’1’ % B=0
B = B0*[ 0 0 0 ];
case ’2’
% Bz
B = B0*[ 0 0 bfieldstrength ];
case ’3’
% dBz/dy
B = B0*[ 0 0 1+bfieldstrength*y ];
case ’4’
% dBz/dz cusp
B = B0*[ -0.5*x (-1)*bfieldstrength*y bfieldstrength*z];
case ’5’
% B toroidal
% Verwendung von Toroidal Koordinaten
R = sqrt(x^2 +y^2 + 1.e-10); cosphi = x/R; sinphi = y/R;
B = B0*(R0/R)*[ -sinphi cosphi 0 ]*bfieldstrength;
case ’6’
% B tokamak
% Verwendung von Toroidal Koordinaten
R = sqrt(x^2 +y^2 + 1.e-10); cosphi = x/R; sinphi = y/R;
rho = sqrt((R-R0)^2 + z^2+1.e-10); sint = -z/rho; cost = (R-R0)/rho;
Bt = B0 * (R0/R)*[-sinphi cosphi 0];
Bp = B0 *(rho/R)*[-cosphi*sint -sinphi*sint -cost ];
B = Bt +Bp/q;
case ’7’
% B dipole
r2 = x^2 +y^2 +z^2;
B = m0*r2^(-2.5)*[ 3*x*z 3*y*z 3*z^2-r2]*bfieldstrength;
case ’8’
% user defined
%
% ==> Hier kann ein beliebiges Vektorfeld eingefuegt werden
%
otherwise
% Hier wird der Fehler, der schon im GUI abgefangen wird nochmals zur
% Sicherheit ueberprueft
msgbox(’Diese Konfiguration existiert nicht’,’ERROR’,’error’);
end
Btot = norm(B); OmegaC=charge/mass*norm(B);
% E-FELD: Definition der einzelnen Feldstrukturen
switch num2str(eFieldTyp)
case ’1’ % E=0
E = E0*[ 0 0 0 ];
case ’2’
% Ex
E = E0*[ efieldstrength 0 0 ];
case ’3’
% Ez
E = E0*[ 0 0 efieldstrength ];
case ’4’
% dEy/dt
E = E0*[ 0 t*charge/mass*Btot 0]*efieldstrength;
case ’5’
% user defined
%
% ==> Hier kann ein beliebiges Vektorfeld eingefuegt werden
%
otherwise
% Hier wird der Fehler, der schon im GUI abgefangen wird nochmals zur
% Sicherheit ueberprueft
msgbox(’Diese Konfiguration existiert nicht’,’ERROR’,’error’);
end
% diese Funktion setzt die initialisierenden Werte fuer die Berechnung
function [tspan,y0,options] = init
global charge mass B0 E0 R0 m0 nLarmor velscala
global xposv yposv zposv velx vely velz
72