Аппроксимация трехпараметрического множества неопределенности при помощи эллипсоида | Статья в журнале «Молодой ученый»

Библиографическое описание:

Архипова А. Б., Лузин А. А. Аппроксимация трехпараметрического множества неопределенности при помощи эллипсоида // Молодой ученый. — 2010. — №8. Т. 1. — С. 13-16. — URL https://moluch.ru/archive/19/1908/ (дата обращения: 21.09.2018).

Для установления закономерностей каких-либо явлений проводятся экспериментальные исследования, в ходе которых измеряют значение тех или иных физико-химических величин. При проведении любого эксперимента используются приборы различной степени точности. Поэтому результаты любого измерения всегда содержат ошибки и погрешности.

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

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

Задача принадлежит к классу аппроксимации с целью построения одной общей кривой, проходящей через все заданные точки. Существует два случая:

·         Аппроксимирующая кривая проходит через точки заданной функции в узлах

·         Кривая может не проходить через точки, заданные таблицей, но отклонение от них минимально.

В обоих случаях пытаются минимизировать отклонение. Трудность заключается в том, что значения одной из координат точек заданы в виде интервалов. Чтобы обойти подобное затруднение используют метод центра неопределенностей, где фигурой, оценивающей разброс параметров аппроксимирующих функций, является эллипсоид.

Рассмотрим задачу построения зависимости вида  с интервальными значениями y.

Пусть в результате m измерений получены экспериментальные данные, содержащие точные значения входных и интервальные значения выходных переменных:

,

,

........................

.

Ограничения входных и выходных данных:

1. ;

2. ;

3. Данные должны подчиняться определенному закону, а не быть произвольными.

Согласно принципам интервально-статистического анализа запишем систему неравенств:

,.

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

Таким образом, искомая область неопределенности заключена в пересечении данных плоскостей и представляет из себя шестигранник. Для того, чтобы получить грубое приближение, нужно найти пересечения для всех троек непараллельных плоскостей путем решения восьми систем линейных уравнений с тремя переменными (для первых трех уравнений):

Решив каждую из этих систем, получим координаты соответствующей вершины. Обозначим для удобства искомые переменные как a,b,c, а коэффициенты как kij. Тогда система запишется в виде:

 

В результате преобразований получим оценки для a, b, c:

;

Далее зная координаты точек, можно построить эллипсоид, описанный вокруг шестигранника, который будет начальной грубой оценкой. Так как шестигранник образован парами параллельных плоскостей, то он является параллелепипедом. Его центр находится на пересечении диагоналей, центр эллипсоида будет совпадать с центром параллелепипеда.

Для того чтобы натянуть на параллелепипед эллипсоид минимального объема, найдем матрицу аффинного преобразования с помощью которой данный параллелепипед получается из куба. Выберем новую систему координат, оси которой будут начинаться в центре параллелепипеда и заканчиваться на центрах трех смежных граней. Обозначим  - координаты центра параллелепипеда и эллипсоида. , где k=1..3, — координаты центров трех смежных граней. Переход от одних  координат к другим в матричном виде запишется так: , где X – старые координаты, X' – новые, А – матрица перехода. В наших обозначениях, матрица перехода имеет вид:

 

В новых координатах эллипсоид представляет из себя сферу. При подстановке в формулу сферы координат X'=A-1X, получим явную формулу эллипсоида. Таким образом, получили грубую оценку трехпараметрического множества неопределенности при помощи эллипсоида.

При наличии дополнительных экспериментальных данных, зная первоначальное приближение, необходимо провести уточнение эллипсоида неопределенности. Данные нового измерения в пространстве представляют из себя полосу, ограниченную двумя параллельными плоскостями, которые при пересечении с эллипсоидом грубой оценки дадут новое множество неопределенности. Опишем его новым эллипсоидом минимального объема.

Пусть в R3 заданы эллипсоид E и полоса R, определяемые неравенствами:

Построим эллипсоид минимального объема, содержащий тело . Множество R можно задать следующим образом: , если обозначить

В силу новых обозначений запишем уравнение множества R в эквивалентном виде:

Имеет место нормировочное соотношение: .

Пусть вещественные числа ω и u удовлетворяют условиям:

.

Комбинируя неравенства для эллипсоида и полосы получим однопараметрическое семейство эллипсоидов, содержащих тело :

Обозначим . Тогда уравнение примет вид:

Обозначив , получим выражение для однопараметрического семейства эллипсоидов:

.

Объем такого эллипсоида определяется соотношением:

.

Определитель матрицы В пропорционален определителю матрицы A и определяется по формуле:

.

Коэффициент пропорциональности зависит от u. Обозначим этот коэффициент:

Далее найдем максимум функции . Для этого введем новую функцию  , прологарифмируем ее:

.

Функции q(u) и q1(u) достигают максимума в одной точке. Исследуя на экстремум функцию q1(u), получим два корня:

Подставляя найденное значение в соответствующие выражения, получаем матрицу коэффициентов B, и объем искомого эллипсоида.

 

Литература

1.    Белов, В.М. Оценивание параметров эмпирических зависимостей методом центра неопределенности/ В.М. Белов [и др]. – Новосибирск: Наука,2001. – 176 с.

Основные термины (генерируются автоматически): эллипсоид, минимальный объем, координата, данные, центр параллелепипеда, матрица перехода, однопараметрическое семейство эллипсоидов, грубая оценка.


Похожие статьи

Сечение поверхностей 2-го порядка общего вида по эллипсу...

