blob: 681106b4bacf050718e3855cd0e06671e936d127 (
plain)
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
88
89
|
#N canvas 1 53 991 458 10;
#X obj 21 12 inlet;
#X obj 260 12 inlet;
#X obj 260 66 purepd/float_argument \$1 0.15;
#X obj 271 39 loadbang;
#X obj 59 421 outlet;
#X obj 21 57 trigger float bang;
#X obj 126 153 float;
#X text 162 153 St1;
#X obj 260 102 clip 0 1;
#X obj 175 421 outlet;
#X obj 450 12 inlet;
#X obj 461 39 loadbang;
#X obj 450 102 clip 0 1;
#X text 428 265 bt = y * (St - St1) + (1 - y) * bt1;
#X text 509 102 y;
#X text 319 102 a;
#X obj 176 265 expr $f4 * ($f1 - $f2) + (1 - $f4) * $f3;
#X obj 294 154 float;
#X text 332 153 bt1;
#X obj 450 66 purepd/float_argument \$2 0.3;
#X obj 59 207 expr $f4*$f1+(1-$f4) * ($f2 + $f3);
#X obj 302 391 pddp/pddplink http://www.itl.nist.gov/div898/handbook/pmc/section4/pmc435.htm
;
#X obj 245 421 outlet;
#X obj 630 12 inlet;
#X obj 641 39 loadbang;
#X obj 630 102 clip 0 1;
#X text 689 102 B;
#X text 274 208 St = a*(Yt / It1+(1-a)*St1;
#X text 524 223 we are only using 2L as the "season";
#X obj 245 349 expr $f4 * ($f2/$f1) + (1 - $f4) * $f3;
#X text 482 351 It = B * (Yt/St) + (1 - B) * It1;
#X text 467 311 It;
#X obj 630 66 purepd/float_argument \$3 0.6;
#X obj 810 12 inlet;
#X obj 821 39 loadbang;
#X obj 810 102 int;
#X obj 810 127 min 1;
#X text 857 99 how many periods;
#X obj 810 66 purepd/float_argument \$4 1;
#N canvas 198 494 450 300 keep 0;
#X obj 313 36 inlet;
#X obj 214 255 outlet;
#X obj 71 34 inlet;
#X obj 71 92 list append;
#X obj 186 165 list split;
#X connect 0 0 4 1;
#X connect 2 0 3 0;
#X connect 3 0 4 0;
#X connect 4 0 3 1;
#X connect 4 1 1 0;
#X connect 4 2 3 1;
#X restore 334 309 pd keep old It values;
#X connect 0 0 5 0;
#X connect 1 0 2 0;
#X connect 2 0 8 0;
#X connect 3 0 2 0;
#X connect 5 0 20 0;
#X connect 5 1 6 0;
#X connect 5 1 17 0;
#X connect 6 0 16 1;
#X connect 6 0 20 1;
#X connect 8 0 20 3;
#X connect 10 0 19 0;
#X connect 11 0 19 0;
#X connect 12 0 16 3;
#X connect 16 0 17 1;
#X connect 16 0 9 0;
#X connect 16 0 29 1;
#X connect 17 0 16 2;
#X connect 17 0 20 2;
#X connect 19 0 12 0;
#X connect 20 0 4 0;
#X connect 20 0 6 1;
#X connect 20 0 16 0;
#X connect 20 0 29 0;
#X connect 23 0 32 0;
#X connect 24 0 32 0;
#X connect 25 0 29 3;
#X connect 29 0 22 0;
#X connect 29 0 39 0;
#X connect 32 0 25 0;
#X connect 33 0 38 0;
#X connect 34 0 38 0;
#X connect 35 0 36 0;
#X connect 36 0 39 1;
#X connect 38 0 35 0;
#X connect 39 0 29 2;
|