Discussion :

[Stalinou]

Bonjour,

Je suis étudiant de mathématiques et je souhaite calculer une intégrale faisant intervenir une indicatrice. Plus exactement, je souhaite intégrer la fonction
f(u,v)=u si v<u et 0 sinon (dans la fenêtre [0,1]x[0,1] par exemple). J'ai d'abord programmé la fonction en faisant une boucle en if. Mais lorsque je cherche à l'intégrer et que j'utilise la commande:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

dblquad(f,0,1,0,1)

, il m'est affiché:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
Error in ==> quad at 71
if ~isfinite(y(7)).

Que dois-je faire?

Merci pour vos réponses!

Stalinou.

[Invité]

Bonjour,

Et si tu commençais par nous montrer ton code?
Et particulièrement celui de ta fonction f?
Tu peux remarquer que la fonction quad utilisée par dblquad fait rentrer au début un vecteur u de 7 valeurs:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
x = [a a+h a+2*h (a+b)/2 b-2*h b-h b];
y = f(x, varargin{:});

et donc voudrait bien obtenir un vecteur de même longueur aussi en sortie (et non un simple 0).
Enfin je pense que le problème vient de là.

[Stalinou]

Merci pour ta réponse. Voici mon programme:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
function z=f(u,v)
if u<v
    z=u
else
    z=0
end

[Invité]

Bonsoir,

C'est bien ce que je pensais, quand u est un vecteur, tu retournes seulement 0.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
function z=f(u,v)
 
z = double(u<v).*u;

ou 'sans utiliser de fonction':

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
f = @(u,v) double(u<v).*u;
dblquad(f,0,1,0,1);

Par contre je remarque que tu as inversé la condition (u<v) par rapport à ton premier post...

PS: Pense à l'indentation de ton code Ctrl+A , Ctrl+I

[Stalinou]

Merci pour ta réponse et tes conseils qui m'ont bien aidé.
A bientôt.
Stalinou.