21.03.2011, Elektrotechnika I stopień PWSZ Leszno, SEMESTR II, Metody Numeryczne, 2. 21.03.2011


0x08 graphic
Sprawozdanie

laboratoria Metod Numerycznych 21.03.2011

*Metoda Eulera - jest metodą iteracyjnego wyznaczania zagadnienia początkowego

*Algorytm z odwracaniem- jest związany z twierdzeniem o próbkowaniu. Częstotliwość próbkowania musi być co najmniej dwukrotnie większa od częstotliwości sygnału próbkowanego.

*Wzór EULERA:

0x01 graphic

I . Metoda „Eulera w przód” (forword)

  1. Implementacja:

Function eulerFordword(K,M,B,x0,t,F,deltat)

%B współczynnik tarcia

%M - masa klocka

%x0- punkt startowy

%t- czas symulacji

%deltat - krok czasowy

x= x0;

N=t/deltat;

A=[0 1; -(K/M) -(B/M)]

B=[0; (1/M)]

for i=1: N

x-x+deltat*(A*x+B*F);

x1_w(i) = x(1);

x2_w(i) = x(2);

t_w(i)= deltat*I;

end;

  1. Sprawdzenie sposobu działania dla danych:

A=[0 1; -(K/M) -(B/M)]

B=[0; (1/M)]

  1. Symulacja układu dla danych:

M= 1 ;

B = 1;

x0=[0;0];

t=0x01 graphic

F= 1;

deltat= 0x01 graphic


0x01 graphic

K= 100;

M= 1 ;

B = 1;

x0=[0;0];

t=0x01 graphic

F= 1;

deltat= 0x01 graphic

0x01 graphic

K= 100;

M= 1 ;

B = 1;

x0=[0;0];

t=10;

F= 5;

deltat= 0x01 graphic

0x01 graphic

II . Metoda „Eulera w tył” (backword)

  1. Function eulerBackword(K,M,B,x0,t,F,deltat)

%B współczynnik tarcia

%M - masa klocka

%x0- punkt startowy

%t- czas symulacji

%deltat - krok czasowy

x= x0;

N=t/deltat;

I=[10;01]

A=[0 1; -(K/M) -(B/M)]

B=[0; (1/M)]

for i=1 :N

x=inv(I-dt*A)*x+inv(I-dt*A)*dt*B*F;

x1_w(i) = x(1);

x2_w(i) = x(2);

t_w(i)= dt*i;

end;

  1. Sprawdzenie sposobu działania dla danych:

A=[0 1; -(K/M) -(B/M)]

B=[0; (1/M)]

  1. Symulacja układu dla danych:

M= 1 ;

B = 1;

x0=[0;0];

t=0x01 graphic

F= 5;

deltat= 0x01 graphic

0x01 graphic

III. Wnioski:

- są to metody jedne z prostszych lecz nie zawsze wychodzi prawidłowy wynik

Wykonali : Krzysztof Maćkowski, Daniel Szepe



Wyszukiwarka