  We present a new method for multi-component power spectra estimation in
multi-frequency observations of the CMB. Our method is based on matching a
model to the cross and auto power spectra of observed maps. 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
simulations containing five astrophysical components and white noise. The beam
smoothing effect is taken into account.
