Не ботай

Метод наименьших квадратов: прямая через облако точек

В лабораторной точки никогда не лягут на прямую идеально. МНК находит прямую, которая ошибается меньше всех остальных, — и это две формулы плюс аккуратная таблица.

Формулы

Для nn точек (xi; yi)(x_i;\, y_i) коэффициенты прямой y=ax+by = ax + b:

a=n∑xiyi−∑xi∑yin∑xi2−(∑xi)2,b=∑yi−a∑xina = \frac{n\sum x_i y_i - \sum x_i \sum y_i}{n\sum x_i^2 - \left(\sum x_i\right)^2}, \qquad b = \frac{\sum y_i - a\sum x_i}{n}

Выглядит громоздко, но на практике это одна таблица с четырьмя суммами: ∑x\sum x, ∑y\sum y, ∑xy\sum xy, ∑x2\sum x^2.

Полный разбор

Экспериментальные точки: (1; 2)(1;\, 2), (2; 3)(2;\, 3), (3; 5)(3;\, 5).

Шаг 1. Таблица сумм:

xix_iyiy_ixiyix_i y_ixi2x_i^2
11221
22364
335159
∑\sum6102314

Шаг 2. Наклон (n=3n = 3):

a=3⋅23−6⋅103⋅14−62=69−6042−36=96=1,5a = \frac{3 \cdot 23 - 6 \cdot 10}{3 \cdot 14 - 6^2} = \frac{69 - 60}{42 - 36} = \frac{9}{6} = 1{,}5

Шаг 3. Сдвиг:

b=10−1,5⋅63=13≈0,33b = \frac{10 - 1{,}5 \cdot 6}{3} = \frac{1}{3} \approx 0{,}33

Ответ: y=1,5x+0,33y = 1{,}5x + 0{,}33.

Проверка через точку средних: xˉ=2\bar{x} = 2, yˉ=103≈3,33\bar{y} = \tfrac{10}{3} \approx 3{,}33; подставляем: 1,5⋅2+0,33=3,331{,}5 \cdot 2 + 0{,}33 = 3{,}33 ✓. Прямая прошла через центр облака — коэффициенты верны.

Суть

Смысл метода в одной фразе: из всех прямых мира выбирается та, у которой сумма квадратов вертикальных расстояний до точек минимальна. Форма ∑(yi−axi−b)2→min⁡\sum (y_i - ax_i - b)^2 \to \min — и формулы выше получаются из равенства нулю частных производных по aa и bb.

Типичные ошибки

Частые вопросы

Почему именно квадраты отклонений?

Квадрат убивает знак: отклонения вверх и вниз не сокращаются, а честно складываются, и большие промахи штрафуются сильнее маленьких. Минимум суммы квадратов даёт единственную «лучшую» прямую.

Что означают a и b в уравнении y = ax + b?

a — наклон: на сколько в среднем растёт y при росте x на единицу. b — сдвиг: значение y при x = 0. В физической лабе наклон обычно и есть искомая величина (сопротивление, ускорение).

Как проверить, что прямая посчитана верно?

Два быстрых теста: прямая обязана проходить через точку средних (x̄; ȳ), а подставленные в неё иксы должны давать значения, близкие к экспериментальным y — без систематического перекоса в одну сторону.

А если зависимость не линейная?

МНК работает и для парабол, и для экспонент — меняется только набор формул. Классический трюк лабораторных: прологарифмировать данные, чтобы экспонента стала прямой, и применить обычный линейный МНК.

Связанные темы

Идёшь по программе с нуля? Маршрут раздела «Математическая статистика» подскажет, что учить следующим.