
- •5. Преобразование непрерывных сигналов в дискретные
- •5.1. Преимущества цифровой формы пердставления сигналов
- •5.2. Общая постановка задачи дискретизации
- •5.3. Способы восстановления непрерывного сигнала
- •5.4. Критерии качества восстановления
- •5.5. Методы дискретизации посредством выборок
- •5.6. Равномерная дискретизация. Теорема Котельникова
- •5.6.1. Дискретизация по частотному критерию
- •5.6.2. Теорема Котельникова
- •5.7. Теоретические и практические аспекты использования теоремы Котельникова
- •5.8. Дискретизация по критерию наибольшего отклонения
- •5.9. Дискретизация с использованием экстраполирующих многочленов Тейлора
- •5.10. Адаптивная дискретизация
- •5.11. Квантование сигналов
- •5.12. Квантование сигнала при наличии помех
- •Список литературы
- •Оглавление
5.6.2. Теорема Котельникова
Теорема устанавливает принципиальную возможность полного восстановления детерминированной функции с ограниченным спектром по ее отсчетам и указывает предельное значение интервала времени между отсчетами, при которой такое восстановление еще возможно. Она формулируется следующим образом: функция u(t),допускающая преобразование Фурье и имеющая непрерывный спектр, ограниченный полосой частот от 0 доFc=C/(2), полностью определяется дискретным рядом своих мгновенных значений, отсчитанных через интервалы времени
Физическая основа теоремы выявляется при рассмотрении связи между формой функции и шириной ее спектра. Только в случае, когда спектр функции безграничен, ее значения в сколь угодно близкие моменты времени могут изменяться произвольно (корреляционная связь между ними отсутствует). Сокращение высокочастотной части спектра до граничной частоты 1равнозначно устранению из временной функции выбросов, которые могли быть сформированы этими высокочастотными составляющими (рис5.2,а). При меньших граничных частотах2(рис.5.2,б) и3(рис.5.2,в) имеем более сглаженные функции времени. Поскольку значения этих функций в моменты времениu(t1) иu(t1+t) в пределах некоторого интервалаtне могут изменяться существенно, можно ограничиться значениями функции, взятыми через интервалыt(отсчетами).
Рис.5.2.
В результате сложных вычислений, окончательно получим
(147)
Итак, функция u(t) выражена через ее дискретные значения, взятые в моменты времениtn=nΔt=nπ/ωc.
Так как при любых целых k и п справедливы соотношения
ωc(kΔt-nΔt)=(k-n)ωcΔt=(k-n)π ,
то
(148)
Благодаря этому свойству значения функции u(k) в моменты времениtn= nΔtпредставляют собой не что иное, как ее отсчеты.
Представление функции u(t) в виде ряда (147) (ряда Котельникова) является частным случаем разложения
Роль коэффициентов Сkвыполняют отсчетыu(nΔt) функцииu(t).Базисными являются функции вида
(149)
Они называются функциями отсчетов.
Графики
этих функций прип=0 ип=1 приведены
на рис.5.3. Каждая функцияψn(t)
имеет неограниченную протяженность во
времени и достигает своего наибольшего
значения, равного единице, в момент
времениt=nπ/ωc;относительно этого момента времени она
симметрична. В моменты времениt=kπ/ωc,гдеkп,
функция обращается в нуль. Все функции
ортогональны между собой на бесконечно
большом промежутке времени, что легко
проверяется путем вычисления интеграла:
(150)
Рис.5.3. Графики функций отсчетов
Каждую функцию отсчета можно рассматривать как реакцию (отклик) идеального фильтра нижних частот с граничной частотой fcна дельта-импульс, приходящий в момент времениtn=nΔtи имеющий площадь, равнуюu(nΔt).
Теорема Котельникова распространяется на непрерывный в среднеквадратическом смысле стационарный случайный процесс с ограниченным энергетическим спектром (Sn(ω)=0 при|ω|>ωП=2πFп).
Такой процесс представляется суммой квазидетерминированных процессов, где роль ортогональных детерминированных функций выполняют функции отсчета, а случайных коэффициентов — величины выборок:
(151)
где Δt=π/ωП=1/(2FП).
Таким образом, при указанных ограничениях случайный процесс полностью определяется счетным множеством случайных величин — координат процесса.
Пример.Определить по теореме Котельникова шаг дискретизации Δtдля детерминированной функции
ориентируясь на практическую ширину спектра
с η=0,95.
По формуле
(*)
находим спектральную характеристику
откуда
Практическую ширину спектра определяем, пользуясь соотношением (*):
Поскольку
имеем
По таблице значений тангенсов получаем
ωп=13,1 1/с.
Следовательно, Δt=π/ωп=0,24 c.