Does Matlab's integral function pass the integration variable to the function handle one scalar at a time or as a vector?
32 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
yuval dahan
el 24 de Jun. de 2024 a las 7:29
Comentada: yuval dahan
el 24 de Jun. de 2024 a las 14:54
Hello,
My integrand function contains a sum that depends on the integration variable, as follows:
tau = 1;
T = 20;
q = integral(@(t) f(t,tau),0,T)
function y = f(t,tau)
N = 0:floor(t,tau)
for i = 1:length(N)
n = N(i);
y = y + g(n,mod(t,tau) + (n+1)*tau) - g(n,mod(t,tau) + n*tau);
end
Where g is some function.
I'm not sure if the integral function is handling it correctly, because if it passes a vector t, the sum will not iterate over the correct range.
I would really appreciate any insights!
Thanks!
0 comentarios
Respuesta aceptada
Stephen23
el 24 de Jun. de 2024 a las 7:42
Editada: Stephen23
el 24 de Jun. de 2024 a las 7:53
"Does Matlab's integral function pass the integration variable to the function handle one scalar at a time or as a vector?"
By default it assumes that the function is vectorized, i.e. accepts a vector and returns a vector of the same size. This is explained in the documentation: "for scalar-valued problems, the function y = fun(x) must accept a vector argument, x, and return a vector result, y."
If your function only accepts a scalar, then set the 'ArrayValued' option to true:
q = integral(@(t) f(t,tau), 0,T, 'ArrayValued',true)
3 comentarios
John D'Errico
el 24 de Jun. de 2024 a las 9:24
All that has been said is true. HOWEVER, I would add a caveat. If you are having problems with the integral, then there is some potential the problem lies in your function g, and in what you do. The use of mod makes your function highly non-smooth.
As such, I would strongly recommend you plot the function you are trying to integrate here. I would do so, but you don't supply g, so we cannot help you more.
Más respuestas (0)
Ver también
Categorías
Más información sobre Numerical Integration and Differentiation 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!