We present the results of nine years of the blazar observing programme at the RATAN-600 radio telescope (2005-2014). The data were obtained at six frequency bands (1.1, 2.3, 4.8, 7.7, 11.2, 21.7 GHz) for 290 blazars, mostly BL Lacs. In addition, we used data at 37 GHz obtained quasi-simultaneously with the Metsahovi radio observatory for some sources. The sample includes blazars of three types: high-synchrotron peaked (HSP), low-synchrotron peaked (LSP), and intermediate-synchrotron peaked (ISP). We present several epochs of flux density measurements, simultaneous radio spectra, spectral indices and properties of their variability. The analysis of the radio properties of different classes of blazars showed that LSP and HSP BL Lac blazars are quite different objects on average. LSPs have higher flux densities, flatter spectra and their variability increases as higher frequencies are considered. On the other hand, HSPs are very faint in radio domain, tend to have steep low frequency spectra, and they are less variable than LSPs at all frequencies. Another result is spectral flattening above 7.7 GHz detected in HSPs, while an average LSP spectrum typically remains flat at both the low and high frequency ranges we considered.

We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for six very low synchrotron peaked (VLSP) blazars. These objects were selected as VLSP candidates (with the $\nu_{peak}^{s} \leq 10^{13}$ Hz) from the archival data. We have build spectral energy distribution using quasi-simultaneous observations at the Zeiss-1000 and RATAN-600 telescopes of the Special astrophysical observatory of RAS and made an estimation of the $\nu_{peak}^{s}$. We confirmed classification as a VLSP for two sources ([HB89] 1308+326 and 3C 345), for four other blazars we have calculated $\nu_{peak}^{s}>10^{13}$ Hz.

We study the synchrotron component of the spectral energy distribution (SED) on the sample of 877 blazars using the ASDC SED Builder Tool with available broadband data from the literature. Our sample includes 423 flat-spectrum radio sources (FSRQs), 361 BL Lac objects and candidates, and 93 blazars of uncertain type. We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for the 875 objects and further classified them as high, intermediate and low synchrotron peaked sources (HSPs/ISPs/LSPs). There are 42 HSPs with $\nu_{peak}^{s} > 10^{16.5}$ Hz, 222 ISPs with $10^{14.5} < \nu_{peak}^{s} < 10^{16.5}$ Hz, and 611 LSPs with $\nu_{peak}^{s} < 10^{14.5}$ Hz in our sample. We have calculated an average value of $\nu_{peak}^{s}$ to be $10^{13.4 \pm 1.0}$ Hz for FSRQs and $10^{14.6 \pm 1.4}$ Hz for BL Lacs. We found out that $\nu_{peak}^{s}$ and the flux density at 4.8 GHz have a different distribution (as indicated by Kolmogorov-Smirnov test at significance level 0.05) for the FSRQ and BL Lac blazars, and for the RBL and XBL types of BL Lacs. Distribution of $\nu_{peak}^{s}$ values is broader for BL Lacs, than for FSRQs. There are no ultra-high energy peaked objects (with $\nu_{peak}^{s} > 10^{19}$ Hz) in our BL Lac sample according to our estimations. The significant part of FSRQs (41%) and small part of BL Lacs (9%) in our sample could be considered as candidates to the very-low synchrotron peaked blazars (with $\nu_{peak}^{s} < 10^{13}$ Hz). Our foundations confirm results of the previous studies made on samples with significantly smaller number of objects.

We study the correlation between gamma-ray and radio band radiation for 123 blazars, using the Fermi-LAT first source catalog (1FGL) and the RATAN-600 data obtained at the same period of time (within a few months). We found an apparent positive correlation for BL Lac and flat-spectrum radio quasar (FSRQ) sources from our sample through testing the value of the Pearson product-moment correlation coefficient. The BL Lac objects show higher values of the correlation coefficient than FSRQs at all frequencies, except 21.7 GHz, and at all bands, except $10-100$ GeV, typically at high confidence level (> 99%). At higher gamma-ray energies the correlation weakens and even becomes negative for BL Lacs and FSRQs. For BL Lac blazars, the correlation of the fluxes appeared to be more sensitive to the considered gamma-ray energy band, than to the frequency, while for FSRQ sources the correlation changed notably both with the considered radio frequency and gamma-ray energy band. We used a data randomization method to quantify the significance of the computed correlation coefficients. We find that the statistical significance of the correlations we obtained between the flux densities at all frequencies and the photon flux in all gamma-ray bands below 3 GeV is high for BL Lacs (chance probability $\sim 10^{-3} - 10^{-7}$). The correlation coefficient is high and significant for the $0.1-0.3$ GeV band and low and insignificant for the $10-100$ GeV band for both types of blazars for all considered frequencies.

We present a new catalogue of the RATAN-600 multi-frequency measurements for BL Lac objects. The purpose of this catalogue is to compile the BL Lac multi-frequency data that is acquired with the RATAN-600 simultaneously at several frequencies. The BL Lac objects emit a strongly variable and polarized non-thermal radiation across the entire electromagnetic spectrum from radio to gamma-rays and represent about one percent of known AGNs. They belong to the blazar population and differ from other blazars featureless optical spectrum, which sometimes have absorption lines, or have weak and narrow emission lines. One of the most effective ways of studying the physics of BL Lacs is the use of simultaneous multi-frequency data. The multi-frequency broadband radio spectrum was obtained simultaneously with an accuracy of up to 1-2 minutes for four to six frequencies: 1.1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz. The catalogue is based on the RATAN-600 observations and on the data from: equatorial coordinate and redshift, R-band magnitude, synchrotron peak frequency, SED classes and object type, literature. The present version of the catalogue contains RATAN-600 flux densities measurements over nine years (2006-2014), radio spectra at different epochs, and their parameters for more than 300 BL Lacs objects and candidates. The BL Lacs list is constantly updated with new observational data of RATAN-600.