Ортогональный базис

@STV_AL · Регистрация: 04.11.2018

Студворк — интернет-сервис помощи студентам

Помогите решить вопрос и исправить выпадающую ошибку.
Напишите программу, которая находит наилучшее решение системы линейных алгебраических уравнений методом наименьших квадратов.

Формат входных данных:

В первой строке задаются два числа: количество уравнений n и количество неизвестных m. Количество уравнений не меньше количества неизвестных. Далее идут n строк, каждая из которых содержит m+1 число. Первые m чисел -- это коэффициенты i-го уравнения системы, а последнее, (m+1)-е число -- коэффициент bi, стоящий в правой части i-го уравнения.
Формат выходных данных:
В качестве результата следует вывести решение системы в виде m чисел, разделенных пробелом.
Sample Input:

4 2
4 2 8
5 2 4
2 6 2
3 0 8

Sample Output:

1.6531165311653115 -0.30894308943089427

Написал:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
import numpy.linalg as lg
import numpy as np
 
n, m = [int(i) for i in input().split(" ")]
 
A = []
b = []
i = 0
 
while i != n:
    s = [int(i) for i in input().split()]
    A.append(s[:len(s)-1]), b.append(s[-1])
    i+=1
 
A = np.matrix(A)
b = np.array(b)
 
m, c = lg.lstsq(A, b, rcond=None)[0]
 
print(m, c)

Ввожу:
4 2
4 2 8
5 2 4
2 6 2
3 0 8

Получаю:

Traceback (most recent call last):
File "jailed_code", line 19, in <module>
m, c = lg.lstsq(A, b, rcond=None)[0]
ValueError: not enough values to unpack (expected 2, got 1)

ПОМОГИТЕ НАЙТИ ОШИБКУ!!!!

@u235 · 19.01.2020, 08:53

Сообщение от STV_AL

Напишите программу, которая находит наилучшее решение системы линейных алгебраических уравнений методом наименьших квадратов.

Если не путаю, то это делается через псевдообращение матрицы.
numpy.linalg.pinv
Т.е. находите псевдообратную матрицу коэффициентов и умножаете ее на вектор b.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .