To overcome the limitation of Krylov subspace techniques using standard preconditioners on speed and memory requirement in harmonic balance simulations of the steady-state response under multi-tones drive

a new sparse-subsecional matrix as the preconditioner is developed in this paper.This preconditioner

different form the ordinary block-diagonal matrix

not only is a sufficient approximation of Jacobian matrix

but also exploits the techniques of sparseness and subsecional condensation.With this preconditioner

it could be faster by 50% and ucenwry requirement is reduced by 60% than that of block-diagond matrix.Some numeric experiments illustrating it are presented.