矩形矩阵的QR分解及其在线性最后平方中的应用

Question

我想为 “矩形矩阵的QR分解及其在线性最后平方中的应用” 我尝试了这段代码

    import numpy as np
    import scipy.linalg as linalg
    A = np.matrix('[1, 2, 3; 4, 5, 6]')
    Q,R = linalg.qr(A)
    b = np.matrix('[7; 8]')
    x = linalg.solve(A,b)
    y = dot(Q.T,b)
    xQR = linalg.solve(R,y) 

    print("\nSolution compared")
    print(x.T,'Ax=b')
    print(xQR.T,'Rx=y')

但出现错误

    ValueError: Input a needs to be a square matrix.

我该如何解决？ 谢谢