Linear Grey-Box system identification: idgrey()

9 Ansichten (letzte 30 Tage)
wei zhang
wei zhang am 3 Aug. 2019
Kommentiert: wei zhang am 5 Aug. 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?

Akzeptierte Antwort

Rajiv Singh
Rajiv Singh am 5 Aug. 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 Kommentar
wei zhang
wei zhang am 5 Aug. 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

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Linear Model Identification finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by