Documentation of timeseries

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

r = rand(1, 100000)-0.5;
r = r./std(r); r = r*5;

Cross-Reference Information

Listing of script timeseries

x = zeros(1, 100000);
y = 0;
nday = 15;
a = exp(-(7200./(nday*24*3600)))
b = sqrt((1-a.^2))
r=out7;

for i = 1:100000;
  x(i) = a*y + b*(r(i));
  y = x(i);
end

cd /home/disk/hayes2/dvimont/ocean/coup
load taux;
x = taux;

[p, f] = spectrum(x, 2048, 1024);
f2 = f*6;
subplot(2,1,1);
  plot(log10(f2), f2.*p(:,1), 'k');
  xt = [128 64 32 16 8 4 2 1 0.5 0.25 1/6];
  set(gca, 'XTick', log10(1./xt), 'XTickLabel', xt);
subplot(2,1,2);
  plot([1:1000]/6, x(1:1000), 'k');