I assume that you can calculate ω(k) where k is a vector in reciprocal space. The magnitude of the k vector is 2π/λ where λ is the wavelength of the phonon and the direction is the direction of propagation of the phonon. For a solid of dimension Lx,Ly,Lz the allowable kx values are 0, 2π/Lx, 4π/Lx, ... Similar relation holds for ky and kz. Thus, the number of states is proportional to the volume in reciprocal space.
Therefore, to find phonon spectrum, you have to construct a constant ω surface in reciprocal space for all the values of ω, calculate the volume enclosed by the volume to get Vk(ω). Once you've done it, differentiate it with respect to ω and you have your spectrum.
To illustrate the point, let's consider anisotropic solid and low frequency phonons. The phonon frequency ω is related to the magnitude of the k vector by the relation
ω = c|k|
where c is the sound velocity. Constant ω surfaces will be spheres in reciprocal space and the volume given by
Vk(ω) = (4π/3) k3 = (4π/3c3)ω3
Now, it's easy to differentiate wrt ω to get
g(ω) = (4π/3) k3 = (4π/c3)ω2
Of course, you would need to normalize the expression, but these are details.
Two things to note:
g(ω) is proportional to ω2
and inversely proportional to c3