s4.m

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

add_block('built-in/Inport',[sys,'/','osc/wr//wb'])
set_param([sys,'/','osc/wr//wb'],...
		'position',[30,45,50,65])

add_block('built-in/Integrator',[sys,'/','osc/sin'])
set_param([sys,'/','osc/sin'],...
		'orientation',2,...
		'position',[265,123,295,157])

add_block('built-in/Gain',[sys,'/','osc/Gain1'])
set_param([sys,'/','osc/Gain1'],...
		'orientation',2,...
		'Gain','wb',...
		'position',[340,120,385,160])

add_block('built-in/Product',[sys,'/','osc/Product1'])
set_param([sys,'/','osc/Product1'],...
		'orientation',2,...
		'position',[140,128,160,172])

add_block('built-in/Outport',[sys,'/','osc/cos_thetar'])
set_param([sys,'/','osc/cos_thetar'],...
		'position',[470,105,490,125])

add_block('built-in/Outport',[sys,'/','osc/sin_thetar'])
set_param([sys,'/','osc/sin_thetar'],...
		'Port','2',...
		'position',[470,160,490,180])
add_line([sys,'/','osc'],[260,140;165,140])
add_line([sys,'/','osc'],[250,80;270,80;270,115;465,115])
add_line([sys,'/','osc'],[395,70;415,70;415,140;390,140])
add_line([sys,'/','osc'],[55,55;365,55])
add_line([sys,'/','osc'],[250,80;365,80])
add_line([sys,'/','osc'],[185,80;210,80])
add_line([sys,'/','osc'],[135,150;115,150;115,80;130,80])
add_line([sys,'/','osc'],[335,140;300,140])
add_line([sys,'/','osc'],[55,55;90,55;90,190;200,190;200,160;165,160])
add_line([sys,'/','osc'],[260,140;230,140;230,170;465,170])


%     Finished composite block 'osc'.

set_param([sys,'/','osc'],...
		'position',[450,205,475,240])


%     Subsystem  'qdr2abc'.

new_system([sys,'/','qdr2abc'])
set_param([sys,'/','qdr2abc'],'Location',[56,346,800,637])

add_block('built-in/Mux',[sys,'/','qdr2abc/Mux'])
set_param([sys,'/','qdr2abc/Mux'],...
		'position',[125,34,155,196])

add_block('built-in/Inport',[sys,'/','qdr2abc/id'])
set_param([sys,'/','qdr2abc/id'],...
		'Port','2',...
		'position',[60,165,80,185])

add_block('built-in/Inport',[sys,'/','qdr2abc/iq'])
set_param([sys,'/','qdr2abc/iq'],...
		'position',[60,125,80,145])

add_block('built-in/Inport',[sys,'/','qdr2abc/sin_thetar'])
set_param([sys,'/','qdr2abc/sin_thetar'],...
		'Port','5',...
		'position',[60,85,80,105])

add_block('built-in/Inport',[sys,'/','qdr2abc/cos_thetar'])
set_param([sys,'/','qdr2abc/cos_thetar'],...
		'Port','4',...
		'position',[60,45,80,65])

add_block('built-in/Fcn',[sys,'/','qdr2abc/qdr2qs'])
set_param([sys,'/','qdr2abc/qdr2qs'],...
		'Expr','u[1]*u[3] + u[2]*u[4]',...
		'position',[190,64,310,86])

add_block('built-in/Fcn',[sys,'/','qdr2abc/qdr2ds'])
set_param([sys,'/','qdr2abc/qdr2ds'],...
		'Expr','u[1]*u[4] - u[2]*u[3]',...
		'position',[190,140,310,160])

add_block('built-in/Mux',[sys,'/','qdr2abc/Mux1'])
set_param([sys,'/','qdr2abc/Mux1'],...
		'inputs','3',...
		'position',[375,39,405,261])

