Laboratorium metod numerycznych
Laboratorium 6
ALBERT GAŁACH ANDRZEJ GREGUŁA	Grupa ED 3.1	Data wykonania ćw. 5.01.2009	Ocena

STANY NIEUSTALONE METODĄ KLASYCZNĄ

clear;//czyszczenie pamiąci

xdel(winsid());//czyszczenie okien graficznych

clc;//czyszczenie konsoli

function[pochodne]=stany(t,x);//definiuje funkcje opisujaca uklad rownan rozniczkowych

pochodne(1)=(E-x(1)*((L/R2)+R1*C)-x(2)*((R1/R2)+1))/(L*C);

pochodne(2)=x(3)/C;

pochodne(3)=C*x(1);

endfunction

R1=input("podaj wartosc rezystora 1: ");//podajemy wartosci układu

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

L=input("podaj wartosc cewki: ");

C=input("podaj wartosc kondensatora: ");

E=input("podaj wartosc zasilania: ");

duC0=0;//definiuje parametry poczatkowe zdefiniowanych zmiennych stanu i czasu

uC0=0;

i0=0;

t0=0;

t=0:0.001:0.1;//definiuje wektor czasu

roz=ode([duC0;uC0;i0],t0,t,stany);//polecenie rozwiazania ukladu

i3=roz(3,:);//przebiegi

uC=roz(2,:);

i2=uC/R2;

i1=i2+i3;

uR1=R1*i1;

uR2=R2*i2

uC=uR2

uL=E-uR1-uR2

Ez=uL+uR1+uR2

subplot(221);//dzieli okno graficzne i rysuje przebiegi

plot2d(t,i1,style=[9],leg='i1(t)',axesflag=5);

plot2d(t,[(5*i1)' i2' (1000*i3)'], style=[1,2,5], leg='i1(t)@i2(t)@i3(t)', axesflag=5);

subplot(212);

plot2d(t,[uC' uR1' uR2' uL' uC' Ez'], style=[3,2,5,4,7,9], leg='uC(t)@uR1(t)@uR2(t)@uL(t)@uC(t)@E', axesflag=5);

SCILAB

podaj wartosc rezystora 1: -->2

podaj wartosc rezystora 2: -->1

podaj wartosc cewki: -->0.2

podaj wartosc kondensatora: -->0.000005

podaj wartosc zasilania: -->5

uR2 =

column 1 to 8

0. 0.0000124 0.0000497 0.0001121 0.0001995 0.0003118 0.0004492 0.0006116

column 9 to 15

0.0007989 0.0010112 0.0012485 0.0015108 0.0017981 0.0021103 0.0024475

uC =

column 1 to 8

0. 0.0000124 0.0000497 0.0001121 0.0001995 0.0003118 0.0004492 0.0006116

column 9 to 15

0.0007989 0.0010112 0.0012485 0.0015108 0.0017981 0.0021103 0.0024475

uL =

column 1 to 8

5. 4.9999626 4.9998503 4.9996629 4.9994005 4.9990632 4.9986509 4.9981636

column 9 to 15

4.9976013 4.9969641 4.9962519 4.9954648 4.9946028 4.9936658 4.992654

Ez =

column 1 to 16

5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5.

column 17 to 32

5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5. 5.

WYKRES

---