Usando Polinômios de Darboux na Obtenção de Integrais
Primeiras Liouvillianas de Equações Diferenciais Ordinárias
Racionais de Segunda Ordem
André Claudino
Como este trabalho está organizado?
Abordagem Darbouxiana
Dado uma 1DO racional
\dfrac{dy}{dx} = \dfrac{M(x,y)}{N(x,y)}
Existe um fator integrante...
R = \prod_i p_i(x,y)^{n_i}
n_i \in \mathbb{Q}
Que satisfaz...
\dfrac{\partial (RM)}{\partial y} = - \dfrac{\partial (RN)}{\partial x}
com...
\mathbb{D} \equiv N(x,y) \partial_x + M(x,y) \partial_y
\dfrac{\mathbb{D}[R]}{R} = -(N_x + M_y)
teremos...
\mathbb{D} [p_i] = p_i q_i
o problema se torna...
R = \prod_i p_i(x,y)^{n_i}
onde
para uma solução liouviliana...
\mathbb{D}[B_i] = B_i \lambda_i
portanto...
R = e^{\frac{A}{B}} \prod_i p_i^ {c_i}
\dfrac{\mathbb{D} [B]}{B} = \sum_i m_i \dfrac{\mathbb{D}[B_i]}{B_i}
\mathbb{D}[A]-A \sum_i m_i \lambda_i = -\prod_i B_i ^{m_i} \left( {\sum_i c_i\tau_i + M_y + N_x} \right)
\mathbb{D}[p_i] = p_i \tau_i
e o problema se torna...
\mathbb{D}[B_i] = \lambda_i B_i
\mathbb{D} [p_i] = p_i q_i
Observações
-
não sabemos os graus máximos dos polinômios de Darboux
-
O número de indeterminados cresce rápido com o grau dos polinômios de Darboux
-
resolver o sistema de equações para cada i é complicado, mesmo para graus baixos
Método da expansão
-
Gera um sistema de equações em formato de cadeira regular
- A complexidade do sistema cresce mais devagar
- Permite solucionar para graus maiores de polinômios de Darboux
- O algoritmo é mais comportado do ponto de vista de consumo de recursos
Ideia
\dfrac{dy}{dx} = \dfrac{M}{N}
M = \sum_{i} M(y)_{i} x^i \\
N = \sum_{i} N(y)_{i} x^i
P = \sum P(y)_{i} x^i \\
Q = \sum Q(y)_{i} x^i
M_0 \dfrac{d P_0}{dy} + N_0 P_1 = Q_0 P_0
M_0 \dfrac{d P_1}{dy} + M_1 \dfrac{d P_0}{dy} + 2 N_0 P_2 + N_1 P_1 = Q_0 P_1 + Q_1 P_0
\cdots
Exemplo
\dfrac{dy}{dx} = -\dfrac{3 + y^2 x^4}{x^4}
\mathbb{D} = x^4\partial_x - (3+y^2 x^4)\partial_y
P = P_0(y) + P_1(y) x + P_2(y) x^2 + P_3(y) x^3 + P_4(y) x^4
Q = Q_0(y) + Q_1(y) x + Q_2(y) x^2 + Q_3(y) x^3 + Q_4(y) x^4
g(P) \equiv 6
g_x(Q) \equiv max(g_x(M), g_x(N)-1) = max(4, 4 -1) = max(4, 3) = 4
g_x(p) \equiv 4
g(Q) \equiv max(g(M), g(N)) = 5
-3 \frac{d P_0}{dy} - P_0 Q_0 = 0 \\
-P_0 Q_1 - P_1 Q_0 -3 \frac{dP_1}{dy} = 0 \\
-3 \frac{dP_4}{dy} - y^2\frac{dP_0}{dy} + P_1 - P_0 Q_4 - P_1 Q_3 - P_2 Q_2 - P_3 Q_1 - P_4 Q_0 = 0 \\
-P_0 Q_3 - P_1 Q_2 - P_2 Q_1 -P_3Q_0 - 3 \frac{dP_3}{dy} = 0 \\
-y^2 \frac{dP_4}{dy} - P_4 Q_4 = 0 \\
2 P_2 - P_1Q_4 -P_2Q_3 -P_3 Q_2 - P_4 Q_1 - y^2 \frac{dP_1}{dy} = 0 \\
3 P_3 -P_2 Q_4 -P_3Q_3 - P_4Q_2 - y^2\frac{dP_2}{dy} = 0 \\
4 P_4 - P_3Q_4 - P_4Q_3 -y^2\frac{dP_3}{dy} = 0 \\
-P_0 Q_2 - P_1 Q_1 - P_2 Q_0 - 3 \frac{dP_2}{dy} = 0
-y^2 \frac{dP_4}{dy} - P_4 Q_4 = 0
Resovendo
P_4 = C e^{\int{-\frac{Q_4}{y^2}}dy}
temos
\int{-\frac{Q_4}{y^2}}dy
deve ser logaritmo ou zero
max(g(P_4)) = max(g_y(P_4)) = 2
g(Q_4) = g_y (Q_4) = 1 \Rightarrow Q_4 = k y
P_4 = C e^{\int{-\frac{k y}{y^2}}dy} = C e^{- k \int{\frac{y}{y^2}}dy} = C e^{- k \int{\frac{dy}{y}}} = C e^{-kln(y)} = C y^{-k}
C \equiv 1 \\
k \in {0, -1, -2}
P_4 = y^{2}
Q_4 = -2y
Resolvendo
Substituindo
4 P_4 - P_3Q_4 - P_4Q_3 -y^2\frac{dP_3}{dy} = 0
4 y^2 + 2P_3y - y^2Q_3 - y^2\frac{dP_3}{dy} = 0
\Rightarrow P_3 = \left( \int{\frac{4-Q_3}{y^2}}dy + C_1 \right)y^2
Q_3 \equiv a_0 + a_1y + a_2y^2
usando
P_3 = -a_1 y^2 ln(y) -a_2 y^3 + y^2 K_1 + (a_0 - 4)y
P_3 = -a_2 y^3 + y^2 K_1 + (a_0 - 4)y
P_3
é polinômio
\Rightarrow a_1 \equiv 0
Q_3 = a_0 + a_2 y^2
3P_3 - P_2 Q_4 -P_3 Q_3 -P_4 Q_2 - y^2 \frac{d P_2}{dy} = 0
Resovendo
Temos
P_2 = y^2 \int {\dfrac{a_2^2 y^4 -a_2 K_1 y^3 + a_2 y^2 + (3K_1 -Q_2 -a_0 k_1)y -12 + 7a_0 -a_0^2}{y^3}} dy + K_2 y^2
Q_2 \equiv b_0 + b_1 y + b_2 y^2 +b_3 y^3 \Rightarrow
P_2 = (a_2 - b_1)y^2 ln(y) + \left( \frac{a_2^2 -b^3}{2} \right) y^4 - (b^2 + a_2 k_1)y^3 + y^2 K_2 + (b_0 + a_0 K_1 -3K_1)y
+ \frac{a_0^2-7a_0}{2} + 6
b_1 = a_2
(polinômio)
P_2 = \frac{a_3^2-b^3}{2} y^4 - (b_2 + a_2 K_1)y^3 + K_2 y^2 + (K_1(a_0 - 3) + b_0 )y + \frac{a_0^2 - 7a_0}{2} + 6 \\
Q_2 = b_0 + a_2 y + b_2 y^2 + b_3 y^3
\Rightarrow
2 P_2 - P_1 Q_4 - P_2 Q_3 - P_3 Q_2 - P_4 Q_1 - y^2\frac{dP_1}{dy} = 0
Resovendo
Temos
P_1 = y^2 \int{\dfrac{a_2(a_2 -3 b_3)y^6 -2(a_2(K_1+2b_2) -K_1b_3)y^5 \\
+ (a_2^2(a_0-4) + 2 a_2K_2 + b_3(a_0 -6) + 2 b_2 K_1)y^4 \\
- 4b_2 y^3}{-2y^4}
}dy + K_3y^2
Q_1 \equiv c_0 + c_1 y + c_2 y^2 +c_3 y^3 + c_4 y^4 \Rightarrow
P_1 = - \frac{a_2^2 - 3 a_2 b_3 + 2c_4}{6}y^5
+ \frac{K_1 a_2^2 + 2a_2 b_2 - K_1 b_3 -c_3}{2}y^4
+ y^2 \int{\dfrac{
(a_2(a_0^2 - 5 a_0 + 4) + 2 b_0K_1 + 2K_2(a_0 -2) + 2 Q_1)y^2 \\
+ 2(2 K_1 a_0^2 + a_0(2 b_0 - 5K_1) + 6K_1 6 b_0) y \\
a_0^3 - 9 a_0^2 + 26 a_0 -24)}{-2y^4}
}dy + K_3y^2
\begin{cases}
2b_2 - c_1 = 0 \Rightarrow c_1 = 2b_2 \\
a_0^3 - 9 a_0^2 + 26 a_0 - 24 = 0 \Rightarrow a_0 \in \{ 2, 3, 4 \}
\end{cases}
+ (a_2^2 - a_2 K_2 + 2b_3 -2b_2 K_1 -c_2)y^3
+ K_2 y^2 + (b_0 K_1 -a_2 + c_0)y -b_0
Q_1 = c_0 + 2 b_2 y + c_2 y^2 + c_3 y^3 + c_4y^4
-3 \dfrac{d P_0}{dy} - P_0 Q_0 = 0
Resovendo
Temos
P_0 = K_4 e^{-\int\frac{Q_0}{3}dy}
{-\int\frac{Q_0}{3}dy}
Deve ser logaritmo ou zero
-\int\frac{Q_0}{3}dy = 0 \Rightarrow
\begin{cases}
Q_0 = 0 \\
P_0 = K_4
\end{cases}
continuando, teremos...
-6 {K_3}-2 {K_4} {b_2}=0\\-3 {b_3}+3 {
c_2}=0\\-1/2 {a_2} {b_3} {c_4}+1/3 {{c_4}}^{2}+1/6 {{
a_2}}^{3}{c_4}=0\\
-6-{{b_0}}^{2}+{c_0}+2 {K_4}=0\\\,\,\,
6+{{b_0}}^{2}-2 {K_4}-{c_0}=0\\3 {K_1}+{b_0} {c_0}-3 {
b_0}-{K_4} {b_0}=0\\
-{K_4} {c_0}+3 {a_2}-3 {
c_0}-3 {b_0} {K_1}=0\\-{K_4} {c_4}+5/2 {{a_2}}^
{3}-15/2 {a_2} {b_3}+5 {c_4}=0\\
1/6 {{a_2}}^{3}{b_3}-1
/2 {{a_2}}^{2}{c_4}-1/2 {a_2} {{b_3}}^{2}+5/6 {b_3}
{c_4}=0\\4/3 {c_4} {a_2}+1/6 {{a_2}}^{4}+1/2 {{b_3}}
^{2}-{{a_2}}^{2}{b_3}=0\\
6 {c_3}+6 {b_3} {K_1}-6 {
K_1} {{a_2}}^{2}-12 {a_2} {b_2}-{K_4} {c_3}=0\\
-{K_1} {c_0}-{K_3}-{b_0} {K_2}+3 {b_2}+{
a_2} {K_1}=0\\
-2 {b_2}+2 {a_2} {b_0}-6 {K_1}+{
K_1} {c_0}-{{b_0}}^{2}{K_1}-2 {b_0} {c_0}=0\\ -{
K_1} {c_4}+3/2 {b_3} {K_1} {a_2}-3/2 {b_2}
{{a_2}}^{2}-1/2 {K_1} {{a_2}}^{3}+3/2 {c_3} {a_2}+3/2 {b_2} {b_3}=0\\ -{K_4} {a_2}-6 {K_2}+2 {
b_0} {b_2}+{a_2} {c_0}-{{c_0}}^{2}-{b_0} {K_1} {c_0}=0\\ -{K_4} {c_2}+9 {c_2}-18 {b_3}+9 {
a_2} {K_2}-9 {{a_2}}^{2}+9 {b_2} {K_1}=0\\ 1/2 {
b_3} {K_1} {c_4}-1/2 {K_1} {{a_2}}^{2}{c_4}+
1/6 {{a_2}}^{3}{c_3}-{a_2} {b_2} {c_4}-1/2 {a_2}
{b_3} {c_3}+5/6 {c_4} {c_3}=0\\ 9 {b_2}-{K_4} {
b_2}-{K_3} {c_0}+9 {a_2} {K_1}-2 {b_0} {
K_1} {b_2}+{b_0} {c_2}+2 {a_2} {b_2}-2 {c_0} {b_2}=0\\ 1/6 {{a_2}}^{3}{b_2}-1/2 {K_1} {{a_2}}^
{2}{b_3}+4/3 {c_4} {b_2}+{a_2} {K_1} {c_4}-1/2
{{a_2}}^{2}{c_3}-3/2 {a_2} {b_3} {b_2}+{b_3} {
c_3}+1/2 {{b_3}}^{2}{K_1}=0\\ -{K_1} {c_3}+2 {b_2
} {K_1} {a_2}+{K_2} {{a_2}}^{2}+2 {c_2} {a_2}-2 {a_2} {b_3}+7/3 {c_4}-{b_3} {K_2}-4/3 {{
a_2}}^{3}+{{b_2}}^{2}=0\\ -{K_3} {a_2}+2 {a_2} {b_2
}+1/2 {K_1} {{a_2}}^{2}-1/2 {b_0} {b_3}-{b_2} {
K_2}-{K_1} {c_2}+5/2 {c_3}-1/2 {{a_2}}^{2}{b_0}+3/2 {b_3} {K_1}=0\\ {a_2} {K_2} {b_0}+{a_2}
{b_2}-{b_0} {b_3}-2 {b_2} {K_2}-{c_3}-{{a_2}}^{2}{b_0}-{K_3} {a_2}+{K_1} {c_2}+{a_2} {
K_1} {c_0}=0\\ 9 {a_2}-{K_2} {c_0}-{c_2}-{b_0}
{b_2}-{K_4} {a_2}-{a_2} {c_0}-{b_0} {K_1
} {a_2}-{K_3} {b_0}+{{a_2}}^{2}+2 {b_2} {K_1
}=0\\ -6 {{a_2}}^{2}-{{a_2}}^{2}{c_0}-2 {b_3} {c_0}+{
b_2} {K_1} {c_0}+{b_0} {c_3}+{a_2} {K_2}
{c_0}+6 {b_3}+{c_2} {a_2}-2 {K_3} {b_2}-{
b_0} {K_1} {c_2}-{K_4} {b_3}=0\\ -{{a_2}}^{2}{
c_4}+{b_2} {K_1} {c_4}-{a_2} {b_2} {c_3}+{
a_2} {K_2} {c_4}-1/2 {a_2} {b_3} {c_2}+1/2 {
{c_3}}^{2}-1/2 {K_1} {{a_2}}^{2}{c_3}+ +1/2 {b_3} {
K_1} {c_3}+1/6 {{a_2}}^{3}{c_2}+4/3 {c_4} {c_2
}-2 {b_3} {c_4}=0\\ -{{a_2}}^{2}{c_3}-{a_2} {b_3} {
b_2}+2/3 {c_4} {b_2}+{a_2} {K_2} {c_3}-2 {
b_3} {c_3}-{a_2} {b_2} {c_2}+{b_2} {K_1} {
c_3}+ +1/2 {b_3} {K_1} {c_2}+3/2 {c_3} {c_2}-1/
2 {K_1} {{a_2}}^{2}{c_2}-{K_3} {c_4}+1/3 {{a_2}}^{3}{b_2}=0\\ -2 {b_2} {{a_2}}^{2}-{b_0} {K_1} {
c_3}-4 {b_2} {b_3}+2 {{b_2}}^{2}{K_1}-{a_2} {
b_2} {c_0}-1/2 {c_3} {c_0}+{c_3} {a_2}+ +{b_0}
{c_4}+2 {a_2} {K_2} {b_2}-1/2 {K_1} {{a_2
}}^{2}{c_0}+2 {c_2} {b_2}-{K_3} {c_2}+1/2 {b_3
} {K_1} {c_0}=0\\ 3/2 {b_3} {K_1} {b_2}-2 {{b_3}}^{2}-{a_2} {{b_2}}^{2}-1/2 {{a_2}}^{2}{c_2}+{a_2}
{K_1} {c_3}+3/2 {b_3} {c_2}+1/6 {{a_2}}^{4}+ +1/3 {c_4} {a_2}+3/2 {c_3} {b_2}-1/2 {K_1} {{a_2
}}^{2}{b_2}+{a_2} {K_2} {b_3}-{K_2} {c_4}-3/2
{{a_2}}^{2}{b_3}=0\\ -2/3 {b_0} {c_4}-{K_2} {c_3}
+1/6 {{a_2}}^{3}{b_0}-{b_2} {b_3}-1/2 {K_1} {{
a_2}}^{3}+1/2 {c_3} {a_2}-1/2 {a_2} {b_3} {b_0}
+ +{{b_2}}^{2}{K_1}+{a_2} {K_2} {b_2}+1/2 {b_3}
{K_1} {a_2}+2 {c_2} {b_2}+{K_1} {c_4}-{
K_3} {b_3}+{a_2} {K_1} {c_2}-3 {b_2} {{a_2}}^{2}=0\\ -2 {b_3} {c_2}+{b_2} {K_1} {c_2}-{b_0
} {K_1} {c_4}-{K_1} {{a_2}}^{2}{b_2}+{a_2} {
K_2} {c_2}+{c_3} {b_2}-1/2 {a_2} {b_3} {c_0}+ +1/6 {{a_2}}^{3}{c_0}+{b_3} {K_1} {b_2}-2/3 {
c_4} {c_0}-{{a_2}}^{2}{c_2}+{{c_2}}^{2}-{K_3} {
c_3}+{c_4} {a_2}-2 {a_2} {{b_2}}^{2}=0\\ {c_2} {a_2}-{a_2} {b_3}-1/2 {{a_2}}^{2}{c_0}-{c_4}-{{a_2}}
^{3}+2 {{b_2}}^{2}-{K_2} {c_2}-{a_2} {b_2} {b_0}+{K_1} {c_3}+ -1/2 {b_0} {c_3}-{K_3} {b_2}-
1/2 {K_1} {{a_2}}^{2}{b_0}+3 {b_2} {K_1} {a_2}-1/2 {b_3} {c_0}-1/2 {b_3} {K_1} {b_0}+{K_2} {{a_2}}^{2}=0
\{K_1, K_2, K_3, K_4, a_2, b_0, \\b_2, b_3, c_0, c_2, c_3, c_4\}
K_4=3 \\
K_1=K_2=K_3=0\\
a_2=0\\
b_0=b_2=b_3=0\\
c_0=c_2=c_3=c_4=0
P = 3 + x^2 - 2 x^3 y + x^4 y^2 \\
Q = 2 x^3 - 2 y x^4
\Rightarrow
Método da expansão
Simultâneo
\dfrac{dy}{dx} = \dfrac{M}{N}
M = \sum_{i} M(y)_{i} x^i \\
N = \sum_{i} N(y)_{i} x^i
p = \sum p(y)_{i} x^i \\
q = \sum q(y)_{i} x^i
M = \sum_{i} M(x)_{j} y^j \\
N = \sum_{i} N(x)_{j} y^j
p = \sum p(x)_{j} y^j \\
q = \sum q(x)_{j} y^j
Expande em
x
Expande em
y
Procedimento
- Expande em ambas as variáveis
- Resolve o sistema de menor complexidade (sistema A)
- Resolva o sistema restante (sistema B) usando os resultados encontrados no sistema A
- Caso haja indeterminados, resolva no sistema remanescente
Exemplo
\dfrac{dy}{dx} = -\dfrac{3 + y^2 x^4}{x^4}
\mathbb{D} = x^4\partial_x - (3+y^2 x^4)\partial_y
P = p_0(x) + p_1(x) y + p_2(x) y^2
Q = q_0(x) + q_1(x) y
g(P) \equiv 6
g_y(Q) \equiv max(g_y(M) - 1, g_y(N)) = max(2-1, 0) = max(1, 0) = 1
g_y(p) \equiv 2
g(Q) \equiv max(g(M), g(N)) = 5
-3 p_1(x) -p_0(x) q_0(x) +x^4\frac{d}{dx}p_0(x) = 0\\
-6 p_2(x) + x^4\frac{d}{dx}p_1(x) - p_0(x) q_1(x) - p_1(x) q_0(x) = 0 \\
-x^4 p_1(x) + x^4 \frac{d}{dx}p_2(x) - p_1(x) q_1(x) -p_2(x) q_0(x) = 0 \\
-2 x^4 p_2(x) -p_2(x) q_1(x) = 0
Resulta no sistema B
-2 x^4 p_2(x) -p_2(x) q_1(x) = 0
Resolvendo
q_1(x) = -2 x^4
temos...
Com a informação do sistema A
q_1(x) = -2 x^4
P_4 = y^{-k} \\
Q_4 = k y^1
E do sistema B
k = -2 =>
= y^2
Q = q_0(x) -2 x^4 y
\begin{cases}
P_4 = y^2 \\
p_2(x) \equiv a_0 + a_1 x + a_2 x^2 + a_3 x^3 + a_4 x^4
\end{cases}
\Rightarrow a_4 = 1
p_2(x) = a_0 + a_1 x + a_2 x^2 + a_3 x^3 + x^4
P_3 = \left( \int{\frac{4-Q_3}{y^2}}dy + C_1 \right)y^2
g(Q_3) = 1 \Rightarrow Q_3 \equiv A_0 + A_1y
P_3 = (K_1-A_1 ln(y))y^2 + (A_0 - 4)y
P_3 = K_1 y^2 + (A_0 -4)y
P_3
é polinômio
\Rightarrow A_1 \equiv 0
Q_3 = A_0
P_2 = (B_1 ln(y) + K_2) y^2 + (B_0 + A_0 K_1 -3K_1) + \frac{A_0^2 - 7A_0}{2} + 6
B_1 = 0
P_2 = K_2 y^2 + (B_0 + A_0 K_1 -3K_1) + \frac{A_0^2 - 7A_0}{2} + 6
Q_2 = B_0
Q_2 = B_0 + B_1 y \Rightarrow
P_1 = ( -C_1\,\ln(y)+K_3) {y}^{2}+ ( -2\,K_2+A_0\,K_2+B_0\,K_1+C_0) y+3\,K_1-5/2\,A_0\,K_1 \\
-3\,B_0+1/2\,{A_0}^{2} K_1 + A_0\,B_0+{\frac{-4-3/2\,{A_0}^{2}+1/6\,{A_0}^{3}+13/3\,A_0}{y}}
P_1 = K_3 {y} ^ {2} + (-2 \, K_2 + A_0 \, K_2 + B_0 \, K_1 + C_0) y + 3 \, K_1-5 / 2 \, A_0 \, K_1-3 \, B_0 + 1/2 \, {A_0} ^ {2} K_1 + A_0 \, B_0
Q_1 = C_0
Q_1 = C_0 + C_1 y \Rightarrow
\begin{cases}
C_1 = 0 \\
{\frac {-4-3 / 2 \, {A_0} ^ {2} +1/6 \, {A_0} ^ {3} + 13 / 3 \, A_0} {y}} = 0
\end{cases}
\Rightarrow \begin{cases}
C_1 = 0 \\
A_0 \in \{ 2, 3, 4\}
\end{cases}
P_0 = K_4 e^{\int-\frac{Q_0}{3} dy}
P = {x}^{4}{y}^{2}+ (-4+A_0)y{x}^{3}+K_1{y}^{2}{x}^{3}+ \frac{6-7A_0+{A_0}^{2}}{2} {x}^{2}+(-3K_1+\\
A_0K_1+B_0)y{x}^{2}+K_2{y}^{2}{x}^{2}+\frac{6K_1-5A_0K_1-6B_0+{A_0}^{2}K_1+2A_0B_0}{2}\,x+ \\ (-2\,K_2+A_0\,K_2+B_0\,K_1+C_0)yx+K_3\,{y}^{2}x+D_0,
Q = C_0\,x+B_0\,{x}^{2}+A_0\,{x}^{3}-2\,{x}^{4}y
\Rightarrow \begin{cases}
Q_0 = 0 \\
P_0 = K_4
\end{cases}
Substituindo
-x^4 p_1(x) + x^4 \frac{d}{dx}p_2(x) - p_1(x) q_1(x) -p_2(x) q_0(x) = 0
{x}^{4}p_1(x) +{x}^{4} (a_1+2\,a_2\,x+3\,a_3\,{x}^{2}+4\,{x}^{3}) - (a_1\,x+a_2\,{x}^{2}+a_3\,{x}^{3}+{x}^{4})q_0(x) = 0
Obtemos
p_2(x) = a_0 + a_1 x + a_2 x^2 + a_3 x^3 + x^4
em
Q = C_0\,x+B_0\,{x}^{2}+A_0\,{x}^{3}-2\,{x}^{4}y \Rightarrow
q_0(x) \equiv b_0+b_1\,x+b_2\,{x}^{2}+b_3\,{x}^{3}
q_0(x)
q_0(x)
teria grau 5, porém:
p_1(x) = (b_3 -4)x^3 + (b_2 + a_3 b_3 -3 a_3)x^2 + (b_1 -2 a_2 + a_2 b_3 a_3 b_2)x + b_0 + a_2b_2 + a_1 b_3 + a_3 + a_3 b_1 -a_1 + \\
\frac{a_1 b_2 + a_2 b_1 + a_3 b_0}{x} + \frac{a_1b_1 + a_2a_2 b_0}{x^2} + \frac{a_1 b_0}{x^3}
continuando, teremos...
\begin{cases}
\frac{9{A_0}^{2}-{A_0}^{3}}{2} +12 -13 A_0 = 0 \\
\frac{15 A_0 B_0 - 3 A_0^2B_0}{2} -9 B_0 = 0 \\
-\frac{3 A_0 B_0^2}{2} -6 A_0 + 3 B_0^2 + 12 = 0 \\
\frac{B_0^3}{2}+6B_0 = 0
\end{cases}
P = x^4 y^2 - 2yx^3 + 3 + x^2
Q = 2 x^3 - 2 x^4 y
Método da expansão para segunda ordem
Cadeia regular em segunda ordem
y'' = z' = \dfrac{M(x,y,z)}{N(x,y,z)}
(z \equiv y')
\mathbb{D} \equiv N\partial_x + zN \partial_y + M\partial_z
\mathbb{D} = (N_0 + N_1 z + N_2 z^2) \partial_x + z((N_0 + N_1 z + N_2 z^2)\partial_y + (M_0 + M_1 z + M_2z^2)\partial_z
P = P_0(x,y) + P_1(x,y)z + P_2(x,y)z^2
N \equiv N_0(x,y) + N_1(x,y) z + N_2(x,y) z^2
M \equiv M_0(x,y) + M_1(x,y) z + M_2(x,y) z^2
Q = Q_0(x,y) + Q_1(x,y)z + Q_2(x,y)z^2 + Q_3(x,y) z^3
g_z(Q) = max(g_z(M)-1, g_z(N)+1)
= max(1, 3) = 3
Teremos o sistema...
M_{{0}} P_{{1}} +N_{{0}}
{\frac {\partial }{\partial x}}P_{{0}}
-P_{{0}} Q_{{0}}
= 0 \\
2M_{{0}} P_{{2}} +M_{{1}}
P_{{1}} +N_{{0}}
{\frac {\partial }{\partial x}}P_{{1}} +
N_{{0}} {\frac {\partial }{\partial y}}P_{{0}}
+N_{{1}} {\frac {\partial }{
\partial x}}P_{{0}} -P_{{1}} Q_{
{0}} -
P_{{0}} Q_{{1}}
= 0 \\
2M_{{1}} P_{{2}} +M_{{2}}
P_{{1}} +N_{{0}}
{\frac {\partial }{\partial x}}P_{{2}} +
N_{{0}} {\frac {\partial }{\partial y}}P_{{1}}
+N_{{1}} {\frac {\partial }{
\partial x}}P_{{1}} +N_{{1}} {\frac {\partial }{\partial y}}P_{{0}} +
N_{{2}} {\frac {\partial }{\partial x}}P_{{0}}
-P_{{1}} Q_{{1}} -P_{{2
}} Q_{{0}} -
P_{{0}} Q_{{2}} = 0 \\
2M_{{2}} P_{{2}} +N_{{0}}
{\frac {\partial }{\partial y}}P_{{2}}
+N_{{1}} {\frac {\partial }{\partial x}}P_{{2}} +
N_{{1}} {\frac {\partial }{\partial y}}P_{{1}} +N_{{2}} {
\frac {\partial }{\partial x}}P_{{1}} +N_{{2}}
{\frac {\partial }{\partial y}}P_{{0}} -
P_{{1}} Q_{{2}} -P_{{2
}} Q_{{1}} -P_{{0
}} Q_{{3}} = 0 \\
N_{{1}} {\frac {\partial }{\partial y}}P_{{2}}
+N_{{2}} {\frac {\partial}{
\partial x}}P_{{2}} +N_{{2}} {
\frac {\partial }{\partial y}}P_{{1}} -
P_{{1}} Q_{{3}} - P_{{2}} Q_{{2}} = 0 \\
N_{{2}} {\frac {\partial }{\partial y}}P_{{2}}
- P_{{2}} Q_{{3}} =0.
z'=-{\frac {{x}^{2}{y}^{2}+{x}^{2}yz+2\,xy{z}^{2}+x{z}^{3}+{z}^{4}-{
x}^{2}z+x{y}^{2}+y{z}^{2}-xy}{ 2\,\left( xy+{z}^{2}-x \right) z}}
Exemplo
D =(2\, \left( xy+{z}^{2}-x \right) z)\,\partial_x + z\,(2\, \left( xy+{z}^{2}-x \right) z)\,\partial_y +
(-{x}^{2}{y}^{2}-{x}^{2}yz-2\,xy{z}^{2}-x{z}^{3}-{z}^{4}+{x}^{2}z-x{y}^{2}-y{z}^{2}+xy)\,\partial_z
P = P_{{0}} \left( x,y \right) +P_{{1}} \left( x,y \right) z+P_{{2}} \left( x,y \right) {z}^{2}
Q = Q_{{0}} \left( x,y \right) +Q_{{1}} \left( x,y \right) z+Q_{{2}} \left( x,y \right) {z}^{2}+Q_{{3}} \left( x,y \right) {z}^{3}+Q_{{4}}
\left( x,y \right) {z}^{4}
Exemplo
\left( -{y}^{2}{x}^{2}-{y}^{2}x+xy \right) P_{{1}} \left( x,y
\right) -P_{{0}} \left( x,y \right) Q_{{0}} \left( x,y \right) = 0 \\
2\, \left( xy-x \right) {\frac {\partial }{\partial x}}P_{{0}} \left(x,y \right) -P_{{1}} \left( x,y \right) Q_{{0}} \left( x,y \right) -P_
{{0}} \left( x,y \right) Q_{{1}} \left( x,y \right) + 2\, \left( -{y}^{2}{x}^{2}-{y}^{2}x+xy \right) P_{{2}} \left( x,y \right) + \left( -y{x
}^{2}+{x}^{2} \right) P_{{1}} \left( x,y \right) = 0 \\
\\
2\, \left( -2\,xy-y \right) P_{{2}} \left( x,y \right) -xP_{{1}}
\left( x,y \right) +2\, \left( xy-x \right) {\frac {\partial }{
\partial x}}P_{{2}} \left( x,y \right) + 2\,{\frac {\partial }{\partial x}}P_{{0}} \left( x,y \right) \\
-P_{{1}} \left( x,y \right) Q_{
{2}} \left( x,y \right) -P_{{2}} \left( x,y \right) Q_{{1}} \left( x,y
\right) -P_{{0}} \left( x,y \right) Q_{{3}} \left( x,y \right) + 2\, \left( xy-x \right) {\frac {\partial }{\partial y}}P_{{1}} \left( x,y
\right) = 0 \\
-P_{{1}} \left( x,y \right) Q_{{4}} \left( x,y \right) -P_{{2}}
\left( x,y \right) Q_{{3}} \left( x,y \right) +2\,{\frac {\partial }{
\partial x}}P_{{2}} \left( x,y \right) + 2\,{\frac {\partial }{\partial y}}P_{{1}} \left( x,y \right) -2\,P_{{2}} \left( x,y \right)
=0 \\
-P_{{2}} \left( x,y \right) Q_{{4}} \left( x,y \right) +2\,{\frac {
\partial }{\partial y}}P_{{2}} \left( x,y \right)
= 0
2\,{\frac {\partial }{\partial x}}P_{{1}} \left( x,y \right) +2\,
\left( xy-x \right) {\frac {\partial }{\partial y}}P_{{2}} \left( x,y
\right) +2\,{\frac {\partial }{\partial y}}P_{{0}} \left( x,y
\right) - P_{{1}} \left( x,y \right) Q_{{3}} \left( x,y \right) \\
-P_{{2}} \left( x,y \right) Q_{{2}} \left( x,y \right) -P_{{0}} \left( x,y
\right) Q_{{4}} \left( x,y \right) - 2\,xP_{{2}} \left( x,y \right) -P_{{1}} \left( x,y \right)
= 0
2\, \left( xy-x \right) {\frac {\partial }{\partial x}}P_{{1}} \left(
x,y \right) +2\, \left( xy-x \right) {\frac {\partial }{\partial y}}P_
{{0}} \left( x,y \right) -P_{{1}} \left( x,y \right) Q_{{1}} \left( x,
y \right) \\
- P_{{2}} \left( x,y \right) Q_{{0}} \left( x,y \right) -P_{{0
}} \left( x,y \right) Q_{{2}} \left( x,y \right) +2\, \left( -y{x}^{2}
+{x}^{2} \right) P_{{2}} \left( x,y \right) + \left( -2\,xy-y \right) P_{{1}} \left( x,y \right)
= 0
Resolvendo
-P_{{2}} \left( x,y \right) Q_{{4}} \left( x,y \right) +2\,{\frac {
\partial }{\partial y}}P_{{2}} \left( x,y \right) = 0
P_2 = F_0(x) e^{\int{\frac{Q_4}{2}}dy}
\int{\frac{Q_4}{2}}dy
deve ser logaritmo ou zero
P_2 = F_0(x) e^0 = F_0(x)
Q_4 = 0
temos um grau de liberdade e P tem grau 2
\Rightarrow F_0(x) \equiv 1 \therefore
P_2 = 1
Resolvendo
-P_{{1}} \left( x,y \right) Q_{{4}} \left( x,y \right) -P_{{2}}
\left( x,y \right) Q_{{3}} \left( x,y \right) +2\,{\frac {\partial }{
\partial x}}P_{{2}} \left( x,y \right) + 2\,{\frac {\partial }{\partial y}}P_{{1}} \left( x,y \right) -2\,P_{{2}} \left( x,y \right)
=0 \\
P_1 = \frac{1}{2}\int{Q_3(x,y)}dy + y + F_1(x)
\begin{cases}
Q_3 = a_0 \\
F_1 = f_0 + f_1 x
\end{cases}
\Rightarrow P_1 = \dfrac{1}{2}a_0 y + y + f_0 + f_1 x
Resolvendo
P_0 = \int{(\frac{(a_0^2 + 3a_0 + 2)y + 2Q_2}{4})}dy + ((a_0 + 1)(f_0 + f_1 x) - f_1)y + F_1(x)
\begin{cases}
Q_2 = b_0 + b_1x + b_2y \\
F_2 = g_0 + g_1 x + g_2 x^2
\end{cases} \Rightarrow
2\,{\frac {\partial }{\partial x}}P_{{1}} \left( x,y \right) +2\,
\left( xy-x \right) {\frac {\partial }{\partial y}}P_{{2}} \left( x,y
\right) +2\,{\frac {\partial }{\partial y}}P_{{0}} \left( x,y
\right) - P_{{1}} \left( x,y \right) Q_{{3}} \left( x,y \right) \\
-P_{{2}} \left( x,y \right) Q_{{2}} \left( x,y \right) -P_{{0}} \left( x,y
\right) Q_{{4}} \left( x,y \right) - 2\,xP_{{2}} \left( x,y \right) -P_{{1}} \left( x,y \right)
= 0
P_1 =
g_0 + g_1 x + g_2 x^2
+ \frac{b_0 + f_0 + a_0 f_0 - 2 f_1}{2}y
+ \frac{2 + b_1 + a_0 f_1 + f_1}{2} xy
+ \frac{2 + 2 b_2 + a_0^2 + 3 a_0}{8} y^2
Resolvendo
2\, \left( -2\,xy-y \right) P_{{2}} \left( x,y \right) -xP_{{1}}
\left( x,y \right) +2\, \left( xy-x \right) {\frac {\partial }{
\partial x}}P_{{2}} \left( x,y \right) + 2\,{\frac {\partial }{\partial x}}P_{{0}} \left( x,y \right) \\
-P_{{1}} \left( x,y \right) Q_{
{2}} \left( x,y \right) -P_{{2}} \left( x,y \right) Q_{{1}} \left( x,y
\right) -P_{{0}} \left( x,y \right) Q_{{3}} \left( x,y \right) + 2\, \left( xy-x \right) {\frac {\partial }{\partial y}}P_{{1}} \left( x,y
\right) = 0 \\
Apenas ainda é desconhecido
Q_1(x,y)
P = xy + z^2
Q = -2 xyz - 2z^3
Função S
Na busca de polinômios de Darboux
O que é
z' = \dfrac{dz}{dx} = \phi = \dfrac{M(x,y,z)}{N(x,y,z)}
z \equiv y'
\dfrac{dz}{dy} = -S(x,y,z)
S(x,y,z) = \dfrac{I_y}{I_z}
S \equiv \dfrac{P(x,y,z)}{N(x,y,z)} \Rightarrow
\begin{cases}
\frac{dz}{dy} = -\frac{P}{N} \\
\frac{D_A[R]}{R} = P_z - N_y
\end{cases}
D_A \equiv N\partial_y - P\partial_z
Nos casos
Dado que...
\dfrac{dz}{dx} = \dfrac{M(x,y,z)}{N(x,y,z)}
\dfrac{dz}{dy} = -S(x,y,z)
\}
Tem
R
como fator integrante
\sum_i n_i \dfrac{D_A[p_i]}{p_i} = P_z - N_y
\Downarrow
Cenário 1
- Solução Liouviliana não elementar
- Só existe um polinômio no fator integrante
R = e^{\frac{1}{p_1}}p_1^{n_1}
\dfrac{D_A [R]}{R} = n_1 \dfrac{D_A[p_1]}{p_1} - \dfrac{D_A[p_1]}{p_1^2} = P_z - N_y
n_1 D_A[p_1] - \dfrac{D_A[p_1]}{p_1} = (P_z - N_y) p_1
\Downarrow
q_{A1} \equiv \dfrac{D_A[p_1]}{p_1} \Rightarrow
n_1 D_A[p_1] - q_{A1} = (P_z - N_y) p_1
\Downarrow
Cenário 2
- Solução elementar
- Existem vários polinômios no fator integrante, mas só um tem grau elevado
R = \prod_i p_i ^{n_i}
\dfrac{D_A [R]}{R} = n_1 \dfrac{D_A[p_1]}{p_1} + \sum^{k}_{i=2} n_i \dfrac{D_A[p_i]}{p_i} = P_z - N_y
n_1 q_{A1} + \sum_{i=2}^k n_i q_{Ai} = (P_z - N_y)
\Downarrow
D_A[p_1] = q_{A1} p_1
\Downarrow
com o cofator determinado
Problema dos graus máximos
- Não existe limite superior teórico para os graus dos polinômios de Darboux
- Todos os métodos são semi-algorítmicos
Problema
Objetivo
Estimar os graus em que existam polinômios de Darboux
Modelo semântico
Semântica
- Representação vetorial de objetos com métrica que faça sentido para comparações no problema, baseada nas relaçõe
Frases e Palavras
- m1: Rock and Roll Music in the 1960’s
- m2: Recipe for a Drum Roll, a Demonstration of Techniques
- m3: Drum and Bass Composition
- m4: A Perspective of Rock Music in the 90’s
- m5: Music and Composition of Popular Bands
- b1: How to make Bread and Rolls, a Demonstration
- b2: Ingredients for Crescent Roll Dough
- b3: A Recipe for Sourdough Bread
- b4: A Quick Recipe for Pizza Dough with Natural Ingredients
L = \left[\begin{matrix}0.0 & 1.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0 & 0.0\\1.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0\\0.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 1.0\\0.0 & 1.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0 & 1.0\\0.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0 & 1.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 1.0 & 0.0 & 1.0 & 0.0\\1.0 & 0.0 & 0.0 & 1.0 & 0.0 & 1.0 & 0.0 & 0.0 & 0.0\\1.0 & 0.0 & 0.0 & 0.0 & 0.0 & 1.0 & 0.0 & 1.0 & 0.0\end{matrix}\right]
- m1: Rock and Roll Music in the 1960’s
- m2: Recipe for a Drum Roll, a Demonstration of Techniques
- m3: Drum and Bass Composition
- m4: A Perspective of Rock Music in the 90’s
- m5: Music and Composition of Popular Bands
- b1: How to make Bread and Rolls, a Demonstration
- b2: Ingredients for Crescent Roll Dough
- b3: A Recipe for Sourdough Bread
- b4: A Quick Recipe for Pizza Dough with Natural Ingredients
L = U_nWV_n^T
Análise de componentes principais...
U_8 =
V_8 =
W = \displaystyle \left[\begin{matrix}2.6 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 2.19 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 1.91 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 1.59 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 1.26 & 0.0 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.99 & 0.0 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.72 & 0.0 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.24 & 0.0\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\end{matrix}\right]
Vetores semânticos
U_4 = \displaystyle \left[\begin{matrix}0.0 & -0.59 & 0.16 & -0.21\\-0.38 & -0.07 & -0.63 & 0.0\\-0.08 & -0.19 & -0.4 & 0.59\\0.0 & -0.59 & 0.16 & -0.21\\-0.02 & -0.49 & -0.05 & 0.31\\-0.19 & -0.01 & -0.35 & -0.52\\-0.36 & 0.04 & 0.39 & 0.23\\-0.56 & 0.03 & 0.03 & -0.28\\-0.59 & 0.04 & 0.31 & 0.2\end{matrix}\right]
Semântica de polinômios
p_1 = a_1 x^2 y^1 + a_2 x^2 y^2
p_2 = b_1 x y^1 + b_2 x^2 y^2
\Rightarrow
0 x y^1 + a_1 x^2 y^1 + a_2 x^2 y^2
b_1 x y^1 + 0 x^2 y^1 + b_2 x^2 y^2
\begin{pmatrix}
0 & a_1 & a_2 \\
b_1 & 0 & a_2
\end{pmatrix}
x y^1
x^2 y^1
x^2 y^2
p_1
p_2
Preparação dos dados
- Separar numeradores e denominadores das equações
- Cria semânticas separadas para numeradores e denominadores
- Normaliza os vetores
- Entradas numerador e denominador
- Reúne ambos no bloco racional
- Codificação com redução
- Prediz grau
Taxa de acerto em torno de 80%
Problemas
- É preciso treinar um modelo para cada grau
- Adicionar um novo grau implica num modelo treinado do zero
- Treinamento muito demorado
- Necessita de muitas amostras
- Semântica evolui para novos graus, deve ser reconstruída.
Modelo de transferência
Rede triplice siamesa
- Ancora e Positivo têm o mesmo grau máximo no PD, negativo têm grau diferente
- f(E) é a equação que determina o modelo
- S é o a função objetivo a ser minimizada
- Classes se baseiam no maior polinômio conhecido que compõe o fator integrante
- Redução de dimensionalidade se dá com base na distância entre classes
- Aprende a diferença entre as equações antes de aprender o padrão comum
|| f(E_a) - f(E_p)|| < ||f(E_a) - f(E_n)||
S(E_a, E_p, E_n) \equiv max(|| f(E_a) - f(E_p)|| - ||f(E_a) - f(E_n)|| + \alpha, 0)
S(E_a, E_p, E_n) \equiv || f(E_a) - f(E_p)|| - ||f(E_a) - f(E_n)|| + \alpha
ou
\Downarrow
Modelo de classificação
Resultados gerais
