Error Conversion to double from sym is not possible
显示 更早的评论
My code is defined as follows:
U = sym('u', [101 1]);
syms u P Q
P(1,1) = 1;
Q(1,1) = 0;
for k = 1:100
P(k+1,1) = .05*(-.3*P(k,1) +.65*Q(k,1)+U(k,1))+P(k,1);
Q(k+1,1) = .05*(-.65*Q(k,1) + .3*P(k,1)-U(k,1))+Q(k,1);
end
plot(P,U)
I know that the subs() function may be a solution to the conversion for symbolic functions so I also tried a "plot(P, double(subs(U,P)))" command instead, but this also pops up as an error. Does anyone know of another way for me to possibly plot this symbolic matrix?
采纳的回答
Walter Roberson
2013-12-17
Change
P(1,1) = 1;
Q(1,1) = 0;
to
P(1,1) = sym(1);
Q(1,1) = sym(0);
otherwise the arrays will be numeric and not able to store symbolic data.
I see that you tried subs(U,P) but unfortunately for you, all of your P and Q are defined in terms of the undefined symbols u1, u2, and so on. Your P(2) is defined in terms of U(1), your P(3) is defined in terms of P(2) and U(2) and so P(3) is defined in terms of U(1) and U(2), and so on, with your final P(101) being defined in terms of all of the U(1:100).
Perhaps your intent is that U(1) will gain a numeric value, and that U(2) will gain a numeric value because U(1) has gained a numeric value, and U(3) would gain a numeric value because now U(2) has gained a numeric value, and so on, but subs() does a simultaneous substitution, so at best after one step, U(1) would be given a numeric value and U(2) would still be in terms of the symbol U(1) and so on. If you are trying to define U iteratively in such a way that by the end each U will have a numeric value, you cannot do it through simultaneous substitution with subs(): you would need to do the assignment iteratively just like you do for P and Q.
更多回答（1 个）
Azzi Abdelmalek
2013-12-16
0 个投票
What are values affected to U?
2 个评论
Azzi Abdelmalek
2013-12-16
Mary commented
U is a control input on the probability P, so it ranges from 0 to 1.
Azzi Abdelmalek
2013-12-16
How will you plot (P,U) if you don't assign values to U?
类别
在 帮助中心 和 File Exchange 中查找有关 Common Operations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
