leeern by Yuri

Notas de YuriRod


Página tipo blog en el que voy a publicar mis notas de aprendizaje, en especial de temas como matemáticas, física y quizá algo de programación

Redes neuronales
Redes neuronales
Redes neuronales
Redes neuronales

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))=0F(x,y,y',y'',…,y^{(n)})=0
  • xx es la variable independiente.
  • y=y(x)y=y(x) es la función desconocida o variable dependiente.
  • y,y,...,y(n)y',y'', ..., y^{(n)} son las derivadas de yy respecto a xx hasta el 𝑛-ésimo orden.

EDO de primer orden

M(x,y)dx+N(x,y)dy=0M(x,y)dx+N(x,y)dy=0

M(x,y)M(x,y) 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)dx+N(y)dy=0M(x)dx+N(y)dy=0N(y)dy=M(x)dxN(y)dy = -M(x)dxN(y)dy=M(x)dx\int N(y)dy = -\int M(x)dxn(y)=m(x)n(y)=m(x)y=n1(m(x))y=n^{-1}(m(x))

Convertir una EDO de primer orden a una EDO de variables separables:


Factorización
M(x)dx+N(y)dy=0M(x)dx+N(y)dy=0

M(x)dxN(y)dy=1\frac{M(x)dx}{N(y)dy}=-1dxdy=N(y)M(x)\frac{dx}{dy}=-\frac{N(y)}{M(x)}

dxdy=f(x,y)\frac{dx}{dy}=f(x,y)

Se debe verificar que f(x,y)f(x,y) se pueda factorizar f(x,y)=g(x)h(y)f(x,y)=g(x)h(y)

Si f(x,y)f(x,y) es factorizable entonces de cumple que:

x0,y0/f(x0,y0)0    f(x0,y0)f(x,y)=f(x0,y)f(x,y0)\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)

Resolviendo:

dxdy=f(x,y)\frac{dx}{dy}=f(x,y)dxdy=g(x)h(y)\frac{dx}{dy}=g(x)h(y)g(x)1dx=h(y)dyg(x)^{-1}dx=h(y)dy
Funciones homogéneas
f(λx,λy)=λkf(x,y)f(\lambda x,\lambda y)=\lambda^{k}f(x,y)

f(x,y)f(x,y) es una función homogénea de grado kk

λ=1x\lambda=\frac{1}{x}f(xx,yx)=1xkf(x,y)f(\frac{x}{x}, \frac{y}{x})=\frac{1}{x}^{k}f(x,y)f(1,yx)=xkf(x,y)f(1, \frac{y}{x})=x^{-k}f(x,y)


xkf(1,yx)=f(x,y)x^{k} f(1, \frac{y}{x})=f(x,y)ykf(xy,1)=f(x,y)y^{k} f( \frac{x}{y},1)=f(x,y)
M(x,y)dx+N(x,y)dy=0M(x,y)dx+N(x,y)dy=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
ykM(xy,1)dx+ykN(xy,1)dy=0y^{k} M(\frac{x}{y},1)dx+y^{k} N(\frac{x}{y},1)dy=0M(xy,1)dx+N(xy,1)dy=0M(\frac{x}{y},1)dx+ N(\frac{x}{y},1)dy=0v=xyv=\frac{x}{y}yv=xyv=xydv+vdy=dxydv + vdy= dxM(v,1)(ydv+vdy)+N(v,1)dy=0M(v,1)(ydv + vdy)+ N(v,1)dy=0
N más fácil que M
xkM(1,yx)dx+xkN(1,yx)dy=0x^{k} M(1,\frac{y}{x})dx+x^{k} N(1,\frac{y}{x})dy=0M(1,yx)dx+N(1,yx)dy=0 M(1,\frac{y}{x})dx+ N(1,\frac{y}{x})dy=0v=yxv=\frac{y}{x}yv=yyv=yxdv+vdx=dyxdv+vdx=dyM(1,v)dx+N(1,v)(xdv+vdx)=0 M(1,v)dx+ N(1,v)(xdv+vdx)=0
Ecuaciones con coeficientes lineales
(Ax+By+CLinealL1:)dx+(Dx+Ey+FLinealL2:)dy=0(\overset{L1:}{\underset{Lineal}{\underbrace{Ax+By+C}}})dx + (\overset{L2:}{\underset{Lineal}{\underbrace{Dx+Ey+F}}})dy =0Se 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

