I'm not sure exactly how to compute this efficiently, but I know that a result for the general form of the solution. A solution only exists if the range of $WB$ is contained in the range of $A$, so the first step is to pick any $W$ such that the range of $WB$ is contained in the range of $A$. Clearly this is possible with basically no computation, but perhaps there are some choices that would be numerically wise - I can't help you there. Then it is known that $V$ will exist and the general solution for $V$ will be of the form
$$ V = A^{\dagger}WB + (I - A^{\dagger}A)Y$$
where $A^{\dagger}$ is the Moore-Penrose pseudoinverse of $A$, and $Y$ is any matrix where the number of rows of $Y$ is equal to the number of columns of $A$, and the number of columns of $Y$ is equal to the number of columns of $B$. So you can basically just pick any $Y$.