How can i calculate e^A*t
조회 수: 233 (최근 30일)
이전 댓글 표시
How can i calculate e^A*t without using Markov Chain?
Where e=exp , A is a square matrix, and t is a variable
댓글 수: 10
채택된 답변
추가 답변 (5개)
Kye Taylor
2012년 5월 31일
Use the expm function for computing a matrix exponential
댓글 수: 4
KJ N
2017년 11월 9일
exp() only does computes the exponential of A element-by-element, as shown above like this: >> a = [1 2 3 ; 2 5 2; 1 4 3]
a =
1 2 3
2 5 2
1 4 3
>> syms t
>> exp(a*t)
ans =
[ exp(t), exp(2*t), exp(3*t)]
[ exp(2*t), exp(5*t), exp(2*t)]
[ exp(t), exp(4*t), exp(3*t)]
If that's what you're going for, that's great, but not terribly difficult to compute by hand for even somewhat large n x n matrices with integer elements. However, the original poster said they wanted to avoid using the markov chain (a somewhat onerous process, especially when done by hand for large matrices, even with simple integer values as the elements), leading me to understand they were referring to the matrix exponential, not the element-by-element exponential, hence the correct answer in this case would be to use expm(). I had been looking for the same answer, and Kye Taylor was the only post saying use expm instead of exp, so I thought I would try to ensure those in the future looking for the same answer as myself would be helped by a clarification.
Walter Roberson
2017년 11월 9일
We tried a number of times to get the original poster to clarify, but all we got was that they want the exp() solution and that they are looking for a "deeper reason" for something. The poster effectively defined the exp() solution as being the correct one.
Your analysis might well be what the poster really needed, but it is contrary to what little they defined as being correct for their needs.
Shenhai
2017년 1월 20일
편집: Shenhai
2017년 1월 20일
I guess it is not always possible to get the close form solution of exp(At)...
Sometimes I can get result with: exp(At) = iL(sI-A)^-1, where iL is the inverse Laplace transformation, like:
syms s t
A = [0 1;0 0];
expAt = ilaplace(inv(s*eye(size(A,1))-A),s,t);
This will give the result as: [1 t;0 1]
Any other ideas?
댓글 수: 0
Shahram Bekhrad
2012년 6월 8일
As far as I'm aware you probably need it for finding the answer of a state space equation. I myself couldn't find any good function or command yet, so you might have to write a Script file (m-file) and find it. you can use about 3 or 4 way of calculating the said statement. These things are taught in courses like modern control theory. I used the following expression but still have some difficulties. exp(A.t)=I+At+ (At)^2/2! + (At)^3/3!+ (At)^4/4!+. . .
댓글 수: 0
ABCD
2016년 9월 29일
Dear Nick, do you mean this?
>> a = [1 2 3 ; 2 5 2; 1 4 3]
a =
1 2 3
2 5 2
1 4 3
>> syms t >> exp(a*t)
ans =
[ exp(t), exp(2*t), exp(3*t)] [ exp(2*t), exp(5*t), exp(2*t)] [ exp(t), exp(4*t), exp(3*t)]
댓글 수: 1
ABCD
2016년 9월 29일
>> a = [1 2 3 ; 2 5 2; 1 4 3]
a =
1 2 3
2 5 2
1 4 3
>> syms t
>> exp(a*t)
ans =
[ exp(t), exp(2*t), exp(3*t)]
[ exp(2*t), exp(5*t), exp(2*t)]
[ exp(t), exp(4*t), exp(3*t)]
참고 항목
카테고리
Help Center 및 File Exchange에서 Linear Algebra에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!