newtonraphson.m

来自「对平稳时间序列中的参数估计方法（包括矩估计」· M 代码 · 共 22 行

q=input('输入滑动平均阶数：')
max=input('输入所要求的精度：')
V=input('输入样本方差：');
for j=1:q+1
    E(j)=input('输入各阶样本自协方差函数：');           %E(j)为j-1阶样本自协方差函数，注意这里是j-1阶，其中第一个就是样本方差
end
E=E';
A=zeros(q+1,1);
A(1,1)=sqrt(V);
while 1
    AA=A-inv(jz1(A,q))*jz2(A,E,q);                      %迭代得到的值
    if abs(AA-A)<max                                    %迭代精度
        break
    else
        A=AA;
    end
end
for j=2:q+1                      
    F(j)=-A(j)/A(1);                                     %F(j)的值为滑动平均参数的估计值
end
F(1)=A(1)*A(1);                                          %F(1)的值为噪声方差的估计值
F=F'