Construction of orthogonal multiscaling functions and multiwavelets with higher approximation order based on the matrix extension algorithm