Sprawozdanie
Laboratorium z metod numerycznych
Ćwiczenie nr 4 : Metoda oczkowa
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

B=[1,-1,-1,1,0,0,0,0,0; 0,0,0,-1,-1,1,0,0,0; 0,0,1,0,1,0,1,1,0; 0,1,0,0,0,0,0,-1,1] //macierz incydencji

Zd=diag([Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8,Z9]) //macierz diagonalna impedancjji gałęziowych

Z=B*Zd*B' //macierz impedancji

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

Ioczkowe=inv(Z)*B*(Eg-Zd*Izr) //prąd oczkowy

Ig=B'*Ioczkowe+Izr //wektor prądów gałęziowych

abs(Ig) // moduł 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

B =

1. - 1. - 1. 1. 0. 0. 0. 0. 0.

0. 0. 0. - 1. - 1. 1. 0. 0. 0.

0. 0. 1. 0. 1. 0. 1. 1. 0.

0. 1. 0. 0. 0. 0. 0. - 1. 1.

Zd =

column 1 to 5

46. + 188.49556i 0 0 0 0

0 0 0 0 0

0 0 47. 0 0

0 0 0 204.20352i 0

0 0 0 0 48.

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

column 6 to 9

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

49. + 219.91149i 0 0 0

0 235.61945i 0 0

0 0 0 0

0 0 0 50. + 251.32741i

Z =

column 1 to 3

93. + 392.69908i - 204.20352i - 47.

- 204.20352i 97. + 424.11501i - 48.

- 47. - 48. 95. + 235.61945i

0 0 0

column 4

0

0

0

50. + 251.32741i

Eg =

250.

300.

0.

0.

0.

0.

0.

350.

0.

Izr =

0.

0.

0.

0.

0.

0.

0.

0.

0.

Ioczkowe =

- 0.2505661 - 0.0526897i

- 0.2279284 - 0.1284495i

0.4501853 - 1.2075201i

- 0.0380718 + 0.1913695i

Ig =

- 0.2505661 - 0.0526897i

0.2124943 + 0.2440593i

0.7007514 - 1.1548304i

- 0.0226377 + 0.0757598i

0.6781137 - 1.0790706i

- 0.2279284 - 0.1284495i

0.4501853 - 1.2075201i

0.4882571 - 1.3988897i

- 0.0380718 + 0.1913695i

ans =

0.2560460

0.3236028

1.3508093

0.0790697

1.2744534

0.2616307

1.2887093

1.4816502

0.1951199