Основные термины (генерируются автоматически): повернутое положение, трехосный эллипсоид, заданная площадь, фронтальный очерк, эллипс, поверхность, общее положение

Оценка погрешности кубатурных формул общего вида над фактор-пространством Соболева.

Выражение объемов n-мерного симплекса и n-мерного...

Для невырожденной квадратной матрицы порядка выполняется следующее равенство.

. Подставив координаты вершин , где , в известную формулу для объема n-мерного симплекса.

. Объем этого параллелепипеда равен

Алгоритм интервального оценивания параметров нелинейных...

В частном случае, когда линейна, доверительная область есть эллипсоид.

eps1 — минимальный объём ДМ; eps2 — точность определения границ

Требуется найти интервальные МНК-оценки параметров на основании данных табл. 2.

Использование матриц комбинаторного типа для построения...

где , , Xi – матрицы комбинаторного типа, содержащие координаты i-го объекта класса K1, Yj – аналогичные матрицы для класса K2.

Обозначив эту последнюю меру близости через j(x) = r(P, P(x)), задачу минимальной матричной коррекции модели (исходной задачи) можно...

Изменение ширины координатной зоны проекции Гаусса...

Это позволяет сократить объем вычислений при преобразовании координат Гаусса — Крюгера из одной зоны в другую.

Таким образом, использование системы координат Гаусса — Крюгера не ограничивается традиционной разбивкой земного эллипсоида на 3-х, или 6-градусные зоны.

Расчет параметров при оценке характеристик комплексированной...

— сжатие эллипсоида. Углы пространственной ориентации (курс, крен, тангаж) получаем

Введем нормальную систему координат OXgYgZg(по ГОСТ 20058–80), центр которой связан с

Матрицу поворота (матрица направляющих косинусов) из связанной системы координат...

Вычисление дисперсии оценки временного положения...

Библиографическое описание: Лишак М. Ю., Рэнцэнмягмар Б. Вычисление дисперсии оценки временного положения радиосигнала на выходе

Теоретическая (минимально достижимая) дисперсия ошибки измерения задержки определяется выражением [1, стр.190, (3.3.38)].

Решение транспортных задач с использованием свойств...

Переменными транспортной задачи являются xij — объемы перевозок от i-го поставщика каждому j-му потребителю.

Таблица 2. Пример транспортной задачи (исходные данные).

Это параллелепипед ABCDA1B1C1D1.

Исследование статической задачи несимметричной теории...

В прямоугольной декартовой системе координат эти уравнения (1.1) и (1.2) перепишутся

, (1.3). , (1.4). . В качестве объема можно выбрать бесконечно малый тетраэдр с тремя

, . (2.4). Если ввести матрицу порядка и N-компонентные векторы и , то система примет матричный вид.

Обсуждение

Социальные комментарии Cackle

Похожие статьи

Сечение поверхностей 2-го порядка общего вида по эллипсу...

Основные термины (генерируются автоматически): повернутое положение, трехосный эллипсоид, заданная площадь, фронтальный очерк, эллипс, поверхность, общее положение

Оценка погрешности кубатурных формул общего вида над фактор-пространством Соболева.

Выражение объемов n-мерного симплекса и n-мерного...

Для невырожденной квадратной матрицы порядка выполняется следующее равенство.

. Подставив координаты вершин , где , в известную формулу для объема n-мерного симплекса.

. Объем этого параллелепипеда равен

Алгоритм интервального оценивания параметров нелинейных...

В частном случае, когда линейна, доверительная область есть эллипсоид.

eps1 — минимальный объём ДМ; eps2 — точность определения границ

Требуется найти интервальные МНК-оценки параметров на основании данных табл. 2.

Использование матриц комбинаторного типа для построения...

где , , Xi – матрицы комбинаторного типа, содержащие координаты i-го объекта класса K1, Yj – аналогичные матрицы для класса K2.

Обозначив эту последнюю меру близости через j(x) = r(P, P(x)), задачу минимальной матричной коррекции модели (исходной задачи) можно...

Изменение ширины координатной зоны проекции Гаусса...

Это позволяет сократить объем вычислений при преобразовании координат Гаусса — Крюгера из одной зоны в другую.

Таким образом, использование системы координат Гаусса — Крюгера не ограничивается традиционной разбивкой земного эллипсоида на 3-х, или 6-градусные зоны.

Расчет параметров при оценке характеристик комплексированной...

— сжатие эллипсоида. Углы пространственной ориентации (курс, крен, тангаж) получаем

Введем нормальную систему координат OXgYgZg(по ГОСТ 20058–80), центр которой связан с

Матрицу поворота (матрица направляющих косинусов) из связанной системы координат...

Вычисление дисперсии оценки временного положения...

Библиографическое описание: Лишак М. Ю., Рэнцэнмягмар Б. Вычисление дисперсии оценки временного положения радиосигнала на выходе

Теоретическая (минимально достижимая) дисперсия ошибки измерения задержки определяется выражением [1, стр.190, (3.3.38)].

Решение транспортных задач с использованием свойств...

Переменными транспортной задачи являются xij — объемы перевозок от i-го поставщика каждому j-му потребителю.

Таблица 2. Пример транспортной задачи (исходные данные).

Это параллелепипед ABCDA1B1C1D1.

Исследование статической задачи несимметричной теории...

В прямоугольной декартовой системе координат эти уравнения (1.1) и (1.2) перепишутся

, (1.3). , (1.4). . В качестве объема можно выбрать бесконечно малый тетраэдр с тремя

, . (2.4). Если ввести матрицу порядка и N-компонентные векторы и , то система примет матричный вид.

Задать вопрос