📄 sample1b.m
字号:
set_param([sys,'/','Scope4'],...
'Vgain','2.000000',...
'Hgain','70.000000',...
'Vmax','4.000000',...
'Hmax','140.000000',...
'Window',[410,361,736,593],...
'position',[190,35,220,65])
add_block('built-in/Scope',[sys,'/','Scope3'])
set_param([sys,'/','Scope3'],...
'Vgain','2.000000',...
'Hgain','70.000000',...
'Vmax','4.000000',...
'Hmax','140.000000',...
'Window',[606,290,910,555])
open_system([sys,'/','Scope3'])
set_param([sys,'/','Scope3'],...
'position',[305,35,335,65])
add_block('built-in/Constant',[sys,'/',['tank max',13,'inflow']])
set_param([sys,'/',['tank max',13,'inflow']],...
'orientation',2,...
'Value','0.5',...
'position',[330,194,375,226])
add_block('built-in/Scope',[sys,'/','Scope2'])
set_param([sys,'/','Scope2'],...
'orientation',2,...
'Vgain','2.000000',...
'Hgain','20.000000',...
'Vmax','4.000000',...
'Hmax','40.000000',...
'Window',[207,553,469,687],...
'position',[410,195,440,225])
% Subsystem 'tank 2'.
new_system([sys,'/','tank 2'])
set_param([sys,'/','tank 2'],'Location',[109,460,655,768])
add_block('built-in/Gain',[sys,'/','tank 2/1//Height'])
set_param([sys,'/','tank 2/1//Height'],...
'Gain','1/ht',...
'position',[185,77,235,133])
add_block('built-in/Outport',[sys,'/',['tank 2/overflow',13,'flag']])
set_param([sys,'/',['tank 2/overflow',13,'flag']],...
'Port','3',...
'position',[495,30,515,50])
add_block('built-in/Outport',[sys,'/',['tank 2/water',13,'level']])
set_param([sys,'/',['tank 2/water',13,'level']],...
'position',[495,100,515,120])
add_block('built-in/Gain',[sys,'/','tank 2/1//area'])
set_param([sys,'/','tank 2/1//area'],...
'Gain','1/area',...
'position',[220,167,270,223])
add_block('built-in/Gain',[sys,'/','tank 2/outletArea'])
set_param([sys,'/','tank 2/outletArea'],...
'Gain','outletArea',...
'position',[410,167,460,223])
add_block('built-in/Fcn',[sys,'/','tank 2/sqrt(2gh)'])
set_param([sys,'/','tank 2/sqrt(2gh)'],...
'Expr','sqrt(2*9.8*u)',...
'position',[300,171,390,219])
% Subsystem ['tank 2/tank',13,'volume'].
new_system([sys,'/',['tank 2/tank',13,'volume']])
set_param([sys,'/',['tank 2/tank',13,'volume']],'Location',[0,0,791,233])
add_block('built-in/Saturation',[sys,'/',['tank 2/tank',13,'volume/Saturation']])
set_param([sys,'/',['tank 2/tank',13,'volume/Saturation']],...
'Lower Limit','lb',...
'Upper Limit','ub',...
'position',[630,121,675,149])
add_block('built-in/Fcn',[sys,'/',['tank 2/tank',13,'volume/Fcn']])
set_param([sys,'/',['tank 2/tank',13,'volume/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,'/',['tank 2/tank',13,'volume/In_1']])
set_param([sys,'/',['tank 2/tank',13,'volume/In_1']],...
'position',[20,130,40,150])
add_block('built-in/Outport',[sys,'/',['tank 2/tank',13,'volume/Out_1']])
set_param([sys,'/',['tank 2/tank',13,'volume/Out_1']],...
'position',[720,125,740,145])
add_block('built-in/Integrator',[sys,'/',['tank 2/tank',13,'volume/Integrator']])
set_param([sys,'/',['tank 2/tank',13,'volume/Integrator']],...
'Initial','xi',...
'position',[540,126,585,144])
add_block('built-in/Mux',[sys,'/',['tank 2/tank',13,'volume/Mux']])
set_param([sys,'/',['tank 2/tank',13,'volume/Mux']],...
'inputs','2',...
'position',[85,116,115,149])
add_line([sys,'/',['tank 2/tank',13,'volume']],[680,135;715,135])
add_line([sys,'/',['tank 2/tank',13,'volume']],[45,140;80,140])
add_line([sys,'/',['tank 2/tank',13,'volume']],[120,135;145,135])
add_line([sys,'/',['tank 2/tank',13,'volume']],[525,135;535,135])
add_line([sys,'/',['tank 2/tank',13,'volume']],[590,135;625,135])
add_line([sys,'/',['tank 2/tank',13,'volume']],[600,135;600,35;55,35;55,125;80,125])
set_param([sys,'/',['tank 2/tank',13,'volume']],...
'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,'/',['tank 2/tank',13,'volume']],...
'Mask Dialogue','Limited integrator|Lower bound:|Upper bound:|Initial condition',...
'Mask Translate','lb = @1; ub = @2; xi = @3;')
set_param([sys,'/',['tank 2/tank',13,'volume']],...
'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,'/',['tank 2/tank',13,'volume']],...
'Mask Entries','0\/ht * area\/initialHt * area\/')
% Finished composite block ['tank 2/tank',13,'volume'].
set_param([sys,'/',['tank 2/tank',13,'volume']],...
'position',[150,174,200,216])
add_block('built-in/Sum',[sys,'/','tank 2/sum'])
set_param([sys,'/','tank 2/sum'],...
'inputs','+-',...
'position',[105,175,125,210])
add_block('built-in/Inport',[sys,'/',['tank 2/flow',13,'in']])
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',[390,116,450,174])
% 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',[325,123,375,167])
add_block('built-in/Derivative',[sys,'/','Derivative'])
set_param([sys,'/','Derivative'],...
'position',[160,140,195,170])
add_block('built-in/Sum',[sys,'/','error'])
set_param([sys,'/','error'],...
'inputs','+-',...
'position',[110,100,130,120])
add_line(sys,[520,120;525,120])
add_line(sys,[455,165;445,210])
add_line(sys,[380,145;385,145])
add_line(sys,[135,110;215,110])
add_line(sys,[325,210;310,210;320,155])
add_line(sys,[455,125;480,125])
add_line(sys,[200,155;215,155])
add_line(sys,[85,105;105,105])
add_line(sys,[95,105;85,105;85,20;465,20;465,110;480,110])
add_line(sys,[285,135;320,135])
add_line(sys,[455,125;465,125;465,260;85,260;85,115;105,115])
add_line(sys,[310,135;290,135;300,50])
add_line(sys,[135,110;175,110;185,50])
add_line(sys,[85,50;95,50])
add_line(sys,[465,240;485,240])
add_line(sys,[145,110;155,155])
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
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -