49881

Cинтез и анализ цифрового фильтра

Курсовая

Коммуникация, связь, радиоэлектроника и цифровые приборы

Дискретная обработка аналогового сигнала. Математическая модель сигнала. Расчёт спектральной плотности сигнала. Дискретизация сигнала.

Русский

2014-09-21

1.09 MB

25 чел.

13

Министерство общего образования Российской Федерации.

ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СИСТЕМ УПРАВЛЕНИЯ И РАДИОЭЛЕКТРОНИКИ (ТУСУР)

Кафедра теоретических основ радиотехники (ТОР)

ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ

Пояснительная записка

к курсовой работе по дисциплине «РТЦ и С»

Студент гр. 160

_______Киселёв В.В.

_________дата               

Руководитель

Доцент кафедры ТОР

_______ Каратаева Н.А.                  

Томск 2002

Содержание:

Введение......................................................................................................................................3

  1.  Исходные данные к курсовой работе...........................................................................4

Дискретная обработка аналогового сигнала..................................................………4

Математическая модель сигнала....................................................……………………..4

Расчёт спектральной плотности сигнала.....................................................................     4

Определение с помощью порогового критерия Котельникова............................... 5

Дискретизация сигнала....................................................................……………………..5

Восстановление аналогового сигнала, используя ряд Котельникова........................…6

Расчёт спектральной плотности дискретной последовательности........................……7

Z-преобразование дискретной последовательности............................................……...7

ДПФ, расчёт комплексных коэффициентов СК . ………………………………………8  

Восстановление аналогового сигнала по ДПФ    ...................................………………8

Вычисление исходной дискретной последовательности, используя ОДПФ к СК.….9

Расчёт цифровых фильтров............................................................................................10

Расчёт передаточной функции цепи.........................................................................…..10

Расчёт импульсной характеристики цепи................................................................…...10

Синтез цифровой цепи методом билинейного Z - преобразования....................…….11

Расчёт ЧХ ЦФ канонического вида......................................................………………...12

Расчёт импульсной характеристики ЦФ.....................................................……………13

Синтез ЦФ с помощью метода инвариантности импульсной характеристики....…..13

Прохождение дискретного сигнала через цифровой фильтр………………………..16

Выводы по проделанной курсовой работе...................................................................18

Список используемой литературы.......………………………………………………...19.


Введение.

Одним из новых и перспективных направлений современной обработки радиосигналов  является цифровая фильтрация. В её основе лежит преобразование аналоговых сигналов в последовательность чисел и обработка этой последовательности в цифровом вычислительном устройстве, роль которого может играть как универсальная ЭВМ, так и специализированный цифровой процессор.

Применение в радиоэлектронике цифровой фильтрации открывает дополнительные возможности при обработке сигналов. В частности, могут быть реализованы сложные алгоритмы фильтрации, которые аналоговыми методами в ряде случаев вообще не удаётся осуществить. С другой стороны, возможен синтез в цифровой форме аналогов известных радиотехнических устройств различного функционального назначения, а именно: фильтров, преобразователей частоты, детекторов.

Настоящая курсовая работа ставит своей целью привить студентам практические навыки в области дискретной и цифровой обработке сигналов на примере решения конкретной задачи, включающей в себя элементы как синтеза, так и анализа цифрового фильтра.

Задание на курсовую работу состоит из двух частей. Первая часть включает подготовку аналогового сигнала к цифровой обработке. Вторая часть содержит синтез и анализ цифрового фильтра.

  1.  Исходные данные к курсовой работе

Примем следующие обозначения для простоты расчётов:E=1;T=1;Tc=3*T

 

Рис.1.1.-исходный аналоговый сигнал S(t)

Рис.1.2.-исходная линейная цепь.

2.Дискретная обработка сигнала.

2.1. Математическая модель сигнала.

S(t)=(E – E*t/T),  

S(t)=(Е-Е*(t-2*T)/T),  

S(t)=0,    2*T<t

  Или через функции Хэвисайда:

