From: simpark park on
Hi.
Anybody know to get 'only' parameter matrix of ARX model, not equation, like arburg of ar model etc ?
Many thanks
From: Wayne King on
"simpark park" <bullet7(a)empal.com> wrote in message <hrqrjo$9jf$1(a)fred.mathworks.com>...
> Hi.
> Anybody know to get 'only' parameter matrix of ARX model, not equation, like arburg of ar model etc ?
> Many thanks

Hi Simpark, invoke the polydata method on the idpoly object.

A = [1 -1.5 0.7]; B = [0 1 0.5];
m0 = idpoly(A,B);
u = iddata([],idinput(300,'rbs'));
e = iddata([],randn(300,1));
y = sim(m0, [u e]);
z = [y,u];
m = arx(z,[2 2 1]);

% Now m is an idpoly object, invoke polydata method

[Afit,Bfit] = polydata(m);


Hope that helps,
Wayne
From: simpark park on
Thanks alot Wayne.
It was really helpful.


"Wayne King" <wmkingty(a)gmail.com> wrote in message <hrrgjg$945$1(a)fred.mathworks.com>...
> "simpark park" <bullet7(a)empal.com> wrote in message <hrqrjo$9jf$1(a)fred.mathworks.com>...
> > Hi.
> > Anybody know to get 'only' parameter matrix of ARX model, not equation, like arburg of ar model etc ?
> > Many thanks
>
> Hi Simpark, invoke the polydata method on the idpoly object.
>
> A = [1 -1.5 0.7]; B = [0 1 0.5];
> m0 = idpoly(A,B);
> u = iddata([],idinput(300,'rbs'));
> e = iddata([],randn(300,1));
> y = sim(m0, [u e]);
> z = [y,u];
> m = arx(z,[2 2 1]);
>
> % Now m is an idpoly object, invoke polydata method
>
> [Afit,Bfit] = polydata(m);
>
>
> Hope that helps,
> Wayne
From: simpark park on
Thanks alot Wayne.
It was really helpful.

"Wayne King" <wmkingty(a)gmail.com> wrote in message <hrrgjg$945$1(a)fred.mathworks.com>...
> "simpark park" <bullet7(a)empal.com> wrote in message <hrqrjo$9jf$1(a)fred.mathworks.com>...
> > Hi.
> > Anybody know to get 'only' parameter matrix of ARX model, not equation, like arburg of ar model etc ?
> > Many thanks
>
> Hi Simpark, invoke the polydata method on the idpoly object.
>
> A = [1 -1.5 0.7]; B = [0 1 0.5];
> m0 = idpoly(A,B);
> u = iddata([],idinput(300,'rbs'));
> e = iddata([],randn(300,1));
> y = sim(m0, [u e]);
> z = [y,u];
> m = arx(z,[2 2 1]);
>
> % Now m is an idpoly object, invoke polydata method
>
> [Afit,Bfit] = polydata(m);
>
>
> Hope that helps,
> Wayne