s2.m

simulink electrical machine(2)

function [ret,x0,str,ts,xts]=s2(t,x,u,flag);
%s2	is the M-file description of the SIMULINK system named s2.
%	s2 has a the following characteristics:
%		14 continuous states
%		0 discrete states
%		0 outputs
%		0 inputs
%		does not have direct feedthrough
%		1 sample times
%
%	The block-diagram can be displayed by typing: s2.
%
%	SYS=s2(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 s2 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 s2 with a FLAG of zero:
%	[SIZES]=s2([],[],[],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',[12,55,770,545])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     'tstop')
set_param(sys,'Min step size', '0.00001')
set_param(sys,'Max step size', '0.04')
set_param(sys,'Relative error','1e-6')
set_param(sys,'Return vars',   '')


%     Subsystem  'network'.

new_system([sys,'/','network'])
set_param([sys,'/','network'],'Location',[415,208,1159,701])

add_block('built-in/Mux',[sys,'/','network/Mux'])
set_param([sys,'/','network/Mux'],...
		'position',[265,62,295,213])

add_block('built-in/Sum',[sys,'/','network/Sum'])
set_param([sys,'/','network/Sum'],...
		'position',[450,126,470,179])

add_block('built-in/State-Space',[sys,'/','network/IZ*d'])
set_param([sys,'/','network/IZ*d'],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/','network/IZ*d'],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','IZ\/',...
		'position',[350,184,410,216])

add_block('built-in/State-Space',[sys,'/','network/RZ*q'])
set_param([sys,'/','network/RZ*q'],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/','network/RZ*q'],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','RZ\/',...
		'position',[350,123,410,157])

add_block('built-in/Demux',[sys,'/','network/q'])
set_param([sys,'/','network/q'],...
		'position',[495,75,535,235])

add_block('built-in/Outport',[sys,'/','network/out_iq2e'])
set_param([sys,'/','network/out_iq2e'],...
		'Port','2',...
		'position',[585,125,605,145])

add_block('built-in/Outport',[sys,'/','network/out_vq4e'])
set_param([sys,'/','network/out_vq4e'],...
		'Port','4',...
		'position',[585,205,605,225])

add_block('built-in/Demux',[sys,'/','network/d'])
set_param([sys,'/','network/d'],...
		'position',[495,253,535,407])

add_block('built-in/State-Space',[sys,'/','network/IZ*q'])
set_param([sys,'/','network/IZ*q'],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/','network/IZ*q'],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','IZ\/',...
		'position',[350,259,410,291])

add_block('built-in/Sum',[sys,'/','network/Sum1'])
set_param([sys,'/','network/Sum1'],...
		'inputs','-+',...
		'position',[445,301,465,354])

add_block('built-in/State-Space',[sys,'/','network/RZ*d'])
set_param([sys,'/','network/RZ*d'],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/','network/RZ*d'],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','RZ\/',...
		'position',[350,324,410,356])

add_block('built-in/Mux',[sys,'/','network/Mux1'])
set_param([sys,'/','network/Mux1'],...
		'position',[265,264,295,411])

add_block('built-in/Outport',[sys,'/','network/out_id2e'])
set_param([sys,'/','network/out_id2e'],...
		'Port','6',...
		'position',[585,300,605,320])

add_block('built-in/Outport',[sys,'/','network/out_vd4e'])
set_param([sys,'/','network/out_vd4e'],...
		'Port','8',...
		'position',[585,380,605,400])

add_block('built-in/Inport',[sys,'/','network/in_Eq1e'])
set_param([sys,'/','network/in_Eq1e'],...
		'position',[135,75,155,95])

add_block('built-in/Inport',[sys,'/','network/in_Eq2e'])
set_param([sys,'/','network/in_Eq2e'],...
		'Port','2',...
		'position',[175,110,195,130])

add_block('built-in/Inport',[sys,'/','network/in_vq3e'])
set_param([sys,'/','network/in_vq3e'],...
		'Port','3',...
		'position',[135,145,155,165])

add_block('built-in/Outport',[sys,'/','network/out_iq1e'])
set_param([sys,'/','network/out_iq1e'],...
		'position',[630,85,650,105])

add_block('built-in/Outport',[sys,'/','network/out_iq3e'])
set_param([sys,'/','network/out_iq3e'],...
		'Port','3',...
		'position',[630,165,650,185])

add_block('built-in/Outport',[sys,'/','network/out_id1e'])
set_param([sys,'/','network/out_id1e'],...
		'Port','5',...
		'position',[630,260,650,280])

add_block('built-in/Outport',[sys,'/','network/out_id3e'])
set_param([sys,'/','network/out_id3e'],...
		'Port','7',...
		'position',[625,340,645,360])

add_block('built-in/Inport',[sys,'/','network/in_iq4e'])
set_param([sys,'/','network/in_iq4e'],...
		'Port','4',...
		'position',[170,180,190,200])

add_block('built-in/Inport',[sys,'/','network/in_Ed1e'])
set_param([sys,'/','network/in_Ed1e'],...
		'Port','5',...
		'position',[130,275,150,295])

add_block('built-in/Inport',[sys,'/','network/in_Ed2e'])
set_param([sys,'/','network/in_Ed2e'],...
		'Port','6',...
		'position',[170,310,190,330])

add_block('built-in/Inport',[sys,'/','network/in_vd3e'])
set_param([sys,'/','network/in_vd3e'],...
		'Port','7',...
		'position',[130,345,150,365])

add_block('built-in/Inport',[sys,'/','network/in_id4e'])
set_param([sys,'/','network/in_id4e'],...
		'Port','8',...
		'position',[170,380,190,400])
add_line([sys,'/','network'],[470,330;490,330])
add_line([sys,'/','network'],[475,155;490,155])
add_line([sys,'/','network'],[415,340;440,340])
add_line([sys,'/','network'],[415,275;430,275;440,315])
add_line([sys,'/','network'],[415,200;425,200;425,165;445,165])
add_line([sys,'/','network'],[415,140;445,140])
add_line([sys,'/','network'],[300,140;310,140;310,275;345,275])
add_line([sys,'/','network'],[300,340;345,340])
add_line([sys,'/','network'],[300,340;330,340;330,200;345,200])
add_line([sys,'/','network'],[300,140;345,140])
add_line([sys,'/','network'],[540,95;625,95])
add_line([sys,'/','network'],[540,175;625,175])
add_line([sys,'/','network'],[540,215;580,215])
add_line([sys,'/','network'],[200,120;260,120])
add_line([sys,'/','network'],[540,135;580,135])
add_line([sys,'/','network'],[540,270;625,270])
add_line([sys,'/','network'],[540,350;620,350])
add_line([sys,'/','network'],[540,390;580,390])
add_line([sys,'/','network'],[540,310;580,310])
add_line([sys,'/','network'],[160,85;260,85])
add_line([sys,'/','network'],[160,155;260,155])
add_line([sys,'/','network'],[195,190;260,190])
add_line([sys,'/','network'],[155,285;260,285])
add_line([sys,'/','network'],[195,320;260,320])
add_line([sys,'/','network'],[155,355;260,355])
add_line([sys,'/','network'],[195,390;260,390])


%     Finished composite block 'network'.

set_param([sys,'/','network'],...
		'position',[615,355,700,615])

add_block('built-in/Constant',[sys,'/','vde3'])
set_param([sys,'/','vde3'],...
		'Value','0',...
		'position',[500,550,520,570])


%     Subsystem  'iqe4'.

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

add_block('built-in/Clock',[sys,'/','iqe4/Clock'])
set_param([sys,'/','iqe4/Clock'],...
		'position',[45,60,65,80])

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

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

add_block('built-in/Outport',[sys,'/','iqe4/out_1'])
set_param([sys,'/','iqe4/out_1'],...
		'position',[285,60,305,80])
add_line([sys,'/','iqe4'],[155,70;205,70])
add_line([sys,'/','iqe4'],[70,70;105,70])
add_line([sys,'/','iqe4'],[255,70;280,70])
set_param([sys,'/','iqe4'],...
		'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,'/','iqe4'],...
		'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_iq4e\/iq_iq4e\/')


%     Finished composite block 'iqe4'.

set_param([sys,'/','iqe4'],...
		'position',[540,458,580,482])


%     Subsystem  'tmodel'.

new_system([sys,'/','tmodel'])
set_param([sys,'/','tmodel'],'Location',[152,74,931,624])

add_block('built-in/Note',[sys,'/','tmodel/Eqp_'])
set_param([sys,'/','tmodel/Eqp_'],...
		'position',[320,150,325,155])

add_block('built-in/Integrator',[sys,'/','tmodel/Eqp'])
set_param([sys,'/','tmodel/Eqp'],...
		'Initial','Eqpo(1)',...
		'position',[255,154,285,186])

add_block('built-in/Sum',[sys,'/','tmodel/sum'])
set_param([sys,'/','tmodel/sum'],...
		'inputs','-+-',...
		'position',[130,143,145,197])

add_block('built-in/Gain',[sys,'/','tmodel/1//Tpdo'])
set_param([sys,'/','tmodel/1//Tpdo'],...
		'Gain','1/Tpdo(1)',...
		'position',[160,149,240,191])

add_block('built-in/Outport',[sys,'/','tmodel/out_Pgen'])
set_param([sys,'/','tmodel/out_Pgen'],...
		'Port','3',...
		'position',[675,175,695,195])

add_block('built-in/Outport',[sys,'/','tmodel/out_|I|'])
set_param([sys,'/','tmodel/out_|I|'],...
		'Port','2',...
		'position',[730,150,750,170])

add_block('built-in/Outport',[sys,'/','tmodel/out_|Vt|'])
set_param([sys,'/','tmodel/out_|Vt|'],...
		'position',[675,125,695,145])

add_block('built-in/Outport',[sys,'/','tmodel/out_Qgen'])
set_param([sys,'/','tmodel/out_Qgen'],...
		'Port','4',...
		'position',[730,200,750,220])

add_block('built-in/Sum',[sys,'/','tmodel/Sum'])
set_param([sys,'/','tmodel/Sum'],...
		'inputs','-+',...
		'position',[145,440,160,475])

add_block('built-in/Gain',[sys,'/','tmodel/Gain1'])
set_param([sys,'/','tmodel/Gain1'],...
		'Gain','1/Tpqo(1)',...
		'position',[180,441,265,479])

add_block('built-in/Integrator',[sys,'/','tmodel/Edp'])
set_param([sys,'/','tmodel/Edp'],...
		'Initial','Edpo(1)',...
		'position',[290,443,320,477])

add_block('built-in/Outport',[sys,'/','tmodel/out_Edpe'])
set_param([sys,'/','tmodel/out_Edpe'],...
		'Port','9',...
		'position',[670,465,690,485])

add_block('built-in/Inport',[sys,'/','tmodel/in_iqe'])