Integral: "Warning: Minimum step size reached"
22 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Hi,
I'm trying to calculate an Integral the following:
xmin=c0/b;
xmax=zn+c0/b;
tw=@(w)(1./(b.*w.*sqrt(1-power(a,2).*power(b,2).*power(w,2))));
r=@(w)(a.*b.*w)./sqrt(1-power(a,2).*power(b,2).*power(w,2));
rx=integral(r,xmin,xmax)
t=integral(tw,xmin,xmax)
For
a=0.00067024
b=-0.06
c0=1480
cn=1492
z0=1000
zn=799.9523
it is working fine, however, for
a=0.00067024
b=0.0133
c0=1480
cn=1492
z0=1000
zn=1902.3
it gives me the following error:
Warning: Minimum step size reached near x = 112181. There may be a singularity, or the
tolerances may be too tight for this problem.
What am I doing wrong?
0 comentarios
Respuestas (3)
Star Strider
el 27 de Sept. de 2015
It is a warning, not an error. It is telling you that the integral function has encountered s singularity and cannot continue.
I would plot the function that is throwing that error to see what the function is doing at that value of the independent variable. You might want to integrate it piecewise below and above the region of that independent variable value.
0 comentarios
Anke Kügler
el 27 de Sept. de 2015
Editada: Walter Roberson
el 27 de Sept. de 2015
1 comentario
Star Strider
el 27 de Sept. de 2015
The ‘i’ means you do have a complex value for the integration. It probably means you are taking the sqrt of a negative value.
Walter Roberson
el 27 de Sept. de 2015
In your formula r, if power(a,2)*power(b,2)*power(w,2) is greater than 1 then 1 minus that expression is negative and you would be taking sqrt() of a negative number.
power(a,2)*power(b,2) is fixed for any one run, so you can calculate the necessary w as
1/sqrt(power(a,2)*power(b,2))
or more simply 1/(a*b) . And if that is in the range xmin to xmax then your integral involves complex numbers.
0 comentarios
Ver también
Categorías
Más información sobre Variables en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!