Уважаемые специалисты, подскажите методику или алгоритм действий, как правильно передискретизировать сигнал типа ОФДМ, чтобы не исказилась форма сигнала и сохранилась ортогональность поднесущих (частотная сетка)?.

dub01 Да была гдето тема, там говорили что по Лагранжу передискрет вполне качествен.

Лагранж не Ланграж - дело десятое и в принципе понятное. Более туманный вопрос для меня это как управлять процессом передискретизации. То что я записываю и затем обрабатываю сигнал с частотой дискретизации, например, в 8000Гц на самом деле не столь очевидно. Это может быть и 8001 и 7992 и т.д. Точно так же и сам источник, по паспорту одно, а в железе может оказаться чуть другое. Следовательно, если в теории необходимо отлагранжить с 8000 до 7200, то на практике эти цифры будут хоть чуть чуть, но другие и если это не учитывать, то будет нарушение ортогональности частотных каналов. Вопрос, по какому критерию можно контролировать процесс передискретизации. Т.е. из каких свойств ОФДМ сигнала можно вытащить сигнал ошибки для управления мю в интерполяции Лагранжем? На кройняк хотелось бы узнать, как правильно по английски данная проблема называется, чтоб порыться в англоязычных источниках.

alex65111 Ну вощем так скажу, нет проблемы 8001 или 7992 у нормальных систем, там жёсткие требования, например у станаго есле не путаю точность тактовой на передаче должна быть +- очень малые процы, поэтому и непонять наших децких проблем профи, там точно извесно что ловится на какой несущей с каким разносом и тактовой какая мода, что касается нас типо любителй так и тут ж просто настоятельно рекомендуют калибровку карт звуковых, ну а что касаетца аналистов то тут и приходится пускатся во все тяжкие, так как 99 процев просто пишут, не зная какой дискрет, не зная ничего собсно про сигналы, да и откуда знать, что к примеру RUS-90 офдм должен быть принят на 7456.8 дискрета и нижний канал к примеру точно 352.89 герц? О чём грить есле образцы на сайтах того ж станага имеют тактовую 2378? Так что это проблема. Но для конкретной записи всё это можно вычеслить и подогнать этим starche занимается, а так как вариантов тьма и непонятно какой пральный, то дело и захлохло, я сдалсо. Но думаетца мне, что вопрос решить мона, но кумекать надо сильно, имхо ессно. А критерий простой, частота дискрета кратна манипуляции и разносу каналов, но как нетрудно видеть для этого их нужно точно измерить, и есле манипуляцию ещё както можно, то разнос хз как, ошибка в доли герца уводит в делёкие дали чень резко.

Напомню, что сигнал m39 Mil. Сигнал записан на fs=8000 (сколько в действительности(очень точно) не ясно). Передискретизирован на fs=7200, потому как Fs/Nfft для частоты 7200 и ФФТ=128 точек получим частотный разнос 56.25 - что удовлетворяет стандарту 188-110B. Длительность символа по стандарту идёт 22.5мс (162отсчёта = 128 + 34), у нас получается Nfft/Fs = 17,8мс (128 символ) + 4,7мс (34/7200 циклический префикс). Пробовал синхронизироваться по сигнгалу через автокорреляционную функцию. Преамбулу, т.е. 4 немодулированные несущие (0.315с) и 3 модулированные несущие (0.18с) вырезал, оставил только чистую прередачу на 39 несущих + 1 доплеровский корректирующий пилот тон, его ещё не отфильтровал - всё для того, чтобы не было грязи при построении АКФ.

В итоге получил такую АКФ:

http://www.marik2sdox.ru/ofdm/akf/akf1.jpg

Если отмасштабировать можно увидеть и измерить пики:

http://www.marik2sdox.ru/ofdm/akf/akf23.jpg

