清空筛选条件
清空筛选条件

Different result command window and guide

1 次查看（过去 30 天）
john
john 2014-3-6
编辑： Mischa Kim 2014-3-7
Hello,
if I write in command window:
t=1
4.5582758767110180218764941687633*sin(314.0*t) - 16.557811525543529374688574837543*cos(314.0*t) + (16.557811525543529374688574837543*( 3.4028e+308*(1375.8856146096578468897575988832*t) + 0.48398685287825251446260421986529* 3.4028e+308*(1375.8856146096578468897575988832*t)))/exp(1485.7142857142857142857142857143*t)
I get NaN,but in guide I get -17.0711.....why?

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

回答（1 个）

Mischa Kim
Mischa Kim 2014-3-6
John, in your expression you have two huge terms, basically,
Inf/Inf
which results in a NaN. You should get the same result for the same expression in your GUI.
  3 个评论
显示 1更早的评论
john
john 2014-3-6
Or how can I plot for any time:
92.9*cos(314.0*t) + 96.3*sin(314.0*t) - (7.94*(cosh(1399.0*t) + 4.06*sinh(1399.0*t)))/exp(1500.0*t)
Mischa Kim
Mischa Kim 2014-3-7
编辑：Mischa Kim 2014-3-7
You can do
syms t
expr = 92.9*cos(314.0*t) + 96.3*sin(314.0*t) - (7.94*(cosh(1399.0*t) + 4.06*sinh(1399.0*t)))/exp(1500.0*t);
expr = subs(expr, sinh(1399.0*t),(exp(1399.0*t) - exp(-1399.0*t))/2);
expr = subs(expr, cosh(1399.0*t),(exp(1399.0*t) + exp(-1399.0*t))/2);
expr = expand(expr,'ArithmeticOnly',true);
ezplot(expr,[0,0.1])
The two subs titutions are hyperbolic trig identities, thereby avoiding the Inf/Inf issue (in combination with the expand command).

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by