#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*--------------------------------------------------------------------
* FlowSheet with simple linear system
*--------------------------------------------------------------------
* Author: Rafael de Pelegrini Soares
* $Id: sample_linear.mso 539 2008-06-16 16:59:56Z rafael $
*--------------------------------------------------------------------*#

using "types";

FlowSheet linear
	PARAMETERS
	nx		as Integer(Default=2);
	ny		as Integer(Default=2);
	nu		as Integer(Default=2);
	
	A(nx,nx) as Real;
	B(nx,nu) as Real;
	C(ny,nx) as Real;
	D(ny,nu) as Real;
	
	VARIABLES
	x(nx)   as Real (Brief="State Variables");
	y(ny)   as Real (Brief="Output Variables");
	u(nu)   as Real (Brief="Control Variables");
	
 	EQUATIONS
	
	diff(x)*'s' = sumt(A*x) + sumt(B*u);
	y           = sumt(C*x) + sumt(D*u);
	
	SPECIFY
	u(1:nu) = sin((time/'s' + [1:nu])*'rad');
	
	SET
	
	nx = 7;
	ny = 3;
	nu = 5;
	
	A = 0;
	B = 0;
	C = 0;
	D = 0;
	
	A(1,1) = 0.5;
	A(2,2) = 0.5;
	A(3,3) = 0.5;
	
	B(1,1) = 0.75;
	B(2,2) = 0.75;
	B(3,3) = 0.75;
	
	C(2,1) = 0.85;
	D = 0.95;
	
  	INITIAL
	x = 0;
 
	OPTIONS
	TimeStep = 0.1;
	TimeEnd = 2;
	TimeUnit = 's';
	SparseAlgebra = true;
	Linearize = true;
end