• Nie Znaleziono Wyników

13. Równania różniczkowe - rozwiązywanie

N/A
N/A
Protected

Academic year: 2021

Share "13. Równania różniczkowe - rozwiązywanie"

Copied!
175
0
0

Pełen tekst

(1)

13. Równania różniczkowe - rozwiązywanie

Grzegorz Kosiorowski

Uniwersytet Ekonomiczny w Krakowie

(2)

1 Podstawowe oznaczenia i definicje

2 Równania różniczkowe: przykłady ekonomiczne

3 Podstawowe typy równań różniczkowych i sposoby ich rozwiązywania

Równania o zmiennych rozdzielonych Równania różniczkowe liniowe

Podstawienia w równaniach różniczkowych

(3)

Wstęp

Często zdarza się, że poszukując zależności pomiędzy wielkościami ekonomicznymi (np. A i B) napotykamy na informacje dotyczące prędkości wzrostu wielkości A w stosunku do wielkości B. Jeśli potrafimy opisać tę prędkość jako funkcję A0(B) to wystarczy ją scałkować, by uzyskać zależność funkcyjną A(B). Jednak zdarza się, że mamy tylko pewne informacje o A0(B) przedstawione w formie równania, a nie sam wzór. Takie równanie, którego niewiadomą jest funkcja i w którym występują pochodne tej funkcji nazywamy równaniem różniczkowym. Czasem takie równania się da rozwiązać i o tego typu równaniach jest ten rozdział.

Dodam jeszcze, że równania różniczkowe mają poważne zastosowania w jakichkolwiek badaniach naukowych, w których występuje ruch lub zmiana warunków (tzw. układ dynamiczny). Na nich opiera się cała klasyczna fizyka, a coraz częściej są konieczne do modelowania zjawisk z dziedzin innych nauk przyrodniczych i społecznych.

(4)

Wstęp

Często zdarza się, że poszukując zależności pomiędzy wielkościami ekonomicznymi (np. A i B) napotykamy na informacje dotyczące prędkości wzrostu wielkości A w stosunku do wielkości B. Jeśli potrafimy opisać tę prędkość jako funkcję A0(B) to wystarczy ją scałkować, by uzyskać zależność funkcyjną A(B). Jednak zdarza się, że mamy tylko pewne informacje o A0(B) przedstawione w formie równania, a nie sam wzór. Takie równanie, którego niewiadomą jest funkcja i w którym występują pochodne tej funkcji nazywamy równaniem różniczkowym. Czasem takie równania się da rozwiązać i o tego typu równaniach jest ten rozdział.

Dodam jeszcze, że równania różniczkowe mają poważne zastosowania w jakichkolwiek badaniach naukowych, w których występuje ruch lub

(5)

Oznaczenia - komentarz

W tym rozdziale rozważamy funkcje różniczkowalne jednej zmiennej y (x ), gdzie y : R ⊃ Dy → R. Pochodną takiej funkcji oznaczaliśmy przez y0, jednak w tym rozdziale wyjątkowo będziemy zazwyczaj używać bardziej fizycznej notacji dydx (ze względów

„mnemotechnicznych” - jak zobaczymy, taki zapis ułatwi nam rozwiązywanie pewnych typów równań).

(6)

Definicja

Równanie różniczkowe zwyczajne

Niech f : Rn+2 ⊃ Df → R. Wtedy równaniem różniczkowym zwyczajnym rzędu n nazywamy równanie

f (x , y (x ), y0(x ), y00(x ), . . . , y(n)(x )) = 0, w którym niewiadomą jest funkcja jednej zmiennej y i w którym występuje pochodna rzędu n tej funkcji (y(n)(x )) i mogą występować pochodne niższych rzędów, sama funkcja y i zmienna niezależna x . Rozwiązaniem (lub całką) takiego równania jest n-krotnie różniczkowalna funkcja y , która je spełnia dla każdego x w swojej dziedzinie.

Słowo „zwyczajny” w definicji odróżna te równania od równań różniczkowych cząstkowych, które odpowiadają za poszukiwanie funkcji wielu zmiennych na podstawie informacji o ich pochodnych cząstkowych. Takie równania są o wiele bardziej skomplikowane i nie będziemy się nimi zajmować.

(7)

Definicja

Równanie różniczkowe zwyczajne

Niech f : Rn+2 ⊃ Df → R. Wtedy równaniem różniczkowym zwyczajnym rzędu n nazywamy równanie

f (x , y (x ), y0(x ), y00(x ), . . . , y(n)(x )) = 0, w którym niewiadomą jest funkcja jednej zmiennej y i w którym występuje pochodna rzędu n tej funkcji (y(n)(x )) i mogą występować pochodne niższych rzędów, sama funkcja y i zmienna niezależna x . Rozwiązaniem (lub całką) takiego równania jest n-krotnie różniczkowalna funkcja y , która je spełnia dla każdego x w swojej dziedzinie.

Słowo „zwyczajny” w definicji odróżna te równania od równań różniczkowych cząstkowych, które odpowiadają za poszukiwanie funkcji wielu zmiennych na podstawie informacji o ich pochodnych cząstkowych. Takie równania są o wiele bardziej skomplikowane i nie

(8)

Przykład - II zasada dynamiki Newtona

Zasadę dynamiki Newtona znają Państwo na pewno w postaci F = am, gdzie F jest siłą działającą na ciało, a - przyspieszeniem, jakiego ciało nabiera pod wpływem tej siły, a m - masą ciała.

