2С. Вычисление реакции ких-фильтра фнч на воздействие в виде периодической последовательности с периодом :

Входными параметрами function-файла являются:

• вектор коэффициентов КИХ-фильтра (см. п. 1C);

• частота дискретизации (она должна совпадать с заданной в требованиях к АЧХ КИХ-фильтра);

• амплитуды гармоник А₁ и А₂;

• частоты гармоник и.

Вводимые значения частот должны быть согласованы с граничными частотами в требованиях к АЧХ КИХ-фильтра. Частота должна быть расположена в ПП, а частота– в ПЗ. При этом отсутствие растекания спектра (см. п. 3С) гарантируется в том случае, если для частотиотношение будет целым числом.

Реакцию КИХ-фильтра вычислить с помощью функцииfilter.

Вывести графики воздействия и реакции КИХ-фильтра.

Выходными параметрами function-файла являются векторы отсчетов воздействия и реакции.

function [x,y] = Untitled2(b,Fs,fk,ft)

DATA=0;

while DATA==0;

A1 = input('A1 = '); %

A2 = input('A2 = '); % Амплитуды гармоник

f1 = input('f1 = '); %

f2 = input('f2 = '); % Частоты гармоник (Гц)

disp('Проверка данных')

disp(' ')

if f1<=ft & f2>=fk

disp('Частоты гармоник введены корректно')

DATA = 1;

else disp ('Частоты гармоник введены некорректно, повторите ввод')

DATA = 0;

end

end

N=64;

n=0:(N-1);

x=A1*cos(2*pi*f1*n/Fs)+A2*cos(2*pi*f2*n/Fs);

y=filter(b,1,x);

subplot(2,1,1), stem(n,x,'fill','MarkerSize',3), grid

xlabel('n'), ylabel('x(n)'), title('Input Signal')

subplot(2,1,2), stem(n,y,'fill','MarkerSize',3), grid

xlabel('n'), ylabel('y(n)'), title('Output Signal ')

end

>> [x,y] = My1(b,Fs,fk,ft)

A1 = 1

A2 = 10

f1 = 162.5

f2 = 2600

Проверка данных

Частоты гармоник введены корректно

ЗС. Вычисление амплитудных спектров воздействия и реакции КИХ-фильтра ФНЧ.

Входными параметрами function-файла являются векторы отсчетов воздействия и реакции КИХ-фильтра.

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

Вывести графики амплитудных спектров.

function [ output_args ] = Untitled3(x,y)

X=fft(x);

Y=fft(y);

N=64;

k=0:(N-1);

subplot(2,1,1), stem(k,X,'fill','MarkerSize',3), grid

xlabel('k'), ylabel('X(k)'), title('Spectrum of Input Signal')

subplot(2,1,2), stem(k,Y,'fill','MarkerSize',3), grid

xlabel('k'), ylabel('Y(k)'), title('Spectrum of Output Signal ')

end

Выполнение работы

Вариант 12.

1. Ввод требований к ачх.

2. Вычисление оценки порядка ких-фильтра, нормированных частот разрыва и параметра β окна Кайзера.

Выведенные значения нормированных частот разрыва (wc) и параметра β (beta) внести в табл. 5.

Таблица 5. Результаты синтеза КИХ-фильтров методом окон

Тип избира–тельности фильтра	Метод окон с окном Кайзера
	порядок фильтра R	вектор нормированных частот разрыва wc	параметр beta
ФНЧ	45	0.32742	3.3953
ФВЧ	46	0.32742	3.3953
ПФ	45	0.22742 ; 0.64677	3.3953
РФ	46	0.22742 ; 0.64677	3.3953