Overview
Comment: | Handles complex components and repetative components. |
---|---|
Downloads: | Tarball | ZIP archive | SQL archive |
Timelines: | family | ancestors | descendants | both | origin/master | trunk |
Files: | files | file ages | folders |
SHA3-256: |
a42aee795cf1e439ce2ae9eeb3190b66 |
User & Date: | gawthrop@users.sourceforge.net on 1996-08-27 08:04:52 |
Other Links: | branch diff | manifest | tags |
Context
1996-08-27
| ||
08:08:44 | Initial revision check-in: a7f1bb8b6c user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
08:04:52 | Handles complex components and repetative components. check-in: a42aee795c user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
1996-08-26
| ||
13:25:07 | Corrected error handling. check-in: 23a7b08827 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Modified mttroot/mtt/bin/trans/m/cbg2ese.m from [b5879f7a9f] to [0dd01471dd].
|
| | > > > > > > > < > | < < < < | > | > > > > > | > > > > > > > > > > > > > > > > > > > > > > > > > > > > > | < | 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 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 | function structure = cbg2ese(system_name, system_type, full_name, ... repetition,... structure, infofile) % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %%%%% Model Transformation Tools %%%%% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Matlab function cbg2ese.m % Acausal bond graph to causal bond graph: mfile format % Structure matrix [states,nonstates,inputs,outputs,zero_outputs] % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ % %% Revision 1.5 1996/08/24 15:02:23 peter % %% Writes `END;' to keep reduce happy. % %% % %% Revision 1.4 1996/08/19 09:03:41 peter % %% Handles repeating components. % %% % %% Revision 1.3 1996/08/18 20:08:02 peter % %% Included additional structure: structure(5) = zero_outputs. % %% % %% Revision 1.2 1996/08/08 18:08:11 peter % %% Sorted out file naming sceme % %% % %% Revision 1.1 1996/08/08 15:53:23 peter % %% Initial revision % %% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% system_name, system_type, full_name, pc = '%'; if nargin<3 eqnfile = 'stdout'; end; if nargin<4 infofile = 'stdout'; end; % Create the (full) system name if length(full_name)==0 full_name = system_name; system_type = system_name; else full_name = [full_name, '_', system_name]; end; full_name_type = [full_name, '_', system_type]; cbg_name = [full_name_type, '_cbg']; abg_name = [system_type, '_abg']; cmp_name = [system_type, '_cmp']; % Return if cbg file doesn't exist if exist(cbg_name)~=2 return end; % Setup file ese_name = sprintf('%s_%1.0f_ese.r', full_name, repetition); filenum = fopen(ese_name, 'w'); % Evaluate the system function to get the bonds eval(['[junk,components]=', abg_name, ';']); eval(['bonds=', cbg_name, ';']); % Find number of bonds [n_bonds,columns] = size(bonds); if (columns ~= 2)&(n_bonds>0) error('Incorrect bonds matrix: must have 2 columns'); end; % Find number of components [n_components,columns] = size(components); n_components = n_components for i = 1:n_components comp = nozeros(components(i,:)); bond_list = abs(comp); direction = sign(comp)'*[1 1]; % Convert from arrow orientated to component orientated causality comp_bonds = bonds(bond_list,:).*direction; eval([ '[comp_type,comp_name,cr,args,repetitions] = ', cmp_name, '(i);' ]); % change name of 0 and 1 components -- matlab doesn't like numbers here if strcmp(comp_type,'0') comp_type = 'zero'; end; if strcmp(comp_type,'1') comp_type = 'one'; end; comp_type = comp_type ports = length(bond_list) repetitions = repetitions if repetitions>1 port_pairs = ports/2; if round(port_pairs)~=port_pairs; mtt_info(['Repeated component ', comp_name, ... ' has an odd number of ports - ignoring repetitions']); repetitions = 1; end; end; if repetitions>1 odd_bonds = bond_list(1:2:ports-1); even_bonds = bond_list(2:2:ports); next_bond = max(max(abs(components)))+1; end; for k = 1:repetitions if repetitions>1 if k==1 bond_list(1:2:ports-1) = odd_bonds; else bond_list(1:2:ports-1) = bond_list(2:2:ports); end; if k==repetitions bond_list(2:2:ports) = even_bonds; else new_bonds = [next_bond:next_bond+port_pairs-1]; next_bond = next_bond+port_pairs; bond_list(2:2:ports) = new_bonds; end; end; % Invoke the appropriate equation-generating procedure name_r = sprintf('%s_%1.0f', full_name, repetition); eqn_name = [comp_type, '_eqn'] if exist(eqn_name)~=2 % Try a compound component disp('---PUSH---'); structure = cbg2ese(comp_name, comp_type, full_name, k, ... structure, infofile); % Link up the bonds name_k = sprintf('%s_%1.0f', full_name, k); name_comp_name = sprintf('%s_%s_%1.0f', full_name, comp_name, k); for port_number=1:length(bond_list) % Effort if comp_bonds(port_number,1)==1 % Source fprintf(filenum, '%s_MTT_inport_%1.0f := %s;\n', ... name_comp_name, port_number, varname(name_r, ... bond_list(port_number),1)); else % sensor fprintf(filenum, '%s := %s_MTT_outport_%1.0f;\n', ... varname(name_r, ... bond_list(port_number),1), name_comp_name, port_number); end; % flow if comp_bonds(port_number,2)==-1 % Source fprintf(filenum, '%s_MTT_inport_%1.0f := %s;\n', ... name_comp_name, port_number, varname(name_r, ... bond_list(port_number),-1)); else % sensor fprintf(filenum, '%s := %s_MTT_outport_%1.0f;\n', ... varname(name_r, ... bond_list(port_number),-1), name_comp_name, port_number); end; end; disp('---POP---'); else % its a simple component eval(['structure = ', ... eqn_name, ... '(name_r,bond_list,comp_bonds,direction,cr,args,structure,filenum);' ]); end; end; end; fclose(filenum); |