Непрерывных сигналов

Дискретизация и восстановление

Как известно [1, 2] детерминированный сигнал s(t), имеющий конечное значение энергии

(1)

может быть представлен в виде весовой суммы элементарных сигналов

(2)

где – система функций, которые обладают свойством ортогональности

при (3)

t₁, t₂ – моменты времени начала и окончания сигналов.

Система чисел с_n называется обобщенным спектром сигнала s(t) в ортогональной системе функций φ_n(t).

При n=k в равенстве (3) интеграл равен квадрату нормы функции φ_n(t):

. (4)

Выбирая специальным образом постоянные коэффициенты в функциях φ_n(t), можно добиться условия нормировки, при котором при любом n. Тогда система функций φ_n(t) называется ортонормированной [1, 2]. При этом спектральные коэффициенты c_n могут быть найдены из выражения:

(5)

которое является скалярным произведением функций s(t) и φ_n(t) [1, 2].

Равенство (5) может быть доказано подстановкой в него разложения (2) с учетом условий ортогональности (3) и нормировки (4) при .

Частным случаем представления (2) является тригонометрический ряд Фурье:

(6)

где ;

;

;

Т – интервал времени, на котором существует процесс s(t), или период сигнала s(t), если он является периодическим.

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

1, cos ωt, sin ωt, cos 2ωt, sin 2ωt, … (7)

где – частота первой гармоники сигнала s(t).

Как указано выше, представление сигнала s(t) тригонометрическим рядом Фурье справедливо в случае существования сигнала s(t) на конечном отрезке времени длительностью Т (тогда его представление рядом (6) справедливо только для значений времени t, находящихся на этом отрезке) или для периодического сигнала с периодом T. Тогда ряд (6) справедлив для любых моментов времени.

Если сигнал s(t) имеет спектральную функцию (спектральную плотность) G_S(ω), отличающуюся от нуля на отрезке частот , то он может быть представлен своими отдельными значениями (отсчетами), следующими через интервал времени, Δt=1/2F_в, где F_в – верхняя частота спектра s(t).

При этом выполняется равенство:

(8)

в котором – отдельные значения (отсчеты) сигнала, играющие роль спектральных коэффициентов с_n в равенствах (2) и (5), а функции образуют систему ортогональных базисных функций.

(9)

Сформулированное выше положение о возможности представления непрерывного сигнала своими отдельными (дискретными) значениями в отечественной литературе часто называется теоремой Котельникова или теоремой отсчетов. Практическое значение ее заключается в том, что для передачи через канал связи непрерывного сигнала (сообщения) s(t) с ограниченной полосой частот достаточно передавать последовательность его дискретных значений , следующих через интервал дискретизации .

Для восстановления s(t) на приемной стороне при этом необходимо сформировать процесс

(10)

где δ(t) – дельта-функция,

и подать его на вход идеального ФНЧ с частотой среза F_в, импульсной реакцией которого является функция:

(11)

где К_ф(0) – значение коэффициента передачи фильтра на нулевой частоте.

Полагая К_ф(0)=1/2F_в, можно получить на выходе фильтра сигнал, описываемый функцией (8), если на вход подать сигнал s₁(t), описываемый выражением (10).

Однако процессы, имеющие спектральную плотность, удовлетворяющую условиям теоремы отсчетов, могут быть предсказаны на сколь угодно большой отрезок времени вперед и, следовательно, не могут нести информацию. Реальные процессы, являющиеся переносчиками информации, могут иметь, спектральную плотность, равную нулю, только в отдельных точках частотной оси. Поэтому их временная дискретизации должна сопровождаться искажением формы восстановленного процесса и, следовательно, потерями информации [2]. Можно показать, что относительная среднеквадратичная ошибка, вызываемая дискретизацией непрерывного процесса x(t), может быть найдена из выражения

(12)

где S_x(t) – спектральная плотность мощности процесса x(t);

F_д=1/∆t – частота дискретизации.

