[問題] 時間序列 ar(1) model已回收

看板MATLAB作者 (期末慘淡中...)時間16年前 (2008/09/06 23:26), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/4 (看更多)
y(t) = beta*y(t-1) + e(t) where e(t)~iid~ N(0,1) y(1)=0 beta可為任意參數值 想要用y(1)帶入此公式 生成y(2), y(3), ... y(200)=>用y(2)得到y(3), 用y(3)得到y(4)... 以下是我的作法 clear; phi1 =eye(300); y = linspace(1,0,300); y(1)=1; y(2:300)=0; yt=y'; for times=1:300 for t=1:300 y1(t) = phi1*yt + randn(300,1); y1(times) = phi1*y1(t) + randn(300,1); end end 請問 迴圈的地方不知該如何寫 才能讓我用y(1)算出y(2)後 能夠再用y(2)帶入該式得到y(3)...呢? -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 220.139.131.180 ※ 編輯: tnsshpig 來自: 220.139.131.180 (09/06 23:45) ※ 編輯: tnsshpig 來自: 220.139.131.180 (09/06 23:59)
文章代碼(AID): #18mg4ynI (MATLAB)
文章代碼(AID): #18mg4ynI (MATLAB)