sample1a.m

模糊神经网络采用matlab编程 o install NEFCON follow these steps: 1. Unpack the tar file NEFCON.TAR into your MA

set_param([sys,'/',['tank 2/flow',13,'in']],...
		'position',[25,175,45,195])

add_block('built-in/Relay',[sys,'/','tank 2/Overflow sensor'])
set_param([sys,'/','tank 2/Overflow sensor'],...
		'On_switch_value','ht - overflowLimit',...
		'Off_switch_value','ht - overflowLimit',...
		'position',[385,18,425,62])

add_block('built-in/Outport',[sys,'/',['tank 2/flow',13,'out']])
set_param([sys,'/',['tank 2/flow',13,'out']],...
		'Port','2',...
		'position',[495,185,515,205])
add_line([sys,'/','tank 2'],[130,195;145,195])
add_line([sys,'/','tank 2'],[50,185;100,185])
add_line([sys,'/','tank 2'],[430,40;490,40])
add_line([sys,'/','tank 2'],[205,195;215,195])
add_line([sys,'/','tank 2'],[275,195;295,195])
add_line([sys,'/','tank 2'],[395,195;405,195])
add_line([sys,'/','tank 2'],[465,195;490,195])
add_line([sys,'/','tank 2'],[475,195;475,260;85,260;85,200;100,200])
add_line([sys,'/','tank 2'],[275,195;285,195;285,40;380,40])
add_line([sys,'/','tank 2'],[285,110;490,110])
set_param([sys,'/','tank 2'],...
		'Mask Display','WATER\nTANK',...
		'Mask Type','Water tank')
set_param([sys,'/','tank 2'],...
		'Mask Dialogue','Cylindrical water tank with overflow signal|Height of tank (m):|Bottom area (m^2):|Out pipe crossection (m^2):|Overflow sensor distance from top (m):|Initial level height (m):')
set_param([sys,'/','tank 2'],...
		'Mask Translate','ht=@1;area=@2;outletArea=@3;overflowLimit=@4;initialHt=@5;',...
		'Mask Entries','2\/1\/0.05\/0\/0.5\/')


%     Finished composite block 'tank 2'.

set_param([sys,'/','tank 2'],...
		'position',[380,106,440,164])


%     Subsystem  'Subsystem'.

new_system([sys,'/','Subsystem'])
set_param([sys,'/','Subsystem'],'Location',[408,498,736,676])

add_block('built-in/Outport',[sys,'/','Subsystem/out flow'])
set_param([sys,'/','Subsystem/out flow'],...
		'position',[250,95,270,115])

add_block('built-in/Product',[sys,'/','Subsystem/Product'])
set_param([sys,'/','Subsystem/Product'],...
		'position',[180,76,215,134])


%     Subsystem  'Subsystem/Limited Integrator'.

new_system([sys,'/','Subsystem/Limited Integrator'])
set_param([sys,'/','Subsystem/Limited Integrator'],'Location',[0,0,791,233])

add_block('built-in/Saturation',[sys,'/','Subsystem/Limited Integrator/Saturation'])
set_param([sys,'/','Subsystem/Limited Integrator/Saturation'],...
		'Lower Limit','lb',...
		'Upper Limit','ub',...
		'position',[630,121,675,149])

add_block('built-in/Fcn',[sys,'/','Subsystem/Limited Integrator/Fcn'])
set_param([sys,'/','Subsystem/Limited Integrator/Fcn'],...
		'Expr','u[2]*((((u[1]>lb)+(u[2]>=0))>0) * (((u[1]<ub)+(u[2]<=0))>0))',...
		'position',[150,116,520,154])

add_block('built-in/Inport',[sys,'/','Subsystem/Limited Integrator/In_1'])
set_param([sys,'/','Subsystem/Limited Integrator/In_1'],...
		'position',[20,130,40,150])

add_block('built-in/Outport',[sys,'/','Subsystem/Limited Integrator/Out_1'])
set_param([sys,'/','Subsystem/Limited Integrator/Out_1'],...
		'position',[720,125,740,145])

add_block('built-in/Integrator',[sys,'/','Subsystem/Limited Integrator/Integrator'])
set_param([sys,'/','Subsystem/Limited Integrator/Integrator'],...
		'Initial','xi',...
		'position',[540,126,585,144])

