How to integrate a vector-valued function?
26 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I am trying to calculate the attached integral. I have data as vectors r and f(r) tht I load as r and fr, respectively.. The vector x should have the same range as r. Is the following correct? I'm not sure how I should incorporate the fact that f(r) is a function of r, or how the integral depends on x as well.
I try
x = linspace(0.05,17.95,180)';
fun = @(x,r) (r.^2).*(fr-1).*sin(x.*r)./(x.*r);
eq1 = integral(@(r) fun(x,r),0,r(end),'ArrayValued',1);
eq2 = cumtrapz(r,(r.^2).*(fr-1).*sin(x.*r)./(x.*r));
But I get two different answers. Again, r and fr (= f(r)) are both vectors of numbers. Is either eq1 or eq2 one correct? Is neither? I have limited understanding of anonymous functions and numerical integration in MATLAB.
0 commentaires
Réponse acceptée
Ameer Hamza
le 30 Avr 2020
You need to integrate this function w.r.t. 'r' so you shouldn't use x as a vector in your code. Both of your approaches seem to be doing what you might not intend to do. Check this code, it shows how to do this with integral and trapz. Both will give similar results.
R = linspace(1, 100, 1000); % example r
FR = exp(-R); % example fr
x = linspace(0.05,17.95,180)';
fr = @(r) interp1(R, FR, r, 'pchip');
integrand = @(r,x) r.^2.*(fr(r)-1).*sin(x.*r)./(x.*r);
fx_int = @(x) integral(@(r) integrand(r,x), 0, max(R));
fx_int_v = zeros(size(x));
for i=1:numel(x)
fx_int_v(i) = fx_int(x(i));
end
fx_trapz = @(x) trapz(R, integrand(R, x));
fx_trapz_v = zeros(size(x));
for i=1:numel(x)
fx_trapz_v(i) = fx_trapz(x(i));
end
4 commentaires
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!