ODE
Ecuaciones diferenciales ordinarias Una Ecuación Diferencial Ordinaria (EDO) es una ecuación que involucra una función desconocida de una sola variable independiente y sus derivadas. Estas ecuaciones representan relaciones entre una función y sus cambios o tasas de cambio, y son fundamentales en el modelado de fenómenos
Clasificación Orden La mayor derivada
Grado El exponente de la mayor derivada
Forma General de una EDO F ( x , y , y ′ , y ′ ′ , … , y ( n ) ) = 0 F(x,y,y',y'',…,y^{(n)})=0 F ( x , y , y ′ , y ′′ , … , y ( n ) ) = 0 x x x es la variable independiente.y = y ( x ) y=y(x) y = y ( x ) es la función desconocida o variable dependiente.y ′ , y ′ ′ , . . . , y ( n ) y',y'', ..., y^{(n)} y ′ , y ′′ , ... , y ( n ) son las derivadas de y y y respecto a x x x hasta el 𝑛-ésimo orden.EDO de primer orden M ( x , y ) d x + N ( x , y ) d y = 0 M(x,y)dx+N(x,y)dy=0 M ( x , y ) d x + N ( x , y ) d y = 0 M ( x , y ) M(x,y) M ( x , y ) y N ( x , y ) N(x,y) N ( x , y ) son llamados coeficientes de la EDO, que pueden depender de 𝑥, de 𝑦, o de ambos.
EDO de variables separables Para resolver una EDO de variables separables:
M ( x ) d x + N ( y ) d y = 0 M(x)dx+N(y)dy=0 M ( x ) d x + N ( y ) d y = 0 N ( y ) d y = − M ( x ) d x N(y)dy = -M(x)dx N ( y ) d y = − M ( x ) d x ∫ N ( y ) d y = − ∫ M ( x ) d x \int N(y)dy = -\int M(x)dx ∫ N ( y ) d y = − ∫ M ( x ) d x n ( y ) = m ( x ) n(y)=m(x) n ( y ) = m ( x ) y = n − 1 ( m ( x ) ) y=n^{-1}(m(x)) y = n − 1 ( m ( x )) Convertir una EDO de primer orden a una EDO de variables separables:
Factorización M ( x ) d x + N ( y ) d y = 0 M(x)dx+N(y)dy=0 M ( x ) d x + N ( y ) d y = 0 M ( x ) d x N ( y ) d y = − 1 \frac{M(x)dx}{N(y)dy}=-1 N ( y ) d y M ( x ) d x = − 1 d x d y = − N ( y ) M ( x ) \frac{dx}{dy}=-\frac{N(y)}{M(x)} d y d x = − M ( x ) N ( y )
d x d y = f ( x , y ) \frac{dx}{dy}=f(x,y) d y d x = f ( x , y ) Se debe verificar que f ( x , y ) f(x,y) f ( x , y ) se pueda factorizar f ( x , y ) = g ( x ) h ( y ) f(x,y)=g(x)h(y) f ( x , y ) = g ( x ) h ( y )
Si f ( x , y ) f(x,y) f ( x , y ) es factorizable entonces de cumple que:
∃ x 0 , y 0 / f ( x 0 , y 0 ) ≠ 0 ⟹ f ( x 0 , y 0 ) f ( x , y ) = f ( x 0 , y ) f ( x , y 0 ) \exists {x_0, y_0} / f(x_0,y_0) \neq 0 \implies f(x_0,y_0)f(x,y)=f(x_0,y)f(x,y_0) ∃ x 0 , y 0 / f ( x 0 , y 0 ) = 0 ⟹ f ( x 0 , y 0 ) f ( x , y ) = f ( x 0 , y ) f ( x , y 0 ) Resolviendo:
d x d y = f ( x , y ) \frac{dx}{dy}=f(x,y) d y d x = f ( x , y ) d x d y = g ( x ) h ( y ) \frac{dx}{dy}=g(x)h(y) d y d x = g ( x ) h ( y ) g ( x ) − 1 d x = h ( y ) d y g(x)^{-1}dx=h(y)dy g ( x ) − 1 d x = h ( y ) d y Funciones homogéneas f ( λ x , λ y ) = λ k f ( x , y ) f(\lambda x,\lambda y)=\lambda^{k}f(x,y) f ( λ x , λ y ) = λ k f ( x , y ) f ( x , y ) f(x,y) f ( x , y ) es una función homogénea de grado k k k
λ = 1 x \lambda=\frac{1}{x} λ = x 1 f ( x x , y x ) = 1 x k f ( x , y ) f(\frac{x}{x}, \frac{y}{x})=\frac{1}{x}^{k}f(x,y) f ( x x , x y ) = x 1 k f ( x , y ) f ( 1 , y x ) = x − k f ( x , y ) f(1, \frac{y}{x})=x^{-k}f(x,y) f ( 1 , x y ) = x − k f ( x , y )
x k f ( 1 , y x ) = f ( x , y ) x^{k} f(1, \frac{y}{x})=f(x,y) x k f ( 1 , x y ) = f ( x , y ) y k f ( x y , 1 ) = f ( x , y ) y^{k} f( \frac{x}{y},1)=f(x,y) y k f ( y x , 1 ) = f ( x , y ) M ( x , y ) d x + N ( x , y ) d y = 0 M(x,y)dx+N(x,y)dy=0 M ( x , y ) d x + N ( x , y ) d y = 0 Será una EDO de variables separables y a la vez una EDO homogénea, si sus coefiecientes son funciones homogéneas y del mismo grado de homogeneidad
M más fácil que N y k M ( x y , 1 ) d x + y k N ( x y , 1 ) d y = 0 y^{k} M(\frac{x}{y},1)dx+y^{k} N(\frac{x}{y},1)dy=0 y k M ( y x , 1 ) d x + y k N ( y x , 1 ) d y = 0 M ( x y , 1 ) d x + N ( x y , 1 ) d y = 0 M(\frac{x}{y},1)dx+ N(\frac{x}{y},1)dy=0 M ( y x , 1 ) d x + N ( y x , 1 ) d y = 0 v = x y v=\frac{x}{y} v = y x y v = x yv=x y v = x y d v + v d y = d x ydv + vdy= dx y d v + v d y = d x M ( v , 1 ) ( y d v + v d y ) + N ( v , 1 ) d y = 0 M(v,1)(ydv + vdy)+ N(v,1)dy=0 M ( v , 1 ) ( y d v + v d y ) + N ( v , 1 ) d y = 0 N más fácil que M x k M ( 1 , y x ) d x + x k N ( 1 , y x ) d y = 0 x^{k} M(1,\frac{y}{x})dx+x^{k} N(1,\frac{y}{x})dy=0 x k M ( 1 , x y ) d x + x k N ( 1 , x y ) d y = 0 M ( 1 , y x ) d x + N ( 1 , y x ) d y = 0 M(1,\frac{y}{x})dx+ N(1,\frac{y}{x})dy=0 M ( 1 , x y ) d x + N ( 1 , x y ) d y = 0 v = y x v=\frac{y}{x} v = x y y v = y yv=y y v = y x d v + v d x = d y xdv+vdx=dy x d v + v d x = d y M ( 1 , v ) d x + N ( 1 , v ) ( x d v + v d x ) = 0 M(1,v)dx+ N(1,v)(xdv+vdx)=0 M ( 1 , v ) d x + N ( 1 , v ) ( x d v + v d x ) = 0 Ecuaciones con coeficientes lineales ( A x + B y + C ⏟ L i n e a l L 1 : ) d x + ( D x + E y + F ⏟ L i n e a l L 2 : ) d y = 0 (\overset{L1:}{\underset{Lineal}{\underbrace{Ax+By+C}}})dx + (\overset{L2:}{\underset{Lineal}{\underbrace{Dx+Ey+F}}})dy =0 ( L in e a l A x + B y + C L 1 : ) d x + ( L in e a l D x + E y + F L 2 : ) d y = 0 Se debe hacer un traslado al origen la intersección de las rectas con un c.v. y dará como resultado una ecuación homogénea
EDO Exacta d F = ∂ F ∂ x d x + ∂ F ∂ y d y = 0 dF = \frac{\partial F}{\partial x} dx + \frac{\partial F}{\partial y} dy=0 d F = ∂ x ∂ F d x + ∂ y ∂ F d y = 0 ∂ ∂ F ∂ x ∂ y = ∂ ∂ F ∂ y ∂ x \frac{\partial \frac{\partial F}{\partial x}}{\partial y} = \frac{\partial \frac{\partial F}{\partial y}}{\partial x} ∂ y ∂ ∂ x ∂ F = ∂ x ∂ ∂ y ∂ F d F = M ( x , y ) d x + N ( x , y ) d y = 0 dF = M(x,y) dx + N(x,y) dy=0 d F = M ( x , y ) d x + N ( x , y ) d y = 0 ∂ M ( x , y ) ∂ y = ∂ N ( x , y ) ∂ x \frac{\partial M(x,y)}{\partial y} = \frac{\partial N(x,y)}{\partial x} ∂ y ∂ M ( x , y ) = ∂ x ∂ N ( x , y ) ∂ F ∂ x = M ( x , y ) \frac{\partial F}{\partial x} = M(x,y) ∂ x ∂ F = M ( x , y ) ∂ F = M ( x , y ) ∂ x \partial F = M(x,y) \partial x ∂ F = M ( x , y ) ∂ x ∫ ∂ F = ∫ M ( x , y ) ∂ x \int \partial F =\int M(x,y) \partial x ∫ ∂ F = ∫ M ( x , y ) ∂ x F = ∫ M ( x , y ) ∂ x F =\int M(x,y) \partial x F = ∫ M ( x , y ) ∂ x ∂ F ∂ y = ∂ ∫ M ( x , y ) ∂ x ∂ y = N ( x , y ) \frac{\partial F}{\partial y} = \frac{\partial \int M(x,y) \partial x}{\partial y} = N(x,y) ∂ y ∂ F = ∂ y ∂ ∫ M ( x , y ) ∂ x = N ( x , y ) Factor integrante M ( x , y ) d x + N ( x , y ) d y = 0 M(x,y) dx + N(x,y) dy=0 M ( x , y ) d x + N ( x , y ) d y = 0 ∂ M ( x , y ) ∂ y ≠ ∂ N ( x , y ) ∂ x \frac{\partial M(x,y)}{\partial y} \neq \frac{\partial N(x,y)}{\partial x} ∂ y ∂ M ( x , y ) = ∂ x ∂ N ( x , y ) u ( x , y ) M ( x , y ) d x + u ( x , y ) N ( x , y ) d y = 0 u(x,y)M(x,y) dx + u(x,y)N(x,y) dy=0 u ( x , y ) M ( x , y ) d x + u ( x , y ) N ( x , y ) d y = 0 ∂ ( u ( x , y ) M ( x , y ) ) ∂ y = ∂ ( u ( x , y ) N ( x , y ) ) ∂ x \frac{\partial(u(x,y)M(x,y))}{\partial y} = \frac{\partial (u(x,y)N(x,y))}{\partial x} ∂ y ∂ ( u ( x , y ) M ( x , y )) = ∂ x ∂ ( u ( x , y ) N ( x , y )) Cómo hallar u(x,y) ∂ ( u M ) ∂ y = ∂ ( u N ) ∂ x \frac{\partial(uM)}{\partial y} = \frac{\partial (uN)}{\partial x} ∂ y ∂ ( u M ) = ∂ x ∂ ( u N ) u ∂ M ∂ y + M ∂ u ∂ y = u ∂ N ∂ x + N ∂ u ∂ x u\frac{\partial M}{\partial y}+M\frac{\partial u}{\partial y} = u \frac{\partial N}{\partial x} + N \frac{\partial u}{\partial x} u ∂ y ∂ M + M ∂ y ∂ u = u ∂ x ∂ N + N ∂ x ∂ u u ( ∂ M ∂ y − ∂ N ∂ x ) = N ∂ u ∂ x − M ∂ u ∂ y u(\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} ) = N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} u ( ∂ y ∂ M − ∂ x ∂ N ) = N ∂ x ∂ u − M ∂ y ∂ u ∂ M ∂ y − ∂ N ∂ x = N ∂ u ∂ x − M ∂ u ∂ y u \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} }{u} ∂ y ∂ M − ∂ x ∂ N = u N ∂ x ∂ u − M ∂ y ∂ u Si u ( x ) u(x) u ( x ) :
∂ M ∂ y − ∂ N ∂ x = N ∂ u ∂ x − M ∂ u ∂ y u \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N \frac{\partial u}{\partial x} -\cancel{M\frac{\partial u}{\partial y}} }{u} ∂ y ∂ M − ∂ x ∂ N = u N ∂ x ∂ u − M ∂ y ∂ u ∂ M ∂ y − ∂ N ∂ x N = ∂ u ∂ x u = ∂ ln ∣ u ∣ ∂ x \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{N} = \frac{ \frac{\partial u}{\partial x} }{u}= \frac{\partial \ln{|u|}}{\partial x} N ∂ y ∂ M − ∂ x ∂ N = u ∂ x ∂ u = ∂ x ∂ ln ∣ u ∣ ∫ ∂ M ∂ y − ∂ N ∂ x N ∂ x = ∫ ∂ ln ∣ u ∣ = ln ∣ u ∣ \int \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{N} \partial x =\int \partial \ln{|u|}= \ln{|u|} ∫ N ∂ y ∂ M − ∂ x ∂ N ∂ x = ∫ ∂ ln ∣ u ∣ = ln ∣ u ∣ e ∫ ∂ M ∂ y − ∂ N ∂ x N ∂ x = u e^{\int \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{N} \partial x }=u e ∫ N ∂ y ∂ M − ∂ x ∂ N ∂ x = u Si u ( y ) u(y) u ( y ) :
∂ M ∂ y − ∂ N ∂ x = N ∂ u ∂ x − M ∂ u ∂ y u \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{\cancel{N \frac{\partial u}{\partial x}} -M\frac{\partial u}{\partial y} }{u} ∂ y ∂ M − ∂ x ∂ N = u N ∂ x ∂ u − M ∂ y ∂ u ∂ M ∂ y − ∂ N ∂ x − M = ∂ u ∂ y u = ∂ ln ∣ u ∣ ∂ y \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{-M} = \frac{ \frac{\partial u}{\partial y} }{u}= \frac{\partial \ln{|u|}}{\partial y} − M ∂ y ∂ M − ∂ x ∂ N = u ∂ y ∂ u = ∂ y ∂ ln ∣ u ∣ Si u ( x , y ) = x m y n u(x,y)= x^m y^n u ( x , y ) = x m y n :
∂ M ∂ y − ∂ N ∂ x = N ∂ u ∂ x − M ∂ u ∂ y u \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} }{u} ∂ y ∂ M − ∂ x ∂ N = u N ∂ x ∂ u − M ∂ y ∂ u ∂ M ∂ y − ∂ N ∂ x = N y n m x m − 1 − M x m n y n − 1 x m y n \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N y^n m x^{m-1} -M x^m n y^{n-1}}{x^m y^n} ∂ y ∂ M − ∂ x ∂ N = x m y n N y n m x m − 1 − M x m n y n − 1 ∂ M ∂ y − ∂ N ∂ x = N m x − 1 − M n y − 1 \frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = N m x^{-1} -M n y^{-1} ∂ y ∂ M − ∂ x ∂ N = N m x − 1 − M n y − 1 Riccati Es una edo de la forma:
d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+ P(x)y = Q(x) d x d y + P ( x ) y = Q ( x ) que se le multiplica por un factor integrante
u = e ∫ P ( x ) ∂ x u = e^{\int P(x) \partial x} u = e ∫ P ( x ) ∂ x d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+ P(x)y = Q(x) d x d y + P ( x ) y = Q ( x ) Obteniendo el factor integrante para Ricatti
d y + P ( x ) y d x = Q ( x ) d x dy+ P(x)ydx = Q(x)dx d y + P ( x ) y d x = Q ( x ) d x d y + ( P ( x ) y − Q ( x ) ) d x = 0 dy+ (P(x)y-Q(x))dx =0 d y + ( P ( x ) y − Q ( x )) d x = 0 M ( x , y ) = P ( x ) y − Q ( x ) M(x,y)= P(x)y-Q(x) M ( x , y ) = P ( x ) y − Q ( x ) M ( x , y ) d x + d y = 0 M(x,y)dx+dy = 0 M ( x , y ) d x + d y = 0 Hallamos u:
d y + M ( x , y ) d x = 0 dy+ M(x,y)dx = 0 d y + M ( x , y ) d x = 0 ∂ ( u 1 ) ∂ x = ∂ ( u M ( x , y ) ) ∂ y \frac{\partial (u 1)}{\partial x} = \frac{\partial(uM(x,y))}{\partial y} ∂ x ∂ ( u 1 ) = ∂ y ∂ ( u M ( x , y )) ∂ u ∂ x = u ∂ ( M ( x , y ) ) ∂ y + M ( x , y ) ∂ ( u ) ∂ y \frac{\partial u}{\partial x} = u\frac{\partial(M(x,y))}{\partial y}+M(x,y)\frac{\partial(u)}{\partial y} ∂ x ∂ u = u ∂ y ∂ ( M ( x , y )) + M ( x , y ) ∂ y ∂ ( u ) ∂ u ∂ x − M ( x , y ) ∂ ( u ) ∂ y = u ∂ ( M ( x , y ) ) ∂ y \frac{\partial u}{\partial x} -M(x,y)\frac{\partial(u)}{\partial y}= u\frac{\partial(M(x,y))}{\partial y} ∂ x ∂ u − M ( x , y ) ∂ y ∂ ( u ) = u ∂ y ∂ ( M ( x , y )) ∂ u ∂ x − M ( x , y ) ∂ ( u ) ∂ y u = ∂ ( M ( x , y ) ) ∂ y \frac{\frac{\partial u}{\partial x} -M(x,y)\frac{\partial(u)}{\partial y}}{u} = \frac{\partial(M(x,y))}{\partial y} u ∂ x ∂ u − M ( x , y ) ∂ y ∂ ( u ) = ∂ y ∂ ( M ( x , y )) ∂ u ∂ x − M ( x , y ) ∂ ( u ) ∂ y u ≡ ∂ ( M ( x , y ) ) ∂ y = ∂ ( P ( x ) y − Q ( x ) ) ∂ y ≡ P ( x ) \frac{\frac{\partial u}{\partial x} -M(x,y)\frac{\partial(u)}{\partial y}}{u} \equiv \frac{\partial(M(x,y))}{\partial y} \ \ \ \ = \ \ \ \ \frac{\partial(P(x)y-Q(x))}{\partial y} \equiv P(x) u ∂ x ∂ u − M ( x , y ) ∂ y ∂ ( u ) ≡ ∂ y ∂ ( M ( x , y )) = ∂ y ∂ ( P ( x ) y − Q ( x )) ≡ P ( x ) Entonces u depende solamente de x:
∂ u ∂ x − M ( x , y ) ∂ ( u ) ∂ y u ≡ ∂ u ∂ x u ≡ ∂ ln ( ∣ u ∣ ) ∂ x = P ( x ) \frac{\frac{\partial u}{\partial x} - \cancel{M(x,y)\frac{\partial(u) }{\partial y}}}{u} \equiv \frac{\frac{\partial u}{\partial x} }{u} \equiv \frac{\partial \ln (|u|)}{\partial x} \ \ \ = \ \ \ P(x) u ∂ x ∂ u − M ( x , y ) ∂ y ∂ ( u ) ≡ u ∂ x ∂ u ≡ ∂ x ∂ ln ( ∣ u ∣ ) = P ( x ) ∫ ∂ ln ( ∣ u ∣ ) ≡ ln ( ∣ u ∣ ) = ∫ P ( x ) ∂ x \int \partial \ln (|u|) \equiv \ln (|u|) \ \ \ = \ \ \ \int P(x) \partial x ∫ ∂ ln ( ∣ u ∣ ) ≡ ln ( ∣ u ∣ ) = ∫ P ( x ) ∂ x u = e ∫ P ( x ) ∂ x u = e^{\int P(x) \partial x} u = e ∫ P ( x ) ∂ x u ( d y d x + P ( x ) y ) = u Q ( x ) u(\frac{dy}{dx}+ P(x)y) = u Q(x) u ( d x d y + P ( x ) y ) = u Q ( x ) e ∫ P ( x ) ∂ x d y d x + e ∫ P ( x ) ∂ x P ( x ) y = e ∫ P ( x ) ∂ x Q ( x ) e^{\int P(x) \partial x} \frac{dy}{dx}+ e^{\int P(x) \partial x} P(x)y = e^{\int P(x) \partial x} Q(x) e ∫ P ( x ) ∂ x d x d y + e ∫ P ( x ) ∂ x P ( x ) y = e ∫ P ( x ) ∂ x Q ( x ) e ∫ P ( x ) ∂ x d y d x + d e ∫ P ( x ) ∂ x d x y = e ∫ P ( x ) ∂ x Q ( x ) e^{\int P(x) \partial x} \frac{dy}{dx}+\frac{d e^{\int P(x) \partial x}}{dx}y = e^{\int P(x) \partial x} Q(x) e ∫ P ( x ) ∂ x d x d y + d x d e ∫ P ( x ) ∂ x y = e ∫ P ( x ) ∂ x Q ( x ) d ( y e ∫ P ( x ) ∂ x ) d x = e ∫ P ( x ) ∂ x Q ( x ) \frac{d (y e^{\int P(x) \partial x})}{dx} = e^{\int P(x) \partial x} Q(x) d x d ( y e ∫ P ( x ) ∂ x ) = e ∫ P ( x ) ∂ x Q ( x ) ∫ d ( y e ∫ P ( x ) ∂ x ) = ∫ e ∫ P ( x ) ∂ x Q ( x ) d x \int d(y e^{\int P(x) \partial x}) =\int e^{\int P(x) \partial x} Q(x) dx ∫ d ( y e ∫ P ( x ) ∂ x ) = ∫ e ∫ P ( x ) ∂ x Q ( x ) d x y e ∫ P ( x ) ∂ x = ∫ e ∫ P ( x ) ∂ x Q ( x ) d x y e^{\int P(x) \partial x} =\int e^{\int P(x) \partial x} Q(x) dx y e ∫ P ( x ) ∂ x = ∫ e ∫ P ( x ) ∂ x Q ( x ) d x y = ∫ e ∫ P ( x ) ∂ x Q ( x ) d x e ∫ P ( x ) ∂ x y =\frac{\int e^{\int P(x) \partial x} Q(x) dx}{e^{\int P(x) \partial x}} y = e ∫ P ( x ) ∂ x ∫ e ∫ P ( x ) ∂ x Q ( x ) d x Bernoulli d y d x + P ( x ) y = Q ( x ) y n \frac{dy}{dx}+ P(x)y = Q(x)y^n d x d y + P ( x ) y = Q ( x ) y n y − n d y d x + P ( x ) y 1 − n = Q ( x ) y^{-n} \frac{dy}{dx}+ P(x)y^{1-n} = Q(x) y − n d x d y + P ( x ) y 1 − n = Q ( x ) v = y 1 − n v=y^{1-n} v = y 1 − n d v d x = ( 1 − n ) y − n d y d x \frac{dv}{dx}=(1-n)y^{-n} \frac{dy}{dx} d x d v = ( 1 − n ) y − n d x d y 1 1 − n d v d x = y − n d y d x \frac{1}{1-n} \frac{dv}{dx}=y^{-n} \frac{dy}{dx} 1 − n 1 d x d v = y − n d x d y 1 1 − n d v d x + P ( x ) v = Q ( x ) \frac{1}{1-n} \frac{dv}{dx}+ P(x)v = Q(x) 1 − n 1 d x d v + P ( x ) v = Q ( x ) Llegado aquí se resuelve por Riccati
EDO de orden superior EDO lineal de orden superior Es una EDO de grado 1, tal que cumple la linealidad
f ( x ) = y = a n ( x ) d n y d x n + a n − 1 ( x ) d n − 1 y d x n − 1 + . . . + a 1 ( x ) d y d x + a 0 ( x ) y f(x)=y=a_n (x) \frac{d^n y}{dx^n}+a_{n-1} (x) \frac{d^{n-1} y}{dx^{n-1}} + ... + a_1 (x) \frac{dy}{dx}+a_0(x)y f ( x ) = y = a n ( x ) d x n d n y + a n − 1 ( x ) d x n − 1 d n − 1 y + ... + a 1 ( x ) d x d y + a 0 ( x ) y g ( x ) = y = sin ( x ) d 3 y d x 3 + 3 d 2 y d x 2 + e x d y d x + x y g(x)=y= \sin(x) \frac{d^3 y}{dx^3} + 3 \frac{d^2 y}{dx^2}+e^x \frac{d y}{dx} + xy g ( x ) = y = sin ( x ) d x 3 d 3 y + 3 d x 2 d 2 y + e x d x d y + x y g ( x ) = y = sin ( x ) y ′ ′ ′ + 3 y ′ ′ + e x y ′ + x y g(x)=y= \sin(x)y''' + 3 y''+e^x y' + xy g ( x ) = y = sin ( x ) y ′′′ + 3 y ′′ + e x y ′ + x y g ( x ) = y = sin ( x ) D 3 [ y ] + 3 y D 2 [ y ] + e x D [ y ] + x y g(x)=y= \sin(x)D^3[y] + 3 yD^2[y]+e^xD[y] + xy g ( x ) = y = sin ( x ) D 3 [ y ] + 3 y D 2 [ y ] + e x D [ y ] + x y Igual que una transformación lineal, por lo que satisface las propiedades de:
Aditividad D [ y 1 + y 2 ] = D [ y 1 ] + D [ y 2 ] D[y_1+y_2]=D[y_1]+D[y_2] D [ y 1 + y 2 ] = D [ y 1 ] + D [ y 2 ] Homogeneidad D [ k y ] = k D [ y ] D[ky]=kD[y] D [ k y ] = k D [ y ] EDO lineal Homogénea Reducción de orden y ′ ′ + p ( x ) y ′ + q ( x ) y = 0 y''+p(x)y'+q(x)y=0 y ′′ + p ( x ) y ′ + q ( x ) y = 0 Solución conocida: y 1 ( x ) y_1(x) y 1 ( x ) Solución a hallar: y 2 ( x ) = v ( x ) y 1 ( x ) y_2(x)=v(x)y_1(x) y 2 ( x ) = v ( x ) y 1 ( x ) y 2 ( x ) = y 1 ( x ) ∫ e ∫ p ( x ) d x y 1 2 d x y_2(x)=y_1(x) \int \frac{e^{\int p(x)dx}}{y_1 ^2} dx y 2 ( x ) = y 1 ( x ) ∫ y 1 2 e ∫ p ( x ) d x d x Coeficientes constantes y n = e r n x \large y_n=e^{r_n x} y n = e r n x Cauchy-Euler y n = x r n = e r n ln ( x ) \large y_n=x^{r_n}=e^{r_n \ln(x)} y n = x r n = e r n l n ( x ) r = r= r = raíz de la ecuación característica e θ i = cos ( θ ) + i sin ( θ ) e^{\theta i}= \cos(\theta)+ i \sin(\theta) e θ i = cos ( θ ) + i sin ( θ ) Para raíces múltiples se deriva con respecto a r r r EDO lineal no Homogénea E D O l = x 3 + sin ( a x ) + e b x + 0 EDOl = x^{3}+\sin(ax)+e^{bx}+0 E D Ol = x 3 + sin ( a x ) + e b x + 0 y 1 = A x 3 + B x 2 + C x + D y_1 = Ax^{3}+Bx^{2}+Cx+D y 1 = A x 3 + B x 2 + C x + D y 2 = A 1 sin ( a x ) + B 1 cos ( b x ) y_2 = A_1 \sin(ax) + B_1 \cos(bx) y 2 = A 1 sin ( a x ) + B 1 cos ( b x ) y 3 = A 2 e a x y_3 = A_2 e^{ax} y 3 = A 2 e a x y h = Soluci o ˊ n homog e ˊ nea y_h = \text{Solución homogénea} y h = Soluci o ˊ n homog e ˊ nea E D O l ( y 1 ) = x 3 EDOl(y_1) = x^{3} E D Ol ( y 1 ) = x 3 E D O l ( y 2 ) = sin ( a x ) EDOl(y_2) = \sin(ax) E D Ol ( y 2 ) = sin ( a x ) E D O l ( y 3 ) = e b x EDOl(y_3) =e^{bx} E D Ol ( y 3 ) = e b x E D O l ( y h ) = 0 EDOl(y_h) = 0 E D Ol ( y h ) = 0 E D O l ( y 1 + y 2 + y 3 + y h ) = x 3 + sin ( a x ) + e b x + 0 EDOl(y_1+y_2+y_3+y_h) = x^{3}+\sin(ax)+e^{bx}+0 E D Ol ( y 1 + y 2 + y 3 + y h ) = x 3 + sin ( a x ) + e b x + 0 Coeficientes indeterminados aaa
Transformada de Laplace La transformada de Laplace de una función f ( t ) \displaystyle f(t) f ( t ) definida para todos los números reales t ≥ 0 t\geq 0 t ≥ 0 , es la función F ( s ) \displaystyle F(s) F ( s ) definida por:
F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t F(s)=\int _{0}^{\infty }e^{-st}f(t)\,dt F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t siempre y cuando la integral esté definida.
Cuando f ( t ) \displaystyle f(t) f ( t ) es una distribución con una singularidad en 0 0 0 entonces la transformada de Laplace se define como:
F ( s ) = lim ε → 0 ∫ − ε ∞ e − s t f ( t ) d t F(s)=\lim _{\varepsilon \rightarrow 0}\int _{-\varepsilon }^{\infty }e^{-st}f(t)\,dt F ( s ) = ε → 0 lim ∫ − ε ∞ e − s t f ( t ) d t f ( t ) f(t) f ( t ) L { f ( t ) } \mathcal{L} \{ f(t) \} L { f ( t )} 1 1 1 1 s \frac{1}{s} s 1 t n , n ∈ N t^n \ , \ n \in \mathbb{N} t n , n ∈ N n ! s n + 1 \frac{n!}{s^{n+1}} s n + 1 n ! e a t e^{at} e a t 1 s − a \frac{1}{s-a} s − a 1 sin ( a t ) \sin(at) sin ( a t ) a s 2 + a 2 \frac{a}{s^2+a^2} s 2 + a 2 a cos ( a t ) \cos(at) cos ( a t ) s s 2 + a 2 \frac{s}{s^2+a^2} s 2 + a 2 s δ ( t ) \delta(t) δ ( t ) 1 1 1 u ( t − a ) u(t-a) u ( t − a ) e − a s s \frac{e^{-as}}{s} s e − a s sinh ( a t ) \sinh(at) sinh ( a t ) a s 2 − a 2 \frac{a}{s^2-a^2} s 2 − a 2 a cosh ( a t ) \cosh(at) cosh ( a t ) s s 2 − a 2 \frac{s}{s^2-a^2} s 2 − a 2 s
L { 1 } = ∫ 0 ∞ e − s t d t \mathcal{L} \{ 1 \}= \int _{0}^{\infty }e^{-st}\,dt L { 1 } = ∫ 0 ∞ e − s t d t L { 1 } = [ e − s t − s ] 0 ∞ = lim t → ∞ e − s t − s − e − s ⋅ 0 − s = 0 + 1 s \mathcal{L} \{ 1 \}= \left[ \frac{e^{-st}}{-s} \right]_0^\infty=\lim_{t \to \infty} \frac{e^{-st}}{-s} - \frac{e^{-s \cdot 0}}{-s}=0+\frac{1}{s} L { 1 } = [ − s e − s t ] 0 ∞ = t → ∞ lim − s e − s t − − s e − s ⋅ 0 = 0 + s 1 L { 1 } = 1 s \mathcal{L} \{ 1 \}= \frac{1}{s} L { 1 } = s 1 L { sin ( a t ) } = ∫ 0 ∞ e − s t sin ( a t ) d t \mathcal{L} \{ \sin(at) \}= \int _{0}^{\infty }e^{-st}\sin(at)\,dt L { sin ( a t )} = ∫ 0 ∞ e − s t sin ( a t ) d t L { 1 } = [ e − s t − s ] 0 ∞ = lim t → ∞ e − s t − s − e − s ⋅ 0 − s = 0 + 1 s \mathcal{L} \{ 1 \}= \left[ \frac{e^{-st}}{-s} \right]_0^\infty=\lim_{t \to \infty} \frac{e^{-st}}{-s} - \frac{e^{-s \cdot 0}}{-s}=0+\frac{1}{s} L { 1 } = [ − s e − s t ] 0 ∞ = t → ∞ lim − s e − s t − − s e − s ⋅ 0 = 0 + s 1 L { 1 } = 1 s \mathcal{L} \{ 1 \}= \frac{1}{s} L { 1 } = s 1 L { sin ( a x ) } \mathcal{L} \{ \sin(ax) \} L { sin ( a x )}