Zauważmy, że jeśli przez x (t) oznaczymy położenie ciała w czasie t to F może być dowolną funkcją zależną od x i t, a(t) jest drugą pochodną funkcji x , a m jest stałą. Dlatego jeśli zdefiniujemy f (t, x (t), x0(t), x00(t)) = F (x (t), t) − m · x00(t), to

f (t, x (t), x0(t), x00(t)) = 0 jest równaniem różniczkowym drugiego stopnia, opisującym II zasadę dynamiki Newtona. Jeśli założymy, że F jest stałe i mF = a ∈ R to równanie różniczkowe F − mx00(t) = 0 jest równoważne równaniu x00(t) = a i jego rozwiązanie opisuje ruch ciała pod wpływem stałej, działającej na nie siły.

(9)

Przykład - II zasada dynamiki Newtona

Zasadę dynamiki Newtona znają Państwo na pewno w postaci F = am, gdzie F jest siłą działającą na ciało, a - przyspieszeniem, jakiego ciało nabiera pod wpływem tej siły, a m - masą ciała.

Zauważmy, że jeśli przez x (t) oznaczymy położenie ciała w czasie t to F może być dowolną funkcją zależną od x i t, a(t) jest drugą pochodną funkcji x , a m jest stałą.

Dlatego jeśli zdefiniujemy f (t, x (t), x0(t), x00(t)) = F (x (t), t) − m · x00(t), to

f (t, x (t), x0(t), x00(t)) = 0 jest równaniem różniczkowym drugiego stopnia, opisującym II zasadę dynamiki Newtona. Jeśli założymy, że F jest stałe i mF = a ∈ R to równanie różniczkowe F − mx00(t) = 0 jest równoważne równaniu x00(t) = a i jego rozwiązanie opisuje ruch ciała pod wpływem stałej, działającej na nie siły.

(10)

Przykład - II zasada dynamiki Newtona

Zasadę dynamiki Newtona znają Państwo na pewno w postaci F = am, gdzie F jest siłą działającą na ciało, a - przyspieszeniem, jakiego ciało nabiera pod wpływem tej siły, a m - masą ciała.

Zauważmy, że jeśli przez x (t) oznaczymy położenie ciała w czasie t to F może być dowolną funkcją zależną od x i t, a(t) jest drugą pochodną funkcji x , a m jest stałą. Dlatego jeśli zdefiniujemy f (t, x (t), x0(t), x00(t)) = F (x (t), t) − m · x00(t), to

f (t, x (t), x0(t), x00(t)) = 0 jest równaniem różniczkowym drugiego stopnia, opisującym II zasadę dynamiki Newtona.

Jeśli założymy, że F jest stałe i mF = a ∈ R to równanie różniczkowe F − mx00(t) = 0 jest równoważne równaniu x00(t) = a i jego rozwiązanie opisuje ruch ciała pod wpływem stałej, działającej na nie siły.

(11)

Przykład - II zasada dynamiki Newtona

Zasadę dynamiki Newtona znają Państwo na pewno w postaci F = am, gdzie F jest siłą działającą na ciało, a - przyspieszeniem, jakiego ciało nabiera pod wpływem tej siły, a m - masą ciała.

Zauważmy, że jeśli przez x (t) oznaczymy położenie ciała w czasie t to F może być dowolną funkcją zależną od x i t, a(t) jest drugą pochodną funkcji x , a m jest stałą. Dlatego jeśli zdefiniujemy f (t, x (t), x0(t), x00(t)) = F (x (t), t) − m · x00(t), to

f (t, x (t), x0(t), x00(t)) = 0 jest równaniem różniczkowym drugiego stopnia, opisującym II zasadę dynamiki Newtona. Jeśli założymy, że F jest stałe i mF = a ∈ R to równanie różniczkowe F − mx00(t) = 0 jest równoważne równaniu x00(t) = a i jego rozwiązanie opisuje ruch ciała pod wpływem stałej, działającej na nie siły.

(12)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą.

W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t). Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała. Jedyne funkcje spełniające tę zależność to wielomiany

pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1. Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t). Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(13)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą. W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t).

Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała. Jedyne funkcje spełniające tę zależność to wielomiany

pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1. Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t). Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(14)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą. W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t). Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała.

Jedyne funkcje spełniające tę zależność to wielomiany pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1. Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t). Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(15)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą. W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t). Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała. Jedyne funkcje spełniające tę zależność to wielomiany

pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1.

Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t). Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(16)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą. W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t). Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała. Jedyne funkcje spełniające tę zależność to wielomiany

pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1. Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t).

Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(17)

Przykład - II zasada dynamiki Newtona

Spróbujmy rozwiązać równanie x00(t) = a, gdzie a ∈ R jest stałą. W tym celu, wprowadźmy funkcję pomocniczą y (t) = x0(t). Naturalnie y0(t) = x00(t) = a. Poszukujemy zatem funkcji, której pochodna jest stała. Jedyne funkcje spełniające tę zależność to wielomiany

pierwszego stopnia, o współczynniku a przy wyrazie w pierwszej potędze. Możemy więc zapisać, że y (t) = at + C1. Stąd

x0(t) = at + C1. Innymi słowy, x jest funkcją pierwotną do at + C1 (względem zmiennej t). Całkując at + C1 otrzymamy, że x musi być postaci a2t2+ C1t + C2, C1, C2 ∈ R.

