s3eig.m

all electric machinery dynamic simulation is realized by this book

function [ret,x0,str,ts,xts]=s3eig(t,x,u,flag);
%s3eig	is the M-file description of the SIMULINK system named s3eig.
%	The block-diagram can be displayed by typing: s3eig.
%
%	SYS=s3eig(T,X,U,FLAG) returns depending on FLAG certain
%	system values given time point, T, current state vector, X,
%	and input vector, U.
%	FLAG is used to indicate the type of output to be returned in SYS.
%
%	Setting FLAG=1 causes s3eig to return state derivatives, FLAG=2
%	discrete states, FLAG=3 system outputs and FLAG=4 next sample
%	time. For more information and other options see SFUNC.
%
%	Calling s3eig with a FLAG of zero:
%	[SIZES]=s3eig([],[],[],0),  returns a vector, SIZES, which
%	contains the sizes of the state vector and other parameters.
%		SIZES(1) number of states
%		SIZES(2) number of discrete states
%		SIZES(3) number of outputs
%		SIZES(4) number of inputs
%		SIZES(5) number of roots (currently unsupported)
%		SIZES(6) direct feedthrough flag
%		SIZES(7) number of sample times
%
%	For the definition of other parameters in SIZES, see SFUNC.
%	See also, TRIM, LINMOD, LINSIM, EULER, RK23, RK45, ADAMS, GEAR.

% Note: This M-file is only used for saving graphical information;
%       after the model is loaded into memory an internal model
%       representation is used.

% the system will take on the name of this mfile:
sys = mfilename;
new_system(sys)
simver(1.3)
if (0 == (nargin + nargout))
     set_param(sys,'Location',[-7,83,784,483])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '999999')
set_param(sys,'Min step size', '0.0001')
set_param(sys,'Max step size', '10')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')

add_block('built-in/Note',[sys,'/','vqr'])
set_param([sys,'/','vqr'],...
		'position',[356,172,361,177])

add_block('built-in/Note',[sys,'/','vdr'])
set_param([sys,'/','vdr'],...
		'position',[358,224,363,229])

add_block('built-in/Inport',[sys,'/','in_Ex'])
set_param([sys,'/','in_Ex'],...
		'Port','3',...
		'position',[335,275,355,295])

add_block('built-in/Outport',[sys,'/','out_Pgen'])
set_param([sys,'/','out_Pgen'],...
		'position',[705,55,725,75])

add_block('built-in/Outport',[sys,'/','out_Qgen'])
set_param([sys,'/','out_Qgen'],...
		'Port','2',...
		'position',[665,70,685,90])


%     Subsystem  'VIPQ'.

new_system([sys,'/','VIPQ'])
set_param([sys,'/','VIPQ'],'Location',[217,186,928,574])

add_block('built-in/Mux',[sys,'/','VIPQ/Mux'])
set_param([sys,'/','VIPQ/Mux'],...
		'position',[180,75,210,260])

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

add_block('built-in/Fcn',[sys,'/','VIPQ/Fcn'])
set_param([sys,'/','VIPQ/Fcn'],...
		'Expr','sqrt(u[1]*u[1] + u[3]*u[3])',...
		'position',[285,75,450,105])

add_block('built-in/Fcn',[sys,'/','VIPQ/P'])
set_param([sys,'/','VIPQ/P'],...
		'Expr','u[1]*u[2] + u[3]*u[4]',...
		'position',[285,180,450,210])

add_block('built-in/Fcn',[sys,'/','VIPQ/Q'])
set_param([sys,'/','VIPQ/Q'],...
		'Expr','u[1]*u[4]-u[2]*u[3]',...
		'position',[285,230,450,260])

add_block('built-in/Outport',[sys,'/','VIPQ/Pgen'])
set_param([sys,'/','VIPQ/Pgen'],...
		'Port','3',...
		'position',[510,185,530,205])

add_block('built-in/Outport',[sys,'/','VIPQ/Out_|Vt|'])
set_param([sys,'/','VIPQ/Out_|Vt|'],...
		'position',[510,80,530,100])

add_block('built-in/Fcn',[sys,'/','VIPQ/Fcn1'])
set_param([sys,'/','VIPQ/Fcn1'],...
		'Expr','sqrt(u[2]*u[2] + u[4]*u[4])',...
		'position',[285,125,450,155])

add_block('built-in/Outport',[sys,'/','VIPQ/Qgen'])
set_param([sys,'/','VIPQ/Qgen'],...
		'Port','4',...
		'position',[510,235,530,255])

add_block('built-in/Outport',[sys,'/','VIPQ/Out_|It|'])
set_param([sys,'/','VIPQ/Out_|It|'],...
		'Port','2',...
		'position',[510,130,530,150])

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

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