dF=Fxdx+Fydy=0dF = \frac{\partial F}{\partial x} dx + \frac{\partial F}{\partial y} dy=0Fxy=Fyx\frac{\partial \frac{\partial F}{\partial x}}{\partial y} = \frac{\partial \frac{\partial F}{\partial y}}{\partial x} dF=M(x,y)dx+N(x,y)dy=0dF = M(x,y) dx + N(x,y) dy=0M(x,y)y=N(x,y)x\frac{\partial M(x,y)}{\partial y} = \frac{\partial N(x,y)}{\partial x} Fx=M(x,y)\frac{\partial F}{\partial x} = M(x,y) F=M(x,y)x\partial F = M(x,y) \partial xF=M(x,y)x\int \partial F =\int M(x,y) \partial xF=M(x,y)xF =\int M(x,y) \partial xFy=M(x,y)xy=N(x,y)\frac{\partial F}{\partial y} = \frac{\partial \int M(x,y) \partial x}{\partial y} = N(x,y)
Factor integrante
M(x,y)dx+N(x,y)dy=0 M(x,y) dx + N(x,y) dy=0M(x,y)yN(x,y)x\frac{\partial M(x,y)}{\partial y} \neq \frac{\partial N(x,y)}{\partial x} u(x,y)M(x,y)dx+u(x,y)N(x,y)dy=0 u(x,y)M(x,y) dx + u(x,y)N(x,y) dy=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}
Cómo hallar u(x,y)
(uM)y=(uN)x\frac{\partial(uM)}{\partial y} = \frac{\partial (uN)}{\partial x} uMy+Muy=uNx+Nuxu\frac{\partial M}{\partial y}+M\frac{\partial u}{\partial y} = u \frac{\partial N}{\partial x} + N \frac{\partial u}{\partial x} u(MyNx)=NuxMuyu(\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} ) = N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} MyNx=NuxMuyu\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} }{u}

Si u(x)u(x):

MyNx=NuxMuyu\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}MyNxN=uxu=lnux\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}MyNxNx=lnu=lnu\int \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{N} \partial x =\int \partial \ln{|u|}= \ln{|u|}eMyNxNx=ue^{\int \frac{\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x}}{N} \partial x }=u

Si u(y)u(y):

MyNx=NuxMuyu\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}MyNxM=uyu=lnuy\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}

Si u(x,y)=xmynu(x,y)= x^m y^n:

MyNx=NuxMuyu\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = \frac{N \frac{\partial u}{\partial x} -M\frac{\partial u}{\partial y} }{u}MyNx=Nynmxm1Mxmnyn1xmyn\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}MyNx=Nmx1Mny1\frac{\partial M}{\partial y} -\frac{\partial N}{\partial x} = N m x^{-1} -M n y^{-1}

Riccati

Es una edo de la forma:

dydx+P(x)y=Q(x)\frac{dy}{dx}+ P(x)y = Q(x)

que se le multiplica por un factor integrante

u=eP(x)xu = e^{\int P(x) \partial x} dydx+P(x)y=Q(x)\frac{dy}{dx}+ P(x)y = Q(x)

Obteniendo el factor integrante para Ricatti

dy+P(x)ydx=Q(x)dxdy+ P(x)ydx = Q(x)dxdy+(P(x)yQ(x))dx=0dy+ (P(x)y-Q(x))dx =0M(x,y)=P(x)yQ(x)M(x,y)= P(x)y-Q(x)M(x,y)dx+dy=0 M(x,y)dx+dy = 0

Hallamos u:

dy+M(x,y)dx=0dy+ M(x,y)dx = 0(u1)x=(uM(x,y))y\frac{\partial (u 1)}{\partial x} = \frac{\partial(uM(x,y))}{\partial y}ux=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}uxM(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}uxM(x,y)(u)yu=(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}uxM(x,y)(u)yu(M(x,y))y    =    (P(x)yQ(x))yP(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)

Entonces u depende solamente de x:

