eval() returns NaN when my matrix gets really big.
显示 更早的评论
Hi
Im making a model of a mass-spring system where I can increase the numbers of masses. I use modal analysis where I transform my eq's into the Q domain and then back again.
My final result is a [n,1] vector where n is the numbers of masses.
When I use:
n = 4 % Number os masses %
T = 1 % Time which is set as "syms T" %
eval(X) % X is my [n,1] vector %
my result is something like:
ans =
0.0158
0.0739
0.1778
0.8150
but when I pick n to be more then 4 i just get "NaN" and I think it is because my vector gets too big for eval()..
When I look at my X vector it looks like this:
(227072026347336273211456024361092029598452603947709420093308928000*sin(T))/3496036292367407553353910200384804229923697463840499539488117510161+..... and it continues likes this for 2 full pages.
What to do?
Thank you!
采纳的回答
Alexander
2012-4-25
Please try the following:
double(vpa(evalin(symengine, ['subs(' char(X) ', T=' num2str(T) ')'])))
This does the following steps:
- Substitute the sym T inside the symbolic engine.
- Evaluate the result using variable-precision arithmetic.
- Convert the result to double.
2 个评论
Walter Roberson
2012-4-25
I think you should probably be able to simplify that to
double(subs(X))
Adilet Nuratay
2021-1-31
Sorry, I have similar problem.
But I cannot understand the given command. Can anyone clarify it ?
Thank you in advance!
更多回答（1 个）
Mads
2012-4-25
0 个投票
4 个评论
Mads
2012-4-25
Alexander
2012-4-25
Strange, I would have expected Walters suggestion to be much faster but give the NaN answer. Good to know, that both works.
Walter Roberson
2012-4-25
If T is a vector then you would need 'T=[' num2str(T) ']'
Alexander
2012-4-25
Yeah, you are right. My suggestion does not work for vectors, I did not think about that.
类别
在 帮助中心 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
