He estado jugando con un argumento de concavidad y el teorema de Jensen. Para $n=2$ Creo que esto funciona, como se muestra a continuación. Desafortunadamente, no funciona para $n > 2$ . De hecho, para $n > 2$ la restricción a cualquier conjunto convexo $x_1 + \dotsc + x_n = \lambda$ se puede demostrar que no cóncavo. Así que, aquí va el resultado parcial.
Dejemos que
$$ f(x_1, x_2) = \frac{a_1^2x_1 + a_2^2x_2}{a_1x_1+a_2x_2}. $$
Entonces $f$ restringido al segmento convexo $I=\{(x_1,x_2) \mid x_1,x_2 \geq 0 \textrm{ and } x_1+x_2=1\}$ es una función cóncava. Esto se puede comprobar de múltiples maneras, por ejemplo, calculando explícitamente la segunda derivada
$$ \frac{\partial^2}{\partial x^2} f(x, 1-x) = \frac{-2a_1a_2(a_1-a_2)^2}{(a_1x+a_2(1-x))^3} $$
y observando que es no positivo para $x \in [0,1]$ . Supongamos que la distribución conjunta de $(Z_1,Z_2)$ induce una distribución de probabilidad en $I$ entonces es simétrico bajo $(x_1,x_2) \mapsto (x_2, x_1)$ . En particular, los valores esperados para las coordenadas en $I$ son iguales y por lo tanto $\mathbb{E}(x_1) = \mathbb{E}(x_2) = \tfrac{1}{2}$ . Entonces, por el teorema de Jensen para funciones cóncavas sobre $I$ :
$$ \mathbb{E}(f(x_1,x_2)) \leq f(\mathbb{E}(x_1),\mathbb{E}(x_2)) = f(\tfrac{1}{2},\tfrac{1}{2}) = \frac{a_1^2+a_2^2}{a_1+a_2}. $$
Ahora $f$ es una función radial y el mismo razonamiento se aplica a todos los segmentos $x_1+x_2 = \lambda$ para $\lambda > 0$ . Y aunque la distribución inducida depende del parámetro $\lambda$ la desigualdad anterior se cumplirá para cada uno de estos segmentos y, por tanto, también para la distribución conjunta $(Z_1,Z_2)$ .