(18)

Przykład - II zasada dynamiki Newtona

x (t) = a2t2+ C1t + C2, C1, C2 ∈ R.

Jeśli założymy, że ciało, którego ruch badamy, zaczynało w chwili t = 0 z położenia x (0) = 0 i było nieruchome tj. miało zerową prędkość początkową

(y (0) = x0(0) = 0), to możemy obliczyć:

0 = y (0) = a·0+C1 ⇒ C1 = 0; 0 = x (0) = a

2·0+C1·0+C2 ⇒ C2 = 0. Zatem ostateczny wzór na ruch pod wpływem stale działającej siły to x (t) = a2t2 - znany ze szkolnej fizyki wzór na pokonaną drogę w ruchu jednostajnie przyspieszonym.

(19)

Przykład - II zasada dynamiki Newtona

x (t) = a2t2+ C1t + C2, C1, C2 ∈ R. Jeśli założymy, że ciało, którego ruch badamy, zaczynało w chwili t = 0 z położenia x (0) = 0 i było nieruchome tj. miało zerową prędkość początkową

(y (0) = x0(0) = 0), to możemy obliczyć:

0 = y (0) = a·0+C1 ⇒ C1 = 0; 0 = x (0) = a

2·0+C1·0+C2 ⇒ C2 = 0. Zatem ostateczny wzór na ruch pod wpływem stale działającej siły to x (t) = a2t2 - znany ze szkolnej fizyki wzór na pokonaną drogę w ruchu jednostajnie przyspieszonym.

(20)

Przykład - II zasada dynamiki Newtona

x (t) = a2t2+ C1t + C2, C1, C2 ∈ R. Jeśli założymy, że ciało, którego ruch badamy, zaczynało w chwili t = 0 z położenia x (0) = 0 i było nieruchome tj. miało zerową prędkość początkową

(y (0) = x0(0) = 0), to możemy obliczyć:

0 = y (0) = a·0+C1 ⇒ C1 = 0;

0 = x (0) = a

2·0+C1·0+C2 ⇒ C2 = 0. Zatem ostateczny wzór na ruch pod wpływem stale działającej siły to x (t) = a2t2 - znany ze szkolnej fizyki wzór na pokonaną drogę w ruchu jednostajnie przyspieszonym.

(21)

Przykład - II zasada dynamiki Newtona

x (t) = a2t2+ C1t + C2, C1, C2 ∈ R. Jeśli założymy, że ciało, którego ruch badamy, zaczynało w chwili t = 0 z położenia x (0) = 0 i było nieruchome tj. miało zerową prędkość początkową

(y (0) = x0(0) = 0), to możemy obliczyć:

0 = y (0) = a·0+C1 ⇒ C1 = 0; 0 = x (0) = a

2·0+C1·0+C2 ⇒ C2 = 0.

Zatem ostateczny wzór na ruch pod wpływem stale działającej siły to x (t) = a2t2 - znany ze szkolnej fizyki wzór na pokonaną drogę w ruchu jednostajnie przyspieszonym.

(22)

Przykład - II zasada dynamiki Newtona

x (t) = a2t2+ C1t + C2, C1, C2 ∈ R. Jeśli założymy, że ciało, którego ruch badamy, zaczynało w chwili t = 0 z położenia x (0) = 0 i było nieruchome tj. miało zerową prędkość początkową

(y (0) = x0(0) = 0), to możemy obliczyć:

0 = y (0) = a·0+C1 ⇒ C1 = 0; 0 = x (0) = a

2·0+C1·0+C2 ⇒ C2 = 0.

Zatem ostateczny wzór na ruch pod wpływem stale działającej siły to x (t) = at2 - znany ze szkolnej fizyki wzór na pokonaną drogę w

(23)

Postać normalna

Teraz, do końca tego rozdziału, skupimy się na równaniach różniczkowych zwyczajnych pierwszego rzędu, czyli zawierających tylko zmienną niezależną, funkcję oraz jej pierwszą pochodną.

Warto zauważyć, że równanie różniczkowe z definicji dane jest w postaci uwikłanej ze względu na wszystkie zmienne. Najczęściej równania te jednak będą w postaci rozwikłanej ze względu na pierwszą pochodną, czyli tzw. normalnej.

Równanie różniczkowe zwyczajne

Postać normalna równania różniczkowego pierwszego rzędu to równanie postaci dydx = g (x , y (x )), gdzie Dg ⊂ R2.

(24)

Postać normalna

Teraz, do końca tego rozdziału, skupimy się na równaniach różniczkowych zwyczajnych pierwszego rzędu, czyli zawierających tylko zmienną niezależną, funkcję oraz jej pierwszą pochodną. Warto zauważyć, że równanie różniczkowe z definicji dane jest w postaci uwikłanej ze względu na wszystkie zmienne.

Najczęściej równania te jednak będą w postaci rozwikłanej ze względu na pierwszą pochodną, czyli tzw. normalnej.

Równanie różniczkowe zwyczajne

Postać normalna równania różniczkowego pierwszego rzędu to równanie postaci dydx = g (x , y (x )), gdzie Dg ⊂ R2.

(25)

Postać normalna

