Python bandpass filter - singular matrix error
I'd suggest checking again your parameter choices for

b,a = fd.iirdesign(Wp,Ws,Rp,As,ftype='butter')


Looking at the resulting filter coefficients, you have pole/zero cancellation at DC, i.e. at $z=1$.