• API
• FAQ
• Tools
• Archive
daily pastebin goal
76%
SHARE
TWEET

# Untitled

a guest Nov 21st, 2017 60 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. B = 0.01;
2. V = 380;
3. X = 15;
4. Inom= 20;
5. w_nom = 3600*pi/30;
6.
7. %% A vazio
8.
9. Eo = V;
10. Io = 0;
11. ko = Eo/w_nom;
12. tetao = 0;
13. phio = acos(0);
14.
15.
17. %|V| = |E|
18.
19. teta_nom = asin((B*X*w_nom^2)/(3*V^2));
20. I_wn = (380 - (ko*w_nom)*exp(j*teta_nom))/(j*X);
21. phi_nom = acos (V*sin(teta_nom)/(X*I_wn));
22.
23. %%
24. %encontrar w_max -> Ia <= 1.5Inom
25. % T = Bw = 3Va*Ea*sin(teta)/(X*w)
26. w = [ 370: 0.01 : 900];
27.  for i=1:1:length(w)
28.     teta(i) = asin((B*X*w(i)^2)/(3*(ko*w(i))^2));
29.     I(i) = (380 - (ko*w(i))*exp(j*teta(i)))/(j*X);
30.  end
31.  indice = find (abs(I)>1.5*Inom, 1); %primeiro valor de velocidade maior que a maxima
32.  wmax = w(indice-1);
33.  teta_wmax = teta(indice-1);
34.  I_wmax = I(indice-1);
35.  phi_wmax = acos(ko*wmax*sin(teta_wmax)/(X*I_wmax));
36.
37.
38.
39. %%
40. %try to plot
41.
42. ph1o = V*exp(i*0);
43. ph2o = Eo*exp(i*tetao);
44. ph3o = Io*exp(i*phio);
45. ph4o = i*X*Io;
46.
47. p1 = [0; ph1o];
48. vecs = cumsum([ph2o; ph4o]);
49. vecs = [p1+vecs ; vecs];
50. p3 = [0; ph3o];
51.
52. figure;
53. plot(real(p1), imag(p1), '->');
54. %xlim([-30 30]);
55. %ylim([-30 30]);
56. xlabel('real part');
57. ylabel('imaginary part');
58. grid on;
59. hold on;
60. plot(real(vecs), imag(vecs), '-+');
61. plot(real(p3), imag(p3), '-*');
62. %plot(real(p4), imag(p4), '--');
63.
64. %%
65. ph1nom = V*exp(i*0);                 %Va
66. ph2nom = ko*w_nom*exp(i*teta_nom);      %Ea
67. ph3nom = I_wn*exp(i*phi_nom);       %Ia
68. ph4nom = i*X*I_wn;                 %jX*Ia
69.
70. p1 = [0; ph1nom];
71. vecs = cumsum([ph2nom; ph4nom]);
72. vecs = [0 ; vecs];
73. p3 = [0; ph3nom];
74. %p4 = [ph2nom; ph4nom];
75.
76. figure;
77. plot(real(p1), imag(p1), '->');
78. %xlim([-30 30]);
79. %ylim([-30 30]);
80. xlabel('real part');
81. ylabel('imaginary part');
82. grid on;
83. hold on;
84. plot(real(vecs), imag(vecs), '-+');
85. plot(real(p3), imag(p3), '-*');
86. %plot(real(p4), imag(p4), '--');
87. %
88. %
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top