The solution set is a vector space and includes all functions of the form $f(x)=ax+b$. But apart from these observations I have nothing to say about the problem. If it helps, restrict the domain of $f$ to the unit interval.
The problem came up when I was trying to find the end behavior of the following martingale. Let $\epsilon_i$ be an IID sequence of fair coinflips. Set $X_0$ to $x_0$ with $0\le x_0\le 1$, and then set $X_{i+1}$ to $X_i^2$ if $\epsilon_i$ lands tails but to $2X_i-X_i^2$ if $\epsilon_i$ lands heads. Then $X_i$ is a martingale and, being $\mathrm{L}^1$-bounded, converges almost surely by a theorem of Doob. Clearly, the only numbers the martingale can converge to are $0$ and $1$. The chances of convergence to $1$ as a function of $x_0$ must satisfy the functional equation above.