Teraz, do końca tego rozdziału, skupimy się na równaniach różniczkowych zwyczajnych pierwszego rzędu, czyli zawierających tylko zmienną niezależną, funkcję oraz jej pierwszą pochodną. Warto zauważyć, że równanie różniczkowe z definicji dane jest w postaci uwikłanej ze względu na wszystkie zmienne. Najczęściej równania te jednak będą w postaci rozwikłanej ze względu na pierwszą pochodną, czyli tzw. normalnej.

Równanie różniczkowe zwyczajne

Postać normalna równania różniczkowego pierwszego rzędu to równanie postaci dydx = g (x , y (x )), gdzie Dg ⊂ R2.

(26)

Równania różniczkowe - jednoznaczność rozwiązań

Rozwiązywanie równania różniczkowego jest w pewnym sensie

uogólnionym całkowaniem - wiedząc coś o pochodnej, szukamy wzoru funkcji. Dlatego równanie różniczkowe najczęściej ma nieskończenie rozwiązań (tak jak całka).

Aby otrzymać pojedyncze rozwiązanie, narzuca się dodatkowe warunki np. wartość w pewnym „punkcie początkowym”. Tak właśnie było w rozwiązanym przed chwilą problemie: na początku dostaliśmy wzór ogólny na x (t), zależny od stałych C1, C2, a dopiero potem, po uzyskaniu informacji o

początkowym położeniu i prędkości badanego obiektu, otrzymaliśmy jednoznaczny wzór funkcji x . Jak uogólnić ten wynik?

(27)

Równania różniczkowe - jednoznaczność rozwiązań

Rozwiązywanie równania różniczkowego jest w pewnym sensie

uogólnionym całkowaniem - wiedząc coś o pochodnej, szukamy wzoru funkcji. Dlatego równanie różniczkowe najczęściej ma nieskończenie rozwiązań (tak jak całka). Aby otrzymać pojedyncze rozwiązanie, narzuca się dodatkowe warunki np. wartość w pewnym „punkcie początkowym”.

Tak właśnie było w rozwiązanym przed chwilą problemie: na początku dostaliśmy wzór ogólny na x (t), zależny od stałych C1, C2, a dopiero potem, po uzyskaniu informacji o

początkowym położeniu i prędkości badanego obiektu, otrzymaliśmy jednoznaczny wzór funkcji x . Jak uogólnić ten wynik?

(28)

Równania różniczkowe - jednoznaczność rozwiązań

Rozwiązywanie równania różniczkowego jest w pewnym sensie

uogólnionym całkowaniem - wiedząc coś o pochodnej, szukamy wzoru funkcji. Dlatego równanie różniczkowe najczęściej ma nieskończenie rozwiązań (tak jak całka). Aby otrzymać pojedyncze rozwiązanie, narzuca się dodatkowe warunki np. wartość w pewnym „punkcie początkowym”. Tak właśnie było w rozwiązanym przed chwilą problemie: na początku dostaliśmy wzór ogólny na x (t), zależny od stałych C1, C2, a dopiero potem, po uzyskaniu informacji o

początkowym położeniu i prędkości badanego obiektu, otrzymaliśmy jednoznaczny wzór funkcji x .

Jak uogólnić ten wynik?

(29)

Równania różniczkowe - jednoznaczność rozwiązań

Rozwiązywanie równania różniczkowego jest w pewnym sensie

uogólnionym całkowaniem - wiedząc coś o pochodnej, szukamy wzoru funkcji. Dlatego równanie różniczkowe najczęściej ma nieskończenie rozwiązań (tak jak całka). Aby otrzymać pojedyncze rozwiązanie, narzuca się dodatkowe warunki np. wartość w pewnym „punkcie początkowym”. Tak właśnie było w rozwiązanym przed chwilą problemie: na początku dostaliśmy wzór ogólny na x (t), zależny od stałych C1, C2, a dopiero potem, po uzyskaniu informacji o

początkowym położeniu i prędkości badanego obiektu, otrzymaliśmy jednoznaczny wzór funkcji x . Jak uogólnić ten wynik?

(30)

Równania różniczkowe - jednoznaczność rozwiązań

Zagadnienie Cauchy’ego

Zagadnienie Cauchy’ego (lub zagadnienie początkowe) to zadanie polegające na wyznaczeniu rozwiązania równania różniczkowego spełniającego warunek:

y0(x ) = g (x , y (x )) y (x0) = y0 .

Twierdzenie Peano-Piccarda

Jeśli g : R2 → R jest funkcją różniczkowalną w pewnym otoczeniu (x0, y0) to zagadnienie Cauchy’ego zadane tak jak w powyższej definicji posiada dokładnie jedno rozwiązanie y w pewnym otoczeniu x0.

(31)

Równania różniczkowe - jednoznaczność rozwiązań

Zagadnienie Cauchy’ego

Zagadnienie Cauchy’ego (lub zagadnienie początkowe) to zadanie polegające na wyznaczeniu rozwiązania równania różniczkowego spełniającego warunek:

y0(x ) = g (x , y (x )) y (x0) = y0 .

Twierdzenie Peano-Piccarda

Jeśli g : R2 → R jest funkcją różniczkowalną w pewnym otoczeniu (x0, y0) to zagadnienie Cauchy’ego zadane tak jak w powyższej definicji posiada dokładnie jedno rozwiązanie y w pewnym otoczeniu x0.

(32)

Równania różniczkowe - jednoznaczność rozwiązań

