1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
-
+
+
+
+
-
-
-
-
-
-
-
+
+
+
+
+
-
+
-
-
+
+
-
-
+
+
-
+
-
+
-
-
+
+
-
+
-
+
-
+
-
+
-
+
-
+
-
-
+
+
-
+
-
+
-
+
-
+
-
|
function makedef(structure,deffile);
function makedef(structure,deffilenum);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.5 1996/11/09 21:05:44 peterg
% %% Only generates MTTIm when at least 2 states!
% %%
% %% Revision 1.4 1996/08/30 19:42:36 peter
% %% Added newline at end of file.
% %%
% %% Revision 1.3 1996/08/24 15:06:22 peter
% %% Write `END;' at end to please reduce.
% %%
% %% Revision 1.2 1996/08/18 20:05:20 peter
% %% Put unded version control
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
filenum = fopen(deffile,'w');
fprintf(filenum, 'MTTNx := %1.0f;\n', states);
fprintf(filenum, 'MTTNz := %1.0f;\n', nonstates);
fprintf(filenum, 'MTTNu := %1.0f;\n', inputs);
fprintf(filenum, 'MTTNy := %1.0f;\n', outputs);
fprintf(filenum, 'MTTNyz := %1.0f;\n', zero_outputs);
fprintf(deffilenum, 'MTTNx := %1.0f;\n', states);
fprintf(deffilenum, 'MTTNz := %1.0f;\n', nonstates);
fprintf(deffilenum, 'MTTNu := %1.0f;\n', inputs);
fprintf(deffilenum, 'MTTNy := %1.0f;\n', outputs);
fprintf(deffilenum, 'MTTNyz := %1.0f;\n', zero_outputs);
% Declare reduce matrices
fprintf(filenum, '%s Declare reduce matrices\n', pc);
fprintf(deffilenum, '%s Declare reduce matrices\n', pc);
if states>0
fprintf(filenum, 'matrix MTTx(%1.0f,1);\n', states);
fprintf(filenum, 'matrix MTTdx(%1.0f,1);\n', states);
fprintf(deffilenum, 'matrix MTTx(%1.0f,1);\n', states);
fprintf(deffilenum, 'matrix MTTdx(%1.0f,1);\n', states);
end;
if nonstates>0
fprintf(filenum, 'matrix MTTz(%1.0f,1);\n', nonstates);
fprintf(filenum, 'matrix MTTdz(%1.0f,1);\n', nonstates);
fprintf(deffilenum, 'matrix MTTz(%1.0f,1);\n', nonstates);
fprintf(deffilenum, 'matrix MTTdz(%1.0f,1);\n', nonstates);
end;
if inputs>0
fprintf(filenum, 'matrix MTTu(%1.0f,1);\n', inputs);
fprintf(deffilenum, 'matrix MTTu(%1.0f,1);\n', inputs);
end;
if outputs>0
fprintf(filenum, 'matrix MTTy(%1.0f,1);\n', outputs);
fprintf(deffilenum, 'matrix MTTy(%1.0f,1);\n', outputs);
end;
if zero_outputs>0
fprintf(filenum, 'matrix MTTyz(%1.0f,1);\n', zero_outputs);
fprintf(filenum, 'matrix MTTui(%1.0f,1);\n', zero_outputs);
fprintf(deffilenum, 'matrix MTTyz(%1.0f,1);\n', zero_outputs);
fprintf(deffilenum, 'matrix MTTui(%1.0f,1);\n', zero_outputs);
end;
% Make an Nx x Nx unit matrix
if states>0
fprintf(filenum, 'matrix MTTI(%1.0f,%1.0f);\n', states,states);
fprintf(deffilenum, 'matrix MTTI(%1.0f,%1.0f);\n', states,states);
for i = 1:states
fprintf(filenum, 'MTTI(%1.0f,%1.0f) := 1;\n', i, i);
fprintf(deffilenum, 'MTTI(%1.0f,%1.0f) := 1;\n', i, i);
end
end;
% Make an Nx/2 x Nx/2 unit matrix
if states>1
fprintf(filenum, 'matrix MTTIm(%1.0f,%1.0f);\n', states/2,states/2);
fprintf(deffilenum, 'matrix MTTIm(%1.0f,%1.0f);\n', states/2,states/2);
for i = 1:states/2
fprintf(filenum, 'MTTIM(%1.0f,%1.0f) := 1;\n', i, i);
fprintf(deffilenum, 'MTTIM(%1.0f,%1.0f) := 1;\n', i, i);
end
end;
% Set the y, yz, u, x and dx matrices
fprintf(filenum, '%s Set the y, yz, u and x matrices\n', pc);
fprintf(deffilenum, '%s Set the y, yz, u and x matrices\n', pc);
for i=1:outputs
fprintf(filenum, 'MTTy(%1.0f,1) := MTTy%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTy(%1.0f,1) := MTTy%1.0f;\n', i, i);
end;
for i=1:zero_outputs
fprintf(filenum, 'MTTyz(%1.0f,1) := MTTyz%1.0f;\n', i, i);
fprintf(filenum, 'MTTui(%1.0f,1) := MTTui%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTyz(%1.0f,1) := MTTyz%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTui(%1.0f,1) := MTTui%1.0f;\n', i, i);
end;
for i=1:inputs
fprintf(filenum, 'MTTu(%1.0f,1) := MTTu%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTu(%1.0f,1) := MTTu%1.0f;\n', i, i);
end;
for i=1:states
fprintf(filenum, 'MTTx(%1.0f,1) := MTTx%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTx(%1.0f,1) := MTTx%1.0f;\n', i, i);
end;
for i=1:nonstates
fprintf(filenum, 'MTTdz(%1.0f,1) := MTTdz%1.0f;\n', i, i);
fprintf(deffilenum, 'MTTdz(%1.0f,1) := MTTdz%1.0f;\n', i, i);
end;
fprintf(filenum, 'END;\n');
fprintf(deffilenum, 'END;\n');
fclose(filenum);
|