AM для теста http://doris.kiev.ua/RS/FREQ/IQ/783k_000.wav
вот результат: http://doris.kiev.ua/RS/FREQ/IQ/783k_000(AM).wav

mikasa76 АМ демодинг конешн не в дугу, нуна переделывать. Вот резалты файла RadioKoteg после SA http://www.radioscanner.ru/uploader/2011/am_k.rar у вас нелинейка явная какая-то, или RadioKoteg так отдемил. :-) Вобще чекить по вещалкам эт фигня. Эт так на грубые косяки, шлифинг делать нужно на спешиал файликах, с фул контролом параметров. Ну да ладно. В архиве кусок файла, весь не лезет в копилку.

Вот в левом после SA АМ демодинг в правом после dspexe:

http://doris.kiev.ua/RS/FREQ/IQ/783k_000(AM)_MIX_AM_K.wav

Mesh RadioKoteg спасибо за резалты, АМ надо доработать конечно, согласен. беру паузу на устранение всех недостатков и постараюсь учесть все пожелания )))))

mikasa76
Да тут спешить некуда, еще мелкий баг нашел. При конвертации Real -> I/Q размер файла на выходе получается меньше, имхо где-то в конце последнего блока нулей надо добавить.
Потом консоль можно сделать, чтобы сразу со всеми параметрами, через запятую. Тогда любой сможет ее из своей программы запускать, чтобы не париться с вызовом функций, ну это так, для ленивых :)

Programmist меньше потому как задержка при свертке, а в конце нулей я не добавляю

задержка при свертке
Странно, должно быть 1:1, я нули в начале и в конце добавляю, но они только для "разгона" функций обработки и "центровки" картинки, хотя тут реализации разные.
У меня все циклы начинаются с отрицательного числа, на размер блока FFT.

Programmist
Странно, должно быть 1:1, я нули в начале и в конце добавляю, но они только для "разгона" функций обработки и "центровки" картинки, хотя тут реализации разные. так я эту задержку отрезаю, потому как это все равно ненужный мусор, а в конце ничего не добавляю, вот поэтому и получается, что реализация будет короче на длину этой задержки. длина задержки определяется параметром, задающем качество обработки, который в свою очередь определяется длиной формирующего КИХ фильтра. лучше даже наверное вначале добавлять.

mikasa76
Да понятно, оно там не всегда ровно выходит, это так, детали, выровнять можно, потом по краям фэйдом пройти.

з.ы.
24 бита представлял как 3 байта из int32, по этой формуле:
i and $FF000000 shr 24 // не используется (символ $, здесь hex)
i and $FF0000 shr 16
i and $FF00 shr 8
i and $FF
потом берется три последних байта (начиная с последнего), не экономно, но ничего лучше не придумал, нет такого типа данных в Delphi.
Обратно так:
int32:= Byte0 + Byte1 * (1 shl 8) + Byte2 * (1 shl 16);
if int32 >= (1 shl 23) then int32:= int32 -(1 shl 24);

mikasa76, поясните, плз, как выбирать параметр Frequency в режиме узкополосного ресемплинга: представлять себе шкалу частот real-сигнала, т.е. диапазон частот от 0 до Fdiscr (где Fdiscr - дискрет в исходном IQ-файле)? Или же представлять себе диапазон от -Fdiscr/2 до +Fdiscr/2 в исходном IQ-файле?

Заранее спасибо.

Sinus если на входе вещественный сигнал, то частота будет определяться из диапазона 0 ... Fdiscr, если же сигнал комплексный, то тогда -Fdiscr/2 до +Fdiscr/2.

Извините за занудство, но чё-то мну тёти всё покоя не дают ))

mikasa76, конверсия iq --> real у Вас в плагине отлично работает (канешн, строго говоря, не мне об этом судить, ибо сравнивал лишь со своим FFT-резалтом, раскритикованным ЦОСшниками-профи..))

А вот полосовую фильтрацию плагином чего-то никак не могу получить... Скорее-то всего, что чё-то не так делаю, но может быть и баг у Вас, т.к. вроде бы никто насчёт фильтрации полос чётко ещё не отписался... (если не баг, то извините, плз, меня за неловкость)

Потому, возвращаюсь к тестовому файлу с тётушками RadioKoteg-а:
http://doris.kiev.ua/RS/FREQ/IQ/_111_000.rar

В нём - два не особо громких FM-сигнала, с центральными частотами приблизительно FC=1MHz и FC=500 kHz после конверсии iq --> real (FFT-способом). Их спектр в FFT видится так:

http://s007.radikal.ru/i302/1112/7f/62c568830323.jpg

Поскольку плагином не сумел селектировать, то выхватил тёть своим FFT-способом (и около одной из них вырезал вредную несучку), и записал их обратно как iq с пониженным дискретом; на всякий случай называю свои поделки не iq, а АВ:

http://www.divshare.com/download/16474267-dd4
band_AB_in_111_000_FFT-select_FC=1MHz_resFd=31234Hz.wav