add_block('built-in/Inport',[sys,'/','qdr2abc/i0'])
set_param([sys,'/','qdr2abc/i0'],...
		'Port','3',...
		'position',[60,215,80,235])

add_block('built-in/Fcn',[sys,'/','qdr2abc/qds2b'])
set_param([sys,'/','qdr2abc/qds2b'],...
		'Expr','-u[1]/2 -sqrt(3)*u[2]/2 + u[3]',...
		'position',[435,139,600,161])

add_block('built-in/Fcn',[sys,'/','qdr2abc/qds2a'])
set_param([sys,'/','qdr2abc/qds2a'],...
		'Expr','u[1] + u[3]',...
		'position',[435,69,605,91])

add_block('built-in/Outport',[sys,'/','qdr2abc/ia'])
set_param([sys,'/','qdr2abc/ia'],...
		'position',[645,69,670,91])

add_block('built-in/Outport',[sys,'/','qdr2abc/ib'])
set_param([sys,'/','qdr2abc/ib'],...
		'Port','2',...
		'position',[645,139,670,161])

add_block('built-in/Note',[sys,'/','qdr2abc/qdr2qds'])
set_param([sys,'/','qdr2abc/qdr2qds'],...
		'position',[255,15,260,20])

add_block('built-in/Note',[sys,'/','qdr2abc/qds2abc'])
set_param([sys,'/','qdr2abc/qds2abc'],...
		'position',[530,20,535,25])

add_block('built-in/Outport',[sys,'/','qdr2abc/ic'])
set_param([sys,'/','qdr2abc/ic'],...
		'Port','3',...
		'position',[640,208,665,232])

add_block('built-in/Fcn',[sys,'/','qdr2abc/qds2c'])
set_param([sys,'/','qdr2abc/qds2c'],...
		'Expr','-u[1]/2 +sqrt(3)*u[2]/2 + u[3]',...
		'position',[435,209,605,231])
add_line([sys,'/','qdr2abc'],[85,55;120,55])
add_line([sys,'/','qdr2abc'],[85,95;120,95])
add_line([sys,'/','qdr2abc'],[85,135;120,135])
add_line([sys,'/','qdr2abc'],[85,175;120,175])
add_line([sys,'/','qdr2abc'],[160,115;165,115;165,75;185,75])
add_line([sys,'/','qdr2abc'],[160,115;165,115;165,150;185,150])
add_line([sys,'/','qdr2abc'],[315,75;370,75])
add_line([sys,'/','qdr2abc'],[315,150;370,150])
add_line([sys,'/','qdr2abc'],[85,225;370,225])
add_line([sys,'/','qdr2abc'],[410,150;430,150])
add_line([sys,'/','qdr2abc'],[410,150;420,150;430,80])
add_line([sys,'/','qdr2abc'],[610,80;640,80])
add_line([sys,'/','qdr2abc'],[605,150;640,150])
add_line([sys,'/','qdr2abc'],[410,150;420,150;430,220])
add_line([sys,'/','qdr2abc'],[610,220;635,220])


%     Finished composite block 'qdr2abc'.

set_param([sys,'/','qdr2abc'],...
		'Font Number',3,...
		'position',[545,157,580,243])


%     Subsystem  'VIPQ'.

new_system([sys,'/','VIPQ'])
set_param([sys,'/','VIPQ'],'Location',[284,425,995,813])

add_block('built-in/Outport',[sys,'/','VIPQ/Texc'])
set_param([sys,'/','VIPQ/Texc'],...
		'Port','6',...
		'position',[480,300,500,320])

add_block('built-in/Fcn',[sys,'/','VIPQ/Texcitation'])
set_param([sys,'/','VIPQ/Texcitation'],...
		'Expr','u[2]*Ipm*xmd',...
		'position',[270,295,435,325])

add_block('built-in/Outport',[sys,'/','VIPQ/Trel'])
set_param([sys,'/','VIPQ/Trel'],...
		'Port','5',...
		'position',[480,250,500,270])