add_block('built-in/Mux',[sys,'/','Subsystem/Limited Integrator/Mux'])
set_param([sys,'/','Subsystem/Limited Integrator/Mux'],...
		'inputs','2',...
		'position',[85,116,115,149])
add_line([sys,'/','Subsystem/Limited Integrator'],[680,135;715,135])
add_line([sys,'/','Subsystem/Limited Integrator'],[45,140;80,140])
add_line([sys,'/','Subsystem/Limited Integrator'],[120,135;145,135])
add_line([sys,'/','Subsystem/Limited Integrator'],[525,135;535,135])
add_line([sys,'/','Subsystem/Limited Integrator'],[590,135;625,135])
add_line([sys,'/','Subsystem/Limited Integrator'],[600,135;600,35;55,35;55,125;80,125])
set_param([sys,'/','Subsystem/Limited Integrator'],...
		'Mask Display','plot(-1,-0.2,3.5,1.2,[0.05,1,2,2.9],[0,0,1,1]); 1/s    ',...
		'Mask Type','Limited integrator.')
set_param([sys,'/','Subsystem/Limited Integrator'],...
		'Mask Dialogue','Limited integrator|Lower bound:|Upper bound:|Initial condition',...
		'Mask Translate','lb = @1; ub = @2; xi = @3;')
set_param([sys,'/','Subsystem/Limited Integrator'],...
		'Mask Help','Implements a limited integrator with \nan expression of the form:\n\nif(x<=lb and u<0) or (x>=ub and u>0)\n\txdot=0\nelse\n\txdot=u.')
set_param([sys,'/','Subsystem/Limited Integrator'],...
		'Mask Entries','0\/1\/InitialPosition\/')


%     Finished composite block 'Subsystem/Limited Integrator'.

set_param([sys,'/','Subsystem/Limited Integrator'],...
		'position',[80,9,130,51])

add_block('built-in/Inport',[sys,'/','Subsystem/source flow'])
set_param([sys,'/','Subsystem/source flow'],...
		'Port','2',...
		'position',[20,110,40,130])

add_block('built-in/Inport',[sys,'/','Subsystem/control'])
set_param([sys,'/','Subsystem/control'],...
		'position',[20,20,40,40])
add_line([sys,'/','Subsystem'],[135,30;160,30;160,90;175,90])
add_line([sys,'/','Subsystem'],[45,30;75,30])
add_line([sys,'/','Subsystem'],[220,105;245,105])
add_line([sys,'/','Subsystem'],[45,120;175,120])
set_param([sys,'/','Subsystem'],...
		'Mask Display','VALVE',...
		'Mask Type','Valve',...
		'Mask Dialogue','Valve|Initial position [0=closed 1=open]:',...
		'Mask Translate','InitialPosition = @1;',...
		'Mask Entries','0\/')


%     Finished composite block 'Subsystem'.

set_param([sys,'/','Subsystem'],...
		'position',[315,113,365,157])

add_block('built-in/Constant',[sys,'/',['tank max',13,'inflow']])
set_param([sys,'/',['tank max',13,'inflow']],...
		'orientation',2,...
		'Value','0.5',...
		'position',[320,194,365,226])


%     Subsystem  ['PID ',13,'Controller'].

new_system([sys,'/',['PID ',13,'Controller']])
set_param([sys,'/',['PID ',13,'Controller']],'Location',[8,55,370,299])

add_block('built-in/Sum',[sys,'/',['PID ',13,'Controller/Sum']])
set_param([sys,'/',['PID ',13,'Controller/Sum']],...
		'inputs','+++',...
		'position',[245,57,265,93])

add_block('built-in/Gain',[sys,'/',['PID ',13,'Controller/D']])
set_param([sys,'/',['PID ',13,'Controller/D']],...
		'Gain','D',...
		'position',[95,129,115,151])

add_block('built-in/Gain',[sys,'/',['PID ',13,'Controller/Proportional']])
set_param([sys,'/',['PID ',13,'Controller/Proportional']],...
		'Gain','P',...
		'position',[120,13,140,37])

add_block('built-in/Transfer Fcn',[sys,'/',['PID ',13,'Controller/Integral']])
set_param([sys,'/',['PID ',13,'Controller/Integral']],...
		'Numerator','[I]',...
		'Denominator','[1 0]',...
		'position',[110,57,145,93])

