Aide pour une optimisation

Version imprimable

bonsoir tout le monde,comment faire pour avoir les valeurs optimisées des variables suivantes:p,L,alp,Ipk,Vpk

et voici le programme avec les valeurs initiales:

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
 
%function A=optim(x);
Vgs=-0.9:0.1:0;
Vds=0:1:6;
 
Vpk=-0.2;
p=[1.561 1.409 5.379];
alp=1.156;
L=0.00739;
Ipk=0.761;
%Vds=0:1:6;
a0=(-p(1)*Vpk)+(p(2)*Vpk.^2)-(p(3)*Vpk.^3);
a1=p(1)-(2*p(2)*Vpk)+(3*p(3)*Vpk.^2);
a2=p(2)-(3*p(3)*Vpk);
a3=p(3);
 
hold on;
for i=1:10,
psi(i)=a0+(a1*Vgs(i))+(a2*Vgs(i)^2)+(a3*Vgs(i)^3);
for j=1:7 , 
    Ids(i,j)=Ipk*(1+tanh(psi(i)))*(1+(Vds(j)*L))*tanh(Vds(j)*alp);
end;
figure(1);
plot(Vds,Ids(i,:));
pause;
end;
hold off;
figure(2);
plot(Vgs,psi);
%%%%%%%%%%%
function G=optimal(L);
G=Ipk*(1+tanh(psi(i)))*(1+(Vds(j)*L))*tanh(Vds(j)*alp);

et merci :)

28/04/2010, 16h13
duf42

Bonjour,

Je comprends pas bien ce que tu veux...

Que représentent ces variables? Quels critères d'optimisation veux-tu utiliser? Que fait (ou est sensé faire) ton programme?

Bonne apm,
Duf

salut à tous,bon j'ai ecrit mon programme d'une autre façon 1 peu clair le voici:

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
35
36
37
38
function diff=opt(x,Vds,y);
A=x(1);
B=x(2);
C=x(3);
diff=A*(1+tanh(psi(i)))*(1+(Vds(j)*B))*tanh(Vds(j)*C) - y(i,j);
 
Vgs=-0.9:0.15:0;
Vds=0:1:6;
 
Vpk=-0.2;
 
a0=(-1.561*Vpk)+(1.409*Vpk.^2)-(5.379*Vpk.^3);
a1=1.561-(2*1.409*Vpk)+(3*5.379*Vpk.^2);
a2=1.409-(3*5.379*Vpk);
a3=5.379;
 
hold on;
for i=1:7,
    psi(i)=a0+(a1*Vgs(i))+(a2*Vgs(i)^2)+(a3*Vgs(i)^3);
for j=1:7 , 
    y(i,j)=0.761*(1+tanh(psi(i)))*(1+(Vds(j)*0.00739))*tanh(Vds(j)*1.156);
end;
figure(1);
plot(Vds,y(i,:))
pause;
end;
 
et dans une autre fenetre sur matlab l'optimisation en utilisant LSQNONLIN:
 
Vds=0:1:6;
 
y(i,j)=0.761*(1+tanh(psi(i)))*(1+(Vds(j)*0.00739))*tanh(rand(size(Vds(j)))*1.156);
 
options = optimset('Largescale','off');
 
X0=[0.761 0.00739 1.156]';
 
x=lsqnonlin(@opt1,X0,[],[],options,Vds(j),y(i,j))

mon but est d'optimiser les valeurs A,B,C par rapport aux valeurs initiales X0 :cry:

12/05/2010, 14h33
Merel

Sinon, ca n'est pas parce que tu nous réécris ton programme d'une autre manière qu'on va pouvoir faire avancer le problème. Si tu ne nous expliques pas clairement ce que tu veux faire (ou ce que ton programme doit faire), ca va être difficile de t'aider. Surtout si ton programme ne renvoit pas d'erreur.

Tu essayes apparemment de faire un régression par la méthode des moindres carré. C'est ca ? Si tu pouvais nous en dire plus...