Can anyone help to write a code for plotting the following equation with time please?

1 次查看（过去 30 天）
Avan Al-Saffar
Avan Al-Saffar 2015-2-9
评论： Torsten 2015-2-10
x = exp( (-B/omega) * cos(omega * t) ) ...
./ ( (B/A)*(integral(exp( (-B/omega)* cos(omega * t) ))))
Where
A= 1;
B= 10;
omega= 1;
x0 =0.1;
t = 0 :0.0001:1000;
  3 个评论
显示 1更早的评论
Roger Stafford
Roger Stafford 2015-2-9
As it stands, the integral in your expression is an indefinite integral and therefore has an arbitrary constant of integration. You need to specify what that constant is in order to successfully plot x as a function of t.

请先登录，再进行评论。

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

采纳的回答

per isakson
per isakson 2015-2-9
编辑：per isakson 2015-2-9
With a little bit of guessing
A= 1;
B= 10;
omega= 1;
x0 =0.1;
t = 0 :0.0001:1000;
fi = @(ti) exp( (-B/omega).*cos(omega*ti) );
fx = @(tj) exp( (-B/omega) .* cos(omega*tj) ) ...
./ ( (B/A).*(integral( fi, 0, tj )));
ezplot( fx, 0:1e-3:12*pi )
produces this
&nbsp
I don't use x0 =0.1; and I get a warning
Warning: Function failed to evaluate on array inputs; [...]
  4 个评论
显示 2更早的评论
Torsten
Torsten 2015-2-10
If lower limit and upper limit of an integral are identical (t=0 in this case), its value is zero - independent of the function to be integrated.
Best wishes
Torsten.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by