#10EDO

Consultas Pre-Primer Parcial: Modelado Continuo

11 min de lectura

Ejercicio 1

{h′=ρ⋅h⋅(1−hk−p)p′=−1ρ⋅p⋅(1−h)\begin{cases} h'=\rho\cdot h\cdot\left( 1- \frac{h}{k}-p \right) \\ p'= - \frac{1}{\rho}\cdot p\cdot(1-h) \end{cases} {H˙=b⋅H⋅(1−Hk)−s⋅P⋅Hp˙=−dP+e⋅s⋅PH\begin{cases} \dot{H}=b\cdot H\cdot\left( 1-\frac{H}{k} \right)-s\cdot P\cdot H \\ \dot{p}=-dP+e\cdot s\cdot PH \end{cases}

Ejercicio 2

a)a) V:Rn→RV:\mathbb{R}^{n}\to \mathbb{R} C1x∗C^{1}\quad x^{*} mínimo de VV X˙=−∇V  ⟹  x∗\dot{X}=-\nabla V\implies x^{*} es equilibrio del sistema asintóticamente estable.

b)b) Estudiar la estabilidad del (0,0)(0,0) para

{x˙=(X−y)2+sin⁡(x)⋅cos⁡(x)y˙=−(x−y)3\begin{cases} \dot{x}=(X-y)^{2} +\sin(x)\cdot \cos(x) \\ \dot{y}=-(x-y)^{3} \end{cases}

Ejercicio 3

{x′=x−y−2x(x2+y2)y′=x+y−y(x2+y2)\begin{cases} x'=x-y-2x(x^{2}+y^{2}) \\ y'=x+y-y(x^{2}+y^{2}) \end{cases}

Probar que el sistema posee un ciclo límite.

Ejercicio 4

x′=−x(4x2+λ2−1)x'=-x(4x^{2}+\lambda^{2}-1)

Esbozar un diagrama de bifurcación.


h=Hesd;p=Psb;k=Kesd;t=bd⋅t;p=bdh=\frac{Hes}{d};\quad p=\frac{Ps}{b};\quad k=\frac{Kes}{d};\quad t=\sqrt{ bd }\cdot t;\quad p=\sqrt{ \frac{b}{d} }

Estabilidad lineal:

DF(p,h)=[ρ(1−2hk−p)−p⋅h1ρ⋅p−1ρ⋅(1−h)]DF(p,h)=\begin{bmatrix} \rho\left( 1- \frac{2h}{k}-p \right) & -p\cdot h \\ \frac{1}{\rho}\cdot p & \frac{-1}{\rho}\cdot(1-h) \end{bmatrix} DF(p1)=DF(0,0)=[ρ00−1ρ]  ⟹  {λ1=ρ>0λ2=−1ρ<0  ⟹  p1 es sillaDF(p_{1})=DF(0,0)=\begin{bmatrix} \rho & 0 \\ 0 & -\frac{1}{\rho} \end{bmatrix} \implies \begin{cases} \lambda_{1}=\rho>0 \\ \lambda_{2}=-\frac{1}{\rho}<0 \end{cases} \implies p_{1} \text{ es silla}

00:00

DF(p2)=DF(k,0)=[−ρ−rhp⋅k0−(1−k)ρ]  ⟹  {λ1=−ρ>0λ2=−(1−k)ρ{<0k<1>0k>1DF(p_{2})=DF(k,0)=\begin{bmatrix} -\rho & -rhp\cdot k \\ 0 & -\frac{(1-k)}{\rho} \end{bmatrix} \implies \begin{cases} \lambda_{1}=-\rho>0 \\ \lambda_{2}=-\frac{(1-k)}{\rho}\begin{cases} <0\quad k<1 \\ >0\quad k>1 \end{cases} \end{cases}

Si k<1  ⟹  p2k<1\implies p_{2} es nodo estable Si k>1  ⟹  p2k>1\implies p_{2} es silla.

k>1k>1

DF(p3)=DF(1,1−1k)=[ρ⋅(1−2k+1+1k)−ρ1ρ(1−1k)0]=[−ρk−ρ1ρ(1−1k)0]DF(p_{3})=DF\left( 1,1-\frac{1}{k} \right)=\begin{bmatrix} \rho\cdot\left( 1-\frac{2}{k}+1+\frac{1}{k} \right) & -\rho \\ \frac{1}{\rho}\left( 1-\frac{1}{k} \right) & 0 \end{bmatrix} =\begin{bmatrix} -\frac{\rho}{k} & -\rho \\ \frac{1}{\rho}\left( 1-\frac{1}{k} \right) & 0 \end{bmatrix} det⁡(DF(p3))=Δ(p3)=1−1k\det(DF(p_{3}))=\Delta(p_{3})=1- \frac{1}{k}

13:39 (Clave esto)

Si la traza me da positiva es inestable. (revisar esto)

traza(DF(p3))=τ(p3)=−ρk<0traza(DF(p_{3}))=\tau(p_{3})=- \frac{\rho}{k}<0

Para matrices 2x2 el polinomio característico es:

pA(λ)=λ2−τ⋅λ+Δp_{A}(\lambda)=\lambda^{2}-\tau\cdot\lambda+\Delta

Veo el discriminante:

λ=τ±τ2−4⋅Δ2\lambda=\frac{\tau\pm \sqrt{ \tau^{2}-4\cdot \Delta }}{2} τ2−4Δ=ρ2k2−4⋅(1−1k)\tau^{2}-4\Delta=\frac{\rho^{2}}{k^{2}}-4\cdot\left( 1-\frac{1}{k} \right) ρ2k2−4⋅(1−1k)>0ρ2>4⋅k−1k⋅k2ρ2≥4(k−1)⋅k\begin{array}{c} \frac{\rho^{2}}{k^{2}}-4\cdot\left( 1-\frac{1}{k} \right) >0\\ \rho^{2}>4\cdot \frac{k-1}{\cancel{ k }}\cdot k^{\cancel{ 2 }} \\ \rho^{2}\geq 4(k-1)\cdot k \end{array}

p3p_{3} es nodo estable.

Si ρ2<4(k−1)k  ⟹  p3\rho^{2}<4(k-1)k\implies p_{3} es foco estable.

La bifurcación es siempre con nodo a nodo.


Ejercicio 2 discutido en 56:00 Ejercicio 4 discutido en 59:14

f(X)=−x⋅(4⋅x2+λ2−1)∼x⋅(1−λ2)f(X)=-x\cdot(4\cdot x^{2}+\lambda^{2}-1)\sim x\cdot(1-\lambda^{2})

Adjuntar foto. Ejercicio 3: 01:04:00

{x=r⋅cos⁡θy=r⋅sin⁡θr2=x2+y2  ⟹  2⋅r⋅r˙=2xx˙+2yy˙\begin{cases} x=r\cdot \cos \theta \\ y=r\cdot \sin \theta \\ r^{2}=x^{2}+y^{2} \end{cases} \implies{2}\cdot r \cdot \dot{r}=2x \dot{x}+2 y \dot{y} =x(x−y−2xr2)+y(x+y−y⋅r2)=x2−xy−2x2r2+xy+y2−x2⋅r2=r2−2x2r2−y2r2=r2⋅(1−x2−x2−y2)=r2⋅(1−r2⋅cos⁡2θ−r2)=r2⋅(1−r2(cos⁡2θ+1))r˙=r⋅(1−r(1−r2(cos⁡2θ+1)))\begin{array}{c} =x(x-y-2xr^{2})+y(x+y-y\cdot r^{2})=x^{2}-xy-2x^{2}r^{2}+xy+y^{2}-x^{2}\cdot r^{2} \\ =r^{2}-2x^{2}r^{2}-y^{2}r^{2} \\ =r^{2}\cdot(1-x^{2}-x^{2}-y^{2}) \\ =r^{2}\cdot(1-r^{2}\cdot \cos ^{2}\theta-r^{2})=r^{2}\cdot(1-r^{2}(\cos ^{2}\theta+1)) \\ \dot{r}=r\cdot(1-r^{}(1-r^{2}(\cos ^{2}\theta+1))) \end{array}

Ahora tengo que darme cuenta que hay un anillo invariante.

r2<r2⋅(cos⁡2θ+1)<2⋅r2r^{2}<r^{2}\cdot(\cos ^{2}\theta+1)<2\cdot r^{2} 1−r2>1−r2(cos⁡2θ+1)>1−2r2>01-r^{2}>1-r^{2}(\cos ^{2}\theta+1)>1-2r^{2}>0

Tomo un radio tal que r<12r<\sqrt{ \frac{1}{2} } . Este es mi radio interior.

Ahora busco otro radio con direcciones para adentro y que no encierra equilibrios. Digamos r>1r>1

Ahora resulta que el anillo resultante es invariante.

  ⟹  A={(x,y)∈R2∣122<x2+y2<2}\implies A=\left\{ (x,y) \in \mathbb{R}^{2} \bigm| \frac{1}{2\sqrt{ 2 }}<\sqrt{ x^{2}+y^{2} }<2 \right\}

Chequear que (0,0)(0,0) es el único equilibrio.

Temas relacionados