Дана несовместная система линейных уравнений y=Ax (y - n-мерный вектор наблюдаемых значений, A - известная матрица размерности n x n, x - неизвестный m-мерный вектор оцениваемых параметров). Решение задачи по методу наименьших квадратов дается m - мерным вектором x r = (ATA) - 1 ATy , где T - обозначает транспонирование.
2. Обратная задача обработки наблюдений
По известному решению x r вычисляем идеальную реализацию в виде n - мерного вектора y r = A x r . Определяем n - мерный вектор ошибок Δ = y - y r . Находим максимальную ошибку e = max (|Δ|1 , ..., |
Δ|n). Строим множество реализаций в виде интервального вектора [y] = [y1 , y2], где y1 = y - εI , y2 = y + εI , I = (1,...,1) T - единичный n-мерный вектор.
3. Интервальное, нечеткие и классические решения
Интервальное решение прямой задачи дается m - мерным вектором [x] = {x |x = (ATA) - 1ATy , y ´ [y]}. Нечеткие решения являются подмножествами [x], которое рассматривается как нечеткое множество. Этим методом решены задачи [1,2], при этом, если e в [1] выбирается по всем измерениям, то в [2] в виде двух чисел и по промежутку стабильности (последние 5 измерений). В классическом случае [3] система несовместных уравнений заменяется на систему y = Ax + v , где v - ошибка, принадлежащая нормальному закону распределения с заданными статистическими характеристиками, что очень трудно проверяется. Решением является случайный вектор, математическое ожидание которого для случая [1] дает тот же закон растворимости NaNO3 в виде y = 67.5 + 0.87z, величина растворимости которого для температуры z=320 будет [3] (стр. 32) доверительным интервалом при надежности 0.9 в виде [94.6 , 96.0]. По рассмотренной методике интервалом для растворимости будет [93.67 , 97.0] (при надежности 1), при этом, вычисления гораздо проще и нагляднее, чем в [3]. Применимость методов классической вероятности [3] для задачи [2] , вообще, является сомнительным, поскольку из 18 лет обработки наблюдений первые 13 лет имеют место большие отклонения от прямой регрессии.
СПИСОК ЛИТЕРАТУРЫ:
- Тарушкин В.Т., Тарушкин П.В., Тарушкина Л.Т. Интервальное решение задачи Д.И. Менделеева - А.А. Маркова - Ю.В. Линника. Электронная конференция РАЕН "Современные проблемы науки и образования", 15 - 20 ноября 2006.
- Тарушкин В.Т., Тарушкин П.В., Тарушкина Л.Т. Интервальная и нечеткая линейная регрессия для ВВП России. Электронная конференция РАЕН "Прикладные исследования и разработки по приоритетным направлениям науки и техники", 15 - 20 января 2007.
- Линник Ю.В. Метод наименьших квадратов. М.: ГИФМЛ, 1958.