Как следует из равенства (12), величина среднеквадратичной ошибки, вызванной дискретизацией, определяется энергией части сигнала х(t), содержащейся в участке спектра, отброшенном предположением о верхнем значении его частоты F_в = 0,5F_д.

Наряду с указанной, существуют дополнительные ошибки, которые вызваны невозможностью формирования сигнала на приемной стороне в полном соответствии с формулой (8), так как на практике невозможно сформировать импульс, подобный δ-функции, и невозможно создать идеальный ФНЧ. Замена δ-функции импульсом конечной амплитуды и конечной длительности, как и замена идеального ФНЧ реальным, частотная характеристика которого не имеет нулевого коэффициента передачи на отрезке частот конечной длительности, приводит к ошибкам, которые в первом приближении могут быть оценены методами, разработанными для оценки значений комбинационных искажений, возникающих при демодуляции сигналов с АИМ фильтром нижних частот.

3. Описание лабораторной установки:

Лабораторная работа выполняется компьютерной модели, структурная схема которой приведена на рис 1. Она составляется из блоков пакета Simulink, оформленных в виде пользовательской библиотеки Student. Часть блоков этой библиотеки составляют специальные блоки, созданные для выполнения функций узлов радиотехнических систем: усилителей, дискретизаторов, безынерционных нелинейных элементов (моделей транзисторов и диодов), RC-цепей и т. п.

Рисунок 1 – Структурная схема модели «Теорема отсчетов»

Генератор импульсов Pulse Generator формирует периодическую последовательность положительных видеоимпульсов и радиоимпульсов со скважностью, равной двум, которую можно наблюдать на осциллографе ContSignal.

Параметры генератора импульсов задаются в специальном окне, которое вызывается на экран двойным нажатием правой клавиши мыши (ПКМ) на изображении блока. Амплитуда импульсов равна 1 вольту (параметр Amplitude), длительность импульсов задаётся присваиванием параметру Period (secs) одного из трёх значений 0.2 мс, 0,4 мс или 0,8 мс, скважность – параметру Pulse Width (% of Period) значения 50. Параметр Phase delay (secs) равен 0.

Аналогичным образом устанавливаются параметры несущей (генератора гармонических колебаний Sine Wave): с помощью параметра Amplitude задаётся амплитуда, равная 1 вольту, постоянная составляющая (параметр Bias) равна 0, частота несущего колебания, равная 10 кГц, задаётся как 2*pi*10e3 (параметр Frequency, rad/secs), а фаза (параметр Phase, rad) – равным pi/2 (число π записывается как pi). Параметром Phase (rad) можно менять начальную фазу несущей и пределах от 0° до 360^о, задавая этот параметр также в радианах.

Сформированный сигнал поступает на устройство дискретизации – блок Discretizator (рисунок 2). Частота дискретизации задается с помощью параметра Frequency digitations, Hz, и может принимать значения 10, 20, 40 и 80 кГц. Последовательность отсчетов может наблюдаться с помощью блока-осциллографа DiscrSignal. При настройке осциллографов следует сделать засечку в пункте Save data to workspace для сохранения результатов моделирования в рабочем пространстве системы MATLAB.

Рисунок 2 – Блок-схема и панель настройки подсистемы Discretizator

Дискретный сигнал через переключатель Switch2 поступает на вход одного из низкочастотных фильтров Filter1 (ФНЧ₁) или Filter2 (ФНЧ₂) с целью восстановления непрерывного сигнала. Оба фильтра имеют частоту среза около 15 кГц, но отличаются скоростью уменьшения коэффициента передачи в области полосы затухания. Для ФНЧ₁ порядок фильтра задаётся равным 4, для ФНЧ₂ порядок фильтра задаётся равным 12. Амплитудно-частотные характеристики ФНЧ₁ и ФНЧ₂ можно определить, подключая на их входы с помощью переключателя Switch2 гармоническое колебание от генератора Sine Wave. Если же подать на входы фильтров последовательность коротких положительных импульсов с выхода блока Discretizator, то можно наблюдать импульсные реакции.