s4.m

Gives all the matlab codes for dynamic simulation of electric machinery by Chee-Mun Ong

add_block('built-in/Constant',[sys,'/','PSS_sw'])
set_param([sys,'/','PSS_sw'],...
		'orientation',2,...
		'Value','pss_sw',...
		'position',[185,385,260,405])

add_block('built-in/Switch',[sys,'/','Switch'])
set_param([sys,'/','Switch'],...
		'orientation',2,...
		'position',[115,329,145,361])

add_block('built-in/Saturation',[sys,'/',['signal',13,'limits']])
set_param([sys,'/',['signal',13,'limits']],...
		'orientation',2,...
		'Lower Limit','-pss_limit',...
		'Upper Limit','pss_limit',...
		'position',[200,316,235,354])

add_block('built-in/Transfer Fcn',[sys,'/','wash-out'])
set_param([sys,'/','wash-out'],...
		'orientation',2,...
		'Numerator','[Ks 0]',...
		'Denominator','[Tw 1]',...
		'position',[260,315,340,355])

add_block('built-in/Transfer Fcn',[sys,'/','comp2'])
set_param([sys,'/','comp2'],...
		'orientation',2,...
		'Numerator','[T3 1]',...
		'Denominator','[T4 1]',...
		'position',[375,315,420,355])

add_block('built-in/Transfer Fcn',[sys,'/','comp1'])
set_param([sys,'/','comp1'],...
		'orientation',2,...
		'Numerator','[T1 1]',...
		'Denominator','[T2 1]',...
		'position',[460,316,510,354])

add_block('built-in/Mux',[sys,'/','Mux'])
set_param([sys,'/','Mux'],...
		'orientation',3,...
		'inputs','5',...
		'position',[337,45,483,70])

add_block('built-in/To Workspace',[sys,'/','To Workspace'])
set_param([sys,'/','To Workspace'],...
		'mat-name','y',...
		'buffer','10000',...
		'position',[540,27,590,43])

add_block('built-in/Clock',[sys,'/','Clock'])
set_param([sys,'/','Clock'],...
		'orientation',2,...
		'position',[555,80,575,100])

add_block('built-in/Sum',[sys,'/','Sum'])
set_param([sys,'/','Sum'],...
		'position',[190,239,210,266])


%     Subsystem  'exciter'.

new_system([sys,'/','exciter'])
set_param([sys,'/','exciter'],'Location',[149,293,965,637])

add_block('built-in/Integrator',[sys,'/','exciter/VR'])
set_param([sys,'/','exciter/VR'],...
		'Initial','VR',...
		'position',[325,89,355,121])

add_block('built-in/Sum',[sys,'/','exciter/Sum1'])
set_param([sys,'/','exciter/Sum1'],...
		'inputs','-+',...
		'position',[460,68,475,117])

add_block('built-in/Saturation',[sys,'/','exciter/VRmax//VRmin'])
set_param([sys,'/','exciter/VRmax//VRmin'],...
		'Lower Limit','VRmin',...
		'Upper Limit','VRmax',...
		'position',[385,82,435,128])

add_block('built-in/Gain',[sys,'/','exciter/Gain2'])
set_param([sys,'/','exciter/Gain2'],...
		'Gain','1/TE',...
		'position',[500,77,565,113])

add_block('built-in/Gain',[sys,'/','exciter/Gain3'])
set_param([sys,'/','exciter/Gain3'],...
		'orientation',2,...
		'Gain','KE/TE',...
		'position',[575,157,685,193])

add_block('built-in/Fcn',[sys,'/','exciter/Se'])
set_param([sys,'/','exciter/Se'],...
		'orientation',2,...
		'Expr','AEx*exp(BEx*u[1])',...
		'position',[480,21,660,49])

add_block('built-in/Outport',[sys,'/','exciter/out_Ef'])
set_param([sys,'/','exciter/out_Ef'],...
		'position',[745,95,765,115])

add_block('built-in/Gain',[sys,'/','exciter/Gain'])
set_param([sys,'/','exciter/Gain'],...
		'Gain','KA/TA',...
		'position',[170,74,265,116])

add_block('built-in/Gain',[sys,'/','exciter/Gain1'])
set_param([sys,'/','exciter/Gain1'],...
		'orientation',2,...
		'Gain','1/TA',...
		'position',[265,158,350,192])

add_block('built-in/Gain',[sys,'/','exciter/Gain5'])
set_param([sys,'/','exciter/Gain5'],...
		'Gain','1/TF',...
		'position',[230,274,295,316])

