s3.m

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

add_line([sys,'/','gen'],[495,410;620,410])
add_line([sys,'/','gen'],[530,40;575,40])
add_line([sys,'/','gen'],[530,90;575,90])
add_line([sys,'/','gen'],[530,115;635,115])
add_line([sys,'/','gen'],[530,65;635,65])
add_line([sys,'/','gen'],[180,250;205,250;205,215;255,215])
add_line([sys,'/','gen'],[265,135;240,135;240,250;255,250])
add_line([sys,'/','gen'],[545,290;620,290])
add_line([sys,'/','gen'],[545,315;690,315])


%     Finished composite block 'gen'.

set_param([sys,'/','gen'],...
		'position',[145,46,215,234])

add_block('built-in/Constant',[sys,'/','Tmech'])
set_param([sys,'/','Tmech'],...
		'Value','Tmech',...
		'position',[45,200,100,220])


%     Subsystem  'Vq4'.

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

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

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

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

add_block('built-in/Clock',[sys,'/','Vq4/Clock'])
set_param([sys,'/','Vq4/Clock'],...
		'position',[45,60,65,80])
add_line([sys,'/','Vq4'],[255,70;280,70])
add_line([sys,'/','Vq4'],[70,70;105,70])
add_line([sys,'/','Vq4'],[155,70;205,70])
set_param([sys,'/','Vq4'],...
		'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,'/','Vq4'],...
		'Mask Translate','period = max(@1); t = @1; y = @2;',...
		'Mask Help','Repeats cycle given in table. Time values should be monotonically increasing.')
set_param([sys,'/','Vq4'],...
		'Mask Entries','[0 0.2 0.2 0.275 0.275 0.7 ]\/[1. 1. 0.5 0.5 1. 1.]*Vq4\/')


%     Finished composite block 'Vq4'.

set_param([sys,'/','Vq4'],...
		'position',[615,145,650,185])

add_block('built-in/Constant',[sys,'/','Vd4'])
set_param([sys,'/','Vd4'],...
		'Value','Vd4',...
		'position',[640,259,675,281])

add_block('built-in/Sum',[sys,'/','Sum4'])
set_param([sys,'/','Sum4'],...
		'position',[715,180,730,215])

add_block('built-in/Sum',[sys,'/','Sum5'])
set_param([sys,'/','Sum5'],...
		'position',[710,242,725,273])


%     Subsystem  'Cap'.

new_system([sys,'/','Cap'])
set_param([sys,'/','Cap'],'Location',[145,152,806,492])

add_block('built-in/Mux',[sys,'/','Cap/Mux'])
set_param([sys,'/','Cap/Mux'],...
		'position',[110,60,130,260])

add_block('built-in/Integrator',[sys,'/','Cap/vq'])
set_param([sys,'/','Cap/vq'],...
		'Initial','VqCs',...
		'position',[500,92,525,118])

add_block('built-in/Fcn',[sys,'/','Cap/Fcn'])
set_param([sys,'/','Cap/Fcn'],...
		'Expr','-wb*xc(1)*u[2] - wb*u[4] ',...
		'position',[190,88,470,122])

add_block('built-in/Outport',[sys,'/','Cap/out_vqa2b'])
set_param([sys,'/','Cap/out_vqa2b'],...
		'position',[595,95,615,115])

add_block('built-in/Integrator',[sys,'/','Cap/vd'])
set_param([sys,'/','Cap/vd'],...
		'Initial','VdCs',...
		'position',[500,192,525,218])

add_block('built-in/Outport',[sys,'/','Cap/out_vda2b'])
set_param([sys,'/','Cap/out_vda2b'],...
		'Port','2',...
		'position',[595,195,615,215])

add_block('built-in/Inport',[sys,'/','Cap/in_iqa2b'])
set_param([sys,'/','Cap/in_iqa2b'],...
		'position',[15,125,35,145])

add_block('built-in/Inport',[sys,'/','Cap/in_ida2b'])
set_param([sys,'/','Cap/in_ida2b'],...
		'Port','2',...
		'position',[50,175,70,195])

add_block('built-in/Fcn',[sys,'/','Cap/Fcn1'])
set_param([sys,'/','Cap/Fcn1'],...
		'Expr','-wb*xc(1)*u[3] + wb*u[1] ',...
		'position',[190,188,470,222])
add_line([sys,'/','Cap'],[135,160;160,160;160,205;185,205])
add_line([sys,'/','Cap'],[475,205;495,205])
add_line([sys,'/','Cap'],[530,105;565,105;565,40;90,40;90,85;105,85])
add_line([sys,'/','Cap'],[475,105;495,105])
add_line([sys,'/','Cap'],[135,160;160,160;160,105;185,105])
add_line([sys,'/','Cap'],[530,105;590,105])
add_line([sys,'/','Cap'],[40,135;105,135])
add_line([sys,'/','Cap'],[75,185;105,185])
add_line([sys,'/','Cap'],[530,205;590,205])
add_line([sys,'/','Cap'],[530,205;560,205;560,285;90,285;90,235;105,235])


