A simple way to do what you want would be by using barycentric coordinates.
One problem that you may be running into is that you're trying to find some way of specifying a point in the plane as a combination of four different points. In fact, three are enough -- this is why the article I linked to talks about barycentric coordinates in the context of a triangle instead of a square. What this means is that you'll always be able to find some such expression, but it won't be unique. For example, your point at the middle of the square could be described as $25\%$ of each corner, or it could be $50\%$ $a$ and $50\%$ $c$, or any number of other things in between. If you want a unique representation for each point, you'll have to impose some other condition.
If you have the coordinates of some point $P=(x_0,y_0)$, you can find how to write it as a combination of your four points by solving a system of linear equations. Let's say we want the fraction of $a$ in your point to be $P_a$, and similarly with $b,c,d$. Then the equations we get are:
$P_a+P_b+P_c+P_d=1$ (that is, the sum of all your proportions should add up to $1$ or $100\%$), and $ x_0=0*P_a+0*P_b+1*P_c+1*P_d \\ y_0=0*P_a+1*P_b+0*P_c+1*P_d $ (that is, whatever proportion you're combining $a,b,c,d$ in to get $P$, you should combine the $x$-coordinates of $a,b,c,d$ in to get $P$'s $x$-coordinate, and similarly with $y$-coordinates).
Any numbers $P_a,P_b,P_c,P_d$ which satisfy these three equations will give you a set of viable proportions in which you can combine the points $a,b,c,d$ to get the point $P$. There are only three equations, and four variables $P_a,P_b,P_c,P_d$ to solve for; this is the source of the non-uniqueness problem I was talking about.