un script sur plusieurs fichiers "functions.m"

Version imprimable

1 pièce(s) jointe(s)

bonjour :king:

donc j'ai programmé mon premier script en Matlab 7.1:mouarf:

mais je veut le diviser sur plusier fichiers ( fonctions) *.m:oops:

bon voila mon script
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 26 27 28 29 30 31 32 33 34 % script calculant l"intégrale d'une fonction en utilisant la méthode de % Gauss Legendre disp('ceci est un programme calculant l"intégrale de la fonction') f = inline('exp(-t)','t') p(1,1)=1; p(2,1:2)=[1 0]; disp('pour continuer la calcule de l"intégrale de cette fonction ') n=input('vous devez specifiez le nombre "n" qui est le degrée du poly de Legendre '); for k=2:n p(k+1,1:k+1)=((2*k-1)*[p(k,1:k) 0]-(k-1)*[0 0 p(k-1,1:k-1)])/k; end sprintf('le polynome de legendre d ordre %d est ',n); pp=p(n+1,:); x=roots(pp); pp1=polyder(pp); for i=1:numel(x) w(i)=2/( (1-x(i)^2)*polyval(pp1,x(i))^2 ); end if sum(w)<=2 for i=1:numel(x) ss(i)=w(i)*feval(f,x(i)) end disp('l"intégrale vaut ') sum(ss) else disp('il ya une erreur d"arrondissage et d"approximation pour les wi') return; end
je veut que
le polynome de legndre ( appele pp ( la dernière ligne de la matrice p)) qui'il soit dans un fichier PLG.m

le calcul des x(i) ( les racines de ces polynomes) et les w(i) (les poids associé à ces racines)) dans un autre fichier XWX.m

puis un petit script qui peut lire une fonction 'inline' pour évaluer son intégrale entre -1 et 1 en utilisant les deux fonctions PLG et XWX

Merci d'avance 8-)
j'ai attaché le script en format txt.

18/12/2007, 16h15
Caro-Line
Et quelle est la question au juste ? :koi:
En gros il suffit de copier-coller la partie qui t'intéresses en réfléchissant aux entrées/sorties.
exemple :
Code:

1 2 3 4 5 6 function pp=PLG(n,p) for k=2:n p(k+1,1:k+1)=((2*k-1)*[p(k,1:k) 0]-(k-1)*[0 0 p(k-1,1:k-1)])/k; end sprintf('le polynome de legendre d ordre %d est ',n); pp=p(n+1,:);
Et dans ton code tu remplaces la partie copiée/collée par l'appel à la fonction :

Code:

pp=PLG(n,p);
18/12/2007, 16h20
tug83

euuh c'est quoi que tu n'arrives pas à faire?, il te suffit d'utiliser le mot function
18/12/2007, 16h22
prologO

Citation:

Envoyé par caro95470

Et quelle est la question au juste ? :koi:
En gros il suffit de copier-coller la partie qui t'intéresses en réfléchissant aux entrées/sorties.