Naturalnym pytaniem jest, czy istnieją zagadnienia Cauchy’ego, które nie mają rozwiązań lub mają więcej niż jedno. Oczywiście, w takiej sytuacji funkcja g z definicji zagadnienia nie może być

różniczkowalna.

Jeśli funkcja g nie jest ciągła w (x0, y0), to zagadnienie może w ogóle nie mieć rozwiązań. Przykładowo:

y0(x ) =

0, gdy x < 0; 1, gdy x ­ 0 y (0) = 0

,

nie ma w ogóle rozwiązań.

(33)

Równania różniczkowe - jednoznaczność rozwiązań

Naturalnym pytaniem jest, czy istnieją zagadnienia Cauchy’ego, które nie mają rozwiązań lub mają więcej niż jedno. Oczywiście, w takiej sytuacji funkcja g z definicji zagadnienia nie może być

różniczkowalna.

Jeśli funkcja g nie jest ciągła w (x0, y0), to zagadnienie może w ogóle nie mieć rozwiązań. Przykładowo:

y0(x ) =

0, gdy x < 0;

1, gdy x ­ 0 y (0) = 0

,

nie ma w ogóle rozwiązań.

(34)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie.

Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem

Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(35)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem

Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(36)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0

(bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem

Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(37)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi

3x2, a prawa również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem

Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(38)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa

również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem

Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(39)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa również 3x2 dla dowolnego x ∈ R.

Nie jest to sprzeczne z twierdzeniem Peano-Piccarda, gdyż funkcja g (x , y ) = 3√3

y2 nie ma pochodnej cząstkowej gy0 w punkcie (0, 0).

(40)

Równania różniczkowe - jednoznaczność rozwiązań

Z kolei, jeśli funkcja g jest ciągła, ale nieróżniczkowalna w (x0, y0), to zagadnienie może mieć dla tego warunku początkowego więcej niż jedno rozwiązanie. Przykładowo rozważmy:

y0(x ) = 33 y2 y (0) = 0 .

Oczywiście rozwiązaniem jest funkcja stała y (x ) = 0 (bo obie strony równania różniczkowego się zerują), ale też funkcja y (x ) = x3, bo lewa strona równania różniczkowego wynosi 3x2, a prawa również 3x2 dla dowolnego x ∈ R. Nie jest to sprzeczne z twierdzeniem√

(41)

Model Friedmana oczekiwań inflacyjnych

W warunkach długotrwałej inflacji ludzie przyjmują pewne

oczekiwania inflacyjne (czyli przewidują przyszłą inflację na pewnym poziomie), co wpływa na ich zachowania rynkowe.

Milton Friedman modelował te oczekiwania za pomocą prostego równania

różniczkowego:

d π

dt = a(p − π),

gdzie a ∈ (0, 1], π(t) oznacza oczekiwaną, a p(t) rzeczywistą stopę inflacji w momencie t. Model powstał na bazie następującej

obserwacji: jeśli w danej chwili p > π to stopa oczekiwana będzie rosnąć (czyli d πdt > 0), zaś jeśli p < π to stopa oczekiwana będzie maleć (czyli d πdt < 0).

(42)

Model Friedmana oczekiwań inflacyjnych

W warunkach długotrwałej inflacji ludzie przyjmują pewne

oczekiwania inflacyjne (czyli przewidują przyszłą inflację na pewnym poziomie), co wpływa na ich zachowania rynkowe. Milton Friedman modelował te oczekiwania za pomocą prostego równania

różniczkowego:

d π

dt = a(p − π),

gdzie a ∈ (0, 1], π(t) oznacza oczekiwaną, a p(t) rzeczywistą stopę inflacji w momencie t.

Model powstał na bazie następującej obserwacji: jeśli w danej chwili p > π to stopa oczekiwana będzie rosnąć (czyli d πdt > 0), zaś jeśli p < π to stopa oczekiwana będzie maleć (czyli d πdt < 0).

(43)

Model Friedmana oczekiwań inflacyjnych

W warunkach długotrwałej inflacji ludzie przyjmują pewne

oczekiwania inflacyjne (czyli przewidują przyszłą inflację na pewnym poziomie), co wpływa na ich zachowania rynkowe. Milton Friedman modelował te oczekiwania za pomocą prostego równania

różniczkowego:

d π

dt = a(p − π),

gdzie a ∈ (0, 1], π(t) oznacza oczekiwaną, a p(t) rzeczywistą stopę inflacji w momencie t. Model powstał na bazie następującej

obserwacji: jeśli w danej chwili p > π to stopa oczekiwana będzie rosnąć (czyli d πdt > 0), zaś jeśli p < π to stopa oczekiwana będzie maleć (czyli d πdt < 0).

(44)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ).

W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ): I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

(45)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ). W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ): I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

(46)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ). W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt.

Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

(47)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ). W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności.

Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

(48)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ). W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t).

Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

(49)

Model wzrostu gospodarczego Domara

Zbudujemy najprostszy model wzrostu gospodarczego, w celu oszacowania wpływu wielkości inwestycji (I ) na dochód narodowy (Y ). W modelu zakładamy, że inwestycje wpływają na dochód w przyszłości, ale zależą od dochodu obecnego.

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać

(50)

Model wzrostu gospodarczego Domara

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

Łącząc te równości otrzymujemy:

= sd Υ

dt = sρdK

dt = ρsI .

(51)

Model wzrostu gospodarczego Domara

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

Łącząc te równości otrzymujemy:

dI

dt = sdY

