Основы ЦОС: 25. Многоскоростные системы (ссылка на скачивание скриптов в описании)

[@Argentum634]

Увидел видос по матлабу - лайк не глядя!

[@rodiontimofeev392]

Здравствуйте. Поясните пожалуйста как, после уменьшения частоты дискретизации, вы определили, что получившийся сигнал равен 1600 герц ?

[@Exponenta-tech]

Здравствуйте. Во-первых, это нам показал спектр результирующего сигнала. Пики на 1000 и 1600 Гц. Во-вторых, это и должно произойти в соответствии с теорией. Зона Найквиста при понижении частоты дискретизации в шесть раз теперь ограничивается 2500 Гц. Сигнал 3400 Гц стоит теперь "по ту сторону" видимой части спектра, и на такой частоте дискретизации мы видим его отсчёты, как будто бы они являются отсчётами его спектральной копии, которая стоит "по эту сторону", и так же отстоит от границы зоны Найквиста на 3400-2500 = 900 Гц. То есть исходный сигнал располагается на частоте 2500+900 Гц, его копия - на частоте 2500-900 Гц. Отсюда и имеем 1600 Гц.

[@rodiontimofeev392]

@@Exponenta-tech Спасибо

[@ЮрийПопов-л6я]

Все красиво написано однако про тонкости автор нислова не разобрался. Есть такая удобная функция изменения частоты дискретизации resample. Вставил p,q, входной сигнал, она сама фильтр рассчитывает и выдает передискретизированный сигнал кажется супер гуд!!!! Но есть как говорится нюанс. Итак берем сигнал вещественный с полосой 4000 Гц (спектр хорошо ограничен по полосе) и частота дискретизации Fd=9600 Гц. Далее преобразовываем его в комплексный функцией hilbert. Ну и теперь бы избавится от половины отсчетов в соответствии с главным законом квантования две точки на период. Сигнал у нас комплексный поэтому можно взять по точке а период из каждой квадратуры. Призываю функцию resample понизить частоту в 2 раза и она жутко портит спектр!!!!! Вот такая особенность. Поборол проблему конечно. Но время на ее обнаружении потратил. Доверял матлабу сильно. А тут оказалось доверяй но проверяй. А как обойти проблему - предлагаю чтобы автор дополнил свой видосик и рассказал сам. А за текущий видосик - лайк однозначно!!!

[@Exponenta-tech]

Юрий, функция resample - весьма высокоуровневая, и поддерживает достаточно большой функционал. Но широкий функционал накладывает определённые ограничения. Конкретно - на вход она может принимать только действительный double, комплексные отсчёты не поддерживаются. Для задачи понижения частоты вдвое в Вашем случае необходимо использовать именно команду downsample. В описании к функции указано, что она поддерживает комплексный вход и типы данных single и double. Сравнение отсчётов действительной и мнимой частей аналитического сигнала до и после понижения частоты, а также анализ спектра сигнала показывают корректный результат операции.

[@ЮрийПопов-л6я]

Примерно так я и сделал Только смею заметить что resample прекрасно работает с комплексным сигналом. Что было многократно проверено и вверх и вниз. и очень удобен тем что сам фильтр считает внутри. т.е меняешь p q а он все сделает остальное. Только вот на такой узкой задаче упоролся. Пришлось лишнее условие добавить чтобы вовремя на downsample переходить. Спасибо за ответ

[@ЮрийПопов-л6я]

Надесь это не последнее видео про multirate. Есть книга М. Чобану. Вот там у него мультирейт многомерный по полной программе. Так что давайте матлаб в контексте многомерного мультирейта рассмотрим