Discussion :

[mimosa13]

Bonjour
j'ai eu un problème en programmant une fonction sur matlab, je dois faire le graphique et je dois chercher le c_min pour que f>=0

voilà le code

Code :

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

1
2
3
4
5
6
7
alpha=2;
a=0.05;
n=3;
f='(alpha/(alpha-1))*((1-a-((n-1)*c)).^(1-(1/alpha))-(c.^(1-(1/alpha))))-((((1/n)*(1-a))-c)*(((n-1)*(1-a-((n-1)*c))).^(-1/alpha))+(c.^(-1/alpha)))';
c=[0:0.01:((1/n)*(1-a))];
y= eval(f);
plot(c,y);

[Jean Dumoncel]

Bonjour,

tu n'es pas parti dans la bonne direction, il faut que tu apprennes à écrire des fonctions : function et Create Functions in Files et que tu l'appliques à ton problème.

[mimosa13]

Bonsoir
j'ai besoin de votre aide, le code ne marche toujours pas, je n'ai pas réussi à avoir l'output que je cherche à savoir le c_min et le graphique
voilà le code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
function [f,c_min]= myfun(alpha, a, n,c)
 
alpha=2;
a=0.05;
n=3;
f='(alpha/(alpha-1))*((1-a-((n-1)*c)).^(1-(1/alpha))-(c.^(1-(1/alpha))))-((((1/n)*(1-a))-c)*(((n-1)*(1-a-((n-1)*c))).^(-1/alpha))+(c.^(-1/alpha)))';
c=[0:((1/n)*(1-a))];
while f>=0
    c_min=c;
end
y= eval(f);
plot(c,f);grid;

[Jean Dumoncel]

Ou as-tu vu qu'il fallait utiliser eval pour calculer la valeur d'une expression? C'est à éviter absolument! de plus ta boucle tant que génère une boucle infinie.

Je ne suis pas sur de ce que tu veux faire mais voici un exemple avec ton code pour évaluer la fonction f. Tu peux définir la fonction myfun ainsi :

Code :

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

1
2
3
function f= myfun(alpha, a, n,c)

f=(alpha/(alpha-1))*((1-a-((n-1)*c)).^(1-(1/alpha))-(c.^(1-(1/alpha))))-((((1/n)*(1-a))-c).*(((n-1)*(1-a-((n-1)*c))).^(-1/alpha))+(c.^(-1/alpha)));

J'ai ajouté un point en rouge car tu multiplies deux tableaux et j'ai supposé qu'il faut les multiplier terme à terme.
Et pour appeler cette fonction :

Code :

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

1
2
3
4
5
6
alpha=2;
a=0.05;
n=3;
c=0:((1/n)*(1-a));
 
f= myfun(alpha, a, n,c);

Mais je te conseille de revoir le vecteur c car il ne contient qu'un seul terme (0).