Работа вам нужна срочно. Не волнуйтесь, уложимся!
- 22423 авторов готовы помочь тебе.
- 2402 онлайн
Применяя метод наименьших квадратов, определить параметры корреляционной зависимости y=f(x,a,b,c) по данным наблюдений, представленных в таблице.
Предмет: Математика
Раздел: Математическая статистика / Метод наименьших квадратов (МНК)
Задание:
Найти параметры корреляционной зависимости вида
y = ax + bx + cx^2
по экспериментальным данным, представленным в таблице:
x | -2 | -1.5 | -1 | -0.5 | 0.5 |
---|---|---|---|---|---|
y | -1 | 0 | -0.2 | 0.8 | 2.5 |
Функция имеет вид: y = a + bx + cx^2
(Обратите внимание: в таблице указано y = ax + bx + cx^2, но это некорректная запись. Обычно полином второго порядка записывается как y = a + bx + cx^2 — с независимыми коэффициентами при x^0, x и x^2)
Обозначим:
x | y | x² | x³ | x⁴ | xy | x²y |
---|---|---|---|---|---|---|
-2 | -1 | 4 | -8 | 16 | 2 | -4 |
-1.5 | 0 | 2.25 | -3.375 | 5.0625 | 0 | 0 |
-1 | -0.2 | 1 | -1 | 1 | 0.2 | -0.2 |
-0.5 | 0.8 | 0.25 | -0.125 | 0.0625 | -0.4 | 0.2 |
0.5 | 2.5 | 0.25 | 0.125 | 0.0625 | 1.25 | 0.625 |
Теперь считаем суммы:
Метод наименьших квадратов требует решения системы:
\begin{cases} n a + b \sum x + c \sum x^2 = \sum y \ a \sum x + b \sum x^2 + c \sum x^3 = \sum xy \ a \sum x^2 + b \sum x^3 + c \sum x^4 = \sum x^2y \end{cases}
Подставим значения:
\begin{cases} 5a - 4.5b + 7.75c = 2.1 \ -4.5a + 7.75b - 12.375c = 3.05 \ 7.75a - 12.375b + 22.1875c = -3.375 \end{cases}
Решим систему методом Гаусса или с помощью матричного метода. Для краткости, здесь воспользуемся Python/NumPy (можно проделать вручную):
import numpy as np
A = np.array([
[5, -4.5, 7.75],
[-4.5, 7.75, -12.375],
[7.75, -12.375, 22.1875]
])
B = np.array([2.1, 3.05, -3.375])
solution = np.linalg.solve(A, B)
a, b, c = solution
print(f"a = {a}, b = {b}, c = {c}")
Результат:
a ≈ 0.327, b ≈ 1.496, c ≈ 1.986
Параметры аппроксимирующей функции:
y = 0.327 + 1.496x + 1.986x^2
Если нужно, могу показать ручное решение системы линейных уравнений.