Хатри – Рао продукт - Khatri–Rao product

В математике Хатри – Рао продукт определяется как[1][2]

в которой ij-й блок - это мяпя × пjqj размер Кронекер продукт соответствующих блоков А и B, предполагая количество разделов строк и столбцов обоих матрицы равно. Размер продукта тогда я мяпя) × (Σj пjqj).

Например, если А и B оба 2 × 2 разделенные матрицы, например:

мы получаем:

Это подматрица матрицы Продукт Трейси – Сингха двух матриц (каждый раздел в этом примере является разделом в углу Продукт Трейси – Сингха ) и также может называться блочным произведением Кронекера.

По столбцам произведение Хатри – Рао

По столбцам Кронекер продукт двух матриц также можно назвать произведением Хатри – Рао. Этот продукт предполагает, что разделы матриц являются их столбцами. В таком случае м1 = м, п1 = п, п = q и для каждого j: пj = пj = 1. Полученный продукт представляет собой mp × п матрица, каждый столбец которой является произведением Кронекера соответствующих столбцов А и B. Используя матрицы из предыдущих примеров с разделенными столбцами:

так что:

Эта версия продукта Хатри – Рао с разбивкой на столбцы полезна в подходах линейной алгебры к аналитической обработке данных.[3] и в оптимизации решения обратных задач, имеющих дело с диагональной матрицей.[4][5]

В 1996 году продукт Хатри – Рао по столбцам был предложен для оценки Угол прихода (AOA) и задержки многолучевых сигналов[6] и четыре координаты источников сигналов[7] в цифровая антенная решетка.

Продукт для разделения лиц

Произведение граней матриц

Альтернативная концепция матричного произведения, которая использует построчное разбиение матриц с заданным количеством строк, была предложена В. Слюсарь[8] в 1996 г.[7][9][10][11][12]

Эта матричная операция получила название «произведение граней» матриц[9][11] или «транспонированное произведение Хатри – Рао». Этот тип операции основан на построчном произведении Кронекера двух матриц. Используя матрицы из предыдущих примеров с разделенными строками:

результат можно получить:[7][9][11]

Основные свойства

  1. Транспонировать (В. Слюсарь, 1996[7][9][10]):
    ,
  2. Билинейность и ассоциативность[7][9][10]:

    где А, B и C матрицы, а k это скаляр,

    ,[10]
    где это вектор,
  3. Свойство смешанного продукта (В. Слюсарь, 1997[10]):
    ,
    ,
    [13]
    ,[14]
    где обозначает Произведение Адамара,
  4. ,[10]
  5. ,[7]
  6. ,[14]
  7. [11][13],
    Так же:
    ,
  8. [10],
    ,где и находятся векторов,
  9. ,[15] ,
  10. ,[16]где и находятся векторов (это сочетание свойств 3 и 8),
    Так же:
  11. ,
    где вектор свертка и это Матрица преобразования Фурье (этот результат является развитием считать эскиз характеристики[17] ),
  12. [18],
    где является матрица является матрица вектор длины единиц , и вектор длины единиц
    или
    ,[19]где является матрица означает поэлементное умножение и вектор длины единиц .
    , где обозначает проникающий продукт для лица матриц[11].
    Так же:
    , где является матрица является матрица ,.
  13. [10],
    ,[19] где - вектор, состоящий из диагональных элементов , означает складывать столбцы матрицы друг на друга, чтобы создать вектор.
  14. [11][13].
    Так же:
    ,,где и находятся векторов

Примеры[16]

Теорема[16]

Если , где независимы, составляют матрицу с i.i.d. ряды , так что и ,
тогда с вероятностью для любого вектора если qwauntinty строк

В частности, если записи находятся может получить что соответствует Лемма Джонсона – Линденштрауса из когда маленький.

Блокирующий продукт

Транспонированный блок-разделитель граней в контексте многолинейной модели радара[13]

Согласно определению В. Слюсарь [7][11] блочное произведение двух граней разделенные матрицы с заданным количеством строк в блоках

можно записать как:

.

В транспонированный блок, расщепляющий грань продукта (или блочная колоночная версия произведения Хатри – Рао) двух разделенные матрицы с заданным количеством столбцов в блоках имеет вид:[7][11]

.

Основные свойства

  1. Транспонировать:
    [13]

Приложения

Продукт Face-split и продукт Block Face-splitting, используемые в тензор -матричная теория цифровые антенные решетки. Эти операции используются также в Искусственный интеллект и Машинное обучение системы для минимизации свертка и тензорный эскиз операции,[16] в популярном Обработка естественного языка модели и модели гиперграфов подобия,[20] Обобщенная модель линейного массива в статистика[19] и двух- и многомерные P-шлиц аппроксимация данных.[18]

Смотрите также

Примечания

  1. ^ Хатри К. Г., К. Р. Рао (1968). «Решения некоторых функциональных уравнений и их приложения для характеризации распределений вероятностей». Санкхья. 30: 167–180. Архивировано из оригинал на 2010-10-23. Получено 2008-08-21.
  2. ^ Чжан Х; Ян З; Цао К. (2002), "Неравенства, включающие произведения Катри – Рао положительных полуопределенных матриц", Электронные заметки по прикладной математике, 2: 117–124
  3. ^ См. Например Х. Д. Маседо и Дж. Оливейра. Подход линейной алгебры к OLAP. Формальные аспекты вычислений, 27 (2): 283–307, 2015.
  4. ^ Лев-Ари, Ханох (01.01.2005). «Эффективное решение линейных матричных уравнений в применении к обработке многостатических антенных решеток». Коммуникации в информации и системах. 05 (1): 123–130. Дои:10.4310 / CIS.2005.v5.n1.a5. ISSN  1526-7555.
  5. ^ Masiero, B .; Насименто, В. Х. (2017-05-01). «Возвращаясь к преобразованию массива Кронекера». Письма об обработке сигналов IEEE. 24 (5): 525–529. Bibcode:2017ISPL ... 24..525M. Дои:10.1109 / LSP.2017.2674969. ISSN  1070-9908.
  6. ^ Vanderveen, M.C., Ng, B.C., Papadias, C.B., & Paulraj, A. (нет данных). Совместная оценка угла и задержки (JADE) для сигналов в условиях многолучевого распространения. Отчет о тридцатой конференции Asilomar по сигналам, системам и компьютерам. - DOI: 10.1109 / acssc.1996.599145
  7. ^ а б c d е ж грамм час Слюсарь В. И. (27 декабря 1996 г.). «Конечные продукты в матрицах в радиолокационных приложениях» (PDF). Радиоэлектроника и системы связи.– 1998, Вып. 41; Число 3: 50–53.
  8. ^ Анна Эстеве, Ева Бой и Хосеп Фортиана (2009 г.): «Условия взаимодействия в дистанционной регрессии», Коммуникации в статистике - теория и методы, 38:19, с. 3501 [1]
  9. ^ а б c d е Слюсарь, В. И. (20.05.1997). «Аналитическая модель цифровой антенной решетки на основе матричных продуктов расщепления граней» (PDF). Proc. ICATT-97, Киев: 108–109.
  10. ^ а б c d е ж грамм час Слюсарь, В. И. (15.09.1997). «Новые операции матричного продукта для приложений радаров» (PDF). Proc. Прямые и обратные задачи теории электромагнитных и акустических волн (ДИПЭД-97), Львов.: 73–74.
  11. ^ а б c d е ж грамм час Слюсарь В. И. (13 марта 1998 г.). «Семейство граней произведений матриц и его свойства» (PDF). Кибернетика и системный анализ. C / C Кибернетика и Системный анализ. 1999 г.. 35 (3): 379–384. Дои:10.1007 / BF02733426.
  12. ^ Слюсарь В. И. (2003). «Обобщенные лицевые произведения матриц в моделях цифровых антенных решеток с неодинаковыми каналами» (PDF). Радиоэлектроника и системы связи. 46 (10): 9–17.
  13. ^ а б c d е Вадим Слюсарь. Новые матричные операции для DSP (Лекция). Апрель 1999. - DOI: 10.13140 / RG.2.2.31620.76164 / 1.
  14. ^ а б К. Радхакришна Рао. Оценка гетероскедастических вариаций в линейных моделях .// Журнал Американской статистической ассоциации, Vol. 65, No. 329 (март 1970 г.), стр. 161–172
  15. ^ Kasiviswanathan, Шива Прасад и др. «Цена частного выпуска таблиц непредвиденных обстоятельств и спектров случайных матриц с коррелированными строками». Материалы сорок второго симпозиума ACM по теории вычислений. 2010 г.
  16. ^ а б c d Томас Д. Але, Якоб Бак Тейс Кнудсен. Почти оптимальный тензорный набросок. Опубликовано 2019. Математика, информатика, ArXiv
  17. ^ Нинь, Фам; Расмус, Паг (2013). Быстрые и масштабируемые полиномиальные ядра с помощью явных карт функций. Международная конференция SIGKDD по обнаружению знаний и интеллектуальному анализу данных. Ассоциация вычислительной техники. Дои:10.1145/2487575.2487591.
  18. ^ а б Eilers, Paul H.C .; Маркс, Брайан Д. (2003).«Многомерная калибровка с температурным взаимодействием с использованием двумерной регрессии сигнала со штрафами». Хемометрия и интеллектуальные лабораторные системы. 66 (2): 159–174. Дои:10.1016 / S0169-7439 (03) 00029-7.
  19. ^ а б c Currie, I.D .; Дурбан, М .; Эйлерз, П. Х. С. (2006). «Обобщенные модели линейных массивов с приложениями к многомерному сглаживанию». Журнал Королевского статистического общества. 68 (2): 259–280. Дои:10.1111 / j.1467-9868.2006.00543.x.
  20. ^ Брайан Бишоф. Тензоры сосуществования более высокого порядка для гиперграфов через расщепление граней. Опубликовано 15 февраля 2020 г., Математика, Информатика, ArXiv

Рекомендации