A pesar de que Erick Wong proporciona una gran solución, voy a publicar uno que se me ocurrió, para la posteridad.
Deje $F_n(x)=\sum_{i=1}^{n}a_i\cos(ix+\phi_i)$$a_i$$[0,1]$$\phi_i$,$\mathbb{R}$. Entonces, para todos $x$, $\frac{F_n(x)+1}{2\sum a_i}$, será en $[0,1]$ desde $\cos$$[-1,1]$. Por último, definir $p(x)$ como la Serie de Taylor para $\frac{F_n(x)+1}{2\sum{a_i}}$,$x=\frac{1}{2}$, con un grado suficientemente alto, de tal manera que $p$ se queda en $[0,1]$.
No es tan limpio como yo esperaba, pero funciona bastante bien y, probablemente, puede ser simplificado. La idea surgió a partir de un polinomio de ajuste aleatorio Serie de Fourier, sin términos en $\sin$ desde que se complicaron las cosas. Los cambios de fase eran necesarias ya que cada una de las $F_n$ estaba saliendo de la $y$ eje en la predicción de puntos.
En el diagrama de abajo, he mostrado, superpuestos sobre $20$ $30$grado $9$ polinomios, el uso de este método, con $n=4$$\phi_i$$[-3,3]$. Claramente que se mantengan dentro de los límites de la unidad de la plaza, rellenar el espacio muy bien y son bastante aleatorios, así que me gustaría considerar la posibilidad de que este método sea un éxito.
Cuando vemos lo que parece, con cerca de $100$ o de $p$, es claro que se agrupan en torno a $y=\frac{1}{2}$, que es lo esperado, por lo que el método no proporciona polinomios con igualdad de oportunidad de estar en cada punto en la unidad de la plaza pero son razonablemente al azar, no obstante.