th_dem_s.m

FISMAT accommodates different arithmetic operators, fuzzification and defuzzification algorithm, imp

		'position',[20,25,40,45])

add_block('built-in/S-function',[sys,'/','Fuzzy Controller/2-D Table Lookup/S-function'])
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup/S-function'],...
		'function name','sftable2',...
		'parameters','xindex, yindex, table',...
		'position',[140,52,190,78])

add_block('built-in/Mux',[sys,'/','Fuzzy Controller/2-D Table Lookup/Mux'])
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup/Mux'],...
		'inputs','2',...
		'position',[90,46,120,79])
add_line([sys,'/','Fuzzy Controller/2-D Table Lookup'],[45,90;65,90;65,70;80,70])
add_line([sys,'/','Fuzzy Controller/2-D Table Lookup'],[195,65;205,65])
add_line([sys,'/','Fuzzy Controller/2-D Table Lookup'],[45,35;65,35;65,55;80,55])
add_line([sys,'/','Fuzzy Controller/2-D Table Lookup'],[125,65;130,65])
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'Mask Display','plot(-10,-10,110,110,[90,50,10],[90,40,30],[90,50,10],[50,26,20],[90,50,10],[22,13,10])',...
		'Mask Type','2-D Table Lookup')
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'Mask Dialogue','Two Dimensional Table Lookup\nThe first input corresponds to X Index and the second input corresponds to the Y Index|X Index|Y Index|Table')
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'Mask Translate','xindex=@1; yindex=@2; table=@3; sftab2chk(xindex,yindex,table);')
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'Mask Help','This block returns a linearly interpolated intersection from the table using the X index (which corresponds to the rows of the table) and the Y index (which corresponds to the columns of the table).  Extrapolation is used.')
set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'Mask Entries','rate\/error\/power\/')


%     Finished composite block 'Fuzzy Controller/2-D Table Lookup'.

set_param([sys,'/','Fuzzy Controller/2-D Table Lookup'],...
		'position',[300,104,345,146])

add_block('built-in/Note',[sys,'/','Fuzzy Controller/Fuzzy Controller'])
set_param([sys,'/','Fuzzy Controller/Fuzzy Controller'],...
		'position',[325,75,326,76])

add_block('built-in/Gain',[sys,'/','Fuzzy Controller/Gain2'])
set_param([sys,'/','Fuzzy Controller/Gain2'],...
		'position',[430,115,450,135])

add_block('built-in/Outport',[sys,'/',['Fuzzy Controller/Blower',13,'switch']])
set_param([sys,'/',['Fuzzy Controller/Blower',13,'switch']],...
		'position',[520,115,540,135])
add_line([sys,'/','Fuzzy Controller'],[455,125;510,125])
add_line([sys,'/','Fuzzy Controller'],[50,160;220,160])
add_line([sys,'/','Fuzzy Controller'],[50,160;50,95;65,95])
add_line([sys,'/','Fuzzy Controller'],[190,105;220,105])
add_line([sys,'/','Fuzzy Controller'],[130,95;155,95])
add_line([sys,'/','Fuzzy Controller'],[130,160;130,110;155,110])
add_line([sys,'/','Fuzzy Controller'],[255,105;290,115])
add_line([sys,'/','Fuzzy Controller'],[255,160;290,135])
add_line([sys,'/','Fuzzy Controller'],[350,125;420,125])
set_param([sys,'/','Fuzzy Controller'],...
		'Mask Display','plot(0,0,100,100,[20,80,80,20,20],[20,20,80,80,20],[70,65,50,35,30],[20,40,50,40,20],[45,40,43],[47,59,45],[30,30],[60,60],[40,40],[65,65],[50,50],[70,70],[60,60],[65,65],[70,70],[60,60])')


%     Finished composite block 'Fuzzy Controller'.

set_param([sys,'/','Fuzzy Controller'],...
		'Drop Shadow',4,...
		'position',[185,51,230,89])

add_block('built-in/Constant',[sys,'/','Set Point'])
set_param([sys,'/','Set Point'],...
		'Value','70',...
		'position',[25,50,45,70])

add_block('built-in/Scope',[sys,'/','Floating Scope'])
set_param([sys,'/','Floating Scope'],...
		'Vgain','90.000000',...
		'Hgain','18000.000000',...
		'Vmax','180.000000',...
		'Hmax','36000.000000',...
		'Window',[0,0,280,220],...
		'position',[480,142,500,168])

add_block('built-in/Note',[sys,'/',['blower',13,'cmd']])
set_param([sys,'/',['blower',13,'cmd']],...
		'position',[255,40,256,41])

add_block('built-in/Note',[sys,'/',['heater',13,'QDot']])
set_param([sys,'/',['heater',13,'QDot']],...
		'position',[310,40,311,41])