http://www.divshare.com/download/16474276-aee
band_AB_in_111_000_FFT-select_FC=500kHz_resFd=31234Hz.wav

Плагин их хорошо FM-демодулирует, проверил с параметрами Frequency = 0, SampleRate = 6000; вот резалты (кое-какие огрехи звучания вызваны, скорее всего, хреновостью FFT-фильтрации, не зря же petr0v её критиковал):

http://www.divshare.com/download/16474277-ca0
mikasa-FM-demod_band_AB-FFT-select_FC=1MHz.wav

http://www.divshare.com/download/16474281-8da
mikasa-FM-demod_band_AB-FFT-select_FC=500kHz.wav

Это для сравнения; если удастся их отфильтровать плагином, то интересно бы сравнить...

Ну лан... чего-то суета уже у нас, и на работе, и дома - времени ни на что толком не хватает. Так что, на всякий случай, - всех с наступающим Новым годом, здоровья, удачи, успехов!

Эти тетушки меня с ума сведут :) Все результаты разные. Здесь маленький архив, где на слух слышно разницу, в том, что Mesh в SA делал и после плагина DSP: http://files.mail.ru/7LCITK, там и файл с вырезанной "несучкой".
Обнулять "несучки" можно только с окном, иначе спектр "растекается" (зачем только это Гиббс придумал) и может быть еще хуже, в некоторых случаях.
До фильтра:

После FFT фильтра:

Все остальное осталось на своих местах, АЧХ и ФЧХ не пострадали.
Конечно, еще нужна доработка, чтобы полностью убрать остатки, но это просто так не получается.
По задумке, дальше в дело вступают кисти и краски, да полная палитра еще не готова. Когда знаешь, что все получится, на первый план выходит лень :)

Sinus конверсия iq --> real у Вас в плагине отлично работает это хорошо )

А вот полосовую фильтрацию плагином чего-то никак не могу получить... Скорее-то всего, что чё-то не так делаю, но может быть и баг у Вас, т.к. вроде бы никто насчёт фильтрации полос чётко ещё не отписался... (если не баг, то извините, плз, меня за неловкость) Потому, возвращаюсь к тестовому файлу с тётушками RadioKoteg-а:
http://doris.kiev.ua/RS/FREQ/IQ/_111_000.rar

По первому сигналу, с ним пришлось повозиться. И так, на первом этапе понижаем частоту дискретизации до промежуточной и переносим в ноль. Парметры настройки указываем следующие:

Mode: 2
Frequency: -500000 (т.к. сигнал комплексный)
SampleRateOut: 50000

промежуточная реализация http://files.mail.ru/X6HYAA

Затем понижаем частоту еще раз, чтобы отфильтровать паразита на частоте примерно +/-1600 Гц, т.к. на детектор надо подать сигнал без посторонних составляющих. Параметры для этого этапа:

Mode: 2
Frequency: 0
SampleRateOut: 3000

следующая промежуточная реализация http://files.mail.ru/3GK2VO

и теперь уже эту реализацию можем спокойно демодулировать ЧМ детектором. Параметры настройки следующие:

Mode: 4
Frequency: 0
SampleRateOut: 3000

вот окончательный результат андреевский спуск

Теперь второй сигнал. Опять же так как входной сигнал комплексный, то частота сигнала будет находится в районе нуля. Выделяем этот элементарный полосовой сигнал и децимируем его. Параметры для этого промежуточного режима:

Mode: 2
Frequency: 0
SampleRateOut: 50000

результат проведенной операции http://files.mail.ru/8KC752

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

Mode: 4
Frequency: 0
SampleRateOut: 22050

вот окончательный результат центральная

Примечание: все режимы обрабатывал с качеством "High"

Встроил вызов плагина в программу (кнопка DSP), но пока это только вызов. Файл должен быть открыт и находится в папке с программой. Имя файла не должно содержать пробелов. Кому интересно, программа (с плагином) здесь: http://files.mail.ru/WDS36V

Programmist работает )))))) правда вызов методов осуществляется через exe-шник, надо будет сделать напрямую, щас подготовлю следующий релиз с учетом предыдущих пожеланий.

mikasa76
Да работать-то он работает, куда денется, только тут надо как-то функции делать, чтобы, к примеру, ввод-вывод мой остался, а обработка Ваша, тогда будет удобно. Сразу не выйдет, еще не понятно, что надо оставить, а что можно выбросить, да и ошибок еще много, как ни включу, все время что-то новое находится.
Для начала, можно поддержку длинных имен файлов (путь) сделать и просто вызов консоли, без exe.
Что-то типа InitDll(FileName: AnsiString).

з.ы.
Если нормально делать, то надо блок-схему рисовать, иначе все это в голове долго не держится.

mikasa76 при сносе с этого файла http://doris.kiev.ua/RS/FREQ/IQ/001/_IF2011-11-2806-00-00.wav
Результат dspexe.exe: http://doris.kiev.ua/RS/FREQ/IQ/001/_IF2011-11-2806-00-00(20111229104101).wav
но в
Programmist можно подготовить, отрезав предварительно фильтром не нужную боковую.