add_block('built-in/Fcn',[sys,'/','VIPQ/Treluctance '])
set_param([sys,'/','VIPQ/Treluctance '],...
		'Expr','u[2]*u[4]*(xd-xq)',...
		'position',[270,244,440,276])

add_block('built-in/Outport',[sys,'/','VIPQ/|It|'])
set_param([sys,'/','VIPQ/|It|'],...
		'Port','2',...
		'position',[480,100,500,120])

add_block('built-in/Fcn',[sys,'/','VIPQ/terminal current mag'])
set_param([sys,'/','VIPQ/terminal current mag'],...
		'Expr','sqrt(u[2]*u[2] + u[4]*u[4])',...
		'position',[265,97,440,123])

add_block('built-in/Outport',[sys,'/','VIPQ/Pm'])
set_param([sys,'/','VIPQ/Pm'],...
		'Port','3',...
		'position',[480,155,500,175])

add_block('built-in/Fcn',[sys,'/','VIPQ/Pmotor'])
set_param([sys,'/','VIPQ/Pmotor'],...
		'Expr','u[1]*u[2] + u[3]*u[4]',...
		'position',[265,150,440,180])

add_block('built-in/Outport',[sys,'/','VIPQ/Qm'])
set_param([sys,'/','VIPQ/Qm'],...
		'Port','4',...
		'position',[480,205,500,225])

add_block('built-in/Fcn',[sys,'/','VIPQ/Qmotor'])
set_param([sys,'/','VIPQ/Qmotor'],...
		'Expr','u[1]*u[4]-u[2]*u[3]',...
		'position',[270,201,440,229])

add_block('built-in/Fcn',[sys,'/','VIPQ/terminal voltage mag'])
set_param([sys,'/','VIPQ/terminal voltage mag'],...
		'Expr','sqrt(u[1]*u[1] + u[3]*u[3])',...
		'position',[270,47,440,73])

add_block('built-in/Outport',[sys,'/','VIPQ/|Vt|'])
set_param([sys,'/','VIPQ/|Vt|'],...
		'position',[480,50,500,70])

add_block('built-in/Inport',[sys,'/','VIPQ/id'])
set_param([sys,'/','VIPQ/id'],...
		'Port','4',...
		'position',[115,245,135,265])

add_block('built-in/Inport',[sys,'/','VIPQ/vd'])
set_param([sys,'/','VIPQ/vd'],...
		'Port','3',...
		'position',[115,200,135,220])

add_block('built-in/Inport',[sys,'/','VIPQ/iq'])
set_param([sys,'/','VIPQ/iq'],...
		'Port','2',...
		'position',[115,155,135,175])

add_block('built-in/Inport',[sys,'/','VIPQ/vq'])
set_param([sys,'/','VIPQ/vq'],...
		'position',[115,110,135,130])

add_block('built-in/Mux',[sys,'/','VIPQ/Mux'])
set_param([sys,'/','VIPQ/Mux'],...
		'position',[185,96,215,279])
add_line([sys,'/','VIPQ'],[140,120;180,120])
add_line([sys,'/','VIPQ'],[140,165;180,165])
add_line([sys,'/','VIPQ'],[140,210;180,210])
add_line([sys,'/','VIPQ'],[140,255;180,255])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,165;260,165])
add_line([sys,'/','VIPQ'],[445,165;475,165])
add_line([sys,'/','VIPQ'],[445,215;475,215])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,215;265,215])
add_line([sys,'/','VIPQ'],[445,60;475,60])
add_line([sys,'/','VIPQ'],[445,110;475,110])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,60;265,60])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,110;260,110])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,260;265,260])
add_line([sys,'/','VIPQ'],[445,260;475,260])
add_line([sys,'/','VIPQ'],[440,310;475,310])
add_line([sys,'/','VIPQ'],[220,190;240,190;240,310;265,310])


%     Finished composite block 'VIPQ'.

