Function to find ar coefficients using Least squres

1 次查看（过去 30 天）
dav
dav 2013-2-24
Hi,
Can someone help me write a FUNCTION to estimate AR(p) coefficients using LEAST SQUARES?
Thanks
  4 个评论
显示 2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Azzi Abdelmalek
Azzi Abdelmalek 2013-2-24
编辑：Azzi Abdelmalek 2013-2-24
Construct your measure matrix Fi,
yourcoeff=FI\y
Example
% y(n)=a1 y(n-1)+a2 y(n-2)+ a3 y(n-3) your system
% if y=[1 2 3 4 5 6 7 8 9 10]
y= [4;5;6;7] % Using 4 samples
FI=[3 2 1 % measure matrix
4 3 2
5 4 3
6 5 4]
% coeff=[a1;a2;a3] coefficients to find
coeff=FI\y
  3 个评论
显示 1更早的评论
dav
dav 2013-2-27
Here is the code that generates some data. I need to fit an arma model to yt1 and get the coef. estimates with that condition that paremeter estimates are positive.
lc; clear;
m=160; n=10;
ftyPI=zeros(10,100); fsigPI=zeros(10,100);
Lbyth=zeros(10,100); Lbsigh=zeros(10,100);
mccounty =zeros(10,100); mccountsig =zeros(10,100);
for mc= 1:100
a0 = 0.05; a1 = 0.1; b1 = 0.85;
epsi=zeros(3000,1); simsig=zeros(3000,1);
for i = 1:3000 if (i==1) simsig(i) = a0 ; s=(simsig(i))^0.5; epsi(i) = normrnd(0,1) * s; else simsig(i) = a0+ a1*(epsi(i-1))^2+ b1*simsig(i-1); s=(simsig(i))^0.5; epsi(i) = normrnd(0,1) * s; end end
yt1=epsi.^2;

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear and Nonlinear Regression 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by