%     Finished composite block 'Cap'.

set_param([sys,'/','Cap'],...
		'position',[560,183,600,272])

add_block('built-in/Constant',[sys,'/','Ef'])
set_param([sys,'/','Ef'],...
		'Value','Efo',...
		'position',[50,154,90,176])


%     Subsystem  'ShuntCap'.

new_system([sys,'/','ShuntCap'])
set_param([sys,'/','ShuntCap'],'Location',[145,152,806,492])

add_block('built-in/Mux',[sys,'/','ShuntCap/Mux'])
set_param([sys,'/','ShuntCap/Mux'],...
		'position',[110,60,130,260])

add_block('built-in/Integrator',[sys,'/','ShuntCap/vq'])
set_param([sys,'/','ShuntCap/vq'],...
		'Initial','Vq1',...
		'position',[500,92,525,118])

add_block('built-in/Fcn',[sys,'/','ShuntCap/Fcn'])
set_param([sys,'/','ShuntCap/Fcn'],...
		'Expr','-wb*xc(2)*u[2] - wb*u[4] ',...
		'position',[190,88,470,122])

add_block('built-in/Outport',[sys,'/','ShuntCap/out_vqa2b'])
set_param([sys,'/','ShuntCap/out_vqa2b'],...
		'position',[595,95,615,115])

add_block('built-in/Integrator',[sys,'/','ShuntCap/vd'])
set_param([sys,'/','ShuntCap/vd'],...
		'Initial','Vd1',...
		'position',[500,192,525,218])

add_block('built-in/Outport',[sys,'/','ShuntCap/out_vda2b'])
set_param([sys,'/','ShuntCap/out_vda2b'],...
		'Port','2',...
		'position',[595,195,615,215])

add_block('built-in/Inport',[sys,'/','ShuntCap/in_iqa2b'])
set_param([sys,'/','ShuntCap/in_iqa2b'],...
		'position',[15,125,35,145])

add_block('built-in/Inport',[sys,'/','ShuntCap/in_ida2b'])
set_param([sys,'/','ShuntCap/in_ida2b'],...
		'Port','2',...
		'position',[50,175,70,195])

add_block('built-in/Fcn',[sys,'/','ShuntCap/Fcn1'])
set_param([sys,'/','ShuntCap/Fcn1'],...
		'Expr','-wb*xc(2)*u[3] + wb*u[1] ',...
		'position',[190,188,470,222])
add_line([sys,'/','ShuntCap'],[135,160;160,160;160,205;185,205])
add_line([sys,'/','ShuntCap'],[475,205;495,205])
add_line([sys,'/','ShuntCap'],[530,105;565,105;565,40;90,40;90,85;105,85])
add_line([sys,'/','ShuntCap'],[475,105;495,105])
add_line([sys,'/','ShuntCap'],[135,160;160,160;160,105;185,105])
add_line([sys,'/','ShuntCap'],[530,105;590,105])
add_line([sys,'/','ShuntCap'],[40,135;105,135])
add_line([sys,'/','ShuntCap'],[75,185;105,185])
add_line([sys,'/','ShuntCap'],[530,205;590,205])
add_line([sys,'/','ShuntCap'],[530,205;560,205;560,285;90,285;90,235;105,235])


%     Finished composite block 'ShuntCap'.

set_param([sys,'/','ShuntCap'],...
		'position',[355,175,390,250])


%     Subsystem  'seriesRL'.

new_system([sys,'/','seriesRL'])
set_param([sys,'/','seriesRL'],'Location',[191,509,852,849])

add_block('built-in/Mux',[sys,'/','seriesRL/Mux'])
set_param([sys,'/','seriesRL/Mux'],...
		'inputs','6',...
		'position',[110,58,130,257])

add_block('built-in/Inport',[sys,'/','seriesRL/in_vqa'])
set_param([sys,'/','seriesRL/in_vqa'],...
		'position',[15,95,35,115])

add_block('built-in/Inport',[sys,'/','seriesRL/in_vqb'])
set_param([sys,'/','seriesRL/in_vqb'],...
		'Port','3',...
		'position',[15,165,35,185])

add_block('built-in/Inport',[sys,'/','seriesRL/in_vda'])
set_param([sys,'/','seriesRL/in_vda'],...
		'Port','2',...
		'position',[50,130,70,150])

add_block('built-in/Inport',[sys,'/','seriesRL/in_vdb'])
set_param([sys,'/','seriesRL/in_vdb'],...
		'Port','4',...
		'position',[45,200,65,220])

add_block('built-in/Integrator',[sys,'/','seriesRL/iq'])
set_param([sys,'/','seriesRL/iq'],...
		'Initial','Iq',...
		'position',[500,92,525,118])

add_block('built-in/Fcn',[sys,'/','seriesRL/Fcn'])
set_param([sys,'/','seriesRL/Fcn'],...
		'Expr','(wb/xl(1))*(u[2]-u[4] - r(1)*u[1] -xl(1)*u[6])',...
		'position',[190,88,470,122])

