Linear Grey-Box system identification: idgrey()

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?

 채택된 답변

Rajiv Singh
Rajiv Singh 2019년 8월 5일

0 개 추천

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

wei zhang
wei zhang 2019년 8월 5일
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

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Grey-Box Model Estimation에 대해 자세히 알아보기

질문:

2019년 8월 3일

댓글:

2019년 8월 5일

Community Treasure Hunt

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

Start Hunting!

Translated by