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

Programmist, отлично спектр сдвинулся. Вот wav-аудио результат FM-демода из Вашего 32-битного iq-файла, уже после ресамплинга к Fd=12 кГц. (тоже потом попробую удалить из копилки). Всё отлично - никаких хвостов!

Теперь давайте iq после ФНЧ со срезом 10 кГц (на прежней высокой Fdiscr).

С ФНЧ у меня нет, надо разобраться, каким-то образом до того сделал IQ, который сразу демодулировался без хвоста.
Вот он. Только там дискрет выставлен неверно. Такому сигналу фильтрация не нужна, по краям все чисто. Не могу вернуть то, каким образом он был получен, ладно, бог с ним. Если конвертация из Real в IQ получилась правильная, значит можно делать ФНЧ для комплексного сигнала.

Хорошо, тогда пока отложим... (до завтра/послезавтра) Пока что, "комплексным сдвигом спектра" Вы получили хороший iq-сигнал. Теперь надо реализовать одинаковую НЧ-фильтрацию этих обоих каналов (можно их по очереди фильтровать как обычный моно-сигнал, лишь бы совершенно одинаково) - это и будет iq-сигнал, как бы извлечённый из исходного, в котором в общем случае могло быть много "радиостанций".

Вот он.
Скачал, удаляйте, завтра посмотрю после работы (И всем спок. ночи. :)

Да, надо отложить, а то так сразу сложно все делать в онлайн. Думаю, что разберусь потихоньку. Примеры теперь есть, главное ничего не потерять.
Сигнал удаляю, спасибо Вам. Спокойной ночи.

Вот он - с тем дискретом (8 кГц) в нём после демода нет частот выше ~ 1 кГц. А какой в нём дискрет надо выставить - не ясно. Да, не будем торопиться, откладываем.

Sinus
На вторые сутки начинает доходить. Файл, который я Вам прислал (IQ_privet.wav), это не IQ. Все, уто угодно, но только не IQ, так как содержит в себе данных в два раза больше, чем нужно. Конечно, его можно отфильтровать, сделать децимацию и т. д. Мне кажется, что конвертер из Real в IQ нужно делать по принципу, как у нас сделан IQ2R, только в обратную сторону. Тогда все получится хорошо, за исключением небольших потерь по краям и в середине спектра.
Спешить здесь не стоит. Пока буду думать. Хранение данных в IQ формате дает огромное преимущество в плане обработки, с другой стороны требует больших затрат на получение конечного или промежуточного результата (просто послушать сигнал).
А ошибка в моей демодуляции понятна. Не понятно только почему вообще она получилась не из IQ, а из двух идентичных реальных сигналов, где данные взяты через один отсчет.

В качестве иллюстрации к умозаключениям здесь положил небольшой архив, где видно, что дала формула:
A[n] := R[n] cos(2пи f n)
B[n] := R[n] sin(2пи f n)
И что хотелось-бы получить с минимальными затратами по времени на обработку. Файл IQ в архиве получен путем инверсии и тройной конвертации IQ2R.

з.ы.
С другой стороны, если посмотреть внимательно, то потери при таком преобразовании все равно есть, имхо теряется фаза сигнала:



Данный вопрос требует более внимательного изучения, т. к. эти потери могут быть и при отображении картинки спектра.
Есть и некоторые приятные моменты. В конструкторе сигналов, который пока в зачаточном состоянии, можно нарисованный сигнал сразу сохранить в IQ, далее получить из него FM или PM, по желанию, но до этого еще далеко. :)

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



Вот и не могу решить, какое из этих двух зол наименьшее. Что-то тут не так, как должно было-бы быть.

Ну, вот и сделал я демодулятор. Версия 55.
В конце концов, пользователю совсем не обязательно знать, что там внутри этой программы происходит. Указанная выше проблема с размерностью IQ массива осталась.
На кнопках |< и >|под плеером работает правая кнопка мыши, чтобы позиционировать метки времени в начало или конец записи.

Programmist


вот запись: http://doris.kiev.ua/RS/P-303/Example2/RF_201302252003278.wav
После ЧМ должно быть такое, по крайней мере это HDSDR: http://doris.kiev.ua/RS/P-303/Example2/IF_201302252003278.wav

У версии 55 получилось подобное, но дискрет вышел вместо 96 аж 121 кГц, и некоторые не стыковки небольшие в спектре в виде отличия. Запись дать не могу так как всплыла ошибка с записью, то есть я пишу в файл а при следующей загрузке вместо нормального спектра разваленный напрочь в виде шума. Могу дать скриншоты но думаю что вы сами получите спектра. Все таки есть лишний шумовой окрас в вашем демодуляторе ЧМ, тяжело определить но он есть.
Вот ЧМ была в COOL+ http://doris.kiev.ua/RS/P-303/Example2/AF_COOL_.wav
а эта в HDSDR http://doris.kiev.ua/RS/P-303/Example2/AF_HDSDR_.wav

RadioKoteg
Дискрет на выходе зависит от изначальной ширины прямоугольника, которым выделяем сигнал. Если вышло на выходе слишком много, то его можно еще раз уменьшить кнопкой "CSR" (не "ASR"), но лучше сразу выделять точнее. А с этой записью я вообще не понял, зачем ее целиком в ЧМ декодировать, там ниже еще какие-то тетки-практикантки в SSB сидят.

Programmist Там многоканалка. После ЧМ http://doris.kiev.ua/RS/P-303/p-303.jpg

