Проектирование цифрового фильтра

   Cибирский Государственный Университет

   телекоммуникаций  и информатики

Кафедра МЭС

Курсовая  работа

на тему:

«Проектирование цифрового фильтра» 
 
 
 
 
 

Выполнила: студентка

группы ЗМ-71

Киселева Ю. А.

Проверил:

Чухров А. С. 
 
 
 
 
 
 

Новосибирск 2010 

СОДЕРЖАНИЕ

ВВЕДЕНИЕ………………………………………………………………………3

      1. Структурная схема цифрового фильтра………………………………………….4

      2. Расчет устойчивости фильтра........................................................................................5

      3. Расчет X(jkw₁) и H(jkw₁) с помощью БПФ…………………………………….7 
 

9

      Расчет  мощности собственных шумов синтезируемого фильтра …..…11

ЗАКЛЮЧЕНИЕ…………………………………………………………………………………….13

СПИСОК ЛИТЕРАТУРЫ…………………………………………………………………….14 

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

ВВЕДЕНИЕ 

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

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

     Под цифровым фильтром понимают дискретную систему, которая описывается уравнением:

 

и реализованную  программным путем на цифровой ЭВМ  или аппаратным путем в виде специализированного  цифрового вычислительного устройства.

     Сигналы на входе x(nT) и выходе y(nT) являются цифровыми, представленными в виде двоичного кода, и цифровом фильтре в соответствии с заданными алгоритмами выполняются операции пересылки, сложения, умножения кодов.

     В курсовой работе приведены расчеты  фильтра во временной и частотной  областях при помощи быстрого дискретного преобразования Фурье (БПФ) и обратного быстрого преобразования Фурье (ОБПФ), а также расчет выходного сигнала.

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

 

1.Синтезировать структурную схему  цифрового фильтра. 

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

,

     К разностному уравнению можно  перейти, зная передаточную характеристику фильтра H(Z), характеризующую цепь в частотной области. Передаточная характеристика фильтра в общем виде:

     Подставляем в общую формулу заданные коэффициенты и получаем передаточную характеристику проектируемого цифрового фильтра:

   

На основании  передаточной функции определяем выходной сигнал:

     Далее переходим к оригиналам и записываем разностное уравнение:

     По  разностному уравнению видим, что  значение выходной величины в любой  момент времени определяется не только значением входной величины, но и предыдущим значением выходной величины. Следовательно, проектируемый фильтр является рекурсивным.

     Наиболее  часто используют структурные схемы  рекурсивных фильтров прямой формы  и прямой канонической формы.

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

     Структурная схема проектируемого фильтра, приведенная  на рис.1.1, содержит два сумматора, умножители и всего три элемента задержки (минимальное число).   

                     

                                                                      0,95  

                                                                                 

                                                                                                                         

                                               0,75                                                    -0,18 
 
 
 

                                              

                                              0,39                                                      0,3 
 
 

                                               0,48                                                   0,42    
 
 
 
 
 

Рис 1. Структурная схема  рекурсивного фильтра   
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

2. Определение  устойчивости цепи.

     При анализе работы любой цифровой структуры  важное место приобретает вопрос устойчивости. Если по каким-либо причинам цепь оказывается неустойчивой, то вместо желаемого фильтра получают генератор.

     Для того чтобы фильтр был устойчивым, полюсы его передаточной функции H(Z) должны располагаться внутри единичного круга плоскости Z.  

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

Приравниваем  знаменатель к нулю и ищем корни полученного уравнения.

1+0,75Z^-1+0,39Z^-2+0,48Z^-3=0 (*Z³) 

Z³ + 0.75Z² + 0.39Z +0.48 = 0 

Z³+0.39Z=-0.75Z²-0.48 

Y1=Z³+0.39Z           Y2= -0.75Z²-0.48 

Z	-1	-0,75	-0,5	-0,25	0	0,25	0,5	0,75	1
Y1	-1,39	-0,714	-0,32	-0,113	0	0,1131	0,32	0,7144	1,39
Y2	-1,23	-0,902	-0,668	-0,527	-0,48	-0,527	-0,668	-0,902	-1,23

 

Z₁=-0,85

     Поделим полином на полином для того, чтобы  узнать полюсы и проверить цепь на устойчивость (рекурсивная или не рекурсивная)

  
 
 

Z3	+0.75Z2	+0.39Z	+0.48	Z	+0.85
Z3	+0.85Z2			Z2	-0.12Z	+0.305
	-0.1Z2	+0.39Z	+0.48
	-0.1Z2	+0.085Z
		0.305Z	+0.48
		0.305Z	+0.259
		0	+0.22

 

Z²-0,12Z+0,305=0

D=(-0,12)² – 4*0,305=0,0144 – 1,22= - 1,2056

Расположение  полюсов на комплексной плоскости показано на рис.3.

                                                         Z2

                                                Z1

 

                                                           Z3 
 
 

Рис.3. Расположение полюсов  H(Z) на комплексной плоскости 

Так как все  три полюса передаточной функции  находятся внутри единичного круга Z-плоскости, фильтр является устойчивым (рекурсивным). 
 
 
 
 
 
 
 
 
 
 
 

3. Расчет X(jkw₁) и H(jkw₁) с помощью БПФ.

     Быстрое преобразование Фурье (БПФ) устанавливает связь между отсчетами во временной и частотной областях. Формула БПФ для входного сигнала:

,

где  N – количество отсчетов во временной и частотной областях; 
  - весовая функция. 

Заданная последовательность:

X(nT)={-0,33; -0,99; -0,45; 0; -0,27; -0,19; 0,36;- 0,39 }. 

     Так как количество отсчетов N=8, БПФ производится в три этапа. Определяем для каждого из этапов значения весовых функций.

     Этап 1.  Количество взаимодействующих элементов – 2,  
 W⁰= = 1 .

     Этап 2.  Количество взаимодействующих элементов – 4,  
 W⁰= = 1 ;  W¹= = = -j .

          Этап 3.  Количество взаимодействующих  элементов – 8,  
 W⁰= = 1 ;  W¹= = = 0,7071 - j 0,7071 ; 
 W²= = = -j; W³= = = - 0,7071-j0,7071 

Рассчитываем  X(jkw₁) в программе, которая называется BPF. И получаем: