Revision 615, 0.5 kB
(checked in by smidl, 15 years ago)
|
tutorial examples
|
Rev | Line | |
---|
[615] | 1 | % name random variables |
---|
| 2 | y = RV({'y'}); |
---|
| 3 | u = RV({'u'}); |
---|
| 4 | |
---|
| 5 | % create f(y_t| y_{t-3}, u_{t-1}) |
---|
| 6 | fy.class = 'mlnorm<ldmat>'; |
---|
| 7 | fy.rv = y; |
---|
| 8 | fy.rvc = RV({'y','u'}, [1 1], [-3, -1]); |
---|
| 9 | fy.A = [0.5, -0.9]; |
---|
| 10 | fy.const = 0; |
---|
| 11 | fy.R = 1e-100; |
---|
| 12 | |
---|
| 13 | % create f(u_t| ) |
---|
| 14 | fu.class = 'enorm<ldmat>'; |
---|
| 15 | fu.rv = u; |
---|
| 16 | fu.mu = 0; |
---|
| 17 | fu.R = 1e-100; |
---|
| 18 | |
---|
| 19 | % create DS |
---|
| 20 | DS.class = 'MpdfDS'; |
---|
| 21 | DS.mpdf.class = 'mprod'; |
---|
| 22 | DS.mpdf.mpdfs = {fy, epdf2mpdf(fu)}; |
---|
| 23 | DS.init_rv = RV({'y','y','y'}, [1,1,1], [-1,-2,-3]); |
---|
| 24 | DS.init_values = [0.1, 0.2, 0.3]; |
---|
| 25 | |
---|
| 26 | experiment.ndat=100; |
---|
| 27 | |
---|
| 28 | M=estimator(DS,{},experiment); |
---|