S(t)= (E – E*t/T) * +E*+E*(t-2*T)/T*

2.2. Расчёт спектральной плотности исходного сигнала.

Спектральную плотность можно вычислить применяя непосредственно интеграл Фурье:

S(w)= 

Для построения частотных характеристик спектральную плотность полезно нормировать относительно интервала описания сигнала ТС=3*Т

Рис.2.1.- Спектральная плотность исходного сигнала.

Рис..2.2.График аргумента комплексной спектральной плотности исходного сигнала.

2.3. Определение wВс помощью порогового критерия Котельникова.

0.1*Smax=S(wВ)

wВ=

2.4. Дискретизация  исходного сигнала.

Зная wВ определим  wД=2* wВ 

wД =

ТД==

Теперь вычислим число степеней свободы N

N=Tc/TД=24

Построим дискретизированный сигнал и произведём выборку:

 Рис.2.3. Отсчёты исходного аналогового сигнала.

Выборка:

E*{1;0.873;0.746;0.619;0.492;0.365;0.238;0.111;0.984;0.857;0.73;0.603;0.476;

0.349;0.222;0.095;0;0;0;0;0;0;0;0;}

2.5. Восстановление аналогового сигнала,используя ряд Котельникова.

 Условия теоремы Котельникова выполняются для полученного дискретного сигнала, поэтому воспользуемся рядом Котельникова, чтобы восстановить из дискретного сигнала аналоговый.

 

Рис.2.4. Аналаговый сигнал, восстановленный с помощью ряда Котельникова

 

   Как видно из рисунка, восстановленный сигнал остался апериодическим и его форма напоминает исходную, однако значения сигнала  и восстановленного точно совпадают лишьв точках отсчёта.

2.6. Расчёт спектральной плотности дискретной последовательности.

Для расчёта спектральной плотности дискретной последовательности применим прямое преобразование Фурье и выполним нормировку относительно N:

 

Рис.2.5. Спектральная плотность дискретизированного и исходного сигналов

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

    

2.7 Z-преобразование дискретной последовательности.

S(z)=

Для нашей дискретной последовательности Z-преобразование будет иметь вид:

S(z)=1+0.873*z-1+0.746*z-2+0.619*z-3+0.492*z-4+0.365*z-5+0.238*z-6+0.111*z-7+

0.984*z-8+0.857*z-9+0.73*z-10+0.603*z-11+0.476*z-12+0.349*z-13+0.222*z-14+0.095*z-15

  1.  ДПФ, Расчёт коэффициентов Ск.

Ск рассчитаем по формуле:

, где  ,

2.6. Графическое представление модулей комплексных коэффициентов Фурье.

На графиках мы видим, что спектр дискретной последовательности точно проходит через отсчёты коэффициентов СК.

Также следует заметить,что ,  номера которых симметричны относительно , образуют комплексно-сопряженные пары .

  1.  Восстановление аналогового сигнала по ДПФ.

На основании найденных коэффициентов восстановим исходный сигнал

S(t)=C0+

Рис.2.7. Аналоговый переодический сигнал, восстановленный по коэффициентам ДПФ.

Восстановленный сигнал является периодической функцией времени, он точно проходит по отсчётам выборки. Очевидно, что при большем значении N восстановление будет точнее.

2.10 Восстановление исходной дискретной последоватнльности используя ОДПФ к СК.

Восстановим коэффициенты  по формуле

,где  ,

Рис.2.8.Восстановленная дискретная последовательность.
3.Расчёт цифровых фильтров
.

3.1.Расчёт передаточной функции цепи.

Расчёт передаточной функции проведём по известной формуле:

Конечное выражение для К(р):

К(р)=, где

Учтём, что , подставим это в выражение для  К(р).Запишем выражение модуля  К(w):

 

Найдём  с помощью порогового критерия:

, получим .

Построим АЧХ передаточной функции цепи:

Рис.3.1.АЧХ аналогового фильтра-прототипа.

3.2.Расчёт импульсной характеристики.