set_param([sys,'/','VIPQ'],...
		'orientation',3,...
		'position',[251,115,354,135])

add_block('built-in/Mux',[sys,'/','Mux'])
set_param([sys,'/','Mux'],...
		'orientation',3,...
		'inputs','9',...
		'position',[223,60,537,80])

add_block('built-in/Note',[sys,'/','id'])
set_param([sys,'/','id'],...
		'position',[360,265,365,270])

add_block('built-in/Note',[sys,'/','delta'])
set_param([sys,'/','delta'],...
		'position',[365,185,370,190])

add_block('built-in/Note',[sys,'/','iq'])
set_param([sys,'/','iq'],...
		'position',[360,165,365,170])

add_block('built-in/Note',[sys,'/','ia'])
set_param([sys,'/','ia'],...
		'position',[595,150,600,155])

add_block('built-in/Constant',[sys,'/','Ipm'])
set_param([sys,'/','Ipm'],...
		'Value','Ipm',...
		'position',[235,240,265,260])

add_block('built-in/Clock',[sys,'/','Clock'])
set_param([sys,'/','Clock'],...
		'position',[175,100,195,120])


%     Subsystem  'm4'.

new_system([sys,'/','m4'])
set_param([sys,'/','m4'],'Location',[200,4653208,260,4653268])

add_block('built-in/Note',[sys,'/',['m4/Masked block of m4.m to',13,'initialize s4.m and plot results']])
set_param([sys,'/',['m4/Masked block of m4.m to',13,'initialize s4.m and plot results']],...
		'position',[10,45,15,50])
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',[70,53,159,95])

add_block('built-in/To Workspace',[sys,'/','Workspace'])
set_param([sys,'/','Workspace'],...
		'orientation',2,...
		'mat-name','y',...
		'buffer','100000',...
		'position',[310,15,330,35])

add_block('built-in/Scope',[sys,'/','sm1'])
set_param([sys,'/','sm1'],...
		'Vgain','3.000000',...
		'Hgain','1.500000',...
		'Vmax','6.000000',...
		'Hmax','3.000000',...
		'Window',[22,407,780,565])
open_system([sys,'/','sm1'])
set_param([sys,'/','sm1'],...
		'position',[430,13,450,37])
add_line(sys,[480,215;540,215])
add_line(sys,[480,230;540,230])
add_line(sys,[205,220;285,220])
add_line(sys,[380,55;380,25;335,25])
add_line(sys,[205,190;285,190])
add_line(sys,[95,260;115,260;115,220;155,220])
add_line(sys,[95,220;105,220;105,200;155,200])
add_line(sys,[95,180;155,180])
add_line(sys,[480,215;480,325;125,325;125,240;155,240])
add_line(sys,[480,230;490,230;490,315;140,315;140,260;155,260])
add_line(sys,[205,190;265,190;265,140])
add_line(sys,[350,185;355,185;355,170;540,170])
add_line(sys,[350,185;355,185;355,160;290,160;290,140])
add_line(sys,[350,225;445,225])
add_line(sys,[350,285;400,285;400,185;540,185])
add_line(sys,[350,285;400,285;400,155;340,155;340,140])
add_line(sys,[240,280;285,280])
add_line(sys,[205,220;235,220;235,165;315,165;315,140])
add_line(sys,[380,55;380,25;425,25])
add_line(sys,[295,110;295,105;275,105;275,85])
add_line(sys,[310,110;310,85])
add_line(sys,[200,110;240,110;240,85])
add_line(sys,[585,170;610,170;610,130;520,130;520,85])
add_line(sys,[350,245;435,245;435,155;485,155;485,85])
add_line(sys,[350,225;425,225;425,145;450,145;450,85])
add_line(sys,[350,205;415,205;415,85])
add_line(sys,[340,110;340,105;380,105;380,85])
add_line(sys,[325,110;325,100;345,100;345,85])
add_line(sys,[270,250;285,250])

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