This paper presents a simulation method to
compute noise in autonomous and forced circuits
using the Harmonic Balance (HB) formulation.
Thanks to iterative linear solvers and multidimension
conversion matrices, it can handle
large multi-tone circuits to calculate phase noise
spectra near and far from the carrier frequency.
Its efficiency is illustrated with two real circuits.