Analytical formulas for the static multipole polarizabilities of hydrogen-like ions are derived by using the analytical wave functions and the reduced Green function and by applying a numerical fitting procedure. Our results are then applied to the studies of blackbody radiation shifts to atomic energy levels at different temperatures. Our analytical results can be served as a benchmark for other theoretical methods.