Intergration d'une fonction à variable muette

Bonjour,

Je suis en train d'essayer d’intégrer une fonction dont la variable est muette mais je me retrouve bloquer. Ma fonction semble marcher mais j'ai comme réponse de Matlab Ans = []. C'est mon premier problème...

Le second est: lorsque je lance mon programme pour l'intégration de cette fonction ( dans un autre fichier donc) Matlab me renvoie cette erreur:

"??? Attempted to access y(1); index out of
bounds because numel(y)=0.
Error in ==> quadl at 74
if ~isfinite(y(1))
Error in ==> Temperature at 10
 Tac=quadl(@(x) Fction,0,5);"

:weird:


Code de la fonction :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

function T=Fction(D,f,K,Po,z,r) N=350; set(0,'RecursionLimit',N) persistent x D=0.1*10^(-6); f=0.05; w=2*pi*f; K=0.3; Po=10^12; Sigma=sqrt(i*(w/D)); z=0; r =-0.0024:0.0001:0.0024; for n=1:1:length(r); T=(Po/(4*pi*K))*((exp((x^2+Sigma^2)*z)/(x^2+Sigma^2)).*((besselj(0,x*r(n)) end

---

Dans mon autre fichier j'ai:

Code:

---

1 2 3 4 5 6 7

a1=0; b=100; r =-0.0024:0.0001:0.0024; Tac=quadl(@(x) Fction,0,5); figure(1) plot(r,log(abs(Tac.*r)),'-k',r,unwrap(angle(Tac)),'-r')

---

Pour information, je n'ai aucune idée de la valeur du x, cette variable me sers simplement à intégrer. De plus j'ai essayé d'intégrer grâce aux séries mais matlab ne supporte pas....
Avez vous une solution svp?

Merci d'avance!

Rebonjour,

voila j'ai un peu modifié mon code

Code:

---

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

function T=Fction(D,f,K,Po,z,r) N=360; set(0,'RecursionLimit',N) persistent x D=0.1*10^(-6); f=0.05; w=2*pi*f; K=0.3; Po=10^12; Sigma=sqrt(i*(w/D)); z=0; r =-0.0024:0.0001:0.0024; % A=repmat(r1,49,1); % r=A.*eye(49) for n=1:1:length(r); T=@(x)(Po/(4*pi*K))*((exp((x^2+Sigma^2)*z)/(x^2+Sigma^2)).*((besselj(1,x*r(n))))*x); end Tac=quadl(Fction,0,5); figure(1) plot(r,log(abs(T.*r)),'-k',r,unwrap(angle(T)),'-r')

---

Message d'erreur:

??? Maximum recursion limit of 360 reached. Use
set(0,'RecursionLimit',N)
to change the limit.  Be aware that exceeding
your available stack space can
crash MATLAB and/or your computer.

Error in ==>
Fction>create@(x)(Po/(4*pi*K))*((exp((x^2+Sigma^2)*z)/(x^2+Sigma^2)).*((besselj(1,x*r(n))))*x)

(j'ai essayé d'autres techniques entre temps mais rien de marche)
Maintenant ce code me dit qu'il faut plus d'itération donc j'ai augmenté mon N mais au dessus de 360 Matlab crash (version 2008a)
Je pense ne plus être loin de la solution mais je suis bloqué...
Auriez vous des idées svp?
Merci d'avance

Ta Fction est infiniment récursive puisqu'elle s'appelle elle-même à chaque fois. N'est-ce pas plutôt @T que tu voulais intrégrer ??

Merci pour ta réponse.
En effet merci c'est plutôt cette fonction que je voudrais intégrer.
Mais dois-je le faire dans un autre fichier ou dans le même Matlab peut accepter?

J'ai donc fais des tests il semble plus logique de le faire dans un autre fichier (je suis pas forcement un bon utilisateur de Matlab ça se voit...)

voici mon nouveau script:

Code:

---

1 2 3 4 5 6

K=0.3; Po=10^12; a1=0; b=100; r =-0.0024:0.0001:0.0024; Tac=quadl(@T,0,100)

---

Mais voici le message d'erreur:

??? Error using ==> feval
Undefined function or method 'T' for input
arguments of type 'double'.

Error in ==> quadl at 70
y = feval(f,x,varargin{:}); y = y(:).';

Error in ==> Temperature at 9
    Tac=quadl(@T,0,100)

J'avoue ne pas forcement comprendre comment T est reconnu par Matlab... En fait je cherche à intégrer T entre 0 et +infini et je pense que quadl est le meilleure fonction pour ça. (J'ai essayé de passer par les séries mais la fonction trouvée n'est pas la bonne)
Je continue mes recherches mais avez-vous des commentaires sur le code?
Merci d'avance

Citation:

---

Undefined function or method 'T' for input

---

Refais le tout dans un même fichier, la ou est définie T, sinon Matlab ne la trouvera évidement pas. Remplace juste dans ton premier exemple

Code:

---

1 2 3 4

Tac=quadl(Fction,0,5); % par Tac=quadl(T,0,5);

---

Merci pour l'aide,

Nouveau message:

??? Error using ==> mpower
Matrix must be square.

Error in ==>
Fction>@(x)(Po/(4*pi*K))*((exp((x^2+Sigma^2)*z)/(x^2+Sigma^2)).*((besselj(1,x*r(n))))*x)
at 21
    T=@(x)(Po/(4*pi*K))*((exp((x^2+Sigma^2)*z)/(x^2+Sigma^2)).*((besselj(1,x*r(n))))*x);
    
Error in ==> quadl at 70
y = feval(f,x,varargin{:}); y = y(:).';

Error in ==> Fction at 31
Tac=quadl(T,0,100);

Je comprends pas pourquoi il pense que c'est des matrices...
J'ai fait une boucle pour passer r de vecteur lignes à matrice alors j'ai voulu vérifier certaine chose:
size Sigma = 1 5 :weird:
Alors qu'à la vue c'est un scalaire j'ai vérifié: (il me donne un message d'erreur alors que ce matin il me donnait un imaginaire pur.... )

??? Error using ==> sigma at 83
Not enough input arguments.

Bref, je suis désolé de solliciter pour ce qui pourrait être des broutilles mais je suis un peu dans une impasse depuis un moment....

Ce ne sont pas des matrices mais des vecteurs. Il faut utiliser .^2 (carré terme à terme) et non ^2 (carré matriciel, qui lui nécessite des matrices carrés)

Excuse moi, je voulais dire de vecteur ligne à scalaire...
Mais par définition ça devrait être des scalaires.. Sigma devrait avoir une seule composante (complexe)

L'erreur sur sigma vient sûrement du fait que la variable n'était pas définie, et c'est la fonction MATLAB du même nom (qui elle requiert ou moins un paramètre d'entrée) qui était utilisée.
Sinon aucun problème à ce que sigma soit un scalaire je te rassure. C'est juste x qui élevé au carré sans le point pour faire du membre à membre qui posait problème ;)

Merci de ta réponse.
J'ai donc changé le nom de la variable en Sigma1 (ce que j'ai compris c'est que Sigma est déjà utilisée par Matlab) Mais toujours le même message d'erreur.
Mais maintenant le size donne : 1 6 (6 est le nombre de caractères de Sigma1 et 5 celui de Sigma (comme ci dessus) Coïncidence?)
Et il ne reconnait toujours pas la variable....

Citation:

---

Envoyé par Benmu

Mais toujours le même message d'erreur.

---

Lequel ?
Pourrais-tu nous montrer la nouvelle version de ton code ?

Attention que la syntaxe size Sigma correspond à size('Sigma'). Il faut le mettre entre parenthèses :

Code:

---

size(Sigma)

---

Voir la FAQ Qu'est-ce qu'un appel fonctionnel ?

Ah ouai merci pour le size! Mais il ne trouve pas la variable par contre...

Citation:

---

Envoyé par Winjerome

Lequel ?
Pourrais-tu nous montrer la nouvelle version de ton code ?

---

Le message d'erreur est:

??? Error using ==> mpower
Matrix must be square.

Error in ==>
Fction>@(x)(Po/(4*pi*K)).*((exp((x.^2+Sigma1.^2)*z)./(x^2+Sigma1.^2)).*((besselj(0,x*r(n)))).*x)
at 21
    T=@(x)(Po/(4*pi*K)).*((exp((x.^2+Sigma1.^2)*z)./(x^2+Sigma1.^2)).*((besselj(0,x*r(n)))).*x);
    
Error in ==> quadl at 70
y = feval(f,x,varargin{:}); y = y(:).';

Error in ==> Fction at 31
Tac=quadl(T,0,100);

Voila mon code:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

function T=Fction(D,f,K,Po,z,r) N=360; set(0,'RecursionLimit',N) persistent x D=0.1*10^(-6); f=0.05; w=2*pi*f; K=0.3; Po=10^12; Sigma1=sqrt(i*(w/D)); z=0; r =-0.0024:0.0001:0.0024; for n=1:1:length(r); T=@(x)(Po/(4*pi*K)).*((exp((x.^2+Sigma1.^2)*z)./(x^2+Sigma1.^2)).*((besselj(0,x*r(n)))).*x); end Tac=quadl(T,0,100); figure(1) plot(r,log(abs(T.*r)),'-k',r,unwrap(angle(T)),'-r') end

---

J'utilise la boucle for pour passer le "r" de matrice à scalaire mais je ne comprends pas Matlab ne reconnait pas les variables quand je les rentre dans la console pour en voir les valeurs..

Il reste du x^2. Je te conseille aussi de changer les * et / en .* et ./ ; pour les mêmes raisons.

Il te manque un dernier point :

Code:

---

T=@(x)(Po/(4*pi*K)).*((exp((x.^2+Sigma1.^2)*z)./(x^2+Sigma1.^2)).*((besselj(0,x*r(n)))).*x);

---

Par contre de cette boucle, tu ne gardes que la valeurs de la dernière itération, je doute que ce soit ce que tu veuilles. :?
Il faudrait peut-être voir à inclure l'appel à quadl dans la boucle.

À noter que ces lignes

Code:

---

1 2	N=360; set(0,'RecursionLimit',N)

---

ne sont plus nécessaires.

Nickel merci!

Ca marche même si j'obtiens pas ce aue je pensais avoir.... Je vais voir un n-ieme fois si il n'y a pas d'erreur entre ma formule et mon problème mais au niveau code le problème est résolu!

Merci encore.