uxM(x,y)(u)yuuxuln(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) ln(u)ln(u)   =   P(x)x\int \partial \ln (|u|) \equiv \ln (|u|) \ \ \ = \ \ \ \int P(x) \partial x u=eP(x)xu = e^{\int P(x) \partial x}
u(dydx+P(x)y)=uQ(x)u(\frac{dy}{dx}+ P(x)y) = u Q(x)eP(x)xdydx+eP(x)xP(x)y=eP(x)xQ(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)eP(x)xdydx+deP(x)xdxy=eP(x)xQ(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)d(yeP(x)x)dx=eP(x)xQ(x) \frac{d (y e^{\int P(x) \partial x})}{dx} = e^{\int P(x) \partial x} Q(x)d(yeP(x)x)=eP(x)xQ(x)dx\int d(y e^{\int P(x) \partial x}) =\int e^{\int P(x) \partial x} Q(x) dxyeP(x)x=eP(x)xQ(x)dxy e^{\int P(x) \partial x} =\int e^{\int P(x) \partial x} Q(x) dxy=eP(x)xQ(x)dxeP(x)xy =\frac{\int e^{\int P(x) \partial x} Q(x) dx}{e^{\int P(x) \partial x}}
Bernoulli
dydx+P(x)y=Q(x)yn\frac{dy}{dx}+ P(x)y = Q(x)y^nyndydx+P(x)y1n=Q(x)y^{-n} \frac{dy}{dx}+ P(x)y^{1-n} = Q(x)v=y1nv=y^{1-n}dvdx=(1n)yndydx\frac{dv}{dx}=(1-n)y^{-n} \frac{dy}{dx}11ndvdx=yndydx\frac{1}{1-n} \frac{dv}{dx}=y^{-n} \frac{dy}{dx}11ndvdx+P(x)v=Q(x)\frac{1}{1-n} \frac{dv}{dx}+ 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=an(x)dnydxn+an1(x)dn1ydxn1+...+a1(x)dydx+a0(x)yf(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)yg(x)=y=sin(x)d3ydx3+3d2ydx2+exdydx+xyg(x)=y= \sin(x) \frac{d^3 y}{dx^3} + 3 \frac{d^2 y}{dx^2}+e^x \frac{d y}{dx} + xyg(x)=y=sin(x)y+3y+exy+xyg(x)=y= \sin(x)y''' + 3 y''+e^x y' + xyg(x)=y=sin(x)D3[y]+3yD2[y]+exD[y]+xyg(x)=y= \sin(x)D^3[y] + 3 yD^2[y]+e^xD[y] + xy

Igual que una transformación lineal, por lo que satisface las propiedades de:

  • Aditividad
  • D[y1+y2]=D[y1]+D[y2]D[y_1+y_2]=D[y_1]+D[y_2]
  • Homogeneidad
  • D[ky]=kD[y]D[ky]=kD[y]

EDO lineal Homogénea

Reducción de orden
y+p(x)y+q(x)y=0y''+p(x)y'+q(x)y=0
  • Solución conocida: y1(x)y_1(x)
  • Solución a hallar: y2(x)=v(x)y1(x)y_2(x)=v(x)y_1(x)
y2(x)=y1(x)ep(x)dxy12dxy_2(x)=y_1(x) \int \frac{e^{\int p(x)dx}}{y_1 ^2} dx
Coeficientes constantes
yn=ernx\large y_n=e^{r_n x}
Cauchy-Euler
yn=xrn=ernln(x)\large y_n=x^{r_n}=e^{r_n \ln(x)}
  • r=r= raíz de la ecuación característica
  • eθi=cos(θ)+isin(θ)e^{\theta i}= \cos(\theta)+ i \sin(\theta)
  • Para raíces múltiples se deriva con respecto a rr

EDO lineal no Homogénea

EDOl=x3+sin(ax)+ebx+0EDOl = x^{3}+\sin(ax)+e^{bx}+0 y1=Ax3+Bx2+Cx+Dy_1 = Ax^{3}+Bx^{2}+Cx+D y2=A1sin(ax)+B1cos(bx)y_2 = A_1 \sin(ax) + B_1 \cos(bx) y3=A2eaxy_3 = A_2 e^{ax} yh=Solucioˊn homogeˊneay_h = \text{Solución homogénea} EDOl(y1)=x3EDOl(y_1) = x^{3} EDOl(y2)=sin(ax)EDOl(y_2) = \sin(ax) EDOl(y3)=ebxEDOl(y_3) =e^{bx} EDOl(yh)=0EDOl(y_h) = 0 EDOl(y1+y2+y3+yh)=x3+sin(ax)+ebx+0EDOl(y_1+y_2+y_3+y_h) = x^{3}+\sin(ax)+e^{bx}+0
Coeficientes indeterminados
aaa

Transformada de Laplace

La transformada de Laplace de una función f(t)\displaystyle f(t) definida para todos los números reales t0t\geq 0, es la función F(s)\displaystyle F(s) definida por:

F(s)=0estf(t)dt F(s)=\int _{0}^{\infty }e^{-st}f(t)\,dt

siempre y cuando la integral esté definida.

Cuando f(t)\displaystyle f(t) es una distribución con una singularidad en 00 entonces la transformada de Laplace se define como:

F(s)=limε0εestf(t)dt F(s)=\lim _{\varepsilon \rightarrow 0}\int _{-\varepsilon }^{\infty }e^{-st}f(t)\,dt
f(t)f(t)L{f(t)}\mathcal{L} \{ f(t) \}
111s\frac{1}{s}
tn , nNt^n \ , \ n \in \mathbb{N}n!sn+1\frac{n!}{s^{n+1}}
eate^{at}1sa\frac{1}{s-a}
sin(at)\sin(at)as2+a2\frac{a}{s^2+a^2}
cos(at)\cos(at)ss2+a2\frac{s}{s^2+a^2}
δ(t)\delta(t)11
u(ta)u(t-a)eass\frac{e^{-as}}{s}
sinh(at)\sinh(at)as2a2\frac{a}{s^2-a^2}
cosh(at)\cosh(at)ss2a2\frac{s}{s^2-a^2}
L{1}=0estdt \mathcal{L} \{ 1 \}= \int _{0}^{\infty }e^{-st}\,dt L{1}=[ests]0=limtestses0s=0+1s \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}=1s \mathcal{L} \{ 1 \}= \frac{1}{s}
L{sin(at)}=0estsin(at)dt \mathcal{L} \{ \sin(at) \}= \int _{0}^{\infty }e^{-st}\sin(at)\,dt L{1}=[ests]0=limtestses0s=0+1s \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}=1s \mathcal{L} \{ 1 \}= \frac{1}{s}
L{sin(ax)} \mathcal{L} \{ \sin(ax) \}