Advertisement
Guest User

Untitled

a guest
May 1st, 2016
51
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.87 KB | None | 0 0
  1. %obwód 1
  2. clear;
  3. clc;
  4. close all;
  5. f=3000; %3000/5000/7000/9000/11000
  6. w=2*pi*f;
  7. r11=1000;
  8. r12=1000;
  9. Rw=50;
  10. c11=100*10^(-9);
  11. c12=47*10^(-9);
  12. c13=22*10^(-9);
  13. l11=0.01;
  14. E=0.5;
  15. C11=1/(j*w*c11);
  16. C12=1/(j*w*c12);
  17. C13=1/(j*w*c13);
  18. L11=j*w*l11;
  19.  
  20.  
  21. A=[(1/Rw)+(1/r11)+(1/C12),-(1/r11),-(1/C12),0;
  22.     -(1/r11),(1/r11)+(1/C11),0,-(1/C11);
  23.     -(1/C12),0,(1/r12)+(1/L11)+(1/C12),-(1/r12);
  24.     0,-(1/C11),-(1/r12),(1/r12)+(1/C11)+(1/C13)];
  25. %macierz wynikająca z zastosowania metody potencjałów węzłowych
  26.  
  27. W=[(E/Rw);0;0;0]; %prawa strona układu równań
  28.  
  29. X=inv(A)*W; %obliczanie potencjałów V1/V2/V3/V4
  30.  
  31. M=abs(X), %moduły
  32.  
  33. K=angle(X).*180/pi, %przesunięcia fazowe
  34.  
  35.  
  36. %Prądy
  37. pradzr=(E-X(1))/Rw; %E-V1/Rw
  38.  
  39. mpradzr=abs(pradzr),
  40.  
  41. fpradzr=angle(pradzr)*180/pi,
  42.  
  43. pradc12=(X(1)-X(3))/C12; %(V1-V3)*jwC12
  44.  
  45. mpradc12=abs(pradc12),
  46. fpradc12=angle(pradc12)*180/pi,
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement