A procedure to design orthogonal multiwavelets with good time-frequency resolution is introduced. Formulas to compute the time-durations and the frequency-bandwidths of scaling functions and multiwavelets are obtained. Parameter expressions for the matrix coefficients of the multifilter banks that generate symmetric/antisymmetric scaling functions and multiwavelets supported in [O,N] are presented for N=2,...,6. Orthogonal multiwavelets with optimum time-frequency resolution are constructed, and some optimal multifilter banks are provided.