dt = sd Υ dt =

dK

dt = ρsI .

(52)

Model wzrostu gospodarczego Domara

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

Łącząc te równości otrzymujemy:

ρsI .

(53)

Model wzrostu gospodarczego Domara

Inwestycje (I ) można zdefiniować jako stopę przyrostu kapitału (K ):

I = dKdt. Zakładamy, że potencjał produkcyjny (Υ), czyli największy produkt narodowy możliwy do wytworzenia w chwili t jest

proporcjonalny do posiadanego kapitału: Υ(t) = ρK (t), gdzie ρ jest współczynnikiem proporcjonalności. Jeśli potencjał produkcyjny jest w pełni wykorzystany (wtedy mówi się, że gospodarka jest w stanie równowagi), to możemy założyć, że Y (t) = Υ(t). Z kolei, jeśli s jest krańcową skłonnością do oszczędzania dochodu narodowego (czyli częścią tego dochodu, która zostaje zaoszczędzona i może być zainwestowana: zakładamy, że jest stała) to możemy zapisać sdYdt = dIdt.

Łącząc te równości otrzymujemy:

dI

dt = sdY

dt = sd Υ

dt = sρdK

dt = ρsI .

(54)

Ogólne uwagi o rozwiązywaniu równań różniczkowych

Zasadniczo, przygniatającej większości równań różniczkowych nie da się rozwiązać analitycznie, a przynajmniej przedstawić algorytmu, który w skończonej liczbie kroków dawałby ich rozwiązanie.

W wielu wypadkach również rozwiązania numeryczne (czyli przybliżone, obliczone za pomocą komputerów) nie zdają egzaminu ze względu na tzw. silną zależność od warunków początkowych (popularniejszymi, acz mniej ścisłymi określeniami tego zjawiska są chaos lub efekt motyla), czyli fakt, że mała niedokładność pomiaru stanu

początkowego układu, który dane równanie różniczkowe modeluje, może powodować po pewnym czasie dowolnie dużą różnicę między przewidywaniami modelu a stanem rzeczywistym.

W ramach tego wykładu zajmiemy się jedynie kilkoma najprostszymi rodzajami równań różniczkowych, które akurat da się rozwiązać.

(55)

Ogólne uwagi o rozwiązywaniu równań różniczkowych

Zasadniczo, przygniatającej większości równań różniczkowych nie da się rozwiązać analitycznie, a przynajmniej przedstawić algorytmu, który w skończonej liczbie kroków dawałby ich rozwiązanie. W wielu wypadkach również rozwiązania numeryczne (czyli przybliżone, obliczone za pomocą komputerów) nie zdają egzaminu ze względu na tzw. silną zależność od warunków początkowych (popularniejszymi, acz mniej ścisłymi określeniami tego zjawiska są chaos lub efekt motyla), czyli fakt, że mała niedokładność pomiaru stanu

początkowego układu, który dane równanie różniczkowe modeluje, może powodować po pewnym czasie dowolnie dużą różnicę między przewidywaniami modelu a stanem rzeczywistym.

W ramach tego wykładu zajmiemy się jedynie kilkoma najprostszymi rodzajami równań różniczkowych, które akurat da się rozwiązać.

(56)

Ogólne uwagi o rozwiązywaniu równań różniczkowych

Zasadniczo, przygniatającej większości równań różniczkowych nie da się rozwiązać analitycznie, a przynajmniej przedstawić algorytmu, który w skończonej liczbie kroków dawałby ich rozwiązanie. W wielu wypadkach również rozwiązania numeryczne (czyli przybliżone, obliczone za pomocą komputerów) nie zdają egzaminu ze względu na tzw. silną zależność od warunków początkowych (popularniejszymi, acz mniej ścisłymi określeniami tego zjawiska są chaos lub efekt motyla), czyli fakt, że mała niedokładność pomiaru stanu

początkowego układu, który dane równanie różniczkowe modeluje, może powodować po pewnym czasie dowolnie dużą różnicę między

(57)

Najprostsze równania różniczkowe

Najprostsze równania różniczkowe to równania postaci:

dy

dx = f (x ), gdzie f jest pewną funkcją zależną tylko od zmiennej x .

Ich rozwiązania można uzyskać całkując obie strony względem x i są oczywiście wyrażone wzorem:

y (x ) =

Z

f (x )dx + C .

(58)

Najprostsze równania różniczkowe

Najprostsze równania różniczkowe to równania postaci:

dy

dx = f (x ), gdzie f jest pewną funkcją zależną tylko od zmiennej x . Ich rozwiązania można uzyskać całkując obie strony względem x i są oczywiście wyrażone wzorem:

y (x ) =

Z

f (x )dx + C .

(59)

Najprostsze równania różniczkowe - przykład

Ten sposób rozwiązywania stosowaliśmy, zmagając się z II zasadą dynamiki Newtona. Równanie:

dx

dt = at + C1,

przy danych a, C1 ∈ R jest właśnie równaniem tego typu (prawa strona jest funkcją tylko zmiennej niezależnej t), więc rozwiązaliśmy je całkując obustronnie po t i otrzymując:

x (t) =

a

2t2+ C1t + C2.

(60)

Najprostsze równania różniczkowe - przykład

Ten sposób rozwiązywania stosowaliśmy, zmagając się z II zasadą dynamiki Newtona. Równanie:

dx

dt = at + C1,

