Integral of controllability gramian

48 次查看（过去 30 天）
Meena Farag
Meena Farag 2016-5-31
评论： Meena Farag 2021-10-12
I am having trouble finding a way to integrate the controllability gramian in Matlab. My system is unstable so I can't use the built in function. My system also has eigenvalues on the imaginary axis so I can't use a function I found online for unstable systems. So i'm trying to integrate the controllability gramian for some finite time interval on MATLAB, but it seems impossible. Here's what the integral looks like.
Where A is a 6x6 matrix and B is a 6x1 matrix. The matrix exponential in the equation is what's causing me the most trouble.
Any ideas?

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

采纳的回答

Star Strider
Star Strider 2016-6-1
编辑：Star Strider 2016-6-1
Try this with your matrix and vector:
A = rand(6); % Create Data
B = rand(6, 1); % Create Data
f = @(tau) expm(A*tau)*B*B'*expm(A'*tau); % Integrand
W = @(t) integral(f, 0, t, 'ArrayValued',1); % Controllability Gramian
Wt = W(1)
The integral function was introduced in R2012a. Before that, I believe the appropriate function is quadv.
  4 个评论
显示 2更早的评论
Meena Farag
Meena Farag 2016-6-2
Awesome. This simple problem has been stressing me out for a while, lol. Thanks again!
Star Strider
Star Strider 2016-6-2
As always, my pleasure!
I was surprised that your question hasn’t been asked before.

请先登录，再进行评论。

更多回答（4 个）

Sheng Cheng
Sheng Cheng 2017-2-16
编辑：Sheng Cheng 2017-2-20
I have a different way for computing the controllability gramian matrix based on an early paper, 'Computing integrals involving the matrix exponential', by Charles Van Loan. The paper can be found here: https://www.cs.cornell.edu/cv/ResearchPDF/computing.integrals.involving.Matrix.Exp.pdf
I will skip the mathematical rigorous proof in the paper. In fact, all the result you need to compute the gramian matrix is written in the left column on the first page. Especially, equation (1.2) is the form we are looking for. (Please read the paper for the extremely simple structure of this integral (actually the controllability gramian is indeed an integral involving matrix exponential)).
The code is just in two lines
A = rand(6); % Create Data
B = rand(6, 1); % Create Data
temp = expm([-A B*B';zeros(6,6) A']); % Coming from the first equation below (1.4)
Wc = temp(7:12,7:12)'*temp(1:6,7:12); % Coming from the second equation below (1.4)
Here, you don't need to define a function and an integral like the one suggested by Star Strider. All you need is expm and then some very simple matrix operation.
  4 个评论
显示 2更早的评论
Star Strider
Star Strider 2017-2-20
In all my control courses, I never encountered that. The paper you attached is in my collection.
+1 Vote!
Meena Farag
Meena Farag 2021-10-12
It's been five years since I originally posted this question while in grad school. I can't believe i'm only coming across this alternative method now. What a great paper. Thanks for sharing!

请先登录，再进行评论。

Ahmed Rashid
Ahmed Rashid 2016-5-31
Why don't you check the rank of the controllability matrix?
C = rank([B AB A^2B ... A^(n-1)*B])
If C has full rank, then the system is controllable.
  4 个评论
显示 2更早的评论
Meena Farag
Meena Farag 2016-6-1
It is only equivalent for stable systems, my system is unstable. I found that I have no choice but to compute the integral for some finite time.
Bryan Jevon
Bryan Jevon 2018-10-30
May you help me with source or literature that provide information about how weak rank test by using controllability and observability matrix?

请先登录，再进行评论。

Roger Stafford
Roger Stafford 2016-6-1
In matlab there is a very important difference between e.^((A.’)*τ) and e^((A.’)*τ) (without the dot.) The first of these is an element-wise exponentiation and the second a matrix exponentiation. If you use exp((A.’)*τ), it will produce the element-wise version. I rather suspect you want the element-wise version.
  2 个评论
Meena Farag
Meena Farag 2016-6-1
Yes, I am using the element-wise version (and expm, not exp). I am still not able to integrate. I still get the same error:
Error using .*
Matrix dimensions must agree.
Error in @(t)(expm(A.*t)*B1)*(expm(A.*t)*B1)'
Error in quad (line 76)
y = f(x, varargin{:});
Error in DOC (line 63)
quad(fn(:,1),0,1)
Star Strider
Star Strider 2016-6-1
Actually, I looked this up in my control reference (and Wikipedia Controllability Gramian). It’s matrix exponentiation, expm.

请先登录，再进行评论。

Rajani Metri
Rajani Metri 2018-12-1
How to calculate Minimum control u*(t) required to state transfer from x1(t) to x2(t) and from it the states x1*(t) and x2(*)? also how to Plot them?
Thank you

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Computations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by