source: trunk/sample/miscellaneous/sample_linear.mso @ 623

#*-------------------------------------------------------------------
* 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