przy danych a, C1 ∈ R jest właśnie równaniem tego typu (prawa strona jest funkcją tylko zmiennej niezależnej t), więc rozwiązaliśmy je całkując obustronnie po t i otrzymując:

x (t) = a

2t2+ C1t + C2.

(61)

Równania o zmiennych rozdzielonych

Niemal równie proste są równania o rozdzielonych zmiennych, postaci:

f (y )dy

dx = g (x )

Rozwiązujemy je „mnożąc” obie strony przez dx i następnie całkując: f (y )dy

dx = g (x ) =⇒ [f (y )dy = g (x )dx ] =⇒

Z

f (y )dy =

Z

g (x )dx .

Oczywiście, zapis w nawiasie kwadratowym powyżej jest zupełnie niepoprawny (przynajmniej według naszej wiedzy), bo samo dx bądź dy nie ma zdefiniowanego znaczenia.

(62)

Równania o zmiennych rozdzielonych

Niemal równie proste są równania o rozdzielonych zmiennych, postaci:

f (y )dy

dx = g (x )

Rozwiązujemy je „mnożąc” obie strony przez dx i następnie całkując:

f (y )dy

dx = g (x ) =⇒ [f (y )dy = g (x )dx ] =⇒

Z

f (y )dy =

Z

g (x )dx .

Oczywiście, zapis w nawiasie kwadratowym powyżej jest zupełnie niepoprawny (przynajmniej według naszej wiedzy), bo samo dx bądź dy nie ma zdefiniowanego znaczenia.

(63)

Równania o zmiennych rozdzielonych

Niemal równie proste są równania o rozdzielonych zmiennych, postaci:

f (y )dy

dx = g (x )

Rozwiązujemy je „mnożąc” obie strony przez dx i następnie całkując:

f (y )dy

dx = g (x ) =⇒ [f (y )dy = g (x )dx ] =⇒

Z

f (y )dy =

Z

g (x )dx .

Oczywiście, zapis w nawiasie kwadratowym powyżej jest zupełnie niepoprawny (przynajmniej według naszej wiedzy), bo samo dx bądź dy nie ma zdefiniowanego znaczenia.

(64)

Równania o zmiennych rozdzielonych

f (y )dy

dx = g (x ) =⇒ [f (y )dy = g (x )dx ] =⇒

Z

f (y )dy =

Z

g (x )dx . Oczywiście, zapis w nawiasie kwadratowym powyżej jest zupełnie

niepoprawny (przynajmniej według naszej wiedzy), bo samo dx bądź dy nie ma zdefiniowanego znaczenia. To „mnożenie” jest tylko mnemotechniczną sztuczką, dzięki której możemy łatwo zapamiętać sposób rozwiązywania równania o zmiennych rozdzielonych, a nie jakimś znanym nam działaniem. Dlatego, zanim

ćwiczeniowiec/wykładowca się zorientuje, że robimy z tym równaniem

(65)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

Dla stałego poziomu inflacji p, model Friedmana oczekiwań inflacyjnych jest zadany równaniem różniczkowym o zmiennych rozdzielonych:

d π

dt = a(p − π),

gdzie a, p > 0. By to zauważyć, wystarczy przenieść na lewą stronę wszystko, co jest zależne od π np.:

1 p − π

d π dt = a.

(66)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

Dla stałego poziomu inflacji p, model Friedmana oczekiwań inflacyjnych jest zadany równaniem różniczkowym o zmiennych rozdzielonych:

d π

dt = a(p − π),

gdzie a, p > 0. By to zauważyć, wystarczy przenieść na lewą stronę wszystko, co jest zależne od π np.:

(67)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

1 p − π

d π dt = a.

Zgodnie z zasadą rozwiązywania równań o zmiennych rozdzielonych możemy nieformalnie zapisać postać równoważną jako:

[p−π1 d π = adt], a formalnie jako:

Z 1

p − πd π =

Z

adt ⇒ − ln |p − π| + C1 = at + C2

⇒ ln |p − π| = −at + C3

(68)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

1 p − π

d π dt = a.

Zgodnie z zasadą rozwiązywania równań o zmiennych rozdzielonych możemy nieformalnie zapisać postać równoważną jako:

[p−π1 d π = adt], a formalnie jako:

Z 1

p − πd π =

Z

adt ⇒

− ln |p − π| + C1 = at + C2

⇒ ln |p − π| = −at + C3

(69)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

1 p − π

d π dt = a.

Zgodnie z zasadą rozwiązywania równań o zmiennych rozdzielonych możemy nieformalnie zapisać postać równoważną jako:

[p−π1 d π = adt], a formalnie jako:

Z 1

p − πd π =

Z

adt ⇒ − ln |p − π| + C1 = at + C2

⇒ ln |p − π| = −at + C3

(70)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3

⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim

t→∞π(t) = p, skąd można wyciągnąć wniosek, że jeśli stopa inflacji jest stała przez dłuższy czas, to oczekiwania inflacyjne społeczeństwa zaczną się pokrywać z faktyczną stopą inflacji (o ile model Friedmana jest prawdziwy).

(71)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3 ⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim

t→∞π(t) = p, skąd można wyciągnąć wniosek, że jeśli stopa inflacji jest stała przez dłuższy czas, to oczekiwania inflacyjne społeczeństwa zaczną się pokrywać z faktyczną stopą inflacji (o ile model Friedmana jest prawdziwy).