Для нахождения функции  воспользуемся обратным преобразованием Лапласа:

 

Проведя обратное преобразование Лапласа, получим:

g(t)=

Построим АЧХ полученной импульсной характеристики:

Рис.3.2.Импульсная характеристика.

3.3. Синтез цифровой цепи методом билинейного Z-преобразования.

Рассчитаем wd:

wd=2*30*

Тд=

Проведём расчет ЦФ методом билинейного Z – преобразования:

К(р)=,

Перейдём от аналоговых параметров к цифровым

;

K(z)=

Подставляя вместо значение , получим конечное выражение:

K(z)=

Через K(z) найдём коэффициенты РЦФ:

Таблица 3.1.- Коэффициенты РЦФ

А0

А1

А2

в1

в2

0,045

0

-0,045

1,72     

-0,729

Покажем структуру цифрового фильтра:

А0

В1

 

                                                    В2                                          А2  

Работа рекурсивного фильтра описывается алгоритмом:

Yn=0.045*xn-0.045*xn-2+1.72*yn-1-0.729*yn-2

3.4.Расчёт АЧХ ЦФ канонического вида.

Чтобы от системной функции ЦФ перейти к его амплитудно-частотной характеристике, достаточно сделать замену следующего вида:

Z=ei*w*Td

Изобразим АЧХ ЦФ и для сравнения АЧХ аналогового фильтра-прототипа:

Рис.3.3. АЧХ цифрового фильтра.

Изобразим ФЧХ ЦФ:

3.5.Расчёт импульсной характеристики цифрового фильтра.

Для получения импульсной характеристики ЦФ необходимо произвести обратное Z-преобразование системной функции ЦФ:

Gn=

Gn=Res1+Res2

Посчитав данные вычеты, получим конечное выражение для импульсной характеристики:

Gn=0.093*(0.757)n-1-0.016*(0.963)n-1  , для n1

G0 соответствует а0 из схемы фильтра.

G0=0,045.

Изобразим Gn :

Рис.3.5. Отсчёты импульсной характеристики ЦФ.

  1.  Синтез ЦФ с помощью метода инвариантности импульсной характеристики.

 Для дискретизации импульсной характеристики g(t) необходимо непрерывный аргумент t заменить на дискретный n*Td. Затем пронормировать полученное выражение относительно ТД:

gn=g(n*Td)*Td=

Изобразим дискретизированную импульсную характеристику:

Рис.3.6. Дискретизированная импульсная характеристика.

Системная функция ТЦФ выглядит как :

где an= g(n*Td)*Td..

Трансверсальный ЦФ имеет конечную импульсную характеристику, поэтому для построения ТЦФ нужно взять конечное число отсчётов импульсной характеристики; если брать по пороговому критерию от максимума, то из анализа выпадают отрицательные отсчёты, из-за чего АЧХ ЦФ получается сильно искаженной по отношению к АЧХ фильтра прототипа,  в частности на нулевой частоте, где АЧХ прототипа принимает нулевое значение, АЧХ ЦФ при малом количестве взятых отсчётов достигает достаточно большого значения( при M=5-6 имеет максимум в этой точке.) Также имеются сильные биения АЧХ на всей полосе частот. Поэтому возьмём число отсчётов равное 40.Получим трансверсальный фильтр  с числом коэффициентов an равное 41(считая нулевой).Строить таблицу для всех коэффициентов не имеет смысла, также не имеет смысла рисовать 40 линий задержки на структурной схеме фильтра. Построим таблицу для первых шести коэффициентов an.

 

Таблица 3.2.- Коэффициенты ТЦФ

А0

А1

А2

А3

А4

А5

0,105

0,076

0,054

0,038

0,026

0,016

Изобразим структуру ТЦФ:

                                            

Запишем алгоритм работы ТЦФ:

Yn=0.105*Xn+0.076*Xn-1+…+A40*Xn-40

Для построения АЧХ ТЦФ сделаем в его системной функции замену: Z=ei*w*Td, получим АЧХ.Для сравнения покажем АЧХ аналогового фильтра-прототипа:

Рис.3.7.- АЧХ ТЦФ.

 

Рассчитаем рекурсивный фильтр, который основывается на бесконечной импульсной характеристике.Для этого свернём бесконечно убывающую геометрическую прогрессию( которой является БИХ).Тем самым получим системную функцию РЦФ:

Построим таблицу коэффициентов РЦФ:

Таблица 3.3.-коэффициенты РЦФ.

А0

А1

В1

В2

0,105

-0,105

1,722

-0,741

Структура РЦФ:

Работа рекурсивного БИХ-фильтра описывается алгоритмом:

Yn=0.105*xn-0.105*xn-1+1.722*yn-1-0.741*yn-2

Чтобы от системной функции перейти к АЧХ также сделаем замену следующего вида:

Z=ei*w*Td 

Также для сравнения покажем АЧХ фильтра-прототипа:

 

Рис.3.8.-АЧХ рекурсивного цифрового фильтра.

 Периодичность АЧХ ЦФ стала результатом дискретизации импульсной характеристики аналогового фильтра-прототипаи и представляет собой наложенный вариант частотной характеристики аналогового фильтра.

 Пульсации АЧХ трансверсального фильтра объясняются ограничением импульсной характеристики во времени( конечным числом отсчётов gn).

АЧХ рекурсивного фильтра, полученного методом инвариантности импульсных, более точно повторяет частотную характеристику аналогового фильтра-прототипа.Однако рекурсивный фильтр необходимо проверять на устойчивость.

3.7.Прохождение дискретного сигнала через цифровой фильтр.

Для получения отклика на выходе ЦФ воспользуемся аналогом временной свёртки, то есть дискретной свёрткой, получив при этом значения отсчётов выходного сигнала.

 Дискретная свёртка задаётся следующим соотношением:

Возьмём отсчёты импульсной характеристики  РЦ фильтра, полученного методом билинейного Z-преобразования:

Рис.3.9.-Результат прохождения заданного импульса через РЦФ.

Отклик ЦФ оказывается довольно предсказуемым:

Из полученного графика видно, что аналоговая цепь, как и получившиеся цифровые фильтры являются интегрирующими, а это значит, что они скорее всего ФНЧ. Видим, что огибающая отклика близка по характеру к отклику

аналогового фильтра-прототипа. Это связано с тем, что ТД  достаточно мало

(чем меньше  ТД тем ближе отклик ЦФ к отклику прототипа).


Выводы по проделанной курсовой работе.

В результате проделанной учебной работы были привиты практические навыки в области дискретной и цифровой обработки сигналов на примере решения конкретной задачи, включающей в себя элементы как синтеза, так  и анализа цифрового фильтра.

В отличие от аналоговых сигналов дискретные сигналы описываются последовательностями отсчётных значений в дискретном множестве точек.

Спектр дискретного сигнала состоит из бесконечного числа «копий» спектра исходного аналогового сигнала.

Восстановление исходного сигнала из дискретной последовательности отсчётов неизбежно связано с искажениями.

С точки зрения реализации алгоритма фильтрации принято различать трансверсальные и рекурсивные цифровые фильтры.

При синтезе ЦФ используют импульсную или частотную характеристику аналогового фильтра-прототипа.

Системная функция ЦФ есть Z - преобразование импульсной характеристики и является аналогом операторного коэффициента передачи АЦ, и поэтому также может полностью охарактеризоваться положением своих нулей и полюсов в плоскости комплексной переменной .

Рекурсивный фильтр устойчив в том случае, если все полюсы его системной функции лежат внутри единичного круга с центром в точке ,

а трансверсальный устойчив всегда.

 Выходная последовательность ЦФ есть результат дискретной свёртки входного сигнала и импульсной характеристики фильтра.

 Частотный коэффициент передачи ЦФ является преобразованием Фурье импульсной характеристики фильтра  и представляет собой периодическую функцию частоты с периодом, равным частоте дискретизации.