add_block('built-in/Sine Wave',[sys,'/',['Daily Temp',13,'Variation']])
set_param([sys,'/',['Daily Temp',13,'Variation']],...
		'amplitude','10',...
		'frequency','2*pi/(24*3600)',...
		'position',[165,170,185,190])


%     Subsystem  ['Temp',13,'Convert 2'].

new_system([sys,'/',['Temp',13,'Convert 2']])
set_param([sys,'/',['Temp',13,'Convert 2']],'Location',[577,197,817,328])

add_block('built-in/Gain',[sys,'/',['Temp',13,'Convert 2/Gain']])
set_param([sys,'/',['Temp',13,'Convert 2/Gain']],...
		'Gain','9/5',...
		'position',[60,7,95,43])

add_block('built-in/Outport',[sys,'/',['Temp',13,'Convert 2/Fahrenheit',13,'out']])
set_param([sys,'/',['Temp',13,'Convert 2/Fahrenheit',13,'out']],...
		'position',[180,25,200,45])

add_block('built-in/Inport',[sys,'/',['Temp',13,'Convert 2/Celsius',13,'in']])
set_param([sys,'/',['Temp',13,'Convert 2/Celsius',13,'in']],...
		'position',[15,15,35,35])

add_block('built-in/Constant',[sys,'/',['Temp',13,'Convert 2/Constant']])
set_param([sys,'/',['Temp',13,'Convert 2/Constant']],...
		'Value','32',...
		'position',[65,70,85,90])

add_block('built-in/Sum',[sys,'/',['Temp',13,'Convert 2/Sum']])
set_param([sys,'/',['Temp',13,'Convert 2/Sum']],...
		'position',[130,15,150,50])
add_line([sys,'/',['Temp',13,'Convert 2']],[40,25;50,25])
add_line([sys,'/',['Temp',13,'Convert 2']],[100,25;120,25])
add_line([sys,'/',['Temp',13,'Convert 2']],[155,35;170,35])
add_line([sys,'/',['Temp',13,'Convert 2']],[90,80;105,80;105,40;120,40])
set_param([sys,'/',['Temp',13,'Convert 2']],...
		'Mask Display','C2F')


%     Finished composite block ['Temp',13,'Convert 2'].

set_param([sys,'/',['Temp',13,'Convert 2']],...
		'position',[445,68,475,102])


%     Subsystem  ['Temp',13,'Convert 1'].

new_system([sys,'/',['Temp',13,'Convert 1']])
set_param([sys,'/',['Temp',13,'Convert 1']],'Location',[387,197,620,324])

add_block('built-in/Sum',[sys,'/',['Temp',13,'Convert 1/Sum']])
set_param([sys,'/',['Temp',13,'Convert 1/Sum']],...
		'inputs','+-',...
		'position',[85,15,105,50])

add_block('built-in/Constant',[sys,'/',['Temp',13,'Convert 1/Constant']])
set_param([sys,'/',['Temp',13,'Convert 1/Constant']],...
		'Value','32',...
		'position',[35,70,55,90])

add_block('built-in/Inport',[sys,'/',['Temp',13,'Convert 1/Fahrenheit',13,'in']])
set_param([sys,'/',['Temp',13,'Convert 1/Fahrenheit',13,'in']],...
		'position',[15,15,35,35])

add_block('built-in/Outport',[sys,'/',['Temp',13,'Convert 1/Celsius',13,'out']])
set_param([sys,'/',['Temp',13,'Convert 1/Celsius',13,'out']],...
		'position',[180,25,200,45])

add_block('built-in/Gain',[sys,'/',['Temp',13,'Convert 1/Gain1']])
set_param([sys,'/',['Temp',13,'Convert 1/Gain1']],...
		'Gain','5/9',...
		'position',[125,17,160,53])
add_line([sys,'/',['Temp',13,'Convert 1']],[60,80;70,80;70,40;75,40])
add_line([sys,'/',['Temp',13,'Convert 1']],[40,25;75,25])
add_line([sys,'/',['Temp',13,'Convert 1']],[110,35;115,35])
add_line([sys,'/',['Temp',13,'Convert 1']],[165,35;170,35])
set_param([sys,'/',['Temp',13,'Convert 1']],...
		'Mask Display','F2C')


%     Finished composite block ['Temp',13,'Convert 1'].

set_param([sys,'/',['Temp',13,'Convert 1']],...
		'position',[70,42,100,78])

add_block('built-in/Integrator',[sys,'/','Integrator'])
set_param([sys,'/','Integrator'],...
		'position',[360,15,380,35])

add_block('built-in/Scope',[sys,'/','$$'])
set_param([sys,'/','$$'],...
		'Vgain','3.000000',...
		'Hgain','18000.000000',...
		'Vmax','6.000000',...
		'Hmax','36000.000000',...
		'Window',[0,0,280,220],...
		'position',[475,12,495,38])

