QR-факторизация прямоугольной матрицы и применение для линейных последних квадратов - PullRequest
Новая
       39

QR-факторизация прямоугольной матрицы и применение для линейных последних квадратов

0 голосов
/ 01 декабря 2018

Я хочу написать код для «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.

Как можноЯ это решаю?Спасибо

...