|
|
Если я правильно понял, Вы хотели растянуть/сжать именно (в первую очередь) огибающую спектра: f0 - остается на месте, fs/4 превращается в fs/2, так?
Далее, если я правильно понимаю, то сдвиг периода основного тона меняет только структуру линечатого заполнения, а огибающая остается на месте. Мужские голоса превращаются в женские, детские и т.п. и наоборот. Темп, разборчивость - все как было.
Так неужели получилось то, что Вы спрашивали?
Мне виделся другой путь.
Растяжение спектра: Ваш первый лобовой вариант с выкидыванием нулей. Переход от новой пониженной частоты отсчетов к исходной - через редискретизацию.
Сжатие спектра:
1. Для N отсчетов сигнала на интервале T (частота отсчетов fs=N/T) находим N спектральных к-тов. (первые N/2 комплексных, от 0 до fs/2).
По Вашим данным спектр ограничен и начиная с fs/4 (N/4) все к-ты равны нулю.
2. Расширяем массив, дописывая к этому нулевому хвосту еще нулевой хвостик требуемой длины, например, M нулевых к-тов. Итак, для временного интервала Т теперь имеем не N-точечный спектр, а [2*(M+N/2)]-точечный.
3. Выполняем обратное преобразование. Получаем [N+2M] отсчетов на том же временном интервале. Т.е. новая частота отсчетов fs_new=[N+2M]/T. 4. Если требуется переход к исходной частоте отсчетов fs=N/T, то повышающе/понижающая редискретизация.
Альтернатива:
Можно попробовать растянуть/сжать ИХ синтезирующего БИХ-фильтра, получаемую при линейном предсказании (т.е.растянуть/сжать огибающую спектра). Например, опять таки с помощью редискретизации реакции этого БИХ и использования новых отсчетов ИХ (полученных для новой частоты дискретизации) так, как будто бы они являются отсчетами ИХ со старой частотой дискретизации. Получаем синтезирующий КИХ фильтр, через который пропускаем сигнал возбуждения. Можно с измененным питчем.
===================
И все-таки, что там насчет питч-шифтера?
E-mail: 
info@telesys.ru
