POLITECHHNIKA LUBELSKA W LUBLINIE			Laboratorium z metod numerycznych
						Ćwiczenie numer 3.
Imię: Albert	Nazwisko: Gałach		Grupa: 3.1 a	Semestr: I		Rok akademicki: 2008/2009
Ocena:		Data:			Podpis:

ZAD 2.

clear;//czyszczenie pamieci

xdel;//czyszczenie okien graficznych

clc;//czyszczenie konsoli

a=linspace(0,4*%pi,200);//ustawienie zakresu wykresu

subplot(3,1,1);//dzielenie okna graficznego

y=sin(a);//przebieg sinusoidalny

plot2d(a,y,5);//ustawienie przebiegu fun. i nadanie koloru

xtitle('przeieg sinusodalny','t(s)','y=sin(x)');//opis wykresu

subplot(3,1,2);//2 wykres w 1 oknie graficznym

y=(abs(sin(a))-sin(a+%pi))/2;//funkcja przedstawiajaca wykres 1-polowkowy

plot2d2(a,y,6);//ustawienie przebiegu fun. i nadanie koloru

xtitle('przebieg jednopolowkowy','t(s) ','y=sin(x)');//opis wykresu

subplot(3,1,3);//3 wykres w jednym oknie graficznym

y=abs(sin(a));//funkcja przedstawiajaca wykres 2 polowkowy

plot2d3(a,y,11);//ustawienie przebiegu fun. i nadanie koloru

xtitle('przebieg dwupołówkowy','t(s)','y=sin(x)');//opis wykresu

ZAD 3.

clear;//czyszczenie pamieci

xdel;//czyszcenieokien graficznych

clc;//czyszczenie konsoli

E=input("podaj napiecie zasilające: ");//wart napiecia

f=input("podaj czestotliwosc ukladu: ");//wart Hz

L=input("podaj indukcje cewki: ");//wart indukcj

C=input("podaj pojemnosc kondensatora: ");//wart pojemnosci kondensatora

R=input("podaj wartosc rezystora: ");//wart rezystora

w=2*%pi*f//obliczenie omegi

Zl=w*L*%i//impedancja cewki

Zc=-1*%i/(w*C)//impedancja kondensatora

Zz=Zl+Zc+R//impedancja zastepcza

I=E/Zz;//prad układu

Ia=abs(I)//modul I

Ur=R*I;//napiecie na rez

Ura=abs(Ur)//modul napiecia na rez

Ul=Zl*I;//napiecie na cewce

Ula=abs(Ul)//modul nap na L

Uc=Zc*I;//nap na C

Uca=abs(Uc)//modul nap na C

x=[0;real(E)]

y=[0;imag(E)]

plot2d4(x,y,5,axesflag=0)//wykres nap zasilajacego

x=[0;real(Ur)]

y=[0;imag(Ur)]

plot2d4(x,y,9,axesflag=0)//wykres nap na rezystorze

x=[0;real(Ul)]//

y=[0;imag(Ul)]

plot2d4(x,y,3,axesflag=0)//wykres nap na L

x=[0;real(Uc)]

y=[0;imag(Uc)]

plot2d4(x,y,2,axesflag=0)//wykres nap na C

x=[0;real(I)]

y=[0;imag(I)]

plot2d4(x,y,6,axesflag=0)//wykres pradu zasilajacego

legends(["E","Ur","Ul","Uc","I"],[5,9,3,2,6],opt=2)//legenda wykresu

PRZYKŁAD.

podaj napiecie zasilające: -->123

podaj czestotliwosc ukladu: -->34

podaj indukcje cewki: -->0.31

podaj pojemnosc kondensatora: -->0.000324

podaj wartosc rezystora: -->84

w =

213.6283

Zl =

66.224773i

Zc =

- 14.447616i

Zz =

84. + 51.777157i

Ia =

1.2465088

Ura =

104.70674

Ula =

82.54976

Uca =

18.009081