add_block('built-in/Outport',[sys,'/','seriesRL/out_iqa2b'])
set_param([sys,'/','seriesRL/out_iqa2b'],...
		'position',[595,95,615,115])

add_block('built-in/Integrator',[sys,'/','seriesRL/id'])
set_param([sys,'/','seriesRL/id'],...
		'Initial','Id',...
		'position',[500,192,525,218])

add_block('built-in/Outport',[sys,'/','seriesRL/out_ida2b'])
set_param([sys,'/','seriesRL/out_ida2b'],...
		'Port','2',...
		'position',[595,195,615,215])

add_block('built-in/Fcn',[sys,'/','seriesRL/Fcn1'])
set_param([sys,'/','seriesRL/Fcn1'],...
		'Expr','(wb/xl(1))*(u[3]-u[5] - r(1)*u[6] + xl(1)*u[1])',...
		'position',[190,188,470,222])
add_line([sys,'/','seriesRL'],[530,205;560,205;560,275;90,275;90,245;105,245])
add_line([sys,'/','seriesRL'],[135,160;160,160;160,205;185,205])
add_line([sys,'/','seriesRL'],[475,205;495,205])
add_line([sys,'/','seriesRL'],[530,105;565,105;565,40;90,40;90,70;105,70])
add_line([sys,'/','seriesRL'],[475,105;495,105])
add_line([sys,'/','seriesRL'],[135,160;160,160;160,105;185,105])
add_line([sys,'/','seriesRL'],[530,105;590,105])
add_line([sys,'/','seriesRL'],[40,105;105,105])
add_line([sys,'/','seriesRL'],[75,140;105,140])
add_line([sys,'/','seriesRL'],[40,175;105,175])
add_line([sys,'/','seriesRL'],[70,210;105,210])
add_line([sys,'/','seriesRL'],[530,205;590,205])


%     Finished composite block 'seriesRL'.

set_param([sys,'/','seriesRL'],...
		'position',[480,181,520,269])

add_block('built-in/Sum',[sys,'/','Sum3'])
set_param([sys,'/','Sum3'],...
		'inputs','+-',...
		'position',[315,211,330,244])

add_block('built-in/Sum',[sys,'/','Sum2'])
set_param([sys,'/','Sum2'],...
		'inputs','-+',...
		'position',[285,176,300,209])


%     Subsystem  'm3'.

new_system([sys,'/','m3'])
set_param([sys,'/','m3'],'Location',[600,5177521,792,5177620])

add_block('built-in/Note',[sys,'/',['m3/Masked block of m3.m to initalize',13,'s3.m  in Chapter 10 and plot results ']])
set_param([sys,'/',['m3/Masked block of m3.m to initalize',13,'s3.m  in Chapter 10 and plot results ']],...
		'position',[130,60,135,65])
set_param([sys,'/','m3'],...
		'Mask Display','Initialize\nand plot',...
		'Mask Type','Masked block of m3.m',...
		'Mask Dialogue','eval(''m3'')',...
		'Mask Help','Uses m3.m to initialize and plot')


%     Finished composite block 'm3'.

set_param([sys,'/','m3'],...
		'Drop Shadow',4,...
		'position',[671,62,744,99])
add_line(sys,[525,205;555,205])
add_line(sys,[525,250;555,250])
add_line(sys,[335,230;350,230])
add_line(sys,[305,195;350,195])
add_line(sys,[395,195;475,195])
add_line(sys,[395,230;435,230;435,215;475,215])
add_line(sys,[525,205;540,205;540,170;265,170;265,185;280,185])
add_line(sys,[525,250;540,250;540,290;295,290;295,235;310,235])
add_line(sys,[395,195;400,195;400,35;125,35;125,75;140,75])
add_line(sys,[395,230;405,230;405,280;125,280;125,120;140,120])
add_line(sys,[95,165;140,165])
add_line(sys,[105,210;140,210])
add_line(sys,[735,200;770,200;770,295;450,295;450,235;475,235])
add_line(sys,[730,260;755,260;755,310;465,310;475,255])
add_line(sys,[220,170;255,170;255,200;280,200])
add_line(sys,[220,185;250,185;250,220;310,220])
add_line(sys,[605,205;710,205])
add_line(sys,[680,270;680,265;705,265])
add_line(sys,[605,250;705,250])
add_line(sys,[655,165;690,165;690,190;710,190])
add_line(sys,[220,170;255,170;255,95;440,85])
add_line(sys,[220,185;250,185;250,105;455,105;455,85])
add_line(sys,[605,205;660,205;660,135;470,135;470,85])
add_line(sys,[605,250;665,250;665,130;485,130;485,85])
add_line(sys,[220,140;500,140;500,85])
add_line(sys,[220,200;230,200;230,115;515,115;515,85])
add_line(sys,[220,215;240,215;240,125;530,125;530,85])

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