Результат: http://doris.kiev.ua/RS/FREQ/IQ/001/_IF2011-11-2806-00-00(CE).wav

сори за не попадание в частоту.

Programmist только тут надо как-то функции делать, чтобы, к примеру, ввод-вывод мой остался, а обработка Ваша, тогда будет удобно. Сразу не выйдет, еще не понятно, что надо оставить, а что можно выбросить, да и ошибок еще много, как ни включу, все время что-то новое находится. согласен, надо сначала отладить как следует

RadioKoteg Вы предлагаете сделать просто НЧ фильтрацию с сохранением начального дискрета? я правда не знаю насколько это нужно? но проблемы в этом конечно же нет. )))))

RadioKoteg Вы предлагаете сделать просто НЧ фильтрацию с сохранением начального дискрета? я правда не знаю насколько это нужно? но проблемы в этом конечно же нет. ))))) Я предлагаю сделать снос спектра без заворота, ну то такое как получится.

Децимация нужна хорошая. Есть тут исходник на СИ, если его откомпилировать удастся, может быть и будет работать: http://files.mail.ru/N8N5HB
Совсем простенький такой, для детского сада :) Ну и у меня штук 20 релизов на Delphi, только все не то.
А снос спектра - отдельная тема. Без заворота тоже не выходит.
В общем, чтобы все собрать, "конец света" придется отменить :)
Хорошо сделать можно только разделив все операции, имхо.

RadioKoteg ок, понял, надо подумать как это провернуть

В общем, чтобы все собрать, "конец света" придется отменить :) я тоже за это, отменяем? Кто за?

Кто за?
Отменяем, что делать, хотя зрелище, наверно интересное, но тут тоже не скучно :)

Децимацию добавил, качественная (кнопочка "CSR"). Библиотека чужая, 11 лет чей-то работы:
http://files.mail.ru/ODEQF2
Здесь: http://files.mail.ru/XOPXPM (версия 17).

Programmist Да, это хорошо ресамплинг, но проблема так и не решена сноса спектра. Например как снести спектр шириной допустим 5 кГц на 300 Герц вниз?
Я делаю вашей программой это так: Переношу высоко вверх, низ очищаю, и сношу опять вниз на нужную частоту. Но что делать когда дискрет 10 кГц а спектр 4,5 кГц? Ресамплинг на большую частоту и перенос в вверх потом очистка и вновь в низ.. IMHO оно работает.

RadioKoteg
проблема так и не решена сноса спектра
Да есть модуль на FFT, но там нужна серьезная доработка самого модуля, получается снос за 3-4 шага. Автоматизация еще не готова. Ресемплиг тоже менять надо, не нравится он мне, у настоящего Кула, имхо чем-то лучше будет. Править эту библиотеку - с ума сойдешь, тем более, что она на С++ и CallBack у нее не работает. На больших файлах создается впечатление, что программа зависла.

mikasa76, спасибо, за объяснения; я понял, что неверно задавал частоту, и пытался одновременно полосу выделить (указывая новый дискрет) и демодуляцию произвести; а надо было поэтапно.

Programmist, RadioKoteg, ещё предстоит разные варианты опробовать; всё впереди. В маткаде методом FFT мне удаётся любой снос спектра - с шагом доли герца (шаг в герцах = 1/(длительность сигнала в секундах)). Но этот метод подвергся жёсткой критике профессионалами, поэтому не выкладываю здесь подробностей такого алгоритма - чтобы зря не тратилось на него время и внимание.

Надеюсь, в Новом году совместными усилиями подберём более удачные алгоритмы обработки файлов - произвольной длины, во временной области. Тем более, что mikasa76, насколько я понял, собирался опубликовать подробное описание методов в его плагине, которые, как мы убедились, хорошо работают, - это будет очень интересно. Надеюсь, лень никого из нас не остановит (это главная опасность, по сравнению с ней конец света - чепухня-небылица :))

Надеюсь, в Новом году совместными усилиями подберём более удачные алгоритмы обработки файлов
Думаю, что подберем, тем более, что Новый год начался удачно :) ровно в 23.55 китайское ADSL-чудо салютовало взрывом последнего электролита по питанию, оставив меня без связи. Потребовалась срочная модернизация:


Увеличить
Зато теперь знаю, что взорвалось до того - электролитов было два :)

Вот только от FFT отказываться не стоит, примером может служить последняя встроенная в программу децимация. Ну не дотягивает она по качеству до работы FFT, несмотря на то, что это лучшее, что возможно во временном домене. А задача сноса спектра сейчас на первом плане.

Programmist Ну не дотягивает она по качеству до работы FFT, несмотря на то, что это лучшее, что возможно во временном домене. Однозначна! Всё что вы говорите, чисто правда, вы ж лет поди 11 как в теме? А то и больше. Вам вера безмерная. :-)