Teilchenbewegungen in el./magn. Feldern (Visualisierung)

D.2 eqMotion MATLAB
% digkeit kann bei Aenderung der
% relativen Einheiten auf 1
% gesetzt werden
% diese Variablen werden an dieser Stelle nur fuer eventuelle Erweiterungen
% bzw. Anwendungen auf spezielle Felder eingefuehrt, durch die einzelne
% Feldstaerke, die im GUI eingegeben wird kann die Staerke der Felder direkt
% gesteuert werden
B0 = 1;
E0 = 1;
% globale Variablen
global ode
global bFieldTyp bfieldstrength
global eFieldTyp efieldstrength
global nLarmor
% ODE Funktion
% Variablen fuers B-Feld
% Variablen fuers E-Feld
% Anzahl der Larmorradien
% Kontroll-Funktion des initialisierenden Aufrufs
switch flag
case ’’
% liefert dv/dt = f(t,v) zurueck
varargout{1} = f(t,v);
case ’init’
% liefert [tspan,y0,options] zurueck
[varargout{1:3}] = init;
otherwise
error([’Unknown flag ’’’ flag ’’’.’]);
end
function dydt = f(t,v)
% ===== Bewegungsgleichung =================================================
% setzen der Globalen Variablen
global ode
global bFieldTyp bfieldstrength
global eFieldTyp efieldstrength
global charge mass B0 E0 R0 m0 q
%
[E,B] = feval(’EandB’,v(1),v(2),v(3),t);
% Hier ist die Newton’sche Bewegungsgleichung in vektorieller Form
% dargestellt; die Positionen v(1:3,:) stellen die Position des Teilchens und
% v(4:6,:) die Geschwindigkeit des Teilchens dar
dydt = [ v(4,:)
v(5,:)
v(6,:)
charge/mass*( E(1) +(v(5,:).*B(3) -v(6,:).*B(2)) )
charge/mass*( E(2) +(v(6,:).*B(1) -v(4,:).*B(3)) )
charge/mass*( E(3) +(v(4,:).*B(2) -v(5,:).*B(1)) )
];
% Diese Funktion stellt die einzelnen Feldstrukturen zusammen
function [E,B] = EandB(x,y,z,t)
global charge mass B0 E0 R0 m0
global bFieldTyp bfieldstrength
global eFieldTyp efieldstrength
% B-FELD: Definition der einzelnen Feldstrukturen
switch num2str(bFieldTyp)
71