Como también mostró River, la implicación es equivalente a la siguiente implicación que solo involucra números reales considerando z = x + iy:
(x^2-y^2+1)^2 + 4x^2y^2 \le 4 \Rightarrow (x^3-3xy^2+3x+2)^2 + (3x^2y-y^3+3y)^2 \le 36.
Para demostrar la implicación anterior, River utilizó de manera elegante un método basado en SOSs, pero requiere adivinar varias funciones, lo cual probablemente solo se puede hacer utilizando resolutores específicamente desarrollados para SOSs y conociendo la solución óptima. Aquí, pruebo la implicación anterior utilizando herramientas de análisis elementales.
La implicación anterior puede ser probada si mostramos que el valor óptimo p^* del siguiente problema de optimización es menor o igual a 36:
\max z(x,y)=A(x,y)^2 + B(x,y)^2 sujeto a (x^2-y^2+1)^2 + 4x^2y^2 \le 4 \tag{1}
donde
A(x,y)=x^3-3xy^2+3x+2, B(x,y)=3x^2y-y^3+3y.
Aquí, probamos que p^*=36 utilizando los siguientes pasos:
1) Se puede mostrar que x^2\le 1. De hecho, de
(x^2-y^2+1)^2 + 4x^2y^2 \le 4 \Leftrightarrow (x^2+y^2+1)^2 \le 4(1+y^2), \tag{2}
tenemos x^2 \le 2\sqrt{1+y^2}-(1+y^2) \le \max_{a}(2a-a^2)= 1.
2) La ecuación \nabla Z(x,y)=0 es equivalente a
(x^2-y^2+1)A(x,y)+2xyB(x,y)=0 -2xyA(x,y)+(x^2-y^2+1)B(x,y)=0.
Esto juntos dan ((x^2-y^2+1)^2+4x^2y^2)B(x,y)=0, lo cual tiene una solución solo si B(x,y)=0 \Leftrightarrow y^2=3(1+x^2) \vee y=0. El primer y^2=3(1+x^2) resulta en ninguna solución factible (simplemente ingrésalo en (1)), pero y=0 produce la solución factible (0,0) que no puede ser óptima (compáralo con alguna solución factible trivial como (1,0) o elimínalo al final del paso 3).
3) De paso 2, podemos ver que el valor óptimo p^* debe ser alcanzado en algún punto de la curva (x^2-y^2+1)^2 + 4x^2y^2 = 4, que utilizando (2) puede ser escrito como
\big (x^2+t(x) \big )^2 = 4t(x), t(x)=y^2+1.
Entonces, t(x)=2-x^2+\sqrt{1-x^2} puede ser obtenido explícitamente como una función de x resolviendo la ecuación cuadrática anterior (tiene raíces reales si y solo si x^2\le 1) y notando que t(x)\ge 1. Como la función objetivo z(x,y) depende de y a través de y^2, obtenemos:
f(x)=\big (x^3-3x(t(x)-1)+3x+2\big)^2+(t(x)-1)\big (3x^2-(t(x)-1)+3 \big)^2\\=\left (4x^3-3x\sqrt{1-x^2}+2 \right)^2 +\left (1-x^2+\sqrt{1-x^2}\right)\left (4x^2-1-\sqrt{1-x^2} \right)^2, |x|\le 1.
f(x) [1] es una función univariante de x representada abajo [2]:
\hspace{3cm}![]()
El máximo se alcanza en x=1, lo cual se puede demostrar encontrando las soluciones de f'(x)=0 [3], que son r_1=-0.894803, r_2=-0.163729, y r_3=0.580693, y luego seleccionando la mejor entre -1,r_1,r_2,r_3,1. Por lo tanto, de t(1)=y^2+1, t(1)=1, obtenemos y=0. Por lo tanto, el punto (1,0) con z(1,0)=36, que también domina a (0,0) obtenido en el paso 2, es la solución óptima, y así p^*=36. QED
Siguiendo el comentario de River, quizás podamos demostrar que el máximo ocurre en la frontera sin cálculo o principio del módulo máximo, y el problema es usar la transformación x=\frac{2t}{1+t^2} para deshacernos de \sqrt{1-x^2}. El resultado se presenta aquí, que parece ser complejo. También intenté la transformación x=\sin t, pero el resultado proporcionado aquí todavía es complejo. ¡Cualquier comentario es bienvenido!