How to solve a 2nd degree equation of matrix ?

13 次查看（过去 30 天）
Tahn
Tahn 2013-5-26
Hello,
Pretty new with matlab, I'd like to know how to solve a second-degree equation of matrix : P*B0-P²*C=B1 All of the variables are 7*7 matrix, and I'm looking for P.
I used
temp=solve('P*B0-(P^2)*C=B1')
subs(temp)
I get a result, but it's mainly with NaN and Inf. As it shouldn't be the case, I guess my method isn't the right one, but I've no clue how to do it, and google doesn't really help.
So I'm counting on you !
Thanks !

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

回答（3 个）

Stephan M. H.
Stephan M. H. 2013-5-26
Hi,
Unfortunately you didn't specify what you have and how you defined it but...
the function solve needs symbolic arguments, i.e., you need to create all matrices with the commands sym or syms
The values of B and C matrices should be given to you. The unkown P you could define by (e.g., 2x2 matrix)
syms p1 p2 p3 p4
P_s = sym([p1 p2; p3 p4])
B0_s = sym(B0)
...
Then use
RES = solve(P*B0-P^2*C==B1) % (without '...')
The command subs is only needed if you don't have values for the B and C matrices.
Hope this helps,
Stephan
Tahn
Tahn 2013-5-26
编辑：Tahn 2013-5-26
Thanks for your answer, I'm going to try it right away.
My bad, I didn't specify that I defined B0, B1 and C (Basically, withtout the command sym, just with :
B0 = [1,-0.4,0.44,0,0,0,0;0,1 ...]
So I have to define p1 until p49, and then write
P_s=sym(p1...p49) ?
That seems pretty long though ! But doable, so let's see !
Edit : my bad, I read too fast and didn't understood what sym did. I just have to rewrite P with all the baby variables inside, and change the form of the 3 know matrix with sym.
Tahn
Tahn 2013-5-26
After about 45 minutes of calculation, it still didn't work :( :
>> RES = solve(P_s*B0_s-P_s^2*C_s==B1_s)
Warning: Solutions might be lost. [solvelib::allvalues_warning]
RES =
p1: [1x1 sym]
p2: [1x1 sym]
p3: [1x1 sym]
...
p48: [1x1 sym]
p49: [1x1 sym]
Et si je veux voir ce qu'il y a dans p1 :
>> p1
p1 =
p1
  1 个评论
Stephan M. H.
Stephan M. H. 2013-5-27
I can't help with the first warning, and for the second problem I can only guess what you are writing there in French, but you might want to try
RES.p1
if you used
syms p1 ....
to define everything, p1 etc. still contain only symbolic values. the results is stored in the structure RES
best, Stephan

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Conversion Between Symbolic and Numeric 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by