source: liacs/la/opdr6/t_vul.m@ 285

Last change on this file since 285 was 2, checked in by Rick van der Zwet, 15 years ago

Initial import of data of old repository ('data') worth keeping (e.g. tracking
means of URL access statistics)

File size: 422 bytes
RevLine 
[2]1orig = [ 0, 50, 31;
2 59, 61, 39;
3 120, 79, 57;
4 181, 88, 69;
5 243, 82, 63;
6 304, 62, 40; ];
7
8y_h = orig(:,2);
9y_l = orig(:,3);
10% Design matrix y = B_0 + B_1 cos((2 * pi * t) / 365) +
11% B_2 sin((2 * pi * t) /365)
12
13for i=1:size(orig,1)
14 t = orig(i,1);
15 X(i,1) = 1;
16 X(i,2) = cos((2 * pi * t) / 365);
17 X(i,3) = sin(( 2 * pi * t) / 365);
18end
19
20X
Note: See TracBrowser for help on using the repository browser.