Seguimos el método de este trabajo, utilizando la transformada de Fourier.
Después de su normalizaciones,
$$
\hat f (\omega) = \frac{1}{2\pi} \int_{\mathbb R} f(x) e^{-i \omega x} \, dx, \qquad f(x) = \int_{\mathbb R} \hat f(\omega) e^{i \omega x} \, d\omega,
$$
de los hechos clave que se utilizarán son:
- La identidad de $\cos(bx) = \frac{e^{i b x} + e^{-i b x}}{2}$
- Si $f(x) = e^{ibx}$,$\hat f(\omega) = \delta_b(\omega)$, la delta de Dirac centrado en $\omega = b$.
- La fórmula de la convolución $\widehat{fg} = \hat f * \hat g$.
- La identidad de $\delta_a * \delta_b = \delta_{a+b}$.
La aplicación de estos a su problema,
\begin{align*}
\prod_{k=1}^n \cos\left(\frac{x}{3^k}\right) &= \prod_{k=1}^n \frac{e^{i \frac{x}{3^k}} + e^{-i \frac{x}{3^k}}}{2} \\
\implies \widehat{\prod_{k=1}^n \cos\left(\frac{x}{3^k}\right)} &= \frac{1}{2^n} \left(\delta_{\frac{1}{3^1}} + \delta_{-\frac{1}{3^1}} \right) * \cdots * \left(\delta_{\frac{1}{3^n}} + \delta_{-\frac{1}{3^n}} \right)(\omega) \\
&= \frac{1}{2^n} \sum_{p \in P_n} \delta_p(\omega)
\end{align*}
donde $P_n$ es el conjunto de $2^n$ puntos dado por $P_n = \left\{\sum_{k=1}^n \frac{b_k}{3^k} \;\middle|\; b_k \in \{-1,1\} \; \forall k \right\}$. Tenga en cuenta que los extremos de $\bigcup_{n=1}^\infty P_n$ están delimitadas bien por $\pm \sum_{k=1}^\infty \frac{1}{3^k} = \pm \frac{1/3}{1 - 1/3} = \pm \frac{1}{2}$, y cada punto en $\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]$ es el límite de alguna secuencia $(p_n)_{n=1}^\infty$ donde $p_n \in P_n$. Por lo tanto $\tfrac{1}{2^n} \sum_{p \in P_n} \delta_p(\omega)$ tiende a la densidad uniforme de la masa total $1$$\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]$$n\to\infty$, el cual es dado por la función del indicador de $\chi_{\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]}(\omega)$:
\begin{align*}
\widehat{\prod_{k=1}^\infty \cos\left(\frac{x}{3^k}\right)} &= \chi_{\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]}(\omega).
\end{align*}
Tomando inversa de Fourier, se obtiene
$$
\prod_{k=1}^\infty \cos\left(\frac{x}{3^k}\right) = \int_{\mathbb R} \chi_{\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]}(\omega) e^{i \omega x} \, d\omega = \int_{-1/2}^{1/2} e^{i \omega x} \, d\omega = \frac{\sin\left(\tfrac{x}{2}\right)}{\tfrac{x}{2}}.
$$
ACTUALIZACIÓN: he aquí por qué el anterior es incorrecto. La afirmación de que "cada punto en $\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]$ es el límite de alguna secuencia $(p_n)_{n=1}^\infty$ donde $p_n \in P_n$" es incorrecto; de hecho, los puntos pertenecientes a $P = \bigcup_n P_n$ son los que han equilibrado ternario de expansión consta de los dígitos $1$ $-1$ pero no $0$. No estoy seguro de cómo describir mejor la $P$, pero queremos continuar como antes por la integración de más de $P$ (en lugar de todos los de $\left[ - \tfrac{1}{2}, \tfrac{1}{2} \right]$) para calcular la inversa de la transformada de Fourier dando la respuesta correcta.