I want to construct brownian motion using Haar Wavelets where $\varphi_n, n \in \mathbb{N}$ is an orthonormal basis of $L^2[0,1]$. I take the inner product: $$\langle f,g \rangle = \int_0^1 f(t)g(t)dt$$ Now I take $$\Phi_n(t) = \int_0^t \varphi_n(s) ds$$ This is obviously well defined since $\varphi_n \in L^2[0,1] \subset L[0,1]$.
Now I take a series of i.i.d. random normal variables $Z_n$ to create $$W_N^\varphi(t)= \sum_{n=1}^N Z_n \Phi_n(t)$$
I have shown that this is a martingale (w.r.t $\mathcal{F}_N$) for all $t \in [0,1]$. Furthermore I have shown (by the martingale convergence theorem) that $$W^\varphi(t)=\lim_{N \rightarrow \infty} W_N^\varphi(t) = \sum_{n=1}^\infty Z_n \Phi_n(t)$$ is well defined and Gaussian.
I now want to show that this has independent increments and that for $0\leq s\leq t \leq 1$: $$W^\varphi(t)-W^\varphi(s) \sim N(0,t-s)$$ Since $W^\varphi(t)$ is Gaussian, it is sufficient to show that the increments are uncorrelated, but I also don't know how to do that.
Help would be greatly appreciated!