Discussion :

[raicolin]

Bonjour
Si quelqu'un peut m'aider à implémenter cet algorithmique je suis stagiaire et j'ai pas le temps de lire le langage du début :
On a
un vecteur X à n(fixé) composante connu
on veut résoudre deux equations (en piece jointe ) on doit d'abord resoudre la premiere puis la deuxieme



aprés
on doit calculer le vecteur lambda = (N-(i-1))Φ qui est de longeur n

merci beaucoup

[membreComplexe12]

je pense que si tu vzeux de l'aide il faut que tu fasses au moins quelques efforts en essayant de l'implementer, apres si il y a des erreurs ont pourra te les corriger

A+

[raicolin]

Bonjour
Puisque personne ne veut m'aider, je me lance
donc j'ai fait des recherches sur le net j'ai trouvé la méthode de Newton-Raphson pour résoudre ce genre d'équation
je propose cet algorithme :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
function Y=methnew(solini,erreur)
 
X=0;n=0;sol=0;L=0;
L=input('entrez la taille de X'/n);
for i=1 :L
    r(i)=input('entrez X '/n);
end;
fonc=FonctionN(r,solini); hold on
while abs(fonc)>erreur
    n=n+1; fonc = fonctionN(r,solini);
    dif= FonctionN2(r,solini);
    if dif~=0
        sol=solini -fonc/dif;
    else disp('derivée nulle'); 
    end;
    solini=sol;

toutefois lors de l'execution une erreure "line 4 "end" expected end of input found" que veut dire cela ?

[membreComplexe12]

personnellement j'ai un peu de mal à comprendre le but de ton algo (j'ai l'impression que tu appel des fonctions mais que tu ne nous à pas montré...)

désolé, je ne pense pas pouvoir t'aider je comprends pas vraiment ton probleme...

[raicolin]

FonctionN est la premiere equation d'inconnu N de type f(N)= 0

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
function K = FonctionN(N,X)
A=0;B=0 ;C=0; D=0; E=0;
K=0;
n=length(X);
for i=0:n
    A = A+(1/(N-i+1));
    C = C+1/X(i);
    D = D+((i-1)* X(i));
end
E= C * D ;
B = n/(N-E) ;
K = (A-B) ;
 
 
%return ;

la deuxieme est la dérivée par rapport à N

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
function K = FonctionN2(N,X)
A=0;B=0 ;C=0; D=0; E=0;
K=0;
n=length(X);
for i=0:n
    A = A+(1/(N-i+1).^2);
    C = C+1/X(i);
    D = D+((i-1)* X(i));
end
E= C * D ;
B = n/((N-E)).^2 ;
K = (B-A) ;
 
 
%return ;

l'idée est d'entrer le vecteur x au début aprés apliquer Newton Raphson mais ca marche pas voici les erreurs qu'on me donnent :

Warning: Divide by zero.
> In C:\MATLAB6p1\work\methnew.m at line 6
??? Error using ==> input
Function 'input' is not defined for variables of class 'double'.

Error in ==> C:\MATLAB6p1\work\methnew.m
On line 6  ==>     r(i)=input('entrez X '/n);

[Jerome Briot]

Ton erreur n'a rien à voir avec les mathématiques. C'est une erreur de syntaxe MATLAB. Lis la documentation de la fonction INPUT et tu verras que le "/n" à la fin est faux et inutile.