daily pastebin goal
16%
SHARE
TWEET

Untitled

a guest Nov 21st, 2017 55 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.  
  16. %% velocidade nominal;
  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
Top