In this paper, we solve the problem of designing a robust optimal controller which minimizes the H$\infty$-norm of the mixed sensitivity function matrix for linear multivariable systems. For a given minimized value, ${\gamma}$>o, an algorithm of finding all stabilizing controllers, such that the H$\infty$-norm of the mixed sensitivity function matrix is less than ${\gamma}$, is developed. The proposed algorithm, which is based on the model-matching and the interpolation theory, can be used for the H$\infty$-optimization problem.