# HG changeset patch
# User jwe
# Date 1201849109 0
# Node ID 7d06a404b8f904d77a68a1c0c08555a609be01cf
# Parent  464a55f1a5c2858d4aaa63093cc5beb625afecdd
[project @ 2008-02-01 06:57:50 by jwe]

diff -r 464a55f1a5c2 -r 7d06a404b8f9 scripts/ChangeLog
--- a/scripts/ChangeLog	Fri Feb 01 06:47:48 2008 +0000
+++ b/scripts/ChangeLog	Fri Feb 01 06:58:29 2008 +0000
@@ -1,3 +1,7 @@
+2008-02-01  Dave Goel  <deego3@gmail.com>
+
+	* signal/arch_rnd.m: Correctly index E and Y vectors.
+
 2008-02-01  Bill Denney  <bill@denney.ws>
 
 	* time/weekday.m: Allow vector inputs and speed up.
diff -r 464a55f1a5c2 -r 7d06a404b8f9 scripts/signal/arch_rnd.m
--- a/scripts/signal/arch_rnd.m	Fri Feb 01 06:47:48 2008 +0000
+++ b/scripts/signal/arch_rnd.m	Fri Feb 01 06:58:29 2008 +0000
@@ -63,6 +63,7 @@
     a  = [a, 0];
     la = la + 1;
   endif
+
   lb = length (b);
   b  = reshape (b, 1, lb);
   if (lb == 1)
@@ -79,18 +80,19 @@
   e(1) = sqrt (h(1)) * randn;
   y(1) = b(1) + e(1);
 
-  for t= 2 : M;
+  for t = 2:M
     ta   = min ([t, la]);
-    h(t) = a(1) + a(2:ta) * e(t-1:t-ta+1).^2;
+    h(t) = a(1) + a(2:ta) * e(t-ta+1:t-1).^2;
     e(t) = sqrt (h(t)) * randn;
     tb   = min ([t, lb]);
-    y(t) = b(1) + b(2:tb) * y(t-1:t-tb+1) + e(t);
+    y(t) = b(1) + b(2:tb) * y(t-tb+1:t-1) + e(t);
   endfor
+
   if (T > M)
-    for t = M+1 : T;
-      h(t) = a(1) + a(2:la) * e(t-1:t-la+1).^2;
+    for t = M+1:T
+      h(t) = a(1) + a(2:la) * e(t-la+1:t-1).^2;
       e(t) = sqrt (h(t)) * randn;
-      y(t) = b(1) + b(2:lb) * y(t-1:t-tb+1) + e(t);
+      y(t) = b(1) + b(2:lb) * y(t-tb+1:t-1) + e(t);
     endfor
   endif