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
|
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
|
-
+
+
+
+
+
+
-
+
-
-
-
-
-
+
+
+
+
+
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
|
function write_cbg(system_name,system_type,system);
function write_cbg(system_name,system_type,system,Flipped)
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.2 1998/08/25 06:21:19 peterg
## Just writes additional information; basic info from the abg structure.
##
## Revision 1.1 1998/08/25 05:55:10 peterg
## Initial revision
##
###############################################################
fid=fopen([system_name,"_cbg.m"], "w");# Open file
StatusFormat = " %s.%s.%s.status = %i;\n";
Bformat = " %s.bonds = [\n";
Iformat = " %s.%s.%s.connections(%i) = %i;\n";
fprintf(fid,"function [%s] = %s_cbg\n", system_name, system_name);
fprintf(fid,"# This function is the causal bond graph representation of %s\n",system_name);
fprintf(fid,"# Generated by MTT on %s",ctime(time));
fprintf(fid,"# The file is in Octave format\n");
fprintf(fid,"\n# Acausal bond graph structure\n");
fprintf(fid," [%s] = %s_abg;\n", system_name, system_name);
fprintf(fid," [%s] = %s_abg;\n", system_name, system_type);
fprintf(fid,"\n# Status information\n");
if struct_contains(system,"ports")
for [port,name]=system.ports
fprintf(fid,StatusFormat,system_name,"ports",name,subsystem.status);
endfor;
endif
# if struct_contains(system,"ports")
# for [port,name]=system.ports
# fprintf(fid,StatusFormat,system_name,"ports",name,port.status);
# endfor;
# endif
if struct_contains(system,"subsystems")
for [subsystem,name]=system.subsystems
fprintf(fid,StatusFormat,system_name,"subsystems",name,subsystem.status);
endfor;
endif
[N,M]=size(system.bonds); # Bonds
fprintf(fid,"\n# Causal bond information\n");
fprintf(fid,Bformat,system_name);
for i=1:N
fprintf(fid," ");
for j=1:M
fprintf(fid,"%i ", system.bonds(i,j));
endfor;
fprintf(fid,"\n");
endfor;
fprintf(fid," ];\n");
[N,M]=size(Flipped.ports);
if N>0
fprintf(fid,"\n# Flipped port information\n");
for i=1:N
name=deblank(Flipped.ports(i,:));
eval(["con = system.ports.",name,".connections;"]);
fprintf(fid,Iformat,system_name,"ports",name,1,con);
endfor;
endif;
[N,M]=size(Flipped.subs);
if N>0
for i=1:N
name=deblank(Flipped.subs(i,:));
k=Flipped.cons(i);
eval(["con = system.subsystems.",name,".connections(",num2str(k),");"]);
fprintf(fid,Iformat,system_name,"subsystems",name,k,con);
endfor;
endif;
fclose(fid);
|