Bonjour,

E1 <= 0.01 pour lm=8 sinon , lm= 8*1.1 s non , lm =8*(1.1^2), sinon lm=8*1.1^3 ...etc


Quelqu'un peut-il m'aider. J'ai essayé avec ce programme :


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
donne = [ 7.92    7.85    7.97
    0.001562 0.00194 0.00482
    561      310     78
    400      300     150
    600      400     200
    150      100     50
    0.00003   0.00009  0.00012 ];
 
b= donne (1,: );
c= donne (2,: );
a= donne (3,: );
p0= donne (4,: );
pmax= donne (5,: );
pmin= donne (6,: );
pd=850;
g=donne (7,: );
n=3;
i=1:n;
lm0=8;
%puissance donne pd=(850)
 
 
%resultat
 
 
 
for i=1:3;
 
    for lm=lm0
 
        f(i)= (2*c(i).*p0(i)+ b(i))-(lm*(1-(2*g(i).*p0(i))))
 
        z(i)=2*c(i)+lm*2*g(i);
 
        PL0 = sum ((g(i)).* (p0(i).^2));
 
        p(i)=p0(i)-(f(i)).*(inv(z(i)))
 
        E1= pd- sum(p(i))+PL0
 
        if abs(E1)>=0.01
 
            p0(i)=p(i)
 
        end
 
    end
end