44 votos

El eslabón perdido: una desigualdad

He estado trabajando en un proyecto y resultó un par de resultados relevantes, pero se quedó atascado en un problema difícil:

Conjetura. Si $2\leq n\in\mathbb{N}$ e $0<x<1$ es un número real, entonces $$F_n(x)=\log\left(\frac{(1+x^{4n-1})(1+x^{2n})(1-x^{2n+1})}{(1+x^{2n+1})(1-x^{2n+2})}\right)$$ es una función convexa de $x$.

Tengo una heurística argumento. Se puede ayudar con una rigurosa prueba o herramientas valiosas?

Además de la motivación. Si usted tiene éxito con este, entonces voy a ser un honor contar con usted como un co-autor en este trabajo. El problema en sí mismo puede ser encontrado en la Sección 4.

Nota. Escribir $F=\log\frac{P}Q$,, a continuación, $F''>0$ equivale a la positividad del polinomio $$V:=PQ^2P''+(PQ')^2-P^2QQ''-(P'Q)^2.$$

21voto

psweeney Puntos 16

Fix $n$, y deje $V(x)$ ser el polinomio de la pregunta, para la que queremos mostrar que es positiva en el intervalo abierto $(0,1)$. Conjunto \begin{align*} a(n) &= 1024n^2 - 4096/3n + 320\\ b(n) &= 12288n^3 - 14336n^2 + 4864/3n + 256\\ W(x) &= x^4\left(V(x)-(n+1)^2(2n+1)^2x^{24n-1}(1-x)^4(a(n)x + b(n)(1-x))\right). \end{align*}

Para todos los $n\ge12$ pretendemos que $(1-x)^6$ divide $W(x)$, y que todos los coeficientes de $W(x)/(1-x)^6$ son no negativos. En particular, $W(x)$ es positivo para $x>0$. Como $a(n)$ e $b(n)$ es positivo para todos los $n\ge1$, podemos ver que $V(x)$ es positivo para todos los $0<x<1$, siempre que $n\ge12$. Los casos de $2\le n\le 11$ son fácilmente controlado directamente, por ejemplo, observa que los coeficientes de $(1+x)^{\deg V}V(1/(1+x))$ son positivos.

No estoy seguro acerca de los más fáciles de tipo para demostrar la afirmación acerca de los coeficientes de $W(x)/(1-x)^6$. Se calcula a través de ellos \begin{equation} \frac{W(x)}{(1-x)^6}=W(x)\sum_{k\ge0}\binom{k+5}{5}x^k. \end{equation}

Todo lo que sigue es asistida/confirmado por el Sabio código de abajo.

Los exponentes de la $W(x)$ tiene la forma $2ni+j$ para $0\le i\le 12$ e $0\le j\le 10$. Deje $a_{i,j}$ ser el coeficiente correspondiente. Es un polinomio en $n$.

Tenemos \begin{equation} \frac{W(x)}{(1-x)^6} = \sum_{i,j}a_{i,j}(n)x^{2ni+j}\sum_k\binom{k+5}{5}x^k. \end{equation}

Con $r\ge0$ e $0\le s\le 2n-1$, el coeficiente de $x^{2nr+s}$ en $W(x)\frac{1}{(1-x)^6}$ es la suma de todos los $a_{i,j}(n)\cdot\binom{nr+s-ni-j+y}{5}$ para todos los pares de $(i,j)$ donde $i<r$ e $0\le j\le 10$ o $i=r$ e $j\le s$.

Podemos distinguir los casos de $0\le s\le 9$ de los casos $s\ge10$. En estos últimos casos, escribimos $s=10+y$. Recordemos que $s\le 2n-1$, lo $10+y=2n-1-o$ para un no negativo $o$. Sobre la sustitución de $n$ con $(11+y+1)/2$, obtenemos polinomios en $n$ e $o$.

En estos casos, los coeficientes son polinomios en $n$, que tienen coeficientes positivos sobre la sustitución de $n$ con $n+12$. Así que estos coeficientes son positivos para todos los $n\ge12$.

