bonjours tout le monde la j'ai un code en matlab qui calcul les valeurs g,I1,P2,c2,ni,cfi et me donne les graphs : plot(g,P2);plot(g,c2);plot(g,ni);plot(g,cfi)
en fait ce code n'affiche pas les valeurs de I1 ,P2 ,C2 ,ni , cfi donc j'ai ajouter ce code mais il m'affiche plus de 2000 valeurs et moi j'en veux seulement 20!!!!!!
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
60
61 r1=3.5 r2=2.2 rm=11.7 x1=2.41 x2=3.7 xm=46.22 u1=380 p=2 pmec=20 %pi=4*atan(1) j=1 for g=0.000001:0.0005:1 %j=1:1:1000 %t(j)=g; zm=rm+xm*i; z1=r1+x1*i; z2g(j)=r2/g+x2*i; zeq(j)=((zm*z2g(j))/(zm+z2g(j)))+z1; i1(j)=u1/zeq(j); im(j)=(u1-(i1(j)*z1))/zm; i2(j)=((z1*i1(j))-u1)/z2g(j); I1(j)=abs(i1(j)); I2(j)=abs(i2(j)); Im(j)=abs(im(j)); %e1(j)=z1*i1(j)-u1; %e2(j)=(i2(j)*z2g(j)); omega1=(2*pi*50)/p; n1=60*50/p; %n(j)=(1-g)*n1; omega(j)=(1-g)*omega1; %c(j)=3*p*r2*(u1^2)/(2*pi*50*g*(((r1+r2/g)^2)+((x1+x2)^2))); %p2(j)=c(j)*omega(j); per1(j)=3*r1*(I1(j)^2); per2(j)=3*r2*(I2(j)^2); perfe(j)=3*rm*(Im(j)^2); pem(j)=3*r2*(I2(j)^2)/g; pmec1(j)=pmec*(1-g); P1(j)=per1(j)+perfe(j)+pmec1(j); P2(j)=pem(j)-pmec1(j)-per2(j); c(j)=pmec1(j)/omega1; c2(j)=P2(j)/omega(j); %p1(j)=p2(j)+per1(j)+per2(j)+perfe(j); %pp(j)=p2(j)/p1(j); cfi(j)=P1(j)/(3*u1*I1(j)); ni(j)=P2(j)/P1(j); j=j+1; end g=0.000001:0.0005:1; figure plot(g,I1); I1 figure plot(g,P2); figure plot(g,c2); figure plot(g,ni); figure plot(g,cfi); figure plot(g,c2,'b',g,I1,'r');
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 format short r1=2.2 r2=1.1 rm=1.3 x1=2 x2=0.9 xm=42 u1=220 pmec=30 p=2 pi=3.141592654 for k=1:20 g=[0.0000001 0.015 0.025 0.035 0.045 0.055 0.065 0.075 0.085 0.095 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.9999999]; zm=rm+xm*i; z1=r1+x1*i; z2g(k)=r2/g(k)+x2*i; zeq(k)=((zm*z2g(k))/(zm+z2g(k)))+z1; I1(k)=u1/zeq(k); Im(k)=(u1-(I1(k)*z1))/zm; I2(k)=((z1*I1(k))-u1)/z2g(k); per1(k)=3*r1*(I1(k)^2); per2(k)=3*r2*(I2(k)^2); perfe(k)=3*rm*(Im(k)^2); pem(k)=3*r2*(I2(k)^2)/g(k); pmecG(k)=pmec*(1-g(k)); C(k)=3*p*r2*(u1^2)/(2*pi*50*g(k)*(((r1+r2/g(k))^2)+((x1+x2)^2))); omega1=(2*pi*50)/p; omega(k)=(1-g(k))*omega1; C(k)=pem(k)/omega1; P2(k)=pem(k)*(1-g(k))-pmecG(k); P1(k)=P2(k)+per1(k)+per2(k)+perfe(k); C2(k)=P2(k)/omega(k) End
Partager