Mercurial > hg > octave-nkf

annotate scripts/control/abcddim.m @ 707:d7c4962ec7a0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1994-09-16 13:50:20 by jwe]
author jwe
date Fri, 16 Sep 1994 13:53:18 +0000
parents 16a24e76d6e0
children 3470f1e25a79
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
245
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
1 # Copyright (C) 1993 John W. Eaton
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
2 #
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
3 # This file is part of Octave.
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
4 #
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
5 # Octave is free software; you can redistribute it and/or modify it
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
6 # under the terms of the GNU General Public License as published by the
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
7 # Free Software Foundation; either version 2, or (at your option) any
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
8 # later version.
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
9 #
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
10 # Octave is distributed in the hope that it will be useful, but WITHOUT
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
11 # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
12 # FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
13 # for more details.
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
14 #
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
15 # You should have received a copy of the GNU General Public License
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
16 # along with Octave; see the file COPYING. If not, write to the Free
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
17 # Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 78
diff changeset
18
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
19 function [n, m, p] = abcddim (a, b, c, d)
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
20
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
21 # Usage: [n, m, p] = abcddim (a, b, c, d)
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
22 #
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
23 # Check for compatibility of the dimensions of the matrices defining
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
24 # the linear system (a, b, c, d).
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
25 #
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
26 # Returns n = number of system states,
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
27 # m = number of system inputs,
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
28 # p = number of system outputs.
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
29 #
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
30 # Returns n = m = p = -1 if the system is not compatible.
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
31
78
0fda6e1f90e0 [project @ 1993-08-30 15:40:16 by jwe]
jwe
parents: 70
diff changeset
32 # Written by A. S. Hodel (scotte@eng.auburn.edu) August 1993.
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
33
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
34 if (nargin != 4)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
35 error ("abcddim: illegal number of arguments. need four.")
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
36 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
37
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
38 n = m = p = -1;
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
39
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
40 [an, am] = size(a);
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
41 if (an != am)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
42 fprintf (stderr, "abcddim: a is not square");
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
43 return;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
44 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
45
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
46 [bn, bm] = size(b);
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
47 if (bn != am)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
48 fprintf (stderr, "abcddim: a and b are not compatible");
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
49 return;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
50 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
51
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
52 [cn, cm] = size(c);
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
53 if (cm != an)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
54 fprintf (stderr, "abcddim: a and c are not compatible");
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
55 return;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
56 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
57
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
58 [dn, dm] = size(d);
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
59 if (cn != dn)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
60 fprintf (stderr, "abcddim: c and d are not compatible");
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
61 return;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
62 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
63
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
64 if (bm != dm)
70
a5e8b7c9a3ad [project @ 1993-08-30 14:32:55 by jwe]
jwe
parents: 60
diff changeset
65 fprintf (stderr, "abcddim: b and d are not compatible");
60
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
66 return;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
67 endif
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
68
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
69 n = an;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
70 m = bm;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
71 p = cn;
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
72
671f8bf989d8 [project @ 1993-08-13 21:03:40 by jwe]
jwe
parents:
diff changeset
73 endfunction