%   Script ZUSTTRAJ.M
%   
%   Trajektorie im dreidimensionalen Zustandsraum
%
% J. Lunze  
% 9.10.1995
% für 2. Auflage: 6.11.1998
% für 3. Auflage: 3.4.2001
%
echo off
clear
close all;
A = [-0.5 0 0 ;
     0 -0.3 2;
     0 -2 -0.3];
b = [1;
     1;
     1];
c = [1 1 1];
d=0;
System=ss(A, b, c, d);
x0=[10, 10, 10];
echo on
%
%   Eigenbewegung eines Systems dritter Ordnung
%   dreidimensionale Darstellung im Zustandsraum
%
echo off
figure(1);
T=0:0.01:10;
[Y, T, X]=initial(System, x0, T);
plot3(X(:, 2), X(:, 3), X(:, 1));
hold on
plot3([-10 10], [0 0], [0 0], ':');
plot3([0 0], [-10 10], [0 0], ':');
plot3([0 0], [0 0], [10 0], ':');
xlabel('x_1');
ylabel('x_2');
zlabel('x_3');
grid('off');
title('Trajektorie eines Systems dritter Ordnung im Zustandsraum');
%text(X(1,2), X(1,3), X(1,1), '  x(0)')
text(6, 10, 10, 'x(0)')
axis([-10 10 -10 10 0 10]);
hold off

%
%   .EPS und .PS Datei
%
% epsbild('zusttraj', 'zusttraj');