Rozdíly

Zde můžete vidět rozdíly mezi vybranou verzí a aktuální verzí dané stránky.

Odkaz na výstup diff

temata:17a-matematicka_analyza:diferencialni_rovnice [2011/03/03 11:22]
vagabund
temata:17a-matematicka_analyza:diferencialni_rovnice [2011/04/21 11:46] (aktuální)
vagabund
Řádek 32: Řádek 32:
 **Actual Solution** - GS s dosazením IV **Actual Solution** - GS s dosazením IV
  
-**Implicitní/Explicitní řešení** - y = f(t)/0 = f(y,t)+**Explicitní/Implicitní řešení** - y = f(t)/0 = f(y,t)
  
 </box> </box>
Řádek 40: Řádek 40:
 == Lineární rovnice == == Lineární rovnice ==
  
-<box round 90% green|**Lineární rovnice**>+<box round 90% red|**Lineární rovnice**>
 Tvar: <m>{dy}/{dt} + p(t)y = g(t)</m> Tvar: <m>{dy}/{dt} + p(t)y = g(t)</m>
  
 Řešení: <m>y(t) = {\int{}{}{\mu(t)g(t)dt} + c}/{\mu(t)}</m>, <m>\mu(t) = e^{\int{}{}{p(t)dt}}</m> Řešení: <m>y(t) = {\int{}{}{\mu(t)g(t)dt} + c}/{\mu(t)}</m>, <m>\mu(t) = e^{\int{}{}{p(t)dt}}</m>
 +
 +Řešení2: 
 +
 +1. určíme řešení rovnice <m>y\prime + p(t)y = 0</m> ve tvaru <m>y(t) = cF(t)</m>
 +
 +2. <m>y(t) = cF(t)</m> přepíšeme na <m>y(t) = c(t)F(t)</m> a dosadíme do původní rovnice
 +
 +3. spočíme c(t), dosadíme do <m>y(t) = c(t)F(t)</m>
 +
 +</box>
 +
 +== příklad ==
 +<box round 90% green|**příklad**>
 +
 +<m>v\prime = 9.8 - 0.196v</m>
 +
 +upravíme do správného tvaru
 +
 +<m>v\prime + 0.196v = 9.8</m>
 +
 +<m>\mu(t) = e^{\int{}{}{0.196dt}} = e^{0.196t}</m>
 +
 +<m>y(t) = {\int{}{}{e^{0.196t}.9.8dt} + c}/{e^{0.196t}} = {{9.8}/{0.196}e^{0.196t} + c}/{e^{0.196t}} = 50 + ce^{-0.196t}</m>
 +
 +\\
 +\\
 +
 +<m>v\prime + 0.196v = 0</m>
 +
 +<m>{dv}/{dt} = -0.196v</m>
 +
 +<m>{dv}/{v} = -0.196dt</m>
 +
 +<m>ln|v| = -0.196t + c</m>
 +
 +<m>|v| = e^{-0.196t + c}</m>
 +
 +<m>|v| = e^{-0.196t}e^c</m>
 +
 +<m>|v| = e^{-0.196t}c</m>
 +
 +**pro kladné v**:
 +
 +<m>v = ce^{-0.196t} \doubleright v = c(t)e^{-0.196t}</m>
 +
 +<m>v\prime = c(t)\prime e^{-0.196t} -0.196c(t)e^{-0.196t}</m>
 +
 +dosadíme:
 +
 +<m>c(t)\prime e^{-0.196t} -0.196c(t)e^{-0.196t} + 0.196.c(t)e^{-0.196t} = 9.8</m>
 +
 +<m>c(t)\prime e^{-0.196t} = 9.8</m>
 +
 +<m>c(t)\prime = 9.8e^{+0.196t}</m>
 +
 +<m>c(t) = 50e^{+0.196t} + k</m>
 +
 +<m>v = c(t)e^{-0.196t} = (50e^{+0.196t} + k)e^{-0.196t} = 50 + ke^{-0.196t}</m>
 +
 +**pro záporné v**:
 +
 +<m>-v = ce^{-0.196t} \doubleright v = c(t)(-1)e^{-0.196t}</m>
 +
 +<m>v\prime = c(t)\prime(-1)e^{-0.196t} +0.196c(t)e^{-0.196t}</m>
 +
 +dosadíme:
 +
 +<m>-c(t)\primee^{-0.196t} + 0.196c(t)e^{-0.196t} - 0.196.c(t)e^{-0.196t} = 9.8</m>
 +
 +<m>-c(t)\prime e^{-0.196t} = 9.8</m>
 +
 +<m>c(t)\prime = -9.8e^{+0.196t}</m>
 +
 +<m>c(t) = -50e^{+0.196t} + k</m>
 +
 +<m>v = -c(t)e^{-0.196t} = -(-50e^{+0.196t} - k)e^{-0.196t} = 50 - ke^{-0.196t} = 50 + ke^{-0.196t}</m>
 +
 </box> </box>
  
 == Separabilní dif. rovnice (Separable ...) == == Separabilní dif. rovnice (Separable ...) ==
  
