filter_c = np.fft.fftshift(filter_c)
axs[i // L, i % L].imshow(colorize(filter_c))
axs[i // L, i % L].axis("off")
axs[i // L, i % L].set_title("$j = {}$ \n $\\theta={}".format(i // L, i % L))
i = i+1
After Change
plt.suptitle("The corresponding low-pass filter, also known as scaling function."
"Color saturation and color hue respectively denote complex magnitude and complex phase", fontsize=13)
filter_c = np.abs(filter_c)
plt.imshow(filter_c)