Abbildungen - Kapitel 8

	options = odeset('RelTol',1e-10);
	plot(0,0,'ko','MarkerFaceColor', 'k');	hold on;
	[t,y] = ode45('pendelnl_f',[0 5.3], [3 0],options);
	pfeilportrait(y,'k',[0.722],0.5,0.7);
	[t,y] = ode45('pendelnl_f',[0 2.7], [2 0],options);
	pfeilportrait(y,'k',[0.722],0.5,0.7);
	[t,y] = ode45('pendelnl_f',[0 2.2], [1.1 0],options);
	pfeilportrait(y,'k',[0.74],0.3,0.7);
	[t,y] = ode45('pendelnl_f',[0 2.1], [0.2 0],options);
	pfeilportrait(y,'k',[],0.2,0.7);
	axis([-3.5 3.5 -6.6 6.6])
	
mit der Funktion pfeilportrait.m.
	options = odeset('RelTol',1e-10);
	plot(0,0,'ko','MarkerFaceColor', 'k');	hold on;
	[t,y] = ode45('pendelnlk_f',[0 20], [1 0],options);
	pfeilportrait(y,'k',[0.07 0.15],0.4,0.2);
	axis([-1.2 1.2 -3 3])

	function y = pendelnlk_f(t,x)
	
		g = 9.81;
		k = 0.5;
		y = zeros(2,1);
		y(1) = x(2);
		y(2) = -k.*x(2) - g.*sin(x(1));
					
	options = odeset('RelTol',1e-10);
	plot(0,0,'ko','MarkerFaceColor', 'k');	hold on;
	plot(0.1,0,'k*');
	[t,y] = ode45('pendelks_ohneF_f',[0 0.6], [0.1 0 0 0],options);
	pfeilportrait(y(:,1:2),'--k',[0.49],0.5,0.7);
	[t,y] = ode45('pendelks_mitF_f',[0 10], [0.1 0 0 0],options);
	pfeilportrait(y(:,1:2),'k',[0.2],0.5,0.7);
	axis([-0.1 0.3 -0.2 1])
	
mit den Funktionen pendelks_ohneF_f.m und pendelks_mitF_f.m



→ Zurück zur Homepage www.dgl-buch.de

© 2008 Lars Grüne, Oliver Junge