没有合适的资源?快使用搜索试试~ 我知道了~
arma模型参数
3星 · 超过75%的资源 需积分: 12 12 下载量 106 浏览量
2012-11-16
15:23:34
上传
评论
收藏 855B TXT 举报
温馨提示
试读
2页
arma模型参数的预测
资源推荐
资源详情
资源评论
>>ARMA(1,1)模型的参数估计
0
推荐function [phi1,theta1,delta2]=arma110(Z)
%《随机过程》(汪荣鑫编)P148. ARMA(1,1)模型的参数估计
%调用格式: [phi1,theta1,delta2]=arma110(Z)
%Z为时间序列的向量
%Author: Ji Lin, 2007
W=Z-mean(Z);
n=length(Z);K=15;
r=ones(1,K);p=ones(1,K);
r0=W(1:n)*W(1:n)'/n;
for k=1:K
r(k)=W(1:n-k)*W(1+k:n)'/n;
p(k)=r(k)/r0;
end
q=ones(K,K);Q=ones(1,K);
q(1,1)=p(1);Q(1)=q(1,1);
for k=1:K-1
q(k+1,k+1)=(p(k+1)-p(k:-1:1)*q(k,1:k)')/(1-p(1:k)*q(k,1:k)');
for j=1:k
q(k+1,j)=q(k,j)-q(k+1,k+1)*q(k,k-j+1);
end
Q(k+1)=q(k+1,k+1);
end
phi1=p(2)/p(1);
b=(1-2*p(2)+phi1^2)/(p(1)-phi1);
theta=-(b-sqrt(b^2-4))/2;
theta1=-(b+sqrt(b^2-4))/2;
if(abs(theta)<1)
theta1=theta;
0
推荐function [phi1,theta1,delta2]=arma110(Z)
%《随机过程》(汪荣鑫编)P148. ARMA(1,1)模型的参数估计
%调用格式: [phi1,theta1,delta2]=arma110(Z)
%Z为时间序列的向量
%Author: Ji Lin, 2007
W=Z-mean(Z);
n=length(Z);K=15;
r=ones(1,K);p=ones(1,K);
r0=W(1:n)*W(1:n)'/n;
for k=1:K
r(k)=W(1:n-k)*W(1+k:n)'/n;
p(k)=r(k)/r0;
end
q=ones(K,K);Q=ones(1,K);
q(1,1)=p(1);Q(1)=q(1,1);
for k=1:K-1
q(k+1,k+1)=(p(k+1)-p(k:-1:1)*q(k,1:k)')/(1-p(1:k)*q(k,1:k)');
for j=1:k
q(k+1,j)=q(k,j)-q(k+1,k+1)*q(k,k-j+1);
end
Q(k+1)=q(k+1,k+1);
end
phi1=p(2)/p(1);
b=(1-2*p(2)+phi1^2)/(p(1)-phi1);
theta=-(b-sqrt(b^2-4))/2;
theta1=-(b+sqrt(b^2-4))/2;
if(abs(theta)<1)
theta1=theta;
资源评论
- xiaonan1262013-04-15程序不够好,不权威
- mmao2012-12-06还可以,没看懂
- akliying2014-01-13程序不是很好 理解起来有难度
fzwoniu
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功