source: trunk/sample/estimation/alumina.mso @ 999

Last change on this file since 999 was 574, checked in by Rafael de Pelegrini Soares, 14 years ago

Updated the models to work with some language constraints

File size: 1.3 KB
Line 
1#using "types";
2
3FlowSheet Alumina_Flow
4        VARIABLES
5        y as Real;
6        x1 as Real;
7        x2 as Real;
8
9        PARAMETERS
10        b0 as Real(Default=1);
11        b1 as Real(Default=1);
12        b2 as Real(Default=1);
13
14        SPECIFY
15        x1 = 0.647;
16        x2 = 77.1;
17       
18        EQUATIONS
19        y = b0 + b1*x1 + b2*x2;
20       
21        OPTIONS
22        Dynamic = false;
23        #NLASolver(AbsoluteAccuracy=1e-13);
24end
25
26Estimation Alumina_Est as Alumina_Flow
27
28        ESTIMATE
29        # PAR   START  LOWER  UPPER  UNIT
30        b0     -1.00    -5  5;
31        b1      2.00    -5  5;
32        b2     -0.00    -5  5;
33
34        EXPERIMENTS
35        # FILE                WEIGTH    TYPE
36# This sequence works correctly, but the simulation results go all to the first level of variables (fix me)
37        "alumina_est.dat"               1               "fit";
38        "alumina2_est.dat"              1               "test";
39
40# This sequence does not show all data (fix me)
41#       "alumina2_est.dat"              1               "fit";
42#       "alumina_est.dat"      1                "test";
43
44# This sequence does not compute statistics (fix me)
45#       "alumina2_est.dat"              1               "test";
46#       "alumina_est.dat"               1               "fit";
47
48# This sequence has both problems above (fix me)
49#       "alumina_est.dat"               1               "test";
50#       "alumina2_est.dat"              1               "fit";
51
52        OPTIONS
53        Statistics(
54                Fit=true,
55                Parameter=true,
56                Prediction=false
57        );
58        Significance = 0.95;
59        BiLateral = true;
60        NumJac = true;
61       
62        NLPSolver(
63                AbsoluteAccuracy=1e-9,
64                #File = "complex"
65                File = "ipopt_emso"
66        );
67        Dynamic = false;
68end
Note: See TracBrowser for help on using the repository browser.