formule de Neumann

[invitedce2dcd7]

bonjour tt le monde,
je veux programmer la formule de Neumann pour calculer une inductance mutuelle entre 2 bobine. voici la formule Sans titre.png

et voici mon programme Sans titre2.png. le problème c'est qu'il ne reconnait pas v comme une variable dans l’intégrale. voici le message d'erreur. xyjrthr.png. merci d'avance.
-----

[Antoane]

Bonjour,

Comment Matlab saurait-il que v est un sym dans :
y = sin(v) ;

Une manière facile de calculer l'intégrale est de définir un vecteur v dont les valeurs sont réparties entre 0 et π/2, à calculer l'intégrande sur ces points de v, puis à estimer l'intégrale par application de la méthode des rectangles ou des trapèzes par exemple.
Exemple :

nPts = 1e3 ;
v = linspace(0, pi/2, nPts)
y = (2*sin(v).^2-1)./sqrt(1 - b2*sin(v).^2) ;
dv = diff(v) ;
dv = [dv, dv(end)] ;
integ = sum(y.*dv) ;

Améliorable en choisissant intelligemment les points de v, le calcul du dv, la sommation (e.g. https://fr.mathworks.com/help/matlab/ref/trapz.html)

[invitedce2dcd7]

merci beaucoup, ya des modifs a faire mais ça correspond a ce que je cherche.👍

[Antoane]

Bonjour,

est-ce que ça répond à la partie codage de la question que tu as posée en physique et sur dvp.net ?

[A voir en vidéo sur Futura]

[invitedce2dcd7]

oui, en physique il m'ont donné la direction mathematique a suivre, pour ce que est du code, il m'ont conseillé un forum matlab, c'est ce que je vien de faire et je n'etait pas deçu. merci encore.

[invitedce2dcd7]

bonjour,
je vais vous déranger encore une foi. voila, j'ai essayé de programmer une autre formule inductance.png avec inductance2.png. voici mon programme programme.png. le problème c'est qu'il prend une seule valeur de M (fait pas la boucle), merci.

[Antoane]

Bonsoir,

Il serait plus simple que tu postes le code qu'une copie d'écran.
Le problème est qu'ici, MM est un scalaire dont la valeur change à chaque itération : pour chaque nouvelle valeur du couple (i,j), la ligne "MM=((..." écrase la valeur précédemment calculée.
Finalement, le code est équivalent à celui dans lequel :

"for i=1:N0" est remplacé par "for i=N0"

"for j=1:Nd" est remplacé par "for j=Nd

Suivant ce que tu veux faire (je n'ai pas étudié en détail les calculs), il peut être nécessaire d'écrire "MM(i,j) = ((...", ou "M= M+MM" à chaque itération.

Par ailleurs :

- pourquoi utiliser la fonction @round() ?
- est-il pertinent d'augmenter phi0 et phid à chaque itération ? ne faudrait-il pas remettre phid à zéro "de temps en temps" ?

Plus généralement : je te conseille de commencer par implémenter des intégrales simples et calculables analytiquement (par exemple calculer l'intégrale entre 0 et 10 de x->2*x², puis une fonction de deux variables) pour prendre en main la technique en pouvant vérifier les résultats obtenus.

[JPL]

Il faut encadrer le code par les balises [Code]...[/Code] pour que l’indentation soit bien respectée.

[invitedce2dcd7]

bonjour,
premièrement, j'ai utilisé la fonction @round() car j'ai des problèmes au niveau des sin() et cos() (des valeur de l'ordre de e-5 par foi de e-13).
deuxièmement, phi0 et phid sont des angles d’intégration de 0 a 2*pi (tout le tour de la bobine) avec un pas de calcul de delta_phi0 et delta_phid d'où le delta_phi0=2*pi/N0 (N0 est le nombre de subdivision de la bobine). je sais que mon programme n'est pas totalement faux, il manque juste quelques astuces "Matlabienes" pour que ça fonctionne.

[Antoane]

Bonjour,

Tout n'est effectivement pas à jeter. Le programme n'est cependant pas (seulement) faux car tu manques de compétences en Matlab, il est faux d'un point de vue algorithmique.
C'est pourquoi je te conseille de commencer par coder des intégrales simples (i.e. d'une seule variable et calculables analytiquement, pour vérifier).

premièrement, j'ai utilisé la fonction @round() car j'ai des problèmes au niveau des sin() et cos() (des valeur de l'ordre de e-5 par foi de e-13).

Il y a peut-être une justification autre que tu ne donnes pas, mais ta réponse semble très ad-hoc ! pourquoi ne pas avoir écrit sin(...)+pi/4 au lieu de round(sin(...)) si ce zéro te gène tant ? les deux propositions n'ont aucun sens mathématique.