Выходит что спектра отличаются от того что разные самплрейты.
Странно почему в HDSDR получился 96 а у вас 120 так как что там что сям были сделаны ЧМ с этого файла http://doris.kiev.ua/RS/P-303/Example2/RF_201302252003278.wav

Надо бы придумать методику тестирования точности ЧМ.

Там многоканалка
Понял. :) Да, есть шумовой окрас. Не знаю, вроде все сделано по классическим формулам, скорей всего в HDSDR демодулятор другой и фильтры какие-то установлены.

придумать методику тестирования точности ЧМ
Надо. Пока это так сказать "первый блин", а дальше буду разбираться. Мне все равно нужен ЧМ-модулятор, для рисования сигналов.

Programmist Ясно, значит разные фильтра. Возможно ваш лучше. Это уже надо думать как выяснить что лучше. Может записать WFM вещалку с поднесущими стерео каналов?

RadioKoteg
Боюсь, что стерео не получится, хотя кто его знает, если будет запись, то это уже не плохо, можно будет и ее покрутить, глядишь и выйдет что.

FM Radio http://doris.kiev.ua/RS/FM/HDSDR_20130225_190508Z_18500kHz_RF.wav в количестве 4 .

Попробуйте просто WFM демодулировать, у меня не получилось, шум. Видно что-то не так сделал.

Спасибо, буду разбираться. Слышно отлично, но очень широкая полоса.
А такого понятия, как WFM пока тут нет, просто частотный детектор, на всю ширину, сколько есть.

Слышно отлично, но очень широкая полоса. Бывают и больше. Через пару лет например сделают 128 ядерный 4 гигагерц на ядро процессор, и оцифровка с тактом 1 ГГц уже сейчас реальность. А вы к тому времени поймете что надо пересчет Фурье и всяких нюансов переложить на GPU на NVidea CUDA https://developer.nvidia.com/gpu-accelerated-libraries https://developer.nvidia.com/cufft которое сделает тоже самое что INTELL только в сотни раз быстрее.

Programmistвот http://www.youtube.com/watch?v=ujnkAQ52t0Y FFT , FFT тоже, а вы вместо симуляции океана сигнал загрузите, а в цвет волны от направления так как массив их двух приемников принятое как 3D эфир..

RadioKoteg
NVidea CUDA, В следующей жизни :))
А вещалки прекрасно режутся под дискрет 30КГц, только моно, естественно. Кстати, шума там немного, похоже, что это уже их собственные студийные шумы.

упс, у меня не получается, в результате белый шум но на спектре вроде музыка. Когда делаю запись и потом читаю на экране белый шум.. У нас пробегала эта же ошибка. она была связана с битностью записи и кажеться с неправильными переменными, не помню. Был пункт в меню сохранить в нашел где это.. Во пишу в 32 float все ок с пектром но в наушниках шум белый.

Там надо после демода все, что выше линии отрезать, или сразу CSR на 18КГц.



И амплитуду надо сбросить процентов на 70-80. Что там выше 18-ти, я не знаю, пока не разбирался, может быть информация для стерео сигнала.

Programmist то ж я делаю а в ушах шум...

И амплитуду надо сбросить процентов на 70-80. ага точно.. Все вкурил. Надо понимать этот ньюанс что нет нормализации уровня. Может так и правильно.

Картинку еле загрузил. У меня ее программа просмотра кверху ногами перевернула, ну все, думаю, докрутился спектров, пора завязывать на сегодня. :))
А амплитуду всегда надо нормализовать перед сохранением, особенно в 16 бит, иначе малые уровни потеряться могут.

Что там выше 18-ти, я не знаю, пока не разбирался, может быть информация для стерео сигнала. Да там они, разностный, лефт минус райт. http://en.wikipedia.org/wiki/FM_broadcasting#Stereo_FM а еще выше RDS , вот вам и ответ для чего такие широченные демодуляции ЧМ.

RadioKoteg
Опять куча формул. Тут я без Sinusа не разберусь. Ладно, пойду отдыхать, что-то устал сегодня. Доброй ночи!

Вот сигнал моено пилот стерео и RDS: http://doris.kiev.ua/RS/FM/FULL_201302252213319_201302252140524.wav

Спектр http://doris.kiev.ua/RS/FM/StereoFM.jpg

Programmist
Файл, который я Вам прислал (IQ_privet.wav), это не IQ. Все, что угодно, но только не IQ, так как содержит в себе данных в два раза больше, чем нужно.

Не важно, как его называть; главное - что он годится для проверки формул демодуляции. Целью здесь был сдвиг спектра, а не разработка преобразователя R 2 I/Q. Спектр же сдвигают умножением сигнала на комплексную экспонету. Так вот, полученные по нашим формулам из R[n] два массива A[n] и B[n] - это ни что иное как реальная и мнимая части сигнала R[n], умноженного на комплексную экспоненту. При желании их вполне можно назвать IQ-сигналом, потому что правда заключена вот в этих золотых словах, запечатлённых в славной истории данного топика:

I/Q то они разные бывают
(с) Mesh

I - реальная часть комплексного числа, Q - мнимая, все остальные сущности лишние.
(с) petr0v

Опять куча формул. Формул не надо. Все что выше 38 кГц это левый канал все что ниже правый, правый в правый левый в левый, 19 кГц для работы синхронного SSB детектора.
ВСЕ НЕ ТАК.

http://en.wikipedia.org/wiki/Double-sideband_suppressed-carrier_transmission