Discussion :

[saad.pizzicato]

Bonjour,
J'utilise l'algorithme 'levenberg-marquardt' pour la résolution d'un système d'équations non-linéaires...
J'ai déjà pu avoir de bons résultats avec le même programme, mais à partir d'un certain moment j'ai ce message:

No solution found.

fsolve stopped because the problem appears to be locally singular.

<stopping criteria details>


fsolve stopped because the regularization parameter, 5.000000e+016, exceeds 1e16.

Voici la commande que j'utilise pour Fsolve:

Code :

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

1
2
3
 
 
a= fsolve(@myfun,a0,optimset('Display','iter','TolFun',1e-25,'MaxFunEvals',10000000,'MaxIter',10000000,'Algorithm',{'levenberg-marquardt',.005}, 'TolX',1e-25, 'TolFun',1e-25));

C'est quoi ce paramètre de régulation? et comment changer sa valeur de manière à ne plus avoir cette erreur?

Merci à l'avance.

[Jerome Briot]

Montre nous le code de myfun

[saad.pizzicato]

Le voilà...

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
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
 
 
function f = myfun(a)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
 
load('M');
load('K');
load('B');
 
 
p=9;
 
Amp=0.05;
 
X(1)=Amp;
X(2)=a(2);
X(3)=a(3);
X(4)=a(4);
X(5)=a(5);
X(6)=a(6);
X(7)=a(7);
X(8)=a(8);
X(9)=a(9);
 
 
 
 
 
 
Provisoire=zeros(p,p,p,p);
 
for i=1:p
    for j=1:p
        for k=1:p
            for l=1:p
                Provisoire(i,j,k,l)=X(k)*X(l)*B(i,j,k,l);
            end
        end
    end
end
 
S1=sum(Provisoire,4);
BA=sum(S1,3);
 
 
Cool1= (X*K*X')+ (X*BA*X');
Cool2= X*M*X';
Omeganonlin2=Cool1/Cool2;
 
 
 
SYSTEM=(2*K*X')+(3*BA*X')-(2*Omeganonlin2*M*X');
 
SYSTEM(1)=[];
 
f=SYSTEM;
 
end

[tachmou]

Salut,

Désolé, mais là, je sèche... Il me semble que tu as un problème de précision... Il t'affirme que tu dépasses les 10^16 avec un paramètre de "régularisation". Je ne sais pas à quoi il correspond, mais du coup, en regardant la précision que tu cherches, tu demandes du 10^-25 : je ne sais pas si matlab est capable d'aller jusque là. Essaye avec une précision moins exigeante peut-être ?

Sinon, go potasser la doc

Bon courage

A+

[le fab]

tout comme tachmou, essaye des options moins exigeantes (1e-25 !!)

essayes également d'autres valeurs de λ (0.005) dans ton algo de levenberg-marquardt car a priori vu ton message d'erreur le problème se situe là :

Code :

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

1
2
Trust region radius became too small (trust-region-dogleg algorithm) or
 regularization parameter became too large (levenberg-marquardt algorithm).

ou essayes l'autres algo (trust-region-reflective)