**Bien, viendo que la implementación quizás no es la más adecuada, te propongo la siguiente:**

    def matrix(length, combinationVector, fullMatrix):
        if len(combinationVector) <= length:
            matrixToFillWithOne = combinationVector.copy()
            matrixToFillWithZero = combinationVector.copy()
            matrixToFillWithOne.append(1)
            matrixToFillWithZero.append(0)
            matrix(length,matrixToFillWithOne,fullMatrix)
            matrix(length,matrixToFillWithZero,fullMatrix)
        
        else:
            fullMatrix.append(combinationVector)
            fullMatrix.sort()
        return fullMatrix
    
    
    if __name__ == '__main__':
        matrix(3, [],[])

Pasaremos tres elementos, los cuales serán:

> fullMatrix-> La matriz final que llenaremos con todas las combinaciones.
>
>combinationVector -> Una de las posibles combinaciones [0,0,0][0,1,0].....
>
>length-> Tamaño que queremos.

El funcionamiento es sencillo así que creo que con un esquema quedará más claro tras leer el código:

                                      []
                                    /    \
                                  [0]     [1]
                                  / \     / \
                              [0,0] [0,1] [1,0] [1,1]


Así se irán extendiendo en árbol todas las posibles combinaciones. Una vez la tengamos, la añadiremos a fullMatrix donde serán ordenadas, siendo el resultado final:

    [[0, 0, 0, 0], [0, 0, 0, 1], [0, 0, 1, 0], [0, 0, 1, 1], [0, 1, 0, 0], [0, 1, 0, 1], [0, 1, 1, 0], [0, 1, 1, 1], [1, 0, 0, 0], [1, 0, 0, 1], [1, 0, 1, 0], [1, 0, 1, 1], [1, 1, 0, 0], [1, 1, 0, 1], [1, 1, 1, 0], [1, 1, 1, 1]]