This article pertains to interpolation of Sobolev functions at shrinking lattices hZ d from Lp shift-invariant spaces associated with cardinal functions related to general multiquadrics, φα,c(x) := (|x| 2 + c 2 ) α . The relation between the shift-invariant spaces generated by the cardinal functions and those generated by the multiquadrics themselves is considered. Additionally, Lp error estimates in terms of the dilation h are considered for the associated cardinal interpolation scheme. This analysis expands the range of α values which were previously known to give such convergence rates (i.e. O(h k ) for functions with derivatives of order up to k in Lp, 1 < p < ∞). Additionally, the analysis here demonstrates that some known best approximation rates for multiquadric approximation are obtained by their cardinal interpolants.