Generally when the orthogonal functions are used in system analysis, the time consuming processes of high order matrix inversion for finding the Kronecker products and the truncation errors are occurred. In this paper, a method for the system analysis of bilinear systems via fast walsh transform is devised. This method requires neither the inversion of large matrices nor the cumbersome procedures for finding Kronecker products. Thus, both the computing time and required memory size can be significantly reduced.