Let $x_i$, $i=0,\dots n$ be the partition $P$. Let $h_j(t)=sgn(\alpha(x_{j+1})-\alpha(x_j))1_{(x_i,x_{i+1})}$ and $h=\sum h_j$. $h$ is $\alpha$ integrable since it's simple and $\int_a^b hd\alpha= \sum_1^n |\alpha(x_{i+1})-\alpha(x_i)|$ now we just want to approximate this function. For each $i$ and $\delta>0$ consider the intervals $[x_i+\delta,x_{i+1}-\delta]=J_{i,\delta}$ and consider $f_{i,\delta}=h_i|_{J_{i,\delta}}$. Take $f_{\delta}:\cup J_{i,\delta} \to \mathbb{R}$ as $f_{\delta}=\sum f_{i,\delta}$ and extend linearly (call the extension the same). It's easy to see that the $f_{\delta}$ are continuous, as $\delta \to 0, f_{\delta}(t) \to h(t)$ except possibly at the $x_i$ which is of $\alpha$-measure 0 , and $|f_{\delta}|\leq 1$, so you can apply the dominated convergence theorem to the integrals and conclude (formally you have to apply this last theorem to the integrals associated with the nonnegative measures given by the Jordan decomposition of $\alpha$).
Just to clarify, your $f$ will be a $f_{\delta}$ for $\delta(\varepsilon)$ sufficiently small.