Figure 9.31:
Base addition rate and LC measurement versus time.
Code for Figure 9.31
Text of the GNU GPL.
main.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51 | % Copyright (C) 2001, James B. Rawlings and John G. Ekerdt
%
% This program is free software; you can redistribute it and/or
% modify it under the terms of the GNU General Public License as
% published by the Free Software Foundation; either version 2, or (at
% your option) any later version.
%
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
% General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program; see the file COPYING. If not, write to
% the Free Software Foundation, 59 Temple Place - Suite 330, Boston,
% MA 02111-1307, USA.
% Made this m file in the second printing to get a proper staircase
% plot of the flowrate versus time
% jbr, 12/27/03
%
% This m-file loads data file 'lc.dat'.
% This m-file loads data file 'flow.dat'.
vro = 2370;
teaf = 0.00721;
teaden = 0.728;
flow = load('flow.dat');
lc = load('lc.dat');
tflow = [0. ; flow(:,1)];
flowrate = [0. ; flow(:,2)./teaden];
tlc = lc(:,1);
lcmeas = lc(:,2);
[tflows, flowrates] = stairs (tflow, flowrate);
data = struct();
data.flow = [tflows, flowrates];
data.lc = lc;
gnuplotsave('flowlc.dat', data);
if (~ strcmp (getenv ('OMIT_PLOTS'), 'true')) % PLOTTING
subplot(2,1,1);
plot (tflows, flowrates);
% TITLE flowlc
subplot(2,1,2);
plot (tlc, lcmeas, '+');
% TITLE flowlc
end % PLOTTING
|
lc.dat
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35 | 414 0.17050
424 0.16040
434 0.13250
444 0.10840
493 0.10140
503 0.10420
513 0.10450
523 0.09700
533 0.08240
582 0.06780
592 0.06970
602 0.07020
612 0.06920
621 0.06840
631 0.06750
641 0.06730
651 0.06730
661 0.06540
671 0.05770
681 0.05710
690 0.05630
700 0.04660
710 0.04550
720 0.04490
730 0.04560
740 0.04500
750 0.04340
759 0.04280
769 0.04260
779 0.04090
789 0.03990
799 0.03970
848 0.03940
858 0.03740
868 0.03710
|
flow.dat
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87 | 9.000000 9.7817838E-02 4.786000
19.00000 1.290601 14.21800
29.00000 1.262779 13.99800
39.00000 1.240016 13.81800
49.00000 1.215736 13.62600
59.00000 1.187156 13.40000
69.00000 1.160599 13.19000
79.00000 0.5545961 8.398001
89.00000 1.109508 12.78600
99.00000 1.361925 14.78200
109.0000 1.332333 14.54800
119.0000 1.294648 14.25000
129.0000 1.271885 14.07000
139.0000 1.244569 13.85400
149.0000 1.220794 13.66600
159.0000 1.197778 13.48400
169.0000 1.170463 13.26800
179.0000 1.150988 13.11400
189.0000 1.128983 12.94000
199.0000 1.105968 12.75800
209.0000 1.092563 12.65200
219.0000 1.360155 14.76800
229.0000 1.344473 14.64400
239.0000 1.329045 14.52200
249.0000 1.305524 14.33600
259.0000 1.275426 14.09800
269.0000 1.259239 13.97000
279.0000 1.244063 13.85000
289.0000 1.219783 13.65800
299.0000 1.191961 13.43800
309.0000 1.173751 13.29400
319.0000 1.158322 13.17200
329.0000 0.3163429 6.514000
339.0000 0.5495375 8.358001
349.0000 1.224841 13.69800
359.0000 0.3228684 6.565600
369.0000 9.4845297E-04 4.020000
379.0000 4.6794413E-04 4.016201
389.0000 5.4383278E-04 4.016800
399.0000 4.9326418E-04 4.016400
409.0000 3.6680698E-04 4.015400
419.0000 0.2701340 6.148601
429.0000 1.185891 13.39000
439.0000 0.9299334 11.36600
449.0000 0.5259147 8.171202
459.0000 8.9793204E-04 4.019600
469.0000 8.7256433E-04 4.019400
479.0000 7.9672335E-04 4.018800
489.0000 4.9319270E-04 4.016400
499.0000 4.6799183E-04 4.016200
509.0000 3.6678315E-04 4.015400
519.0000 0.2330809 5.855600
529.0000 0.7784327 10.16800
539.0000 7.5244429E-03 4.072000
549.0000 7.5244186E-03 4.072000
559.0000 5.7539940E-03 4.058000
569.0000 8.0302954E-03 4.076001
579.0000 0.5019882 7.982001
589.0000 -3.0983209E-03 3.988000
599.0000 9.4850064E-04 4.020000
609.0000 -6.3252446E-05 4.012000
619.0000 -1.3278484E-03 4.002000
629.0000 -4.6158792E-03 3.976000
639.0000 -6.3204767E-05 4.012000
649.0000 -2.0865917E-03 3.996000
659.0000 0.1752122 5.398000
669.0000 0.3024322 6.404000
679.0000 4.9952744E-03 4.052001
689.0000 3.9835931E-03 4.044001
699.0000 0.5114223 8.056601
709.0000 3.4151078E-04 4.015200
719.0000 3.1619071E-04 4.015000
729.0000 1.1386871E-04 4.013400
739.0000 6.3300133E-05 4.013000
749.0000 3.4151078E-04 4.015200
759.0000 0.1059619 4.850400
769.0000 6.7019463E-04 4.017800
779.0000 6.8099439E-02 4.551000
789.0000 8.7261200E-04 4.019400
799.0000 9.7372534E-04 4.020200
809.0000 8.7256433E-04 4.019400
819.0000 7.9669955E-04 4.018800
829.0000 7.2083471E-04 4.018200
839.0000 5.6912901E-04 4.017000
849.0000 4.6801567E-04 4.016200
859.0000 5.1856041E-04 4.016600
869.0000 4.4267176E-04 4.016000
|