Таким образом видим, что при нулевой сдвижке (начала АКФ) получаем пик на 2.56*10⁴, ближайший к нему скачок на 2.573*10⁴. Если прикинуть разницу между этими значениями получим 0.013*10⁴, т.е. fft=130 точек. Вопрос почему не 128?Это погрешности измерений или это связано как раз с вопросом о передискретизации? Ведь если исходный сигнал был записан с плавающей частотой дискретизации (8000 по паспорту и в описании файла, но реально не известно сколько, скорее 8000+-какая либо случайная величина), то и когда мы передискретизируем сигнал на 7200 ошибки никуда не исчезают.

Кто желает, кому интересно, выкладываю данный сигнал и программку в MatLab'e для наблюдения АКФ:
С К А Ч А Т Ь

dub01 Ваш сигналик имеет нормальные параметры, у меня получилось следущее тактовая измереная методом SergUA6 44.43 герца, и разнос 56.26 герца. Всё это оченно близко, мои опыты показуют, что точность в 0.1 герца вполне доствточна для вычесления fft блока и других вещей. Чекит по методе описаной в одной из статей, [(56.26-44.43)/44.43]*128=34.081 отсщёта на циклик префикс при длине блока 128 что практически совпадает с реальными даными. Тоесь ни каких неточностей дискрета в вашем случае нету, ну так гдето или неточность впределах разумной погрешности. Всё оствльное типо почему у вас там то, а вы ожидаете другое, незнаю. Сигналик вообще шумный неподелу, имхо принят на какойто девайс с ару глубокой, такой уровень внеполосного шума это не айс.

Всё это оченно близко, мои опыты показуют, что точность в 0.1 герца вполне доствточна для вычесления fft блока и других вещей. Чекит по методе описаной в одной из статей, [(56.26-44.43)/44.43]*128=34.081 отсщёта на циклик префикс при длине блока 128 что практически совпадает с реальными даными. Тоесь ни каких неточностей дискрета в вашем случае нету, ну так гдето или неточность впределах разумной погрешности.

Так вот и не понятна причина даже столь малого расхождения! Кстате вот АКФ чистого сигнала без преамбулы и пилот тона(т.е. отфильтрованного):

http://www.marik2sdox.ru/ofdm/akf/akf-result.jpg

Вот сам - С И Г Н А Л

dub01 Чёто я несовсем пойму, вы хотите на файле с дискретом 8000, получить значение 128 точек fft? Скакой стати есле это не родной дискрет? Ессно вы тут ничего похожего не получите, так как на этой частоте сигнал стакими параметрами не слепить. Имхо, есле сигнал неизвесный, нужно получить как можно точно его скорость и разнос, вычеслить родной дискрет, и размер сетки, передискретить, есле нужно подвинуть и демодулить, всё. То что вы делаете мне непонятно.

Чёто я несовсем пойму, вы хотите на файле с дискретом 8000, получить значение 128 точек fft? Скакой стати есле это не родной дискрет?
Всё в порядке там. Сигнал result.wav действительно 8000Гц , но программка (на Matlab'e), которую я давал выше передискретизирует его на 7200Гц и для него строится АКФ. Т.е. не для 8000, а для 7200, так что тут всё чётко.

dub01 Ну так нечесно, мне откуда знать это? Все файлы по 8000. и так удивился как вы получили 130 точек, ну хз думаю, можт матлаб чудеса творит :( Ну а так проверить матлаб немешает чего и как он передискретит, для спокойствия, а то так долго можно правду искать.

dub01
Чето сдается мне что ты халтуришь. В матлабе посмотри внимательно, чему на самом деле равен Х, измеренный маркером (не то что он отображает - округленно, а на самом деле - более точно).

В матлабе посмотри внимательно, чему на самом деле равен Х, измеренный маркером (не то что он отображает - округленно, а на самом деле - более точно).
Я как раз маркером мерял. Но по Вашему совету увеличил масштаб на каждом из пиков и измерил значения с точностью до 4го знака после запятой. Действительно разница получилась 0.0128*10⁴, т.е. те самые 128 точек)) Спасибо за намёк.

Решил наследить и в Вашем уголке.
Сию разместил очередной свой труд
http://www.radioscanner.ru/info/article286/
Называется "Еще раз об определении
параметров OFDM. "нелегкая это работа"".
Заходите, читайте, комментируйте.
Может пригодится.