You are right of course, and indeed for Sourcegeeks resistor values this will have a large effect. However, if he were to choose his resistor values carefully he could minimise this to a large extent without employing your ideal, yet component requiring, fixed voltage reference.
If the resistor values of R1 and R4 are much larger than R2 and R3, V1 would have a very small effect on V3 and visa versa. If you rearrange the second formula above you get:
Code:
V3 = [V1*R2*R3]/[R1*R2 +R3*(R1+R2)] + [Vdd*R1*R3]/[R1*R2 +R3*(R1+R2)]
If you look at the top line of the second term/fraction you'll see it contains R1, where the top line of the first term contains only R2 and R3. Also note Vdd will always be much larger than V1. Hence, if R1 is much bigger than R2 and R3, the second term will largely dominate, comprising most of the value of V3.
Thats what i meant by saying that decreasing values of R2 and R3 would increase the stability of the virtual ground, or indeed; increasing the values of R1 and R4 proportionally.
Considering the resitor values in a typical feedback opamp circuit should be in the realm of at least 10's of kilo-ohms anyway, this is a viable alternative.