※ 引述《jawychiou ((='.'=)...喵嗚)》之銘言：
: 以下是二階AR model的程式碼
: ---------------------------------------------
: clc
: v = wgn(1000,1,0);
: a = 0.1;
: b = -0.8;
: y = zeros(500,1);
: y(1) = 1;
: y(2) = 1;
: for i = 3:1:500
: y(i) = v(i) - a*y(i-1) - b*y(i-2);
: end
: ----------------------------------------------
: stem 這個指令我試過了
: 完全不像
: 還是有其他的用法
: HELP ME 3Q
我試出來了
以下是程式碼
一些參數做些修改
------------------------------------------------
clc
v = wgn(1000,1,0);
a1 = 0.1;
a2 = -0.8;
r0 = cov (v);
y = zeros(500,1);
y(1) = 1;
y(2) = 1;
for i = 3:1:500
y(i) = v(i) - a1*y(i-1) - a2*y(i-2);
end
r = zeros(100,1);
r(1) = (-a1*r0) / (1 + a2);
r(2) = (-a2 + ((a1^2) / (1 + a2)))*r0;
for m = 3:1:100
r(m) = -a1*r(m-1) - a2*r(m-2);
p = r / r0;
end
ov = r0 + a1*r(1) + a2*r(2)
stem (p);
%plot (y);
--------------------------------------------------
如果版大要放到精華區內，
建議放在DSP吧～
AR MODEL屬適應性濾波器的內容。
--
沒有永遠的贏家和輸家；
從贏家變輸家因為他不懂的改變，
從輸家變贏家因為他懂的改變。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.34.102.155