add_block('built-in/Note',[sys,'/',['House Thermodynamics',13,'(Double click on the "?" for more info)']])
set_param([sys,'/',['House Thermodynamics',13,'(Double click on the "?" for more info)']],...
		'position',[260,240,261,241])


%     Subsystem  'More Info'.

new_system([sys,'/','More Info'])
set_param([sys,'/','More Info'],'Location',[4,42,368,185])

add_block('built-in/Note',[sys,'/',['More Info/This system models the thermodynamics of a house.',13,'The variables are initialized by a file ',13,'called thermodat.m.']])
set_param([sys,'/',['More Info/This system models the thermodynamics of a house.',13,'The variables are initialized by a file ',13,'called thermodat.m.']],...
		'position',[177,25,178,26])

add_block('built-in/Note',[sys,'/',['More Info/Examine the file thermodat.m to see the ',13,'simulation variables you can change at will.',13,'']])
set_param([sys,'/',['More Info/Examine the file thermodat.m to see the ',13,'simulation variables you can change at will.',13,'']],...
		'position',[185,95,186,96])
set_param([sys,'/','More Info'],...
		'Mask Display','?')


%     Finished composite block 'More Info'.

set_param([sys,'/','More Info'],...
		'hide name',0,...
		'Drop Shadow',4,...
		'position',[420,240,447,265])

add_block('built-in/Scope',[sys,'/','Temperature1'])
set_param([sys,'/','Temperature1'],...
		'Vgain','120.000000',...
		'Hgain','18000.000000',...
		'Vmax','240.000000',...
		'Hmax','36000.000000',...
		'Window',[0,0,280,220],...
		'position',[525,12,545,38])


%     Subsystem  'Graph Scope'.

new_system([sys,'/','Graph Scope'])
set_param([sys,'/','Graph Scope'],'Location',[0,0,274,193])

add_block('built-in/S-function',[sys,'/',['Graph Scope/S-function',13,'M-file which plots',13,'lines',13,'']])
set_param([sys,'/',['Graph Scope/S-function',13,'M-file which plots',13,'lines',13,'']],...
		'function name','sfuny',...
		'parameters','ax, color',...
		'position',[130,55,180,75])

add_block('built-in/Inport',[sys,'/','Graph Scope/x'])
set_param([sys,'/','Graph Scope/x'],...
		'position',[65,55,85,75])
add_line([sys,'/','Graph Scope'],[90,65;120,65])
set_param([sys,'/','Graph Scope'],...
		'Mask Display','plot(0,0,100,100,[90,10,10,10,90,90,10],[65,65,90,40,40,90,90],[90,78,69,54,40,31,25,10],[77,60,48,46,56,75,81,84])',...
		'Mask Type','Graph scope.')
set_param([sys,'/','Graph Scope'],...
		'Mask Dialogue','Graph scope using MATLAB graph window.\nEnter plotting ranges and line type.|Time range:|y-min:|y-max:|Line type (rgbw-:*):',...
		'Mask Translate','color = @4; ax = [0, @1, @2, @3];')
set_param([sys,'/','Graph Scope'],...
		'Mask Help','This block plots to the MATLAB graph window and can be used as an improved version of the Scope block. Look at the m-file sfuny.m to see how it works.',...
		'Mask Entries','600\/65\/75\/''y-''\/')


%     Finished composite block 'Graph Scope'.

set_param([sys,'/','Graph Scope'],...
		'Drop Shadow',4,...
		'position',[505,66,535,104])

add_block('built-in/Note',[sys,'/',['Execute th_dat from workspace,',13,'before starting simulation.']])
set_param([sys,'/',['Execute th_dat from workspace,',13,'before starting simulation.']],...
		'position',[490,290,491,291])
add_line(sys,[445,25;465,25])
add_line(sys,[385,25;410,25])
add_line(sys,[315,140;325,140;325,95;340,95])
add_line(sys,[255,140;270,140])
add_line(sys,[105,60;125,60])
add_line(sys,[50,60;60,60])
add_line(sys,[405,85;435,85])
add_line(sys,[415,85;415,225;115,225;115,75;125,75])
add_line(sys,[190,180;210,180;210,145;220,145])
add_line(sys,[190,130;220,130])
add_line(sys,[160,70;175,70])
add_line(sys,[235,70;265,70])
add_line(sys,[480,85;495,85])
add_line(sys,[300,70;340,70])
add_line(sys,[330,70;330,25;350,25])

% Return any arguments.
if (nargin | nargout)
	% Must use feval here to access system in memory
	if (nargin > 3)
		if (flag == 0)
			eval(['[ret,x0,str]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str] = feval(sys);
	end
end