Discussion :

[wolf11]

Bonjour à tous,
j'aimerais modéliser la variation du flux magnétique autour d'une bobine en matlab et, ayant trouvé l'équation de celui-ci, j'ai essayé de résoudre mon intégrale double avec dblquad mais celà ne fonctionne pas...
on m'a parlé de méthode de trapèzes et autres mais je n'y connais absolument rien.. Ce qu'il me faut, c'est juste l'expression du graphique obtenu par la résolution de l'intégrale que voici entre 0 et 2*pi pour les x et entre 0 et entre 0 et 0.0038 pour les y:
f = (0.005-y*cos(x))/((0.005^2+y^2-2*0.005*y*cos(x))^(3/2))y*(2*pi)(2*pi);

En espérant que quelqu'un puisse m'aider dans la réalisation du code. Merci d'avance.

[Invité]

Bonjour,

j'ai essayé de résoudre mon intégrale double avec dblquad mais celà ne fonctionne pas...

Qu'est-ce qui ne fonctionne pas?
Obtiens-tu une erreur?
Montre-nous ce que tu as essayé de coder, même si c'est faux.

Ensuite,

Ce qu'il me faut, c'est juste l'expression du graphique obtenu par la résolution de l'intégrale que voici entre 0 et 2*pi pour les x et entre 0 et entre 0 et 0.0038 pour les y:

La résolution de cette intégrale double te donnera UNE valeur, je ne vois pas ce que tu entends par "l'expression du graphique obtenu" ...

[wolf11]

Effectivement, je me suis mal exprimé. Ce que j'aimerais obtenir en fin de compte,c'est l'expression de la fonction f intégrée selon x allant de 0 à 2*pi. Et donc obtenir une nouvelle fonction variant en fonction de y uniquement..
Voilà ma fonction de départ:
f = (0.005-y*cos(x))/((0.005^2+y^2-2*0.005*y*cos(x))^(3/2))

Le problème, je n'ai encore rien essayé de faire en matlab (mis à part utiliser la fonction 'help'), est que je ne trouve pas de quelle manière, quelle commande utiliser pour trouver l'expression de l'intégrale de f qui ne sera pas un nombre mais une nouvelle fonction de y cette fois (je ne sais pas si je suis très clair...:s).

[Invité]

Dans ce cas, tu peux regarder du côté de la fonction

Code :

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

doc trapz

Qui utilise justement la méthode des trapèzes dont on t'a parlé.

[wolf11]

ok merci beaucoup, je vais regarder à cela.

[wolf11]

je ne comprends pas très bien comment la fonction 'trapz' pourrait me renvoyer une fonction dépendante de y dans ma formule...
Voilà ce que j'ai essayé et le résultat que j'ai obtenu:

Code :

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

1
2
X=0:2*pi/1000:2*pi;
y=(0.005-y*cos(X))/((0.005^2+y^2-2*0.005*y*cos(X))^(3/2));

??? Error using ==> mtimes
Inner matrix dimensions must agree.

Je me doutais bien que ce ne serait pas une bête application de la fonction mais je suis vraiment perdu...