add_block('built-in/Inport',[sys,'/','VIPQ/-idr'])
set_param([sys,'/','VIPQ/-idr'],...
		'Port','4',...
		'position',[115,225,135,245])
add_line([sys,'/','VIPQ'],[140,100;175,100])
add_line([sys,'/','VIPQ'],[140,145;175,145])
add_line([sys,'/','VIPQ'],[140,190;175,190])
add_line([sys,'/','VIPQ'],[140,235;175,235])
add_line([sys,'/','VIPQ'],[215,170;240,170;240,195;280,195])
add_line([sys,'/','VIPQ'],[455,195;505,195])
add_line([sys,'/','VIPQ'],[455,245;505,245])
add_line([sys,'/','VIPQ'],[215,170;240,170;240,245;280,245])
add_line([sys,'/','VIPQ'],[215,170;240,170;240,90;280,90])
add_line([sys,'/','VIPQ'],[215,170;240,170;240,140;280,140])
add_line([sys,'/','VIPQ'],[455,140;505,140])
add_line([sys,'/','VIPQ'],[455,90;505,90])


%     Finished composite block 'VIPQ'.

set_param([sys,'/','VIPQ'],...
		'position',[545,28,580,87])

add_block('built-in/Inport',[sys,'/','in_vdse'])
set_param([sys,'/','in_vdse'],...
		'Port','2',...
		'position',[55,240,75,260])

add_block('built-in/Inport',[sys,'/','in_vqse'])
set_param([sys,'/','in_vqse'],...
		'position',[55,200,75,220])

add_block('built-in/Product',[sys,'/','prod2'])
set_param([sys,'/','prod2'],...
		'orientation',2,...
		'Font Number',3,...
		'position',[420,308,435,337])

add_block('built-in/Inport',[sys,'/','in_Tmech'])
set_param([sys,'/','in_Tmech'],...
		'Port','4',...
		'position',[535,270,555,290])

add_block('built-in/Outport',[sys,'/','out_(wr-we)//wb'])
set_param([sys,'/','out_(wr-we)//wb'],...
		'Port','5',...
		'position',[705,275,725,295])

add_block('built-in/Outport',[sys,'/','out_delta'])
set_param([sys,'/','out_delta'],...
		'Port','3',...
		'position',[705,170,725,190])

add_block('built-in/Outport',[sys,'/','out_Tem'])
set_param([sys,'/','out_Tem'],...
		'Port','4',...
		'position',[705,240,725,260])


%     Subsystem  'q_cct'.

new_system([sys,'/','q_cct'])
set_param([sys,'/','q_cct'],'Location',[168,350,1046,687])

add_block('built-in/Outport',[sys,'/','q_cct/Out_psiq'])
set_param([sys,'/','q_cct/Out_psiq'],...
		'position',[585,45,605,65])

add_block('built-in/Integrator',[sys,'/','q_cct/psiq_'])
set_param([sys,'/','q_cct/psiq_'],...
		'Initial','Psiqo',...
		'position',[390,127,420,153])

add_block('built-in/Fcn',[sys,'/','q_cct/Fcn'])
set_param([sys,'/','q_cct/Fcn'],...
		'Expr','wb*(u[2]-u[3]+(rs/xls)*(u[1]-u[4]))',...
		'position',[180,125,375,155])

add_block('built-in/Inport',[sys,'/','q_cct/In_wrpsid'])
set_param([sys,'/','q_cct/In_wrpsid'],...
		'Port','2',...
		'position',[50,160,70,180])

add_block('built-in/Inport',[sys,'/','q_cct/In_vqr'])
set_param([sys,'/','q_cct/In_vqr'],...
		'position',[50,105,70,125])

add_block('built-in/Mux',[sys,'/','q_cct/Mux'])
set_param([sys,'/','q_cct/Mux'],...
		'position',[135,86,160,189])

add_block('built-in/Mux',[sys,'/','q_cct/Mux3'])
set_param([sys,'/','q_cct/Mux3'],...
		'inputs','2',...
		'position',[515,92,540,163])

add_block('built-in/Fcn',[sys,'/','q_cct/Fcn3'])
set_param([sys,'/','q_cct/Fcn3'],...
		'Expr','xMQ*(u[1]/xls+u[2]/xplkq)',...
		'position',[565,113,715,147])

add_block('built-in/Fcn',[sys,'/','q_cct/Fcn4'])
set_param([sys,'/','q_cct/Fcn4'],...
		'Expr','-(u[1]-u[2])/xls',...
		'position',[565,202,655,228])

