Sprawozdanie
Laboratorium z metod numerycznych
Ćwiczenie nr 5 : Metoda węzłowa
Wykonał: Kilijanek Łukasz

Schemat

Graf zorientowany

Skrypt scilaba

clear;

xdel;

clc;

f=input("podaj wartosc f: ");

R1=input("podaj wartosc R1: ");

R2=input("podaj wartosc R2: ");

R3=input("podaj wartosc R3: ");

R4=input("podaj wartosc R4: ");

R5=input("podaj wartosc R5: ");

L1=input("podaj wartosc L1: ");

L2=input("podaj wartosc L2: ");

L3=input("podaj wartosc L3: ");

L4=input("podaj wartosc L4: ");

L5=input("podaj wartosc L5: ");

C1=input("podaj wartosc C1: ");

C2=input("podaj wartosc C2: ");

C3=input("podaj wartosc C3: ");

E1=input("podaj wartosc E1: ");

E2=input("podaj wrtosc E2: ");

E8=input("podaj wartosc E8: ");

w=2*f*%pi

XL1=(w*L1)*%i

XL2=w*L2*%i

XL3=w*L3*%i

XL4=w*L4*%i

XL5=w*L5*%i

XC1=-1/w*C1*%i

XC2=-1/w*C2*%i

XC3=-1/w*C3*%i

Z1=R1+XL1+XC1

Z2=0

Z3=R2

Z4=XC2+XL2

Z5=R3

Z6=R4+XL3

Z7=XL4

Z8=0

Z9=R5+XL5

Y1=1/Z1

Y2=1/Z2

Y3=1/Z3

Y4=1/Z4

Y5=1/Z5

Y6=1/Z6

Y7=1/Z7

Y8=1/Z8

Y9=1/Z9

A=[1,1,0,0,0,0,0,0,-1; 0,-1,1,0,0,0,0,-1,0; 0,0,-1,-1,1,0,0,0,0; -1,0,0,1,0,1,0,0,0; 0,0,0,0,0,0,-1,1,1] //macierz incydensji

Yd=diag([Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9]) //macierz diagonalna admitancji gałęziowych

Y=A*Yd*A' //macierz admitancji

Eg=[E1;E2;0;0;0;0;0;E8;0] //wektor gałęziowych napięć źródłowych

Izr=[0;0;0;0;0;0;0;0;0] //wektor gałęziowych prądów źródłowych

V=inv(Y)*A*(Izr-Yd*Eg) // wyznaczamypotencjały węzłowe

U=A'*V //obliczamy wektor napięć gałęziowych

Ig=Yd*(U+Eg) //obliczamy prądy gałęziowe

format("v", 5) //ustalenie precyzji wyświetlania liczb

abs(Ig)

Wyniki obliczeń:

podaj wartosc f: -->50

podaj wartosc R1: -->46

podaj wartosc R2: -->47

podaj wartosc R3: -->48

podaj wartosc R4: -->49

podaj wartosc R5: -->50

podaj wartosc L1: -->0.6

podaj wartosc L2: -->0.65

podaj wartosc L3: -->0.7

podaj wartosc L4: -->0.75

podaj wartosc L5: -->0.8

podaj wartosc C1: -->0.0000012

podaj wartosc C2: -->0.0000024

podaj wartosc C3: -->0.0000036

podaj wartosc E1: -->250

podaj wrtosc E2: -->300

podaj wartosc E8: -->350

w =314.15927

XL1 =188.49556i

XL2 =204.20352i

XL3 = 219.91149i

XL4 =235.61945i

XL5 =251.32741i

XC1 =- 3.820D-09i

XC2 = - 7.639D-09i

XC3 = - 1.146D-08i

Z1 =46. + 188.49556i

Z2 =0.

Z3 =47.

Z4 = 204.20352i

Z5 =48.

Z6 =49. + 219.91149i

Z7 =235.61945i

Z8 = 0.

Z9 = 50. + 251.32741i

Y1 =0.0012219 - 0.0050070i

Y2=1/Z2

!--error 27

division by zero...

at line 40 of exec file called by :

exec("C:/Documents and Settings/łukasz/Pulpit/BezNazwa1.sce");

in execstr instruction called by :

ME~1/łukasz/Pulpit/BEZNAZ~2.SCE'));

Wnioski:

Do danego układu scilab ani nikt inny nie da rady obliczyć prądów gałęziowych gdyż dzielenie przez zero jest niewykonalne, a ma ono miejsce przy obliczaniu Y2 i Y8. Gdy wprowadzimy inne Y2 i Y8 prądów w metodach nie da się porównać.