[dblquad] Passage de paramètre variable

Bonsoir,
Je veux savoir comment boucler la fonction fun.m pour différent valeur de P pour calculer le double intégrale Q correspond a chaque valeur de P grâce à la commande matlab de dblquad, puis pour tracer la courbe plot !!!!!

voici le code :
fun.m

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

function z=fun(x,y)     B=3e-6;   k =[1.0000 0 0.2140 0 -0.00623;-0.00623 0 0.0980 0 -0.1035];           W1=1.917877e-8*P; s=0; for i=1:5 for j=1:5 s=s+k(i,j); end end         z=8.86./(B-W1);

---

et

Code:

---

1 2 3

Q=dblquad(@fun,-1.5e-3,1.5e-3);

---

Merci pour toute réponse.
Meilleures salutations.

Bonjour,

>> :faq: Que représente le symbole @?

@ : pour appeler la fonction myfun.m :calim2:
le code global est dans le fichier repcap, le fichier myfun.m est une fonction.
pour calculer le double intégrale en utilisant "dblquad " dans matlab.

Donc mon problème est comment faire calculer le Q (résultat de l’intégrale double) pour différent valeur de P, puisque je veux tracer la courbe de variation de Q en fonction du P.
et merci.

Je te conseille de relire le contenu du lien de la FAQ (et surtout le dernier point dans ton cas: ton P étant a)
Le symbole @ n'est pas fait pour appeler une fonction, mais accéder à son handle.

salut,
Tu peut voire le fichier attachée (page 2, "Then call it in Matlab using @").
P variable compris entre 1 jusqu’au 5 (dans mon cas j'ai fixé P=1).

Dans cette commande :

Code:

---

1 2	Q=dblquad(@fun,-1.5e-3,1.5e-3);

---

j'ai calculé le double intégrale avec a=3e-3

mon problème est : de calculer Q pour différent valeur de P
par exemple : pour
P=1 , Q=26
P=2 , Q=28

P=5 , Q=34
et merci.

Et la solution se trouve dans mon lien :mur:

Code:

---

1 2	P = ... Q=dblquad(@(x,y) myfun(x,y,P) ,-1.5e-3,1.5e-3,-1.5e-3,1.5e-3);

---

Avec:

Code:

---

1 2	function z=myfun(x,y,P) % ...

---

Un effort d'essayer d'appliquer ce qui s'y trouve dedans, est-ce trop demandé? :calim2:
La dernière page de ton pdf te montre même un exemple!

Tu peut m'expliquer encore cette commande

Code:

---

1 2	Q=dblquad(@(x,y) myfun(x,y,P) ,-1.5e-3,1.5e-3,-1.5e-3,1.5e-3);

---

@(x,y) myfun(x,y,P) !!!!!!!!!

merci bien, c'est résolu :ccool: