Given a function f on the positive half-line R + and a sequence (finite or infinite) of points X = {x k } ω k=1 in R n , we define and study matrices S X (f ) = f (|x i − x j |) ω i,j=1 called Schoenberg's matrices. We are primarily interested in those matrices which generate bounded and invertible linear operators S X (f ) on ℓ 2 (N). We provide conditions on X and f for the latter to hold. If f is an ℓ 2 -positive definite function, such conditions are given in terms of the Schoenberg measure σ(f ). We also approach Schoenberg's matrices from the viewpoint of harmonic analysis on R n , wherein the notion of the strong X-positive definiteness plays a key role. In particular, we prove that each radial ℓ 2 -positive definite function is strongly Xpositive definite whenever X is separated. We also implement a "grammization" procedure for certain positive definite Schoenberg's matrices. This leads to Riesz-Fischer and Riesz sequences (Riesz bases in their linear span) of the form F X (f ) = {f (x − x j )} x j ∈X for certain radial functions f ∈ L 2 (R n ). Examples of Schoenberg's operators with various spectral properties are presented.Mathematics Subject Classification (2010). 42A82, 42B10, 33C10, 47B37