Greetings to all, could you help me with this integral please
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Greetings dear friends, I am trying to graph this integral so that I can obtain a graph x =f(t)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/898105/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/898110/image.png)
Thank you!
4 commentaires
John D'Errico
le 17 Fév 2022
Modifié(e) : John D'Errico
le 17 Fév 2022
Plot as a function of WHAT? What will be plotted on each axis? Since you want to integrate over w, that still leaves TWO variables, x and t. And while you seem to think you can plot x on the x axis, and t on the y axis, that makes no sense, since you are asking to plot the INTEGRAL. x = f(t) makes no sense in this context, because there is no functional relationship there, not even an implicit one.
Worse, since there is a complex component in this, what will you plot assuming the result is complex? The real or the imaginary part? And how does the integral factor in?
Réponse acceptée
Torsten
le 17 Fév 2022
Modifié(e) : Torsten
le 17 Fév 2022
t = linspace(0.1,3,29);
x = linspace(-3,3,60);
[T,X] = meshgrid(t,x)
for i = 1:numel(x)
for j = 1:numel(t)
f = @(w) 1/pi*sin(w)./w.*exp(-T(i,j)*w.^2).*cos(w.*X(i,j));
I = integral(f,-Inf,Inf);
F(i,j) = I;
end
end
surf(T,X,F)
You only need to consider the real part of the integrand since the imaginary part is an odd function of w, thus with integral = 0.
4 commentaires
David Goodmanson
le 17 Fév 2022
Hi Luis,
or you could replace the for loops with
F = (1/2)*( erf((X+1)./(2*sqrt(T))) -erf((X-1)./(2*sqrt(T))) );
and get the same result.
Plus de réponses (0)
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!