add_block('built-in/Gain',[sys,'/','exciter/Gain4'])
set_param([sys,'/','exciter/Gain4'],...
		'orientation',2,...
		'Gain','KF/TF',...
		'position',[405,216,520,254])

add_block('built-in/Integrator',[sys,'/','exciter/Ef'])
set_param([sys,'/','exciter/Ef'],...
		'Initial','Efo',...
		'position',[645,90,675,120])

add_block('built-in/Integrator',[sys,'/','exciter/Vs'])
set_param([sys,'/','exciter/Vs'],...
		'Initial','Vs',...
		'position',[345,277,380,313])

add_block('built-in/Note',[sys,'/','exciter/Vstab'])
set_param([sys,'/','exciter/Vstab'],...
		'position',[150,225,155,230])

add_block('built-in/Sum',[sys,'/','exciter/Sum2'])
set_param([sys,'/','exciter/Sum2'],...
		'inputs','+-',...
		'position',[290,84,305,126])

add_block('built-in/Sum',[sys,'/','exciter/Sum3'])
set_param([sys,'/','exciter/Sum3'],...
		'inputs','+-',...
		'position',[600,84,615,126])

add_block('built-in/Sum',[sys,'/','exciter/Sum4'])
set_param([sys,'/','exciter/Sum4'],...
		'orientation',2,...
		'inputs','+-',...
		'position',[305,226,320,264])

add_block('built-in/Sum',[sys,'/','exciter/Sum'])
set_param([sys,'/','exciter/Sum'],...
		'inputs','+---',...
		'position',[135,58,150,127])

add_block('built-in/Inport',[sys,'/','exciter/Vref'])
set_param([sys,'/','exciter/Vref'],...
		'position',[55,45,75,65])

add_block('built-in/Inport',[sys,'/','exciter/Vsupp'])
set_param([sys,'/','exciter/Vsupp'],...
		'Port','2',...
		'position',[55,85,75,105])

add_block('built-in/Inport',[sys,'/','exciter/Vt'])
set_param([sys,'/','exciter/Vt'],...
		'Port','3',...
		'position',[55,130,75,150])

add_block('built-in/Scope',[sys,'/','exciter/Scope'])
set_param([sys,'/','exciter/Scope'],...
		'orientation',2,...
		'Vgain','1.200000',...
		'Hgain','10.000000',...
		'Vmax','2.400000',...
		'Hmax','20.000000',...
		'Window',[177,801,606,971],...
		'position',[315,36,335,64])
add_line([sys,'/','exciter'],[440,105;455,105])
add_line([sys,'/','exciter'],[80,55;100,55;100,70;130,70])
add_line([sys,'/','exciter'],[475,35;440,35;440,80;455,80])
add_line([sys,'/','exciter'],[310,105;320,105])
add_line([sys,'/','exciter'],[155,95;165,95])
add_line([sys,'/','exciter'],[360,105;380,105])
add_line([sys,'/','exciter'],[360,105;365,105;355,175])
add_line([sys,'/','exciter'],[620,105;640,105])
add_line([sys,'/','exciter'],[680,105;710,105;710,175;690,175])
add_line([sys,'/','exciter'],[480,95;495,95])
add_line([sys,'/','exciter'],[570,95;595,95])
add_line([sys,'/','exciter'],[570,175;565,175;565,115;595,115])
add_line([sys,'/','exciter'],[680,105;710,105;710,35;665,35])
add_line([sys,'/','exciter'],[680,105;740,105])
add_line([sys,'/','exciter'],[270,95;285,95])
add_line([sys,'/','exciter'],[260,175;250,175;250,115;285,115])
add_line([sys,'/','exciter'],[400,235;325,235])
add_line([sys,'/','exciter'],[300,295;340,295])
add_line([sys,'/','exciter'],[385,295;405,295;405,255;325,255])
add_line([sys,'/','exciter'],[680,105;710,105;710,235;525,235])
add_line([sys,'/','exciter'],[300,245;205,245;205,295;225,295])
add_line([sys,'/','exciter'],[300,245;110,245;110,115;130,115])
add_line([sys,'/','exciter'],[80,140;100,140;100,100;130,100])
add_line([sys,'/','exciter'],[80,95;100,95;100,85;130,85])


%     Finished composite block 'exciter'.

set_param([sys,'/','exciter'],...
		'position',[135,84,175,156])

