Интервальное распространение

В числовой математике распространение интервалов или распространение ограничений интервалов — это проблема сжатия интервальных доменов, связанных с переменными R, без удаления любого значения, которое согласуется с набором ограничений (т. е. уравнений или неравенств). Его можно использовать для распространения неопределенностей в ситуации, когда ошибки представлены интервалами . [1] Распространение интервалов рассматривает проблему оценки как проблему удовлетворения ограничений .

Атомные подрядчики

Контрактор, связанный с уравнением, включающим переменные x 1 ,..., x n , — это оператор, который сжимает интервалы [ x 1 ],..., [ x n ] (которые должны охватывать x i ), не удаляя никаких значений переменных, которые согласуются с уравнением.

Подрядчик называется атомарным, если он не построен как композиция других подрядчиков. Основная теория, которая используется для построения атомарных подрядчиков, основана на интервальном анализе .

Пример . Рассмотрим, например, уравнение

х 1 + х 2 = х 3 , {\displaystyle x_{1}+x_{2}=x_{3},}

которая включает три переменные x 1 , x 2 и x 3 .

Связанный подрядчик указан в следующих заявлениях

[ х 3 ] := [ х 3 ] ( [ х 1 ] + [ х 2 ] ) {\displaystyle [x_{3}]:=[x_{3}]\cap ([x_{1}]+[x_{2}])}
[ х 1 ] := [ х 1 ] ( [ х 3 ] [ х 2 ] ) {\displaystyle [x_{1}]:=[x_{1}]\cap ([x_{3}]-[x_{2}])}
[ х 2 ] := [ х 2 ] ( [ х 3 ] [ х 1 ] ) {\displaystyle [x_{2}]:=[x_{2}]\cap ([x_{3}]-[x_{1}])}

Например, если

х 1 [ , 5 ] , {\displaystyle x_{1}\in [-\infty ,5],}
х 2 [ , 4 ] , {\displaystyle x_{2}\in [-\infty ,4],}
х 3 [ 6 , ] {\displaystyle x_{3}\in [6,\infty ]}

подрядчик выполняет следующие расчеты

х 3 = х 1 + х 2 х 3 [ 6 , ] ( [ , 5 ] + [ , 4 ] ) = [ 6 , ] [ , 9 ] = [ 6 , 9 ] . {\displaystyle x_{3}=x_{1}+x_{2}\Rightarrow x_{3}\in [6,\infty ]\cap ([-\infty ,5]+[-\infty ,4])=[6,\infty ]\cap [-\infty ,9]=[6,9].}
х 1 = х 3 х 2 х 1 [ , 5 ] ( [ 6 , ] [ , 4 ] ) = [ , 5 ] [ 2 , ] = [ 2 , 5 ] . {\displaystyle x_{1}=x_{3}-x_{2}\Rightarrow x_{1}\in [-\infty ,5]\cap ([6,\infty ]-[-\infty ,4])=[-\infty ,5]\cap [2,\infty ]=[2,5].}
х 2 = х 3 х 1 х 2 [ , 4 ] ( [ 6 , ] [ , 5 ] ) = [ , 4 ] [ 1 , ] = [ 1 , 4 ] . {\displaystyle x_{2}=x_{3}-x_{1}\Rightarrow x_{2}\in [-\infty ,4]\cap ([6,\infty ]-[-\infty ,5])=[-\infty ,4]\cap [1,\infty ]=[1,4].}
Рисунок 1: коробки до сокращения
Рисунок 2: коробки после сжатия

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

х 2 = грех ( х 1 ) , {\displaystyle x_{2}=\sin(x_{1}),}

представлено на рисунках 1 и 2.

Разложение

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

х + грех ( х у ) 0 , {\displaystyle x+\sin(xy)\leq 0,}

можно разложить на

а = х у {\displaystyle а=ху}
б = грех ( а ) {\displaystyle b=\sin(a)}
с = х + б . {\displaystyle c=x+b.}

Интервальные домены, которые должны быть связаны с новыми промежуточными переменными, следующие:

а [ , ] , {\displaystyle a\in [-\infty ,\infty ],}
б [ 1 , 1 ] , {\displaystyle b\in [-1,1],}
с [ , 0 ] . {\displaystyle c\in [-\infty ,0].}

Распространение

Принцип интервального распространения заключается в вызове всех доступных атомарных контракторов до тех пор, пока больше не будет наблюдаться сокращение. [2] В результате теоремы Кнастера-Тарского процедура всегда сходится к интервалам, которые охватывают все возможные значения для переменных. Формализацию интервального распространения можно осуществить благодаря алгебре контракторов . Интервальное распространение быстро сходится к результату и может решать проблемы, включающие несколько сотен переменных. [3]

Пример

Рассмотрим электронную схему, представленную на рисунке 3.

Рисунок 3: Файл:Электронная схема для иллюстрации распространения интервала

Предположим, что из различных измерений мы знаем, что

Э [ 23 В , 26 В ] {\displaystyle E\in [23В,26В]}
я [ 4 А , 8 А ] {\displaystyle I\in [4A,8A]}
У 1 [ 10 В , 11 В ] {\displaystyle U_{1}\in [10В,11В]}
У 2 [ 14 В , 17 В ] {\displaystyle U_{2}\in [14В,17В]}
П [ 124 Вт , 130 Вт ] {\displaystyle P\in [124W,130W]}
Р 1 [ 0 Ω , ] {\displaystyle R_{1}\in [0\Омега ,\infty ]}
Р 2 [ 0 Ω , ] . {\displaystyle R_{2}\in [0\Омега ,\infty ].}

Из схемы имеем следующие уравнения:

П = Э я {\displaystyle P=EI}
У 1 = Р 1 я {\displaystyle U_{1}=R_{1}I}
У 2 = Р 2 я {\displaystyle U_{2}=R_{2}I}
Э = У 1 + У 2 . {\displaystyle E=U_{1}+U_{2}.}

После выполнения интервального распространения получаем

Э [ 24 В , 26 В ] {\displaystyle E\in [24В,26В]}
я [ 4.769 А , 5.417 А ] {\displaystyle I\in [4.769A,5.417A]}
У 1 [ 10 В , 11 В ] {\displaystyle U_{1}\in [10В,11В]}
У 2 [ 14 В , 16 В ] {\displaystyle U_{2}\in [14В,16В]}
П [ 124 Вт , 130 Вт ] {\displaystyle P\in [124W,130W]}
Р 1 [ 1.846 Ω , 2.307 Ω ] {\displaystyle R_{1}\in [1,846\Омега ,2,307\Омега ]}
Р 2 [ 2.584 Ω , 3.355 Ω ] . {\displaystyle R_{2}\in [2,584\Омега ,3,355\Омега ].}

Ссылки

  1. ^ Jaulin, L.; Braems, I.; Walter, E. (2002). Интервальные методы для нелинейной идентификации и надежного управления (PDF) . В трудах 41-й конференции IEEE по принятию решений и управлению (CDC).
  2. ^ Клири, Дж. Л. (1987). Логическая арифметика . Будущие вычислительные системы.
  3. ^ Жолен, Л. (2006). Локализация подводного робота с использованием распространения интервальных ограничений (PDF) . В трудах CP 2006.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Interval_propagation&oldid=1157439535"