source: trunk/sample/estimation/sample_estimation.mso @ 997

Last change on this file since 997 was 847, checked in by Rafael de Pelegrini Soares, 13 years ago

Another estimation sample

File size: 1.1 KB
Line 
1
2#****************************************************************
3# Simple static parameter estimation.
4#
5#
6# Example 4.2 from "Applied Parameter Estimation for
7#                                       Chemical Engineers", Englezos & Kalogerakis
8#
9# Author: Rafael de Pelegrini Soares
10#****************************************************************#
11
12FlowSheet SampleEstimationModel
13        VARIABLES
14        y as Real;
15        x as Real;
16
17        PARAMETERS
18        k1 as Real(Default=1.5);
19        k2 as Real(Default=1.5);
20        k3 as Real(Default=1.5);
21
22        SPECIFY
23        x = 1;
24       
25        EQUATIONS
26        y = k1 + k2 * exp(k3*x);
27       
28        OPTIONS
29        Dynamic = false;
30end
31
32
33#
34# Expected solution is:
35# k = [ 523.305 -156.947 -0.199665 ]
36#
37Estimation SampleEstimation as SampleEstimationModel
38
39        ESTIMATE
40        # PAR   START  LOWER  UPPER  UNIT
41        k1      100    0      1000;
42        k2      -200   -300   0;
43        k3      -1     -10    0;
44
45        EXPERIMENTS
46        # FILE                       WEIGTH
47        "sample_estimation.dat"      1;
48
49        OPTIONS
50        NumJac = false;
51        Statistics (
52        Parameter = true
53        );
54       
55        NLPSolver(
56                File = "ipopt_emso"
57                #File = "nmsearch"
58                #File = "complex", MaxIterations=1000
59        );
60        Dynamic = false;
61end
Note: See TracBrowser for help on using the repository browser.