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
|
// Demonstrate the order of MA(q) as discussed in Example 6.11 on page 182.
// 6.5
exec('plotacf.sci',-1);
exec('label.sci',-1);
xi = 0.1*rand(1,10000,'normal');
a = 1; b = [];
d = [1 1 -0.5];
ar = armac(a,b,d,1,1,1);
v = arsimul(ar,xi);
z = [v xi];
// Plot noise, plant output and ACF
subplot(2,1,1), plot(v(1:500))
label('',4,'','v',4)
subplot(2,1,2), plot(xi(1:500))
label('',4,'n','xi',4)
xset('window',1)
plotacf(v,1,11,1);
|
