Discussion :

[Léti117]

Bonjour tout le monde (ou Bonsoir pour certain)
Voilà encore un pb sous matlab
Alors encore une fois je sollicite votre aide!
Voici le problème:

Je voudrais intégrer (dans une fonction, fichier point m) une expression symbolique et avoir en sortie une valeur exacte, voici mon programme

Je crée la fonctionattention l'expression est un peu longue )

Code :

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

1
2
3
4
5
6
function [gg1] = integ1(T,t)
syms z1
g1=(1/2)*(d*((-4)*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*c1*a+sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*b^2+2*sqrt(-4*c1*a+b^2)*a*cond_fin1*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))))/(a*(sqrt(-4*c1*a+b^2)*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+2*a*cond_fin1*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+b*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))));
int(g1, z1, T, t);
 
end

En vérité j'ai un programme principale et dans se programme je fais appelle à cette fonction donc à part la variable z1, tout est identifié et j'ai toujours ce message d'erreur:

??? Undefined function or variable 'd'.

Error in ==> integ1 at 3

sachant que la variable d est bien définie (elle a une valeur fixe) et idem pour les autres variables.

Comment faire SVP?

[Léti117]

ou plutôt voici une autre version de la fonction:

Code :

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

1
2
3
4
5
6
function [gg1] = integ1(T,t)
syms z1
g1=(1/2)*(d*((-4)*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*c1*a+sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*b^2+2*sqrt(-4*c1*a+b^2)*a*cond_fin1*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))))/(a*(sqrt(-4*c1*a+b^2)*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+2*a*cond_fin1*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+b*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))));
gg1=double(int(g1, z1, T, t));
 
end

[Jerome Briot]

Il faut passer toutes les variables nécessaires en argument d'entrée de la fonction :

Code :

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

1
2
3
4
5
6
7
function gg1 = integ1(T,t,a,b,c1,d)

syms z1
g1=(1/2)*(d*((-4)*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*c1*a+sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))*b^2+2*sqrt(-4*c1*a+b^2)*a*cond_fin1*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))))/(a*(sqrt(-4*c1*a+b^2)*cosh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+2*a*cond_fin1*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))+b*sinh(-0.5*z1*sqrt(-4*c1*a+b^2)+0.5*T*sqrt(-4*c1*a+b^2))));
gg1=double(int(g1, z1, T, t));

end