Laboratorium metod numerycznych |
|||
Laboratorium 6 |
|||
ALBERT GAŁACH |
Grupa ED 3.1 |
Data wykonania ćw. 1.12.2008 12.15-14.00 |
Ocena |
clear;//czyszcenie pamięci
xdel;//czyszczenie okien graficznych
clc;//czyszczenie konsoli
R1=input("podaj wartosc rez 1: ");//podawanie wartości
XC1=input("podaj wartosc indukcji C1: ");
XL1=input("poraj wartość indukcji L1: ");
Z1=(1/(XC1*%i))+(1/(R1))+(1/(XL1))//obliczanie impedancji
Zab=1/Z1
R2=input("podaj wartosc rez 2: ");//podawanie wartości
XL2=input("podaj wartosc indukcji L2: ");
XC2=input("podaj wartosc indukcji C2: ");
R3=input("podaj wartosc rez 3: ");
Zsz=R3+(XC2*%i)//obliczanie impedancji
Z2=(1/(R2))+(1/(XL2*%i))+(1/Zsz)
Zcd=1/Z2
Zglowny=Zab+Zcd
Zabs=abs(Zglowny)
fi_z=atan(imag(Zglowny),real(Zglowny))//obliczanie modułu
Esk=input("definiowanie skutecznej wartosci E: ");//podawanie wartości
fiu=input("definioewanie kąta poczatkowego E: ");
Etl=Esk*cos(fiu)+Esk*%i*sin(fiu)//zamiana na liczbe zespolona
Iglowny=(Esk/Zabs)*cos(fiu-fi_z)+(Esk/Zabs)*sin(fiu-fi_z)*%i
fii=atan(imag(Iglowny),real(Iglowny))//obliczanie modulu
f0=50
Uad=Zglowny*Iglowny//obliczanie napiec i pradow
Uab=Iglowny*Zab
Ic1=Uab/XC1
Ir1=Uab/R1
Il1=Uab/XL1
Ucd=Iglowny*Zcd
Ir2=Ucd/R2
Il2=Ucd/XL2
Isz2=Ucd/Zsz
UXc1=XC1*Ic1
UR1=Ir1*R1
UL1=Il1*XL1
UR2=Ir2*R2
UL2=Il2*XL2
UR3=Isz2*R3
UC2=Isz2*XC2
t=linspace(0,0.1,1000);//cyfry potrzebne do wykresu, od 0 do 1000
it=sqrt(2)*abs(Iglowny)*sin((2*%pi*f0*t)+fii)//podanie wartosci chwilowej pradu
et=Esk*sqrt(2)*sin((2*%pi*f0*t)+fiu)//podanie wrt chwilowej napięcia
fi=fiu-fii//kat nap. - prad.
pt=Esk*abs(Iglowny)*cos(fi)-Esk*abs(Iglowny)*cos(2*(2*%pi*f0*t)+fii+fiu)//moc chwilowa
S=abs(Iglowny)*abs(Esk)//moc poorna
Psk=S*cos(fi)//moc czynna
Qsk=S*sin(fi)//moc bierna
xc1=[0;real(Ic1)]//porawanie wart do wykresow wektorowych
yc1=[0;imag(Ic1)]
xr1=[0;real(Ir1)]
yr1=[0;imag(Ir1)]
xl1=[0;real(Il1)]
yl1=[0;imag(Il1)]
xg=[0;real(Iglowny)]
yg=[0;imag(Iglowny)]
xr2=[real(Iglowny);real(Iglowny)+real(Ir2)]
yr2=[imag(Iglowny);imag(Iglowny)+imag(Ir2)]
xl2=[real(Iglowny);real(Iglowny)+real(Il2)]
yl2=[imag(Iglowny);imag(Iglowny)+imag(Il2)]
xsz2=[real(Iglowny);real(Iglowny)+real(Isz2)]
ysz2=[imag(Iglowny);imag(Iglowny)+imag(Isz2)]
legends(['Prad Ic1 na C1';'Prad Ir1 na R1';'Prad Il1 na L1';'Prad Iglowny';'Prad ir2 na R2';'Prad Il2 na L2';'Prad Isz2 w galezi RC';'Napiecie Uab';'Napiecie Ucd'],[1,2,3,4,5,6,7,9,10],opt=4)//legenda
subplot(1,2,1)//podział okna graficznego
plot2d4(xc1,yc1,1,axesflag=0)//wykresy wektorowe
plot2d4(xr1,yr1,2,axesflag=0)
plot2d4(xl1,yl1,3,axesflag=0)
plot2d4(xg,yg,4,axesflag=0)
plot2d4(xr2,yr2,5,axesflag=0)
plot2d4(xl2,yl2,6,axesflag=0)
plot2d4(xsz2,ysz2,7,axesflag=0)
x1=[0;real(Uab);real(Uab)+real(Ucd);real(Uab)+real(Ucd)-real(Uad)]
y1=[0;imag(Uab);imag(Uab)+imag(Ucd);imag(Uab)+imag(Ucd)-imag(Uad)]
x2=[real(Uab);real(Uab)+real(UC2);real(Uab)+real(UC2)+real(UR3)]
y2=[imag(Uab);imag(Uab)+imag(UC2);imag(Uab)+imag(UC2)+imag(UR3)]
plot2d4(x1,y1,9,axesflag=0)
plot2d4(x2,y2,10,axesflag=0)
subplot(2,2,2)
plot2d1(t,et,5,axesflag=1)//wykresy sinusoidalne
plot2d1(t,it,6,axesflag=1)
plot2d1(t,pt,3,axesflag=5)
PRZYKŁAD
podaj wartosc rez 1: -->1
podaj wartosc indukcji C1: -->2
poraj wartość indukcji L1: -->1
Z1 =
2. - 0.5i
Zab =
0.4705882 + 0.1176471i
podaj wartosc rez 2: -->2
podaj wartosc indukcji L2: -->2
podaj wartosc indukcji C2: -->1
podaj wartosc rez 3: -->1
Zsz =
1. + i
Z2 =
1. - i
Zcd =
0.5 + 0.5i
Zglowny =
0.9705882 + 0.6176471i
Zabs =
1.1504475
fi_z =
0.5667292
definiowanie skutecznej wartosci E: -->5
definioewanie kąta poczatkowego E: -->35
Etl =
- 4.518461 - 2.1409133i
Iglowny =
- 4.312631 + 0.5386120i
fii =
3.0173442
f0 =
50.
Uad =
- 4.518461 - 2.1409133i
Uab =
- 2.0928395 - 0.2539039i
Ic1 =
- 1.0464198 - 0.1269519i
Ir1 =
- 2.0928395 - 0.2539039i
Il1 =
- 2.0928395 - 0.2539039i
Ucd =
- 2.4256215 - 1.8870095i
Ir2 =
- 1.2128108 - 0.9435047i
Il2 =
- 1.2128108 - 0.9435047i
Isz2 =
- 2.1563155 + 0.2693060i
UXc1 =
- 2.0928395 - 0.2539039i
UR1 =
- 2.0928395 - 0.2539039i
UL1 =
- 2.0928395 - 0.2539039i
UR2 =
- 2.4256215 - 1.8870095i
UL2 =
- 2.4256215 - 1.8870095i
UR3 =
- 2.1563155 + 0.2693060i
UC2 =
- 2.1563155 + 0.2693060i
it =
column 1 to 5
0.7617124 0.5695705 0.3768653 0.1837875 - 0.0094721 …
…
S =
21.730675
Psk =
18.333333
Qsk =
11.666667
xc1 =
0.
- 1.0464198
yc1 =
0.
- 0.1269519
xr1 =
0.
- 2.0928395
yr1 =
0.
- 0.2539039
xl1 =
0.
- 2.0928395
yl1 =
0.
- 0.2539039
xg =
0.
- 4.312631
yg =
0.
0.5386120
xr2 =
- 4.312631
- 5.5254417
yr2 =
0.5386120
- 0.4048927
xl2 =
- 4.312631
- 5.5254417
yl2 =
0.5386120
- 0.4048927
xsz2 =
- 4.312631
- 6.4689465
ysz2 =
0.5386120
0.8079180
x1 =
0.
- 2.0928395
- 4.518461
0.
y1 =
0.
- 0.2539039
- 2.1409133
4.441D-16
x2 =
- 2.0928395
- 4.249155
- 6.4054705
y2 =
- 0.2539039
0.0154021
0.2847082