-<box round 90% green|**Separabilní dif. rovnice**>+<box round 90% red|**Separabilní dif. rovnice**>
 Tvar: <m>N(u){dy}/{dx} = M(x)</m> Tvar: <m>N(u){dy}/{dx} = M(x)</m>
  
 Řešení: <m>N(u){dy}/{dx} = M(x) \doubleleftright N(u)dy = M(x)dx \doubleright \int{}{}{N(u)dy} = \int{}{}{M(x)dx}</m> Řešení: <m>N(u){dy}/{dx} = M(x) \doubleleftright N(u)dy = M(x)dx \doubleright \int{}{}{N(u)dy} = \int{}{}{M(x)dx}</m>
 +</box>
 +
 +== příklad ==
 +<box round 90% green|**příklad**>
 +
 +<m>y\prime = {3x^2 + 4x - 4}/{2y - 4}, y(1) = 3</m>
 +
 +<m>{dy}/{dx} = {3x^2 + 4x - 4}/{2y - 4}</m>
 +
 +<m>{2y - 4}{dy} = {3x^2 + 4x - 4}{dx}</m>
 +
 +<m>y^2 - 4y = x^3 + 2x^2 - 4x + k</m>
 +
 +partikularni reseni:
 +
 +<m>9 - 12 = 1 + 2 - 4 + k</m>
 +
 +<m>k = -2</m>
 +
 +<m>y^2 - 4y = x^3 + 2x^2 - 4x - 2</m>
 +
 </box> </box>
  
Řádek 64: Řádek 162:
 </box> </box>
  
-<box round 90% green|**Second Order ODE, Homogenous, constant coeff.**>+<box round 90% red|**Second Order ODE, Homogenous, constant coeff.**>
 Tvar: <m>ay\prime\prime + by\prime + cy = 0</m> Tvar: <m>ay\prime\prime + by\prime + cy = 0</m>
  
Řádek 87: Řádek 185:
 </box> </box>
  
-<box round 90% green|**Second Order ODE, Nonhomogenous, constant coeff.**>+<box round 90% red|**Second Order ODE, Nonhomogenous, constant coeff.**>
  
 Určíme partikulární řešení. Určíme partikulární řešení.
Řádek 98: Řádek 196:
 | <m>asin(\beta t)</m> | <m>Acos(\beta t) + Bsin(\beta t)</m> | | <m>asin(\beta t)</m> | <m>Acos(\beta t) + Bsin(\beta t)</m> |
 | <m>acos(\beta t) + asin(\beta t)</m> | <m>Acos(\beta t) + Bsin(\beta t)</m> | | <m>acos(\beta t) + asin(\beta t)</m> | <m>Acos(\beta t) + Bsin(\beta t)</m> |
 +| <m>3t^2 + 4t + 6</m> | <m>At^2 + Bt + C</m> |
  
 Dopočítáme homogenní a dostaneme: <m>y(t) = Y_P(t) + Y_H(t)</m> Dopočítáme homogenní a dostaneme: <m>y(t) = Y_P(t) + Y_H(t)</m>
