реклама на сайте
подробности

 
 
> Умножение теплицевых матриц
Виктор39
сообщение Apr 10 2015, 05:46
Сообщение #1


Частый гость
**

Группа: Участник
Сообщений: 123
Регистрация: 8-02-13
Из: Минск
Пользователь №: 75 542



Требуется оптимизировать функцию ЦОС, в которой наиболее затратными являются функция умножения комплексных матриц и функция обратной матрицы.

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

{ a+0, a+1, a+2, a+3;
a-1, a+0, a+1, a+2;
a-2, a-1, a+0, a+1 }

{ a+0, a+1, a+2;
a-1, a+0, a+1;
a-2, a-1, a+0;
a-3, a-2, a-1 }

Встречал алгоритмы, которые позволяют работать с теплицевыми матрицами размером [NxN].

Существуют ли оптимизированные алгоритмы, которые позволяют перемножать теплицевые не квадратные матрицы размером [NxM], где N != M ???
Может кто-нибудь встречал, подскажите, пожалуйста...

может быть теплицевыми являются только квадратные матрицы?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
serjj
сообщение Apr 10 2015, 14:23
Сообщение #2


Знающий
****

Группа: Участник
Сообщений: 527
Регистрация: 4-06-14
Из: Санкт-Петербург
Пользователь №: 81 866



Цитата
Не совсем похоже. В OFDM матрицы циркулярные (там халява)

Век живи, век учись). Я думал, что для OFDM Rxx = Rhh = E[h*h^H], где h - оценка канала в частотной области, полученная в пилотных точках. Вектор h мы берем на каждом OFDM символе. Далее w_MMSE = (Rxx + (1/SNR)*I)^-1 * Rxx, h_MMSE = h * w_MMSE^H. Ну вот как-то так. Мне не очень прозрачно как тут переписать матрицу Rhh через матрицу данных, т.к. в самом общем (и худшем) случае h произвольно дышит от блока к блоку.

А для MIMO получаем w_MMSE = (H*conj(H) + (1/SNR)*I)^-1 * conj(H), но тут H - матрица MIMO канала для данной поднесущей, которая в общем (и лучшем) случае произвольная комплексная матрица с линейно-независимыми строками.

ЗЫ: а про какой алгоритм вы говорите? Ниразу не попадалось инфы, чтобы MMSE применялся для SC систем.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Виктор39   Умножение теплицевых матриц   Apr 10 2015, 05:46
- - andyp   Квадратичный алгоритм описан здесь: http://stackov...   Apr 10 2015, 09:19
- - serjj   Речь случайно не о методе наименьших квадратов или...   Apr 10 2015, 09:32
|- - andyp   Цитата(serjj @ Apr 10 2015, 12:32) Речь с...   Apr 10 2015, 11:06
- - Виктор39   ЦитатаКвадратичный алгоритм описан здесь: http://s...   Apr 10 2015, 10:25
- - serjj   Вот книга (11 МБ djvu, не пролазит как прикрепленн...   Apr 10 2015, 11:34
|- - andyp   Цитата(serjj @ Apr 10 2015, 14:34) Ну тог...   Apr 10 2015, 13:47
- - Xenia   Не плохо было бы с помощью SVD посчитать собственн...   Apr 10 2015, 13:09
- - serjj   ЦитатаОн же, кстати, обычно является первой обязат...   Apr 10 2015, 13:33
|- - Xenia   Цитата(serjj @ Apr 10 2015, 16:33) Xenia,...   Apr 10 2015, 14:48
- - serjj   Хы, да это же похоже на MMSE для OFDM/MIMO OFDM ...   Apr 10 2015, 13:55
|- - andyp   Цитата(serjj @ Apr 10 2015, 16:55) Хы, да...   Apr 10 2015, 14:06
|- - andyp   Цитата(serjj @ Apr 10 2015, 17:23) Век жи...   Apr 10 2015, 14:32
- - serjj   Всегда смотрел в сторону хардверной платформы . Д...   Apr 11 2015, 08:04
- - Виктор39   насколько я понял, лучше для оптимизации формулы X...   Apr 13 2015, 09:33
|- - andyp   Цитата(Виктор39 @ Apr 13 2015, 12:33) нас...   Apr 13 2015, 10:52
- - serjj   SVD операция затратная, но: 1) псевдоинверсия вычи...   Apr 13 2015, 10:32
- - serjj   andyp, мне кажется у него другая дилема: считать л...   Apr 13 2015, 10:55
|- - andyp   Цитата(serjj @ Apr 13 2015, 13:55) andyp,...   Apr 13 2015, 11:18
|- - Виктор39   Цитата(andyp @ Apr 13 2015, 14:18) Если э...   Apr 13 2015, 11:58
|- - andyp   Цитата(Виктор39 @ Apr 13 2015, 14:58) В м...   Apr 13 2015, 12:05
- - Виктор39   Цитата(serjj @ Apr 13 2015, 13:55) мне ка...   Apr 13 2015, 11:15
- - serjj   Виктор39, DDE = data directed equalizer? Который N...   Apr 13 2015, 12:25
|- - Виктор39   Цитата(serjj @ Apr 13 2015, 15:25) Виктор...   Apr 13 2015, 12:46
- - serjj   Цитата4) матрица данных у вас должна быть все таки...   Apr 14 2015, 07:19
|- - Виктор39   Цитата(serjj @ Apr 14 2015, 10:19) В доку...   Apr 14 2015, 09:18
- - serjj   Инверсия делается для решения системы линейных ура...   Apr 14 2015, 09:40
|- - andyp   Цитата(serjj @ Apr 14 2015, 12:40) Инверс...   Apr 14 2015, 10:24
- - serjj   ЦитатаМожет статься, будет интересно http://scicom...   Apr 14 2015, 11:05
- - Виктор39   спасибо за помощь   Apr 14 2015, 13:17
- - asoharev   не забываем про супербыстрый алгоритм Шура, описан...   Apr 24 2015, 08:52


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 26th June 2025 - 16:47
Рейтинг@Mail.ru


Страница сгенерированна за 0.01291 секунд с 7
ELECTRONIX ©2004-2016