POLITECHHNIKA LUBELSKA W LUBLINIE	Laboratorium z metod numerycznych
	Ćwiczenie numer 3.
Imię: Łukasz	Nazwisko: Kilijanek
Ocena:	Data:

Zadanie 1.

clear;// kasuje zmienne z pamięci

xdel;// zamka otwarte okna

clc;//czyści pamieć grafik

x=[0:0.03:10];//wprowadzenie zmiennej x

subplot(3,1,1);//wykres 1

a=sin(x);//obliczanie zmiennej a

plot2d(x,a,5);//rysowanie wykresu

xtitle("Sygnał sinusoidalny","x","a=sin(x)")//opis wykresu

subplot(3,1,2);//wykres 2

b=(abs(sin(x))-sin(x+%pi))/2//obliczanie zmiennej b

plot2d(x,b,4)//rysowanie wykresu

xtitle("Sygnał wyprostowany jednopołówkowo","x","b=sin(x)")//opis wykresu

subplot(3,1,3)//wykres 3

c=abs(sin(x))//obliczanie zmiennej c

plot2d(x,c,9)//rysowanie wykresu

xtitle("Sygnał wyprostowany dwupolówkowo","x","c=sin(x)")//opis wykresu

Obliczenia w scilabie

b =

column 1 to 5

- 6.123D-17 0.0299955 0.0599640 0.0898785 0.1197122

column 6 to 10

0.1494381 0.1790296 0.2084599 0.2377026 0.2667314

column 11 to 15

0.2955202 0.3240430 0.3522742 0.3801884 0.4077605

column 16 to 20

0.4349655 0.4617792 0.4881772 0.5141360 0.5396320

column 21 to 25

0.5646425 0.5891448 0.6131169 0.6365372 0.6593847

column 26 to 30

0.6816388 0.7032794 0.7242872 0.7446431 0.7643289

Wykresy

Zadanie 2

clear;//kasuje zmienne z pamięci

xdel();//zamyka otwrte okna

clc;//czyści pamięć grafik

Z=input("podaj impedancje rezystora:");//podajemy zmienną Z

C=input("podaj pojemnosc kondenasatora:");//podajemy zmienną C

L=input("podaj indukcyjnosc cewki:");//podajemy zmienną L

f=input("podaj czestotliwosc:");//podajemy zmienną f

U=input("podaj napiecie na zaciskach:");//podajemy napięcie

i=(U/(Z+2*%i*%pi*f*L-1*%i/(2*%pi*f*C)));//obliczamy prąd

iz=abs(i)

"spadek napiecia na rezystorze"

Ur=i*Z

Urz=abs(Ur)

"spadek napiecia na kondensatorze"

Uc=i*(1/(C*%pi*2*f))

Ucz=abs(Uc)

"spadek napiecia na cewce"

Ul=i*(L*%pi*2*f)

Ulz=abs(Ul);

x=0:real(U):real(U)//wektor U

y=(imag(U)/real(U))*x

plot2d4(x, y, axesflag=0)

x=0:real(i):real(i)//wektor i

y=(imag(i)/real(i))*x

plot2d4(x, y, 2, axesflag=0)

x=0:real(Ur):real(Ur)//wektor Ur

y=(imag(Ur)/real(Ur))*x

plot2d4(x, y, 3, axesflag=0)

x=0:real(Uc):real(Uc)//wektor Uc

y=(imag(Uc)/real(Uc))*x

plot2d4(x, y, 4, axesflag=0)

x=0:real(Ul):real(Ul)//wektor Ul

y=(imag(Ul)/real(Ul))*x

plot2d4(x, y, 5, axesflag=0)

legends(["i", "U", "Ur", "Uc", "Ul"], [1;2;3;4;5], opt=4)

Obliczenia w scilabie

podaj impedancje rezystora:-->50

podaj pojemnosc kondenasatora:-->0.00012

podaj indukcyjnosc cewki:-->0.54

podaj czestotliwosc:-->50

podaj napiecie na zaciskach:-->230

iz =1.5171231

ans =

spadek napiecia na rezystorze

Ur =25.018072 - 71.611819i

Urz =75.856157

ans =

spadek napiecia na kondensatorze

Uc = 13.272499 - 37.99125i

Ucz = 40.242941

ans =

spadek napiecia na cewce

Ul =84.884319 - 242.97318i

Ulz =257.37388

Wykres