We present a new method for multi-component power spectra estimation in multifrequency observations of the CMB. Our method is based on matching the cross and auto power spectra of observation maps to their expected values. All the component power spectra are estimated, as well as their mixing matrix. Noise power spectra are also estimated. The method has been applied to full-sky Planck simulati...