add_block('built-in/Derivative',[sys,'/',['PID ',13,'Controller/Derivative']])
set_param([sys,'/',['PID ',13,'Controller/Derivative']],...
		'position',[150,128,190,152])

add_block('built-in/Outport',[sys,'/',['PID ',13,'Controller/Out_1']])
set_param([sys,'/',['PID ',13,'Controller/Out_1']],...
		'position',[290,65,310,85])

add_block('built-in/Inport',[sys,'/',['PID ',13,'Controller/In_1']])
set_param([sys,'/',['PID ',13,'Controller/In_1']],...
		'position',[25,65,45,85])
add_line([sys,'/',['PID ',13,'Controller']],[120,140;145,140])
add_line([sys,'/',['PID ',13,'Controller']],[270,75;285,75])
add_line([sys,'/',['PID ',13,'Controller']],[50,75;105,75])
add_line([sys,'/',['PID ',13,'Controller']],[65,75;65,140;90,140])
add_line([sys,'/',['PID ',13,'Controller']],[80,75;80,25;115,25])
add_line([sys,'/',['PID ',13,'Controller']],[195,140;215,140;215,85;240,85])
add_line([sys,'/',['PID ',13,'Controller']],[150,75;240,75])
add_line([sys,'/',['PID ',13,'Controller']],[145,25;210,25;210,65;240,65])
set_param([sys,'/',['PID ',13,'Controller']],...
		'Mask Display','PID',...
		'Mask Type','PID Controller',...
		'Mask Dialogue','Enter expressions for proportional, integral, and derivative terms.\nP+I/s+Ds|Proportional:|Integral|Derivative:')
set_param([sys,'/',['PID ',13,'Controller']],...
		'Mask Translate','P=@1; I=@2; D=@3;')
set_param([sys,'/',['PID ',13,'Controller']],...
		'Mask Help','This block implements a PID controller where parameters are entered for the Proportional, Integral and Derivative terms. Unmask this block to see how it works. The derivative term is implemented using a true derivative block.')
set_param([sys,'/',['PID ',13,'Controller']],...
		'Mask Entries','2\/0\/1\/')


%     Finished composite block ['PID ',13,'Controller'].

set_param([sys,'/',['PID ',13,'Controller']],...
		'position',[210,109,250,141])

add_block('built-in/To Workspace',[sys,'/','To Workspace7'])
set_param([sys,'/','To Workspace7'],...
		'hide name',0,...
		'mat-name','u',...
		'buffer','3000',...
		'position',[120,77,170,93])

add_block('built-in/To Workspace',[sys,'/','To Workspace8'])
set_param([sys,'/','To Workspace8'],...
		'hide name',0,...
		'mat-name','y',...
		'buffer','3000',...
		'position',[525,277,575,293])

add_block('built-in/Scope',[sys,'/','Scope1'])
set_param([sys,'/','Scope1'],...
		'Vgain','2.000000',...
		'Hgain','200.000000',...
		'Vmax','4.000000',...
		'Hmax','400.000000',...
		'Window',[732,445,1058,728])
open_system([sys,'/','Scope1'])
set_param([sys,'/','Scope1'],...
		'position',[600,165,630,195])

add_block('built-in/S-Function',[sys,'/','S-Function'])
set_param([sys,'/','S-Function'],...
		'function name','animtank',...
		'parameters','1',...
		'position',[580,86,635,134])

add_block('built-in/Mux',[sys,'/','Mux2'])
set_param([sys,'/','Mux2'],...
		'inputs','2',...
		'position',[525,91,555,124])
add_line(sys,[150,125;205,125])
add_line(sys,[560,110;565,110;565,180;595,180])
add_line(sys,[445,155;435,200])
add_line(sys,[560,110;575,110])
add_line(sys,[370,135;375,135])
add_line(sys,[315,210;300,210;310,145])
add_line(sys,[445,115;485,115;485,305;100,305;100,130;120,130])
add_line(sys,[95,120;120,120])
add_line(sys,[105,120;105,40;485,40;485,100;520,100])
add_line(sys,[105,40;115,85])
add_line(sys,[255,125;310,125])
add_line(sys,[485,285;510,285;520,285])
add_line(sys,[485,115;510,115;520,115])

drawnow

% 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,ts,xts]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str,ts,xts] = feval(sys);
	end
else
	drawnow % Flash up the model and execute load callback
end