(72)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3 ⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim

t→∞π(t) = p, skąd można wyciągnąć wniosek, że jeśli stopa inflacji jest stała przez dłuższy czas, to oczekiwania inflacyjne społeczeństwa zaczną się pokrywać z faktyczną stopą inflacji (o ile model Friedmana jest prawdziwy).

(73)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3 ⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim

t→∞π(t) =

p, skąd można wyciągnąć wniosek, że jeśli stopa inflacji jest stała przez dłuższy czas, to oczekiwania inflacyjne społeczeństwa zaczną się pokrywać z faktyczną stopą inflacji (o ile model Friedmana jest prawdziwy).

(74)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3 ⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim π(t) = p, skąd można wyciągnąć wniosek,

(75)

Równania o zmiennych rozdzielonych - przykład (model Friedmana)

ln |p − π| = −at + C3 ⇒ |p − π| = e−at+C3 = eC3· e−at = Ce−at.

I ostatecznie:

π(t) = p −Ce−at(moduł mogliśmy opuścić, bo C może być ujemne) .

Warto zauważyć, że lim

t→∞π(t) = p, skąd można wyciągnąć wniosek, że jeśli stopa inflacji jest stała przez dłuższy czas, to oczekiwania inflacyjne społeczeństwa zaczną się pokrywać z faktyczną stopą

(76)

Stałe w rozwiązaniach równań różniczkowych - uwaga techniczna

Zauważmy, że w rozwiązanym przykładzie pojawiło się wiele stałych:

C1, C2, C3 itd., a na końcu miała znaczenie tylko stała C , z tych stałych powstająca.

Jakkolwiek ten zapis formalnie jest jedynym poprawnym, to w praktyce równań różniczkowych przyjęło się oznaczać wszystkie stałe jako C np. w zapisie wyniku eC możemy zastąpić przez C bez utraty ogólności rozwiązania. Podobnie, gdy obliczamy całki po obydwu stronach równości (jak w procedurze rozwiązywania równań o zmiennych rozdzielonych), wystarczy by stałą zapisać po jednej stronie (bo i tak można jedną z nich odjąć stronami). Do tej konwencji odtąd będę się stosować.

(77)

Stałe w rozwiązaniach równań różniczkowych - uwaga techniczna

Zauważmy, że w rozwiązanym przykładzie pojawiło się wiele stałych:

C1, C2, C3 itd., a na końcu miała znaczenie tylko stała C , z tych stałych powstająca. Jakkolwiek ten zapis formalnie jest jedynym poprawnym, to w praktyce równań różniczkowych przyjęło się oznaczać wszystkie stałe jako C np. w zapisie wyniku eC możemy zastąpić przez C bez utraty ogólności rozwiązania.

Podobnie, gdy obliczamy całki po obydwu stronach równości (jak w procedurze rozwiązywania równań o zmiennych rozdzielonych), wystarczy by stałą zapisać po jednej stronie (bo i tak można jedną z nich odjąć stronami). Do tej konwencji odtąd będę się stosować.

(78)

Stałe w rozwiązaniach równań różniczkowych - uwaga techniczna

Zauważmy, że w rozwiązanym przykładzie pojawiło się wiele stałych:

C1, C2, C3 itd., a na końcu miała znaczenie tylko stała C , z tych stałych powstająca. Jakkolwiek ten zapis formalnie jest jedynym poprawnym, to w praktyce równań różniczkowych przyjęło się oznaczać wszystkie stałe jako C np. w zapisie wyniku eC możemy zastąpić przez C bez utraty ogólności rozwiązania. Podobnie, gdy obliczamy całki po obydwu stronach równości (jak w procedurze rozwiązywania równań o zmiennych rozdzielonych), wystarczy by stałą zapisać po jednej stronie (bo i tak można jedną z nich odjąć

Do tej konwencji odtąd będę się stosować.

Cytaty

Powiązane dokumenty

Metodę przewidywań możemy stosować w przypadku równań o stałych współczynnikach, gdy wyraz wolny ma jedną z postaci przedstawionych w kolumnie 2 tabeli zamieszczonej w

Podobnie jak w przypadku równań pierwszego i drugiego rzędu, rozwiązywanie równania liniowego niejednorodnego rzędu n-tego polega na wyznaczeniu CORJ, a następnie zastosowaniu

Jeżeli dodatkowo są liniowo niezależne (tworzą układ fundamentalny rozwiązań), to ich kombinacja liniowa jest rozwiązaniem ogólnym układu równań.. Układy

Dla wyznaczenia wartości rozwiązania w punktach odległych od punktu startowego wykorzystuje się procedurę iteracyjną, w której wielokrotnie wykorzystuje się tę samą

Jeżeli powyższy warunek nie jest spełniony to szereg jest rozbieżny. Nie jest on warunkiem wystarczającym zbieżności, tzn. może być spełniony przez szereg rozbieżny..

Szereg funkcyjny jest zbieżny punktowo do funkcji S(x) na zbiorze X, jeżeli ciąg jego sum częściowych jest zbieżny punktowo do funkcji S(x) na tym zbiorze.. Zbiór

Warunki wystarczające na to by suma szeregu Fouriera była równa funkcji, na podstawie której szereg został skonstruowany, nazywane są warunkami Dirichleta..

nazywamy okrąg styczny do krzywej w tym punkcie, leżący po tej samej stronie co krzywa i mający promień równy promieniowi