add_block('built-in/Mux',[sys,'/','q_cct/Mux4'])
set_param([sys,'/','q_cct/Mux4'],...
		'inputs','2',...
		'position',[520,179,540,246])

add_block('built-in/Mux',[sys,'/','q_cct/Mux2'])
set_param([sys,'/','q_cct/Mux2'],...
		'inputs','2',...
		'position',[135,213,160,277])

add_block('built-in/Outport',[sys,'/','q_cct/Out_psimq'])
set_param([sys,'/','q_cct/Out_psimq'],...
		'Port','3',...
		'position',[805,120,825,140])

add_block('built-in/Outport',[sys,'/','q_cct/Out_-iq'])
set_param([sys,'/','q_cct/Out_-iq'],...
		'Port','2',...
		'position',[805,205,825,225])

add_block('built-in/Fcn',[sys,'/','q_cct/Fcn2'])
set_param([sys,'/','q_cct/Fcn2'],...
		'Expr','wb*rpkq*(u[1]-u[2])/xplkq',...
		'position',[185,231,370,259])

add_block('built-in/Integrator',[sys,'/','q_cct/psipkq_'])
set_param([sys,'/','q_cct/psipkq_'],...
		'Initial','Psikqo',...
		'position',[395,232,425,258])

add_block('built-in/Scope',[sys,'/','q_cct/qaxis_Scope'])
set_param([sys,'/','q_cct/qaxis_Scope'],...
		'Vgain','3.000000',...
		'Hgain','0.050000',...
		'Vmax','6.000000',...
		'Hmax','2.000000',...
		'Window',[451,130,627,387],...
		'position',[805,252,825,278])
add_line([sys,'/','q_cct'],[75,170;95,170;95,150;130,150])
add_line([sys,'/','q_cct'],[75,115;95,115;95,125;130,125])
add_line([sys,'/','q_cct'],[165,140;175,140])
add_line([sys,'/','q_cct'],[380,140;385,140])
add_line([sys,'/','q_cct'],[425,140;435,140;435,205;115,205;115,175;130,175])
add_line([sys,'/','q_cct'],[165,245;180,245])
add_line([sys,'/','q_cct'],[375,245;390,245])
add_line([sys,'/','q_cct'],[430,245;460,245;460,305;115,305;115,260;130,260])
add_line([sys,'/','q_cct'],[425,140;435,140;435,110;510,110])
add_line([sys,'/','q_cct'],[430,245;460,245;460,145;510,145])
add_line([sys,'/','q_cct'],[545,130;560,130])
add_line([sys,'/','q_cct'],[545,215;560,215])
add_line([sys,'/','q_cct'],[425,140;470,140;470,195;515,195])
add_line([sys,'/','q_cct'],[720,130;755,130;755,315;495,315;495,230;515,230])
add_line([sys,'/','q_cct'],[720,130;755,130;755,315;90,315;90,230;130,230])
add_line([sys,'/','q_cct'],[720,130;755,130;755,80;115,80;115,100;130,100])
add_line([sys,'/','q_cct'],[720,130;800,130])
add_line([sys,'/','q_cct'],[660,215;800,215])
add_line([sys,'/','q_cct'],[425,140;435,140;435,55;580,55])


%     Finished composite block 'q_cct'.

set_param([sys,'/','q_cct'],...
		'position',[410,154,445,206])


%     Subsystem  'Rotor'.

new_system([sys,'/','Rotor'])
set_param([sys,'/','Rotor'],'Location',[250,375,944,671])

add_block('built-in/Outport',[sys,'/','Rotor/wr//wb '])
set_param([sys,'/','Rotor/wr//wb '],...
		'Port','2',...
		'position',[600,195,620,215])

add_block('built-in/Outport',[sys,'/','Rotor/delta'])
set_param([sys,'/','Rotor/delta'],...
		'position',[595,130,615,150])

add_block('built-in/Outport',[sys,'/','Rotor/(wr-we)//wb'])
set_param([sys,'/','Rotor/(wr-we)//wb'],...
		'Port','4',...
		'position',[595,75,615,95])

add_block('built-in/Inport',[sys,'/','Rotor/Tmech'])
set_param([sys,'/','Rotor/Tmech'],...
		'Port','5',...
		'position',[80,60,100,80])

add_block('built-in/Outport',[sys,'/','Rotor/Tem'])
set_param([sys,'/','Rotor/Tem'],...
		'Port','3',...
		'position',[265,203,290,227])

add_block('built-in/Inport',[sys,'/','Rotor/-idr'])
set_param([sys,'/','Rotor/-idr'],...
		'Port','4',...
		'position',[80,210,100,230])