This paper is concerned with branch and bound technique for the global solution of bilinear matrix inequality (BMI) problem. In this paper, more strict estimation method of lower bound of objective function and new bounding operation are proposed. The effectiveness of proposed method is confirmed in computational experiments.