Abbildungen - Kapitel 7

	hold on;
	h=plot([-pi,0,pi],0,'k.'); 	set(h,'markersize',40);
	pdl = inline('[x(2);-9.81*sin(x(1))]','t','x');
	opt=odeset('RelTol',1e-5);
	for x1=0.4:0.3:3,
	    [t,x]=ode45(pdl,[0 10],[x1;0],opt);
	    h=plot(x(:,1),x(:,2),'k-');
	    set(h,'linewidth',2);
	end
					
	pdl = inline('[x(2);-9.81*sin(x(1))]','t','x');
	opt=odeset('RelTol',1e-5);
	for x1=2.2:0.3:3,
	    [t,x]=ode45(pdl,[0 10],[x1;0],opt);
	    h=plot(x(:,1),x(:,2),'k-');
	    set(h,'linewidth',2);
	end
	axis([2 4 -4 4]);
	
	[t,x]=ode45(pdl,[0 -4],[3.14;0],opt);
	h=plot(x(:,1),x(:,2),'k-');	set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 -4],[3.142;0],opt);
	h=plot(x(:,1),x(:,2),'k-');	set(h,'linewidth',4);
	
	[t,x]=ode45(pdl,[0 4],[3.14;0],opt);
	h=plot(x(:,1),x(:,2),'color',[0.6 0.6 0.6]); set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 4],[3.142;0],opt);
	h=plot(x(:,1),x(:,2),'color',[0.6 0.6 0.6]); set(h,'linewidth',4);
	
	h=plot(pi,0,'k.');
	set(h,'markersize',40);
					
	pdl = inline('[x(2);-9.81*sin(x(1))]','t','x');
	opt=odeset('RelTol',1e-5);
	axis([-4 4 -10 10])
	
	[t,x]=ode45(pdl,[0 -5],[3.14;0],opt);
	h=plot(x(:,1),x(:,2),'k-');	set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 -4],[3.142;0],opt);
	h=plot(x(:,1),x(:,2),'k-');	set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 4],[-3.142;0],opt);
	h=plot(x(:,1),x(:,2),'k-');	set(h,'linewidth',4);
	
	[t,x]=ode45(pdl,[0 5],[3.14;0],opt);
	h=plot(x(:,1),x(:,2),'color',[0.6 0.6 0.6]); set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 4],[3.142;0],opt);
	h=plot(x(:,1),x(:,2),'color',[0.6 0.6 0.6]); set(h,'linewidth',4);
	[t,x]=ode45(pdl,[0 -4],[-3.142;0],opt);
	h=plot(x(:,1),x(:,2),'color',[0.6 0.6 0.6]); set(h,'linewidth',4);
	
	h=plot(pi,0,'k.'); set(h,'markersize',40);
	h=plot(-pi,0,'k.');	set(h,'markersize',40);
					



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

© 2008 Lars Grüne, Oliver Junge