Nun eine passende Integralkurve 3d: octave:89> plot3(circsx,circsy,t); octave:90> plot3(circsx,circsy,t); octave:91> t = linspace(0,4*pi,200)’; octave:92> circsx = cos(t); octave:93> circsy = sin(t); octave:94> plot3(circsx,circsy,t); Voila:
Tag Archives: DGL
Trajektorien einzeichnen
In vorgängig besprochenes Vektorfeld werden Trajektorien zu drei verschiedenen Anfangsbedingungen rot (“cr”, color red) eingezeichnet (x(0)=0,y(0)=5;y(0)=10;y(0)=15): octave:69> t = linspace(0,2*pi,100)’; octave:70> quiver(x,y,-y./z,x./z) octave:71> hold off octave:72> quiver(x,y,-y./z,x./z) octave:73> t = linspace(0,2*pi,100)’; octave:74> circsx = 10.*cos(t); octave:75> circsy = 10.*sin(t); octave:76> hold on octave:77> plot(circsx,circsy,”cr”); octave:78> circsx = 5.*cos(t); octave:79> circsy = 5.*sin(t); octave:80> plot(circsx,circsy,”cr”); octave:81> […]
Normiertes Vektorfeld in octave
In Octave wird ein Vektorfeld für folgendes DGL-System abgebildet: x'(t)=-y(t) y'(t)=x(t) [x, y] = meshgrid (-20:1:20); z=sqrt((-y).^2+x.^2); quiver(x,y,-y./z,x./z); Zunächst werden zwei Matritzen [x,y] generiert. Diese bilden die Grundlage für die einzelnen Punkte im Vektorfeld. Danach wird eine Matrix mit Normierungsfaktoren generiert, indem für jeden Vektor im Vektorfeld der Betrag ausgerechnet wird. (Pythagoras: |v|=sqrt(x^2+y^2)) Schliesslich wird […]