MuPAD unable to get the expected representation

4 次查看（过去 30 天）
Zhu
Zhu 2013-4-19
Hi,
When I open MuPAD Notebook and type the following commands:
Ac := matrix([[0,1,0], [0,0,1], [0,0,0]]):
Dc := matrix([0,0,1]):
expre1 := exp((tk-tau)*Ac)*Dc:
expre2 := expre1*q*expre1':
Q := int(expre2, tau = tk-T..tk):
Simplify(Q)
to calculate the matrix Q, the result returned is an integral representation.
How to get the following expected results by using MuPAD,
Q(1,1) = (T^5*q)/20
Q(1,2) = (T^4*q)/8
Q(1,3) = (T^3*q)/6
Q(2,1) = (T^4*q)/8
Q(2,2) = (T^3*q)/3
Q(2,3) = (T^2*q)/2
Q(3,1) = (T^3*q)/6
Q(3,2) = (T^2*q)/2
Q(3,3) = T*q
Thanks in advance.
  2 个评论
Walter Roberson
Walter Roberson 2013-4-19
Ac is going to be 3 x 3.
Dc is going to be 1 x 3.
For the exp((tk-tau)*Ac) part, are you intending matrix exponentiation or element-by-element exponentiation?
Either way, I think the result will be 3x3 -- and that cannot be right-multiplied by 1 x 3
Zhu
Zhu 2013-4-19
Thanks for reply. As you say, Ac is a 3-by-3 matrix
In MuPAD, Dc := matrix([0,0,1]) is equivalent to Dc := matrix(3,1,[0,0,1]), so Dc is a 3-by-1 column vector
and exp((tk-tau)*Ac) is the matrix exponential of (tk-tau)*Ac
I check the code again. It seems the syntax is valid, but can't get the result as expected

请先登录，再进行评论。

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

采纳的回答

Walter Roberson
Walter Roberson 2013-4-19

更多回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by