I suppose it depends on what you mean by proof here. If you are familiar at all with continuum mechanics you can derive that relationship based on the assumption that you have a Newtonian fluid and a 2-D flow.
To do that, you would start out with Cauchy's equations of motion, which are
[tex]\nabla\cdot\mathbf{T} + \rho \vec{B} = \rho \vec{a}.[/tex]
Here, [itex]\mathbf{T}[/itex] is the stress tensor, [itex]\rho[/itex] is the density, [itex]\vec{B}[/itex] is the body force vector and [itex]\vec{a}[/itex] is the acceleration vector. Basically this is just an expression of Newton's law. It can be rewritten in differential form:
[tex]\dfrac{\partial T_{ij}}{\partial x_j} + \rho B_i = \rho a_i.[/tex]
The hydrostatic pressure in a fluid can be represented by [itex]-p\delta_{ij}[/itex] where [itex]\delta_{ij}[/itex] is the Kronecker delta, and the pressure is always going to be a part of the stress tensor and always in the diagonal terms. If the fluid is undergoing anything other than rigid body motion, there has to be another component of the stress tensor, so the stress tensor is expanded to include a term representing the stresses that will pop up in a fluid due to the rate of deformation (since a fluid is defined as having stress dependent on the rate of deformation rather than the deformation itself).
[tex]T_{ij} = -p\delta_{ij} + T_{ij}^{\prime}[/tex]
where [itex]T_{ij}^{\prime}[/itex] represents the portion of the stress tensor that is a result of the rate of deformation of the fluid under shear. The rate of deformation is
[tex]D_{ij} = \dfrac{1}{2}\left( \dfrac{\partial v_i}{\partial x_j} + \dfrac{\partial v_j}{\partial x_i} \right).[/tex]
By the definition of a Newtonian fluid, the stress due to deformation of the fluid varies linearly with the rate of deformation and the fluid is isotropic. From continuum mechanics, you can show then that
[tex]T_{ij}^{\prime} = \lambda\delta_{ij}D_{kk} + 2\mu D_{ijj}[/tex]
where [itex]\lambda[/itex] and [itex]\mu[/itex] are viscosity coefficients with [itex]\mu[/itex] being the typical dynamic viscosity that you see all over the place. If you go look at the diagonal elements of [itex]T_{ij}^{\prime}[/itex], you get
[tex]T_{ii} = (3\lambda + 2\mu)D_{kk}.[/tex]
Which means that [itex]\lambda + 2\mu/3[/itex] is a proportionality constant known as bulk viscosity relating viscous normal stresses to the rate of dilatation of the fluid. That rate, [itex]D_{kk}[/itex] is zero for an incompressible fluid ([itex]\nabla\cdot\vec{v} = D_{kk} = 0[/itex]), so for an incompressible fluid, you are left simply with a stress tensor of
[tex]T_{ij} = -p\delta_{ij} + 2\mu D_{ij}.[/tex]
If you go and plug this into the Cauchy equations I mentioned before, you get the Navier-Stokes equations. Otherwise, just consider a 2-D, incompressible flow with [itex]v_1 = v_1(x_2)[/itex] and [itex]v_2 = v_3 = 0[/itex]. If you want to look at the shear component [itex]T_{12}[/itex] (the only independent nonzero component), you can simply note that
[tex]\tau = T_{12} = 2\mu D_{12} = 2\mu \left(\dfrac{1}{2}\dfrac{d v_1}{d x_2}\right) = \mu\dfrac{d v_1}{d x_2}.[/tex]
This, in [itex](x,y)[/itex] coordinates, is
[tex]\tau_{xy} = \mu\dfrac{du}{dy}.\hspace{4em}\square[/tex]
Except for your sign convention and your typo where you swapped [itex]u[/itex] with [itex]v[/itex], that is what you have. Is this the sort of thing you were looking for?