Список используемой литературы.

  1.  Баскаков С.И. Радиотехнические цепи и сигналы – 2-е изд. – М.: Высш. шк., 1988. – 448с.
  2.  Каратаева Н.А. Цифровая обработка сигналов – 1-е изд., дораб. и доп.
  •  Т.: ТУСУР, 1998. – 24с.

3.Каратаева Н.А.Дискретная обработка сигналов и цифровая фильтрация-

Методические указания по выполнению курсовой работы: Томск 2002.


C

R

R

C

Z-1

Z-1

Xn

Yn

А1

Аn

Z-1

Z-1

Z-1

А0

А1

В1

В2

Z-1

Z-1

Xn

Yn

Xn

Yn


 

А также другие работы, которые могут Вас заинтересовать

38293. Конституційне право 282 KB
  Декларація про державний суверенітет України від 16 липня 1990 р.Декларація стала основою для прийняття нової Конституції законів України проголосивши три види суверенітету: державний національний і народний. Народ України є єдиним джерелом державної влади в Україні. Від імені всього народу може виступати виключно Верховна Рада України.
38294. Кримінальне право 407 KB
  Часом вчинення злочину визнається час вчинення особою передбаченої законом про кримінальну відповідальність дії або бездіяльності. Питання про КВ дипломатичних представників іноземних держав та інших громадян які за законами України і міжнародними договорами згода на обов'язковість яких надана ВРУ не є підсудні у кримінальних справах судам України в разі вчинення ними злочину на території України вирішується дипломатичним шляхом.8 “Чинність закону про кримінальну відповідальність щодо злочинів вчинених іноземцями або особами без...
38295. Теорія аграрних відносин та аграрна політика 1.08 MB
  Однією із найбільш важливих сфер господарювання людей є відносини які виникають в сфері де головним засобом виробництва є земля. Першою особливістю є специфічний засіб виробництва земля яка характеризується як невіттворюваний елемент виробництва збільшити який практично не можливо. Вплив природних умов викликає в цій галузі особливий технологічний спосіб виробництва. Це вимагає від виробника врахування постійних непередбачуваний змін в технологічному процесі більшу ніж в інших галузях виробництва.
38296. Теорія та метод К.Маркса, А.Маршалла, Дж.М.Кейнса 647 KB
  Навчальний курс «Теорія та метод К. Маркса, А. Маршалла, Дж.М.Кейнса» є своєрідним поглибленням та продовженням політичної економії, мікро- макроекономіки, історії економічних учень та інших наук, які були вивчені студентами раніше, і мали за мету дати майбутнім фахівцям економістам-теоретикам більш чітке й детальне уявлення про історичні умови формування учень, предмет та методологію дослідження наведених вище видатних економістів
38297. Международное право 490.5 KB
  ПОНЯТИЕ МЕЖДУНАРОДНОГО ПРАВА Возникновение и основные этапы развития международного права Международное право: понятие и сущность особенности и функции Соотношение международного и внутригосударственного права: доктрины механизмы воздействия Система международного права. Институты и отрасли международного права Тема 2. ИСТОЧНИКИ МЕЖДУНАРОДНОГО ПРАВА Нормы международного права их особенности и виды.
38298. Види і форми міжнародної відповідальності 198 KB
  Проблема відповідальності в міжнародному праві безпосередньо пов’язана з функціонуванням міжнародного права та забезпеченням миру і правопорядку у світі. Таке розуміння цього міжнародно-правового інституту характерне, як для міжнародної наукової доктрини, так і для вітчизняної.
38299. МІЖНАРОДНО-ПРАВОВЕ РЕГУЛЮВАННЯ ПЕРЕВЕЗЕНЬ МОРСЬКИМИ ВИДАМИ ТРАНСПОРТУ 215 KB
  Перевезення товарів - справа, повязана з ризиком. Перевізники завжди прагнули максимально обмежити свою відповідальність. З іншого боку, вантажовідправники мають договірні зобовязання перед своїми покупцями і теж прагнуть зменшити ризик втрати або псування перевезених товарів.