Řádek 103: Řádek 202:
 </box> </box>
  
 +== příklad ==
 +<box round 90% green|**příklad**>
 +
 +<m>y\prime\prime + 3y\prime - 10y = 0</m>
 +
 +charakteristicka rovnice:
 +
 +<m>{\lambda}^2 + 3{\lambda} - 10 = 0</m>
 +
 +<m>({\lambda} + 5)({\lambda} - 2) = 0</m>
 +
 +<m>{\lambda_1} = 2</m>
 +
 +<m>{\lambda_2} = -5</m>
 +
 +reseni:
 +
 +<m>y = c_1e^{2t} + c_2e^{-5t}</m>
 +
 +//
 +//
 +
 +<m>y\prime\prime - 4y\prime + 4y = 0</m>
 +
 +charakteristicka rovnice:
 +
 +<m>{\lambda}^2 - 4{\lambda} + 4 = 0</m>
 +
 +<m>{({\lambda} -2)}^2 = 0</m>
 +
 +<m>{\lambda} = 2</m>
 +
 +reseni:
 +
 +<m>y = c_1e^{2t} + c_2te^{2t}</m>
 +
 +//
 +//
 +
 +<m>y\prime\prime - 4y\prime + 9y = 0</m>
 +
 +charakteristicka rovnice:
 +
 +<m>{\lambda}^2 - 4{\lambda} + 9 = 0</m>
 +
 +<m>D = 16 - 36 = - 20</m>
 +
 +<m>{\lambda}_{1,2} = {4 \pm 2i\sqrt{5}}/2</m>
 +
 +<m>{\lambda}_{1,2} = 2 \pm i\sqrt{5}</m>
 +
 +reseni:
 +
 +<m>y = c_1e^{2t}sin(\sqrt{5}t) + c_2e^{2t}cos(\sqrt{5}t)</m>
 +</box>
 +
 +<box round 90% green|**příklad**>
 +
 +<m>y\prime\prime - 4y\prime - 12y = 3e^{5t}</m>
 +
 +1.) partikularni reseni hadanim
 +
 +<m>Y_p = Ae^{5t}</m>
 +
 +<m>Y_p\prime = 5Ae^{5t}</m>
 +
 +<m>Y_p\prime\prime = 25Ae^{5t}</m>
 +
 +<m>25Ae^{5t} - 20Ae^{5t} - 12Ae^{5t} = 3e^{5t}</m>
 +
 +<m>-7Ae^{5t} = 3e^{5t}</m>
 +
 +<m>A = -3/7</m>
 +
 +<m>Y_p = -3/7e^{5t}</m>
 +
 +2.) homogenní řešení
 +
 +<m>{\lambda}^2 - 4{\lambda} - 12 = 0</m>
 +
 +<m>(\lambda - 6)(\lambda + 2) = 0</m>
 +
 +<m>\lambda_1 = 6</m>
 +
 +<m>\lambda_2 = -2</m>
 +
 +<m>Y_h = c_1e^{6t} + c_2e(-2t)</m>
 +
 +3.) <m>Y = Y_h + Y_p</m>
 +
 +<m>Y = c_1e^{6t} + c_2e(-2t) + -3/7e^{5t}</m>
 +
 +
 +</box>
 +
 +===== Potvrzení =====
 +
 +<doodle single login|17a3>
 +^ OK ^ !!! ^
 +</doodle>
 +
 +{{tag>vagabund obyc_diferencialni_rovnice linearni_rovnice}}
 +
 +~~DISCUSSION~~
temata/17a-matematicka_analyza/diferencialni_rovnice.1299147751.txt.gz · Poslední úprava: 2011/03/03 11:22 autor: vagabund
Recent changes RSS feed Debian Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki