Demodulat0r O F D M, Demodulat0r O F D M Опции

[Dub01]

Нашёл в сети исходник на MatLab'e по генерации, передачи и получению OFDM сигнала, но у меня задача другая. Имеется сигнал OFDM, записанный с эфира. Надо его прослушать, т.е. демодулировать. Коверкался с исходником, так ничего и не получилось.

Уважаемые эксперты! Помогите, может у кого завалялся демодулятор OFDM сигналов!? Заранее большое спасибо.

[WEST128]

А что известно про источник сигнала ?

[AVR]

Извините за оффтоп, тоже вопросик есть...
Есть запись сигнала, параметры которого полностью известны, но вот демодулировать его у меня не получается... В каждом символе есть небольшая копия конца символа (guard интервал) и сам символ. При движении "окна" DFT по всему этому символу амплитуды ессно остаются почти одинаковыми, но вот фазы крутятся...
Если удачно подогнать смещение от начала, то этот поворот будет таким что его можно скорректировать по pilot-поднесущим (фаза откладывается по вертикали от -180 до + 180):

А если не повезет, то так:

Во втором случае, наверное, ничто не поможет...
Подскажите, пожалуйста, что может быть не так, как быть...
Может частота дискретизации немного не верно выбрана что приводит к такому поведению фаз?..

[alex_os]

В каждом символе есть небольшая копия конца символа (guard интервал) и сам символ. При движении "окна" DFT по всему этому символу амплитуды ессно остаются почти одинаковыми, но вот фазы крутятся...

Не понял, окно DFT вдоль одного символа смещается у Вас? Вы один канал пытаетесь обработать таким образом?

Может частота дискретизации немного не верно выбрана что приводит к такому поведению фаз?..

если частота дискретизации неправильная то все ломается . Многие звуковухи вместо оцифровки на 8000 кГц цифруют на частоте Гц на 50 меньшей. Вообщем осторожно нужно с этим...

[AVR]

Не понял, окно DFT вдоль одного символа смещается у Вас? Вы один канал пытаетесь обработать таким образом?

Запись сигнала в IQ форме, т.е. комплексной, два канала получается. Под окном я подразумеваю то что есть OFDM символ, который состоит из собственно самого информационного символа и из его небольшой части, которая подставлена вперед. Комплексное DFT я беру ессно с точностью на время информационного символа. Начальный отсчет может быть сдвинут "влево" по времени на время длительности этой копии, т.е. это окно у меня движется в пределах одного OFDM символа...

если частота дискретизации неправильная то все ломается . Многие звуковухи вместо оцифровки на 8000 кГц цифруют на частоте Гц на 50 меньшей. Вообщем осторожно нужно с этим...

т.е. как понять "ломается"? Это значит что у меня как раз тот самый случай?

[alex_os]

Запись сигнала в IQ форме, т.е. комплексной, два канала получается. Под окном я подразумеваю то что есть OFDM символ, который состоит из собственно самого информационного символа и из его небольшой части, которая подставлена вперед. Комплексное DFT я беру ессно с точностью на время информационного символа. Начальный отсчет может быть сдвинут "влево" по времени на время длительности этой копии, т.е. это окно у меня движется в пределах одного OFDM символа...
т.е. как понять "ломается"? Это значит что у меня как раз тот самый случай?

А на картинке фаза крутится для одной несущей?
Вообщем, лучше бросить смотреть на непонятную фазу и делать так:
1) Правильно спозиционировать "окно DFT". Для этого считаете корреляцию защитного интервала и последнего куска символа длинной в защитный интервал. Т.е.
Ls - длина символа
Lg - длина защитного интервала
Lu - длина окна
s - Ваш сигнал (вектор-строка)
Ns - длина сигнала в символах

в Матлабовской нотации:

Код

for offset = 1:Ls
    r(offset) = 0;
    index1 = offset;
    index2 = offset+Lu;  
    for i = 1: Ns
        r(offset) = r(offset) + r(index1: index1+Lg-1) * r(index2: index2+Lg-1)';
        index1 = index1+Ls;  index2=index2+Ls;  
    end
end

Идея такая когда защитный интервал и его копия совпадут корреляция будет максимальна. Это широкоизвестный способ символьной синхронизации для ОFDM. Длинна выборки должны быть довольна значительна (порядка сотни символов). График abs( r ) имеет характерную треугольную форму ширина основания треугольника примерно 2*Lg. При отклонении частоты дикретизации от истинной треугольник скругляется, если сильное отклонение вообще исчезает.

2) Если треугольник получился, тогда пол-дела сделано, позиционируется начало вашего окна на offset соответсвующий максимуму(max_offset) r + Lg/2. Считаете ДПФ и пытаетесь увидеть пилоты. Кстати arg(r(max_offset)) величина пропорциональная сдвигу несущей.

Еще по гуглите "OFDM timing recovery" найдете массу ссылок.

ps: а что за сигнал у Вас, если не секрет?