Weighting function for output error IIR filter design
The reason for the worse fit at low frequencies (without additional weighting) is the low density of frequency points at low frequencies due your choice of a logarithmic frequency vector. I'm pretty sure that you will get a better fit at low frequencies if you choose a (more) linear frequency grid.