Teilchenbewegungen in el./magn. Feldern (Visualisierung)

MATLAB
und
ẋ n+1 = ẋ n + h
2m [F (x n+1) + F (x n )]
(C.3)
wobei x die Koordinate und ẋ die Geschwindikeit, m die Masse und F(x) die Kraft,
welche auf die Masse m am Punkt x wirkt.
Im ersten Anschein sieht diese Methode wie eine Prädiktormethode aus, weil hier
der Term x n+1 auftaucht. Dieser wird aber nicht wie bei den Prädiktor Methoden
approximativ bestimmt sondern in Gleichung (C.2) zunächst bestimmt bevor er in
Gleichung (C.3) verwendet wird.
Der Vorteil dieses Algorithmusses liegt darin, dass er über lange Zeiträume integrieren
kann, ohne das der Fehler zu groß wird.
D MATLAB
In diesem Abschnitt sind die MATLAB r○ Dateien abgedruckt. Die Programme sowie
die HTML-Dokumentation befinden sich auf der beigefügten CD-Rom.
D.1 particles
function varargout = particles(varargin)
% !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% BEGIN HEAD (nicht aendern)
% !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Initialisierung des GUI und der mit ihm verbundenen Callbacks
gui_Singleton = 1;
gui_State = struct(’gui_Name’, mfilename, ...
’gui_Singleton’, gui_Singleton, ...
’gui_OpeningFcn’, @particles_OpeningFcn, ...
’gui_OutputFcn’, @particles_OutputFcn, ...
’gui_LayoutFcn’, [] , ...
’gui_Callback’, []);
if nargin & isstr(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% Hier werden Einstellungen getroffen, bevor das GUI sichtbar wird
function particles_OpeningFcn(hObject, eventdata, handles, varargin)
handles.output = hObject;
guidata(hObject, handles);
% Anzeige des Titelbildes und des Logos durch Axis-Objekte
axes(handles.logo)
[a,map]=imread(’logo.png’); image(a); colormap(map); axis equal, axis off
axes(handles.plot_axis)
[a,map]=imread(’title.png’); image(a); colormap(map); axis equal, axis off
function varargout = particles_OutputFcn(hObject, eventdata, handles)
61