Time domain model order reduction using general orthogonal polynomials for K-power bilinear systems.