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

2 Ansichten (letzte 30 Tage)
Avan Al-Saffar
Avan Al-Saffar am 9 Feb. 2015
Kommentiert: Torsten am 10 Feb. 2015
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 Kommentare
1 älteren Kommentar anzeigen
Roger Stafford
Roger Stafford am 9 Feb. 2015
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.

Melden Sie sich an, um zu kommentieren.

Melden Sie sich an, um diese Frage zu beantworten.

Akzeptierte Antwort

per isakson
per isakson am 9 Feb. 2015
Bearbeitet: per isakson am 9 Feb. 2015
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 Kommentare
2 ältere Kommentare anzeigen
Avan Al-Saffar
Avan Al-Saffar am 10 Feb. 2015
But I have cos(omega*t) in the denominator so if t =0, I will get 1 .
Torsten
Torsten am 10 Feb. 2015
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.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Fourier Analysis and Filtering finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by