METODY NUMERYCZNE
Elektrotechnika, rok II, sem. III

ZADANIE PROJEKTOWE

ZESPÓŁ PRZYGOTOWUJĄCY ZADANIE PROJEKTOWE
Lp.	Nazwisko	Imię
I	Madejski	Grzegorz
II	Michalski	Paweł

Dany jest układ automatycznej regulacji:

Wyznaczyć numerycznie rozwiązanie odpowiedzi układu y(t) na wymuszenie: x(t)=1(t) stosując czteropunktową metodę Runge-Kutta. Przyjąć warunki początkowe: y(0)=1; y'(0)=y”(0)=0. Dane:

G₁(S)=
G₂(S)=

Przy czym:

a- liczba liter I imienia

c- liczba liter II imienia

b- suma liter I i II nazwiska

Wyniki obliczeń wykonać w Excelu tabelarycznie i w postaci wykresów.

1.Liczymy transmitancję układu otwartego.

K(S)= G₁(S)*G₂(S)

gdzie: a=8, b=17, c=5

G_z=

G_z=

40*X(S) = Y(S)*(s³+18s²+18s+41) X=1(t)=1

y'''+18y''+18y'+41y = 40

Podstawiając do równania odpowiedni: y'''=x', y''=x, y'=z, y=y

Zamiast jednego równania III rzędu otrzymujemy trzy równania rzędu pierwszego.

Dla odpowiednich pochodnych przyjmujemy współczynniki:

K / y'=z

M / z'=x

N / x'=-18x-18z-41y+40

Następnie liczymy współczynniki:

K₁= h*z_K

M₁= h*x_K

N₁=h*(-18x_k-18z_k-41y_k+40)

K₂ =h*(z_k+M₁/2)

M₂=h*(x_k+N₁/2)

N₂=h*[-18(x_k+N₁/2)-18(z_k+M₁/2)-41(y_k+K₁/2)+40]

K₃= h*(z_k+M₂/2)

M₃= h*(x_k+N₂/2)

N₃= h*[-18(x_k+N₂/2)-18(z_k+M₂/2)-41(y_k+K₂/2)+40]

K₄= h*(z_k+M₃)

M₄= h*(x_k+N₃)

N₄= h*[-18(x_k+N₃)-18(z_k+M₃)-41(y_k+K₃)+40]

I podstawiamy je odpowiednio do wzorów:

y_k+1=y_k+1/6*(K₁+2K₂+2K₃+K₄)

z_k+1=z_k+1/6*(M₁+2M₂+2M₃+M₄)

x_k+1=x_k+1/6*(N₁+2N₂+2N₃+N₄)

Przyjmujemy wartości: h=0,01 oraz k=1,2,3….

gdzie t_k=t₀+k*h

Obliczenia wykonujemy dla 500 pozycji.

---