Linear Grey-Box system identification: idgrey()

2 visualizaciones (últimos 30 días)
wei zhang
wei zhang el 3 de Ag. de 2019
Comentada: wei zhang el 5 de Ag. de 2019
Hello all,
I have a double about the idgrey() function for system identification of linear system.
It is clear that we can either use continuous-time or discrete-time state space description by using idgrey().
The two (continous-time or discrete-time) have different matrix A,B,C,D. for the same dymamic system.
Now, I start with a continous-time state space with idgrey() with 'c' as continous parameter. I can't define the parameter sampling time T>0: this is normal.
My question is: if I use 'cd' as continous-discrete time parameter, and I can assign the sampling time T>0, which kind of state-space description I will get?
Is it a continous-time description with a positive sampling time? as we use the function c2d()?
2.
My intention is to estimate a linear grey-box model with continous-time description, but I have only the sampling data in discrete time.
How can I realize the system identification effectively?
One way, I am using now is idgrey() + idss() + c2d() and then greyest().
is this way correct?

Respuesta aceptada

Rajiv Singh
Rajiv Singh el 5 de Ag. de 2019
The 'cd' option is to support the use case where you want to specify your own discretization formula (not use the standard c2d formulas). In your case, you should declare the file type to be 'c'. The greyest algorithm will take care of any required discretizations internally. Note that the discretization scheme is either zoh or foh dependeing upon your data's InterSample property value.
So the process is quite straightforward:
model = idgrey(file, initial_parameters, 'c');
model2 = greyest(data, model);
  1 comentario
wei zhang
wei zhang el 5 de Ag. de 2019
Hi Rajiv
Thanks for your answer.
In fact, I use the following method through my own discretization formula:
I start from continous time form A, B, C, D; then use the clever trick in
Finally, use greyest(), as you suggest.
Best

Iniciar sesión para comentar.

Más respuestas (0)

Community Treasure Hunt

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

Start Hunting!

Translated by