.

Реализация и анализ ЦФ с КИХ

Язык: русский
Формат: курсова
Тип документа: Word Doc
0 418
Скачать документ

Министерство общего и профессионального образования
Российской Федерации
УРАЛЬСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАФЕДРА ТЕОРЕТИЧЕСКИХ ОСНОВ РАДИОТЕХНИКИ

РЕАЛИЗАЦИЯ И АНАЛИЗ ЦФ С КИХ
Курсовая работа

Руководитель
Коберниченко В.Г.
Студент
Литвинов А.А.
Группа Р-33072

Екатеринбург 2006 г.

1. Задание

Разработать цифровой согласованный фильтр (СФ) с конечной импульсной характеристикой и получить следующие его характеристики:
– спектр входного сигнала;
– спектральную (амплитудно-частотную) характеристику окна;
– АЧХ и ИХ фильтра;
– отклик фильтра на заданный сигнал;
– спектр выходного сигнала.
Проанализировать полученные результаты.
Параметры фильтра (Вариант №16):
Тип фильтра: согласованный с заданным сигналом фильтр;
Тип окна: Ханна;
Тип сигнала: прямоугольный радиоимпульс с несущей частотой, равной fд/4, и внутриимпульсной ЛЧМ (девиация частоты равна fд/4, база сигнала равна 30, скважность – 15).

2. Расчет фильтра

ПРИМЕЧАНИЕ: Все машинные расчеты в данном задании будут проводиться в среде MatLab v 6.5.
Так как в данном задании используется сигнал с B=20, воспользуемся формулами для ЛЧМ-сигналов с большой базой:

, где =dw – частота девиации, а = dw/ti – скорость нарастания частоты импульса.
Аналоговый сигнал имеет вид: при и 0 при .
Импульсная характеристика согласованного фильтра описывается выражением , где k – коэффициент, зависящий от физической реализации устройства (алгоритма), реализующего СФ. Для простоты анализа в дальнейшем амплитуду сигнала включим в k, и приравняем его к 1.
Далее нужно рассчитать, сколько точек необходимо для реализации согласованного фильтра. Сначала сосчитаем, сколько точек нужно для реализации радиоимпульса длиной и.: . Для заданного сигнала
Тогда fд выберем равной 120 Гц, а f0, равную fd/4 – соответственно 30 Гц. В этом случае максимальная частота импульса составит f0+df = 0.25Fd+0.25fd, т.е, ровно половину от частоты дискретизации: 60 Гц, следовательно теорема Котельникова будет выполнена и наложения спектров не наступит. Длительность аналогового импульса равна 1с, дискретного – 120 отсчетов (точек).
Дискретизированный сигнал имеет вид:

Uдискр(n) = Uаналог(n*Tд):

n = 0..Nи-1 = 0..119;

Далее построим выражение для импульсной характеристики фильтра:
Особенностью согласованного фильтра является то, что его импульсная характеристика h(t) является зеркальным отображением сигнала S(t) относительно прямой t=t0/2 (рис.1).

Рисунок 1

Это справедливо и для цифрового согласованного фильтра, поэтому:
Дискретная ИХ СФ:

n=0..Nи-1=0..119;

так как функция cos(t) – 2-периодическая. В MatLab же зеркальное отражение можно осуществить, если инвертировать массив отсчетов дискретизированного импульса, причем n нужно брать не от 0 до Nи-1, а от 1 до Nи, что обусловлено тем, что нумерация элементов в массивах в MatLаb ведется, начиная с единицы.
Полученная импульсная характеристика затем взвешивается окном Ханна:
w(n) = 0.5(1-cos(2π*(n-1)/(Nи-1))) на интервале причем данное окно необходимо сдвинуть вправо на , чтобы перекрывать весь сигнал. В MatLab это окно (уже со сдвигом) строится функцией hann(Nи).
На выходе согласованного фильтра после появления на входе сигнала, с которым он согласован, в момент окончания сигнала и должна появиться автокорреляционная функция(АКФ) этого сигнала. Аппроксимирующее выражение для нормированной АКФ ЛЧМ сигнала имеет вид:

В дальнейшем для определенности, амплитуду и дискретного и наналогового сигнала я беру равной 1.
Теперь приведу необходимые графики(для расчетов использована программа MatLAB):

1) Входной аналоговый и дискретизированный ЛЧМ сигналы S(t) и S(n):

2) Амплитудный спектр (АЧХ) входного сигнала Ws(n):

3) Вид и АЧХ окна:

4) ИХ взвешенного фильтра:

5) АЧХ взвешенного фильтра:

6) Спектр выходного сигнала (подвергшегося согласованной фильтрации). Получается перемножением спектров входного сигнала и спектра взвешенного фильтра:

7) Нормированный отклик фильтра на заданный сигнал – Kssf(tau) /tau – величина сдвига/. Отклик получен сверткой входного сигнала и ИХ взвешенного окном фильтра:

Теперь сравним полученный отклик с аппроксимацией АКФ входгого сигнала (см. следующую страницу):
8) Смещенная на Nи АКФ входного сигнала (отклик СФ без окна на входной сигнал) Kss(tau) /tau – величина сдвига/. Отклик получен сверткой входного сигнала и ИХ ещё не взвешенного окном СФ:

ПОЛУЧЕННЫЕ РЕЗУЛЬТАТЫ

В ходе выполнения настоящего домашнего задания был реализован согласованный фильтр с конечной импульсной характеристикой (КИХ – фильтр).
Коэффициенты этого КИХ-фильтра – это отсчеты его импульсной характеристики, полученной в ходе вычислений. Кроме теоретического алгоритма, приведенного мной в начале задания, разработан непосредственно машинный алгоритм реализации заданного фильтра для выполнения в среде MatLab 6.5.
Листинг программы (с комментариями автора):
i=1; % инициализация программного счетчика
Ni=120; % задание количества отсчётов импульса
h=[1:Ni]; % подготовка массива отсчетов ИХ к последующему заполнению
ti=120; % задание длины импульса в точках
Fd=120; % задание частоты дискретизации
F0=Fd/4; %задание несущей частоты
Wd=2*pi*Fd; % задание угловой частоты дискретизации
W0=2*pi*F0; %задание угловой несущей частоты
dF=Fd/4; % задание девиации частоты
dW=dF*2*pi; % задание угловой девиации частоты
mju=dF/ti; % задание скорости роста частоты (скорость частотной модуляции)
n=0:1/Fd:1-1/Fd; % подготовка дискретов времени с шагом, равным периоду дискретизации
Sn=chirp(n,0,1,dF); % заполнение массива отсчетов дискретизированного сигнала
figure(1)
bar(Sn,0.1); % вывод на экран графика входного сигнала
f=abs(fft(Sn,Ni)); % вычисление спектра входного сигнала с использованием БПФ
figure(2);
stem(f); % постоение графика спектра входного сигнала
w=hann(Ni); % вычисление Ni-точечного окна Ханна
figure(3);
stem(w); % построение графика вычисленного окна
Ww=abs(fft(w,Ni)); % вычисление амплитудного спектра вычисленного окна с помощью БПФ
figure(4);
plot(Ww); % построени найденного спектра
while i

Нашли опечатку? Выделите и нажмите CTRL+Enter

Похожие документы
Обсуждение

Оставить комментарий

avatar
  Подписаться  
Уведомление о
Заказать реферат!
UkrReferat.com. Всі права захищені. 2000-2020