eh oui les E/S que je ne comprend pas en matlab+ la fonction feval (apparement je vais l'utiliser :oops:)
18/12/2007, 16h25
prologO

Citation:

Envoyé par tug83

euuh c'est quoi que tu n'arrives pas à faire?, il te suffit d'utiliser le mot function

j'ai essayé mais il y a tjr des erreurs :ouin: :ouin:
18/12/2007, 16h28
Caro-Line

Citation:

Envoyé par prologO

eh oui les E/S que je ne comprend pas en matlab+ la fonction feval (apparement je vais l'utiliser :oops:)

Montre-nous ce que tu as fait jusqu'à présent.
Je ne vois pas bien où est la difficulté des E/S.
Et pourquoi veux-tu utiliser la fonction feval ? Dans ton cas je n'en vois pas l'intérêt.
Commence déjà par faire une fonction et on te dira ce qui ne va pas.
courage....;)

Code:

1
2
3
4
5
6
7
8
9
10
11
12
function [pp1,x]=plg(n)
p(1,1)=1;
p(2,1:2)=[1 0]; 
 
 
for k=2:n
   p(k+1,1:k+1)=((2*k-1)*[p(k,1:k) 0]-(k-1)*[0 0 p(k-1,1:k-1)])/k;
end
sprintf('le polynome de legendre d ordre  %d est ',n);
pp=p(n+1,:);
x=roots(pp);
pp1=polyder(pp);

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
function xwx(n)
[pp1,x]=plg(n)
for i=1:numel(x)
    w(i)=2/( (1-x(i)^2)*polyval(pp1,x(i))^2 );
 end
if sum(w)<=2 
 
for i=1:numel(x)
    ss(i)=w(i)*feval(f,x(i))
end
disp('l"intégrale vaut ') 
sum(ss)
 
else
    disp('il ya une erreur d"arrondissage et d"approximation pour les wi')
    return;
 end

dans cet fonction
il ya une erreur a la ligne 9

Citation:

??? Undefined function or variable 'f'.
Error in ==> xwx at 9
ss(i)=w(i)*feval(f,x(i))

le script principale

Code:

1
2
3
4
5
6
 
disp('pour continuer la calcule de l"intégrale de cette fonction  ')
n=input('vous devez specifiez le nombre "n" qui est le degrée du poly de Legendre ');
f=inline('exp(-t)')
plg(n)
xwx(n)

18/12/2007, 16h49
Caro-Line

Citation:

Envoyé par prologO

dans cet fonction
il ya une erreur a la ligne 9

Citation:

??? Undefined function or variable 'f'.
Error in ==> xwx at 9
ss(i)=w(i)*feval(f,x(i))

Ben oui, c'est qui f ??? 8O
Plus sérieusement, il faut effectivement que tu mettes la variable f en entrée de ta fonction tout simplement.
18/12/2007, 16h51
prologO

Citation:

Envoyé par caro95470

Ben oui, c'est qui f ??? 8O
Plus sérieusement, il faut effectivement que tu mettes la variable f en entrée de ta fonction tout simplement.

mais j'ai déja introduit la fonction 'f 'dans le script principal 8O
(j'ai réctifier mon dernier msg)
18/12/2007, 17h00
Caro-Line

Citation:

Envoyé par prologO

mais j'ai déja introduit la fonction 'f 'dans le script principal 8O
(j'ai réctifier mon dernier msg)

Les fonctions ne connaissent que les variables qui lui sont passées en entrée et celles définies dans le corps même de la fonction. Donc ta fonction ne connait pas f puisque cette variable est définie dans ton script, il faut donc lui donner en entrée.

Je te conseille cette :faq: : Existe-t-il plusieurs types de définitions de variables ?
18/12/2007, 17h27
prologO

Citation:

Envoyé par caro95470

Les fonctions ne connaissent que les variables qui lui sont passées en entrée et celles définies dans le corps même de la fonction. Donc ta fonction ne connait pas f puisque cette variable est définie dans ton script, il faut donc lui donner en entrée.

Je te conseille cette :faq: : Existe-t-il plusieurs types de définitions de variables ?

j'ai diviser mon script en une fonction et un script
merci à tous
18/12/2007, 17h29
christophe_halgand

Dans ton programme principale :

[sortie que tu souhaites séparer par une virgule]=Nom_de_ta_fonction(entrées dont tu as besoin pour calculer tes sorties séparer par une virgule)

Et dans ton "sous-programme"

function [sortie que tu souhaites séparer par une virgule]=Nom_de_ta_fonction(entrées dont tu as besoin pour calculer tes sorties séparer par une virgule)

et voilà, dans ton sous-programme tu copies puis colles ce que tu avais fait dans ton gros programme en t'arrêtant quand les sorties sont calculées...
18/12/2007, 17h37
prologO

Citation:

Envoyé par christophe_halgand

Dans ton programme principale :

[sortie que tu souhaites séparer par une virgule]=Nom_de_ta_fonction(entrées dont tu as besoin pour calculer tes sorties séparer par une virgule)

Et dans ton "sous-programme"

function [sortie que tu souhaites séparer par une virgule]=Nom_de_ta_fonction(entrées dont tu as besoin pour calculer tes sorties séparer par une virgule)

et voilà, dans ton sous-programme tu copies puis colles ce que tu avais fait dans ton gros programme en t'arrêtant quand les sorties sont calculées...

merci christophe ;)