add_block('built-in/Constant',[sys,'/','Vref*'])
set_param([sys,'/','Vref*'],...
		'Value','Vref',...
		'position',[45,86,75,104])

add_block('built-in/Note',[sys,'/','Vref'])
set_param([sys,'/','Vref'],...
		'position',[100,75,105,80])


%     Subsystem  'DTmech'.

new_system([sys,'/','DTmech'])
set_param([sys,'/','DTmech'],'Location',[5,40,315,196])

add_block('built-in/Outport',[sys,'/','DTmech/out_1'])
set_param([sys,'/','DTmech/out_1'],...
		'position',[285,60,305,80])

add_block('built-in/Look Up Table',[sys,'/','DTmech/Look-Up Table'])
set_param([sys,'/','DTmech/Look-Up Table'],...
		'Input_Values','t',...
		'Output_Values','y',...
		'position',[210,58,250,82])

add_block('built-in/Fcn',[sys,'/','DTmech/Fcn1'])
set_param([sys,'/','DTmech/Fcn1'],...
		'Expr','rem(u[1],period)',...
		'position',[110,60,150,80])

add_block('built-in/Clock',[sys,'/','DTmech/Clock'])
set_param([sys,'/','DTmech/Clock'],...
		'position',[45,60,65,80])
add_line([sys,'/','DTmech'],[255,70;280,70])
add_line([sys,'/','DTmech'],[70,70;105,70])
add_line([sys,'/','DTmech'],[155,70;205,70])
set_param([sys,'/','DTmech'],...
		'Mask Display','plot([t,t+period,t+2*period],[y,y,y])',...
		'Mask Type','Repeating table',...
		'Mask Dialogue','Repeating table.\nEnter values of time and output for first cycle.|Time values:|Output values:')
set_param([sys,'/','DTmech'],...
		'Mask Translate','period = max(@1); t = @1; y = @2;',...
		'Mask Help','Repeats cycle given in table. Time values should be monotonically increasing.',...
		'Mask Entries','time_Dtmech\/tmech_Dtmech\/')


%     Finished composite block 'DTmech'.

set_param([sys,'/','DTmech'],...
		'position',[90,280,145,300])


%     Subsystem  'm4'.

new_system([sys,'/','m4'])
set_param([sys,'/','m4'],'Location',[140,4325536,228,4325622])

add_block('built-in/Note',[sys,'/',['m4/Masked block of m4.m to initialize',13,'s4.m in chapter 10 and to plot results']])
set_param([sys,'/',['m4/Masked block of m4.m to initialize',13,'s4.m in chapter 10 and to plot results']],...
		'position',[55,50,60,55])
set_param([sys,'/','m4'],...
		'Mask Display','Initialize\nand plot',...
		'Mask Type','Masked block of m4.m',...
		'Mask Dialogue','eval(''m4'')',...
		'Mask Help','Uses m4.m to initialize and plot')


%     Finished composite block 'm4'.

set_param([sys,'/','m4'],...
		'Drop Shadow',4,...
		'position',[149,20,220,53])
add_line(sys,[220,165;225,165])
add_line(sys,[220,210;225,210])
add_line(sys,[410,40;410,35;310,35])
add_line(sys,[315,260;495,260])
add_line(sys,[315,235;550,235])
add_line(sys,[315,210;495,210])
add_line(sys,[315,185;550,185])
add_line(sys,[315,160;495,160])
add_line(sys,[315,135;550,135])
add_line(sys,[315,110;495,110])
add_line(sys,[410,40;410,35;535,35])
add_line(sys,[315,110;350,110;350,75])
add_line(sys,[180,120;225,120])
add_line(sys,[315,110;325,110;325,75;35,75;35,145;130,145])
add_line(sys,[370,335;345,335])
add_line(sys,[455,335;425,335])
add_line(sys,[315,235;535,235;535,335;515,335])
add_line(sys,[255,335;240,335])
add_line(sys,[195,335;150,335])
add_line(sys,[110,345;80,345;80,120;130,120])
add_line(sys,[180,395;160,395;150,345])
add_line(sys,[80,95;130,95])
add_line(sys,[155,245;185,245])
add_line(sys,[215,255;225,255])
add_line(sys,[315,160;380,160;380,75])
add_line(sys,[315,185;410,185;410,75])
add_line(sys,[315,210;440,210;440,75])
add_line(sys,[550,90;470,90;470,75])
add_line(sys,[150,290;160,290;160,260;185,260])

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