Does Matlab's integral function pass the integration variable to the function handle one scalar at a time or as a vector?
32 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
yuval dahan
le 24 Juin 2024 à 7:29
Commenté : yuval dahan
le 24 Juin 2024 à 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 commentaires
Réponse acceptée
Stephen23
le 24 Juin 2024 à 7:42
Modifié(e) : Stephen23
le 24 Juin 2024 à 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 commentaires
John D'Errico
le 24 Juin 2024 à 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.
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Numerical Integration and Differentiation dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!