En el último de los casos resulta que en todos los casos menos uno de los coeficientes son positivos. Así que para no negativo $y$ e $o$, los valores son no negativos.

La única excepción que hemos de ver que si multiplicamos con un adecuado polinomio, la resultante de los coeficientes son positivos.

Comentario (en respuesta a su pregunta): es un hecho conocido que un polinomio $V$ es positivo en $(0,1)$ si y sólo si es un no negativo combinación lineal de polinomios $x^i(1-x)^j$. El problema es que puede implicar que $i+j$ es mayor que $\deg V$. (Esto no sucede aquí, sin embargo.) He utilizado un LP de solver para jugar un poco con $V$ e $\frac{V}{(1-x^2)^4x^{2n-2}}$. A partir de que un patrón se mostró que conduzcan a una solución.

Por el camino, $V$ realidad parece ser una combinación lineal no negativa de $x^i(1-x)^{\deg V-i}$. Esto es equivalente a decir que todos los coeficientes de $(1+x)^{\deg V}V(1/(1+x))$ son no negativos. No es difícil calcular las expresiones explícitas para estos coeficientes, pero no lo veo fácil argumentos de por qué no pueden ser negativos.

# Formally compute W
var('X N')
P = (1+X^(4*N-1))*(1+X^(2*N))*(1-X^(2*N+1))
Q = (1+X^(2*N+1))*(1-X^(2*N+2))
V = P*Q^2*P.diff(X,2)+(P*Q.diff(X))^2-P^2*Q*Q.diff(X,2)-(P.diff(X)*Q)^2
a = 1024*N^2 - 4096/3*N + 320
b = 12288*N^3 - 14336*N^2 + 4864/3*N + 256
W = X^4*(V - (2*N + 1)^2 *(N + 1)^2*X^(24*N-1)*(1-X)^4*(a*X + b*(1-X)))

# Check that (1-X)^6 divides W(X)
print all(W.diff(X,i)(X=1).polynomial(QQ) == 0 for i in [0..5])

# Compute the coefficients a_ij for the exponents 2ni+j of W. Somewhat
# clumsy, as I don't know how to deal with polynomials where exponents
# are symbolic expressions.
#
# l is the list of summands of W
l = [z.canonicalize_radical() for z in W.expand().operands()]
K.<n> = QQ[]
aij = {(i,j):K(0) for i in [0..12] for j in [0..10]}
for term in l:
    c = term(X=1)                        # get coefficient of term
    e = (term.diff(X)/c)(X=1)            # get exponent of term
    c = K(c.polynomial(QQ))              # convert c to proper polynomial in n
    i, j = ZZ(e.diff(N))//2, ZZ(e(N=0))  # Clumsy method to compute pairs (i,j)
    aij[i,j] += c

# Check if coefficients aij[i,j] were correctly computed
Wnew = sum(c(n=N)*X^(2*N*i+j) for (i,j),c in aij.items())
print (W-Wnew).canonicalize_radical().is_trivial_zero()

def bino(k): # binomial coefficient binom(-6,k)=binom(k+5,5)
    return prod(k+5-z for z in range(5))/120

# compute the coefficients of W(X)/(1-X)^6
K = K.extend_variables(('y', 'o'))
K.inject_variables()
for r in [0..12]:
    for s in [0..10]:
        d = 1 if s == 10 else 0
        # compute coefficient of X^(2nr+s+d*y)
        f = sum(aij[i,j]*bino(2*n*(r-i)+s+d*y-j) for i in [0..r]
                for j in [0..s if i == r else 10])
        # check non-negativity of the polynomial f
        if d == 0:
        # Checks if the coefficient of X^(2nr+s), which is a
        # polynomial in n, has nonnegative coefficients upon replacing
        # n with n+12
            f = f(n=n+12)
            if min(f.coefficients()+[0]) < 0:
                print "False"
        else:
        # The coefficient of X^(2nr+10+y), which is a polynomial in n
        # and y, has to be nonnegative for 0 <= 10+y <= 2n-1, so 10+y
        # = 2n-1-o for non-negative o. So upon replacing n with
        # (11+y+o)/2, we get a polynomial in y and o which has to be
        # nonnegative for all nonnegative y and o. In all but one
        # case, this holds because the coeffcients are non-negative.
            f = f(n=(11+y+o)/2)
            if min(f.coefficients()+[0]) < 0:
                c = o^2 + 23*o*y + 1360*y^2 + 99*o + 340*y + 1675
                print min(c.coefficients()+(c*f).coefficients()) >= 0

11voto

Iosif Pinelis Puntos 24742

Esta prueba se basa principalmente en algunas de las sugerencias en la primera versión de esta respuesta. Además, estamos construyendo aquí en el comentario de Peter Mueller sobre el polinomio $a$ se define a continuación. Resulta que técnicas similares pueden ser utilizados para todos los demás no trivial de polinomios derivadas de la prueba, la cual nos permite reemplazar a los múltiples usos de Mathematica comando Reducir[] en la última versión de esta respuesta mostrando que los coeficientes de la mencionada relevantes polinomios son todos no negativos. Yo todavía mantienen un cierto uso de Reducir [] - pero sólo para rápidamente hacer la rutina de verificación para $n\le11$.

Tenga en cuenta que \begin{equation} G:= F_n''(x) x^{2-2 n} \left(x^{2 n}+1\right)^2 \left(x^{4 n}+x\right)^2 \left(x^{2 n+1}-1\right)^2 \left(x^{2 n+1}+1\right)^2 \left(x^{2 n+2}-1\right)^2 \end{equation} es un polinomio en $n$, $x$, y $y:=x^{2n}$, de grado $2$ en $n$: \begin{equation} G=a+bn+cn^2, \end{equation} donde $a,b,c$ son ciertos polinomios en $x,y$.

La verificación de $G\ge0$ para $n=2,\dots,11$ es sencillo. Así, supongamos $n\ge12$, de donde $0<y<x^{24}<x^{12}$. Además, siempre se asume que $0<x<1$, $0<y<x^{24}$ , \begin{equation} n=\frac{\ln y}{2\ln x}; \tag{1} \end{equation} esta última relación es simplemente otra forma de la igualdad de $y=x^{2n}$.

Tenga en cuenta que $0<y<1\ \&\ 0<x<1$ implica $x=1/(1+u)\ \&\ y=1/(1+v)$ para algunos $u,v\ge0$. Como sugiere Peter Mueller, resulta que \begin{equation} a\; (1 + u)^{10} (1 + v)^{11}|_{x=1/(1 + u), y=1/(1 + v)} \end{equation} es un polinomio en $u,v$ con todos los coeficientes no negativos. Los detalles de todos los cálculos se puede ver en el cuaderno de Mathematica y su imagen pdf

Así, \begin{equation} a\ge0 \end{equation} para todos los $x,y$ en $(0,1)$.

Siguiente, introducir \begin{equation} x_* := \frac{869}{1000},\quad x_0:=\frac{7655}{10000},\quad y_0:= x_0^{24}\approx0.00164,\quad x_1 := \frac{985}{1000},\quad y_1:= x_1^{24}\approx0.696. \end{equation}

Los siguientes casos/subcases/subsubcases son exhaustivos:

Caso 1: $y\le y_1$. A continuación, la condición de $0<y\le y_1\ \&\ 0<x<1\ \&\ 0<y<x^{24}$ implica $x=x_1/(1+u)\ \&\ y=x^{24}/(1+v)$ o $x=x_1+(1-x_1)/(1+u)\ \&\ y=y_1/(1+v)$ para algunos $u\ge0$ e $v\ge0$, dependiendo de si $x<x_1$ o $x\ge x_1$. Resulta que \begin{equation} c\; (1 + u)^{248} (1 + v)^{10}|_{x=x_1/(1+u)\ \&\ y=x^{24}/(1+v)} \end{equation} y \begin{equation} c\; (1 + u)^{10} (1 + v)^{10}|_{x=x_1+(1-x_1)/(1+u)\ \&\ y=y_1/(1+v)} \end{equation} son cada uno de un polinomio en $u,v$ con todos los coeficientes no negativos. Así, \begin{equation} c\ge0 \quad\text{in Case 1}. \end{equation} Así, para el derivado $G'_n$ de % de $G$ en $n$, tenemos \begin{equation}\tag{2} G'_n=b+2cn\ge b+24c=(G'_n)|_{n=12}. \end{equation}

Ahora tenemos que distinguir dos subcases del Caso 1, con más subsubcases de Subcase 1.1:

Subcase 1.1: $y\le y_1$ e ($y\ge y_0$ o $x\le x_*$).

Subsubcase 1.1.1: $y\le y_1$ e $y\ge y_0$ e $x\ge x_1$. A continuación, $y\le y_1\ \&\ y\ge y_0\ \&\ x\ge x_1 \ \&\ 0<x<1\ \&\ 0<y<x^{24}$ implica $$x=x_{111}:=x_1+(1-x_1)/(1+u)\ \&\ y=y_{111}:=y_0+(y_1-y_0)/(1+v)$$ para algunos $u\ge0$ e $v\ge0$. Resulta que \begin{equation} (G'_n)|_{n=12}\; (1 + u)^{10} (1 + v)^{11}|_{x=x_{111}, y=y_{111}} \end{equation} y \begin{equation} G|_{n=12}\; (1 + u)^{10} (1 + v)^{11}|_{x=x_{111}, y=y_{111}} \end{equation} son cada uno, de nuevo, un polinomio en $u,v$ con todos los coeficientes no negativos.

Subsubcase 1.1.2: $y\le y_1$ e $y\ge y_0$ e $x<x_1$. A continuación, $x=x_{112}:=x_0+(x_1-x_0)/(1+u)\ \&\ y=y_{112}:=y_0+(x^{24}-y_0)/(1+v)$ para algunos $u\ge0$ e $v\ge0$. Resulta que \begin{equation} (G'_n)|_{n=12}\; (1 + u)^{272} (1 + v)^{11}|_{x=x_{112}, y=y_{112}} \end{equation} y \begin{equation} G|_{n=12}\; (1 + u)^{272} (1 + v)^{11}|_{x=x_{112},y=y_{112}} \end{equation} son cada uno, de nuevo, un polinomio en $u,v$ con todos los coeficientes no negativos.

Subsubcase 1.1.3: $y\le y_1$ e $x\le x_*$. A continuación, $x=x_{113}:=x_*/(1+u)\ \&\ y=y_{113}:=x^{24}/(1+v)$ para algunos $u\ge0$ e $v\ge0$. Resulta que \begin{equation} (G'_n)|_{n=12}\; (1 + u)^{272} (1 + v)^{11}|_{x=x_{113},y=y_{113}} \end{equation} y \begin{equation} G|_{n=12}\; (1 + u)^{272} (1 + v)^{11}|_{x=x_{113},y=y_{113}} \end{equation} son cada uno, de nuevo, un polinomio en $u,v$ con todos los coeficientes no negativos.

Por eso, $G|_{n=12}\ge0$ e $(G'_n)|_{n=12}\ge0$ en los tres subsubcases de Subcase 1.1. En vista de $(2)$, tenemos $G\ge0$ en Subcase 1.1.

Subcase 1.2: $y\le y_1$ e $y<y_0$ e $x>x_*$. Luego, por supuesto, la condición de $y\le y_1$ es redundante. Dejamos aquí \begin{equation} \rho(x):=\frac7{2(1-x)}. \end{equation} La condición $x_*<x<1\ \&\ 0<y<y_0$ implica $x=x_{12}:=x_*+(1-x_*)/(1+u)$ e $y=y_{12}:=y_0/(1+v)$ para algunos $u\ge0$ e $v\ge0$. Resulta que \begin{equation} (G'_n)|_{n=\rho(x)/2}\; u\,(1 + u)^{10} (1 + v)^{11}|_{x=x_{12},y=y_{12}} \end{equation} y \begin{equation} G|_{n=\rho(x)/2}\; u^2\,(1 + u)^{10} (1 + v)^{11}|_{x=x_{12},y=y_{12}} \end{equation} son cada uno, de nuevo, un polinomio en $u,v$ con todos los coeficientes no negativos. Así, de manera similar a Subcase 1.1, para todos los $n\ge\rho(x)/2$ obtenemos $G\ge G|_{n=\rho(x)/2}\ge0$. En Subcase 1.2, queda señalar que para nuestro particular $n$, como en $(1)$, que en realidad ha $n\ge\rho(x)/2$. Esto se deduce porque para $y<y_0$ e $x>x_*$ \begin{equation} \ln\frac1y>\ln\frac1{y_0}>1.08\frac7{2(1-x_*)}\,\ln\frac1{x_*} >\frac7{2(1-x)}\,\ln\frac1{x}=\rho(x)\ln\frac1{x}. \end{equation} Por eso, $G\ge0$ en Subcase 1.2 así.

Queda por considerar

Caso 2: $y>y_1$. A continuación, la condición de $y>y_1\ \&\ 0<x<1\ \&\ 0<y<x^{24}[<x^{12}]$ implica $x=x_2:=x_1+(1-x_1)/(1+u),y=y_2:=y_1+(x^{12}-y_1)/(1+v)$ para algunos $u\ge0$ e $v\ge0$. Resulta que \begin{equation} b\;(1 + u)^{140} (1 + v)^{11}|_{x=x_{2},y=y_{2}} \end{equation} también es un polinomio en $u,v$ con todos los coeficientes no negativos. Por lo tanto, $b\ge0$ en el Caso 2. Por lo tanto, si $c\ge0$, entonces trivialmente $G\ge0$ (desde $a\ge0$ siempre). Por lo tanto, sin pérdida de generalidad, $c<0$ y, por tanto, $G$ es cóncava en $n$. Así, es suficiente para soporte de la $n$ como en el $(1)$ entre algunas de las sencillas expresiones racionales $n_1$ e $n_2$ tal que $G|_{n=n_1}\ge0$ e $G|_{n=n_2}\ge0$.

En el Caso 2, $y$ es "cercano" a $1$ y, por lo tanto, en vista de la desigualdad de $y<x^{24}$, por lo que es $x$: $x>y_1^{1/12}=x_1=0.985$. Así que, es un poco más cómodo aquí para cambiar las variables de $x, y$ a "pequeñas variables" $u:=1-x$ e $v:=1-y$. A continuación, $0<v<1-y_1\approx0.304$ e $1-u=x>y^{1/24}=(1-v)^{1/24}>1-v/3$, de donde $u<v/3$. Vamos ahora a $t:=u/v$, por lo que el $0<t<1/3$. Ello se desprende que, en el Caso 2, que $t=(1/3)/(1+r)$ e $v=(1-y_1)/(1+s)$ para algunos $r\ge0$ e $s\ge0$.

Los mencionados soportes de $n_1$ e $n_2$ son, respectivamente, $\ell_1/2$ e $\ell_2/2$, donde
\begin{equation} \ell_1:=\frac1t,\quad\ell_2:=\frac{(2 - v) (1 - t v)}{t (1 - v) (2 - t v)}. \end{equation} No es difícil ver (más detalles sobre esto en el mismo cuaderno de Mathematica), que, de hecho,$n_1 < n < n_2$. Queda por demostrar que $G|_{n=n_1}\ge0$ e $G|_{n=n_2}\ge0$.

Resulta que \begin{equation} (G'_n)|_{n=\ell_1/2;\,x=1-tv,\,y=1-v;\,t=(1/3)/(1+r),\,v=(1-y_1)/(1+s)}\; (1 + r)^{10} (1 + s)^{19} \end{equation} y \begin{equation} (G'_n)|_{n=\ell_2/2;\,x=1-tv,\,y=1-v;\,t=(1/3)/(1+r),\,v=(1-y_1)/(1+s)}\; (1 + r)^{10} (1 + s)^{19}\; q \end{equation} son cada uno de un polinomio en $r,s$ con todos los coeficientes no negativos, donde \begin{multline*} q:= (1 + r)^{10} (1 + s)^{19} \\ \times(11673186598630578538556565100133681446610566511878526881 \\ + 16777216000000000000000000000000000000000000000000000000 s)^2 \\ \times\big(31853088866210192846185521700044560482203522170626175627 \\ + 33554432000000000000000000000000000000000000000000000000 (r + s + r s)\big)^2. \end{multline*}

La prueba está completa.

10voto

Ryan Montgomery Puntos 5153

Yo, el boceto de un método que debe mostrar $F_n$ es convexa para $n$ lo suficientemente grande.

Una expansión de Taylor le da a ese $F_n(x) = x^{2n} (1-x)^2 + O(x^{4n-1})$, y por lo $$F_n''(x) = 2x^{2n-2}[2n^2(1-x)^2 + x^2 - n(1-x)(1+3x)] + O(x^{4n-3}).$$ Esto debería mostrar que $F_n''(x) > 0$ para $n$ grandes, siempre $x < 1-\frac{C}{n}$ para algunos fijos $C>0$.

Por otro lado, de acuerdo a Mathematica, $F_n''(1) = 4n^2 - \frac{16}{3} n + \frac{5}{4}$, que tiene su más grande raíz en $\frac{8 + \sqrt{19}}{12} \approx 1.02991$. (Por cierto, parece que no hay ninguna razón por la que $n$ debe limitarse a ser un número entero, y esto está en línea con algunas parcelas que me indicaba $F_n(x)$ es convexa para $n > 1.03$.) Por otra parte, Mathematica da $$\lim_{n\rightarrow \infty} n^{-2}F_n''(1-\frac{y}{n}) = \frac{16 e^{4y}}{(1+e^{4y})^2} > 0.$$

Poniendo estos dos pedazos juntos, debe ser capaz de demostrar $F_n(x)$ es convexa en $(0,1)$ para $n$ lo suficientemente grande. Desafortunadamente, puede ser muy doloroso trabajo todos los términos de error de manera explícita.

6voto

Peter Puntos 1681

No es una respuesta, solo una gráfica de$F_n(x)$ para$n=2,\ldots,10$:


Fn (x)
Fn (x) 2

5voto

skrebbel Puntos 5183

Esta no es una respuesta, sólo una reformulación de algo que se ve tal vez manejable (a lo largo de la línea del comentario de G. Paseman). Definir, por $m\geq 1$ y $x\in (0,1)$, $$g_m(x):=\frac{d^2}{dx^2}[\log(1+x^m)]=\frac{x^{m-2}(m(m-1)-mx^m)}{(1+x^m)^2}$$ and $$h_m(x):=\frac{d^2}{dx^2}[\log(1-x^m)]=-\frac{x^{m-2}(m(m-1)+mx^m)}{(1-x^m)^2}.$$ Then, for $n\geq 2$, $$F_n''(x)=g_{4n-1}(x)+g_{2n}(x)+h_{2n+1}(x)-g_{2n+1}(x)-h_{2n+2}(x).$$

Dos observaciones

  • La contribución de la primera, segunda, y última términos no es negativo para todos los $n\geq 2$ e $x\in(0,1)$, mientras que en el resto de los términos no es positivo.
  • Por lo tanto, suficiente para mostrar la siguiente: para $n$ e $x$ como en el anterior, $$g_{2n+1}(x)-h_{2n+1}(x)\leq_? g_{4n-1}(x)+g_{2n}(x)-h_{2n+2}(x).$$

Agregado: Desde que me parece que no puede agregar un comentario, sólo quiero señalar que, si bien es, sin duda (siempre) es posible me dan algo, acabo de volver a verificar el álgebra con arce, y creo que las expresiones son correctas (no parece ser un factor que falta de m). Otra nota (además de la simetría observada en el comentario por Paseman) es el (trivial) de identidad $h_{2m}=g_m+h_m$, dando lugar a otra forma de la envolvente: $$g_{2n+1}(x)+g_{n+1}(x)-h_{2n+1}(x)\leq_? g_{4n-1}(x)+g_{2n}(x)-h_{n+1}(x).$$

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X