 |
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
|