1. Айфичер Э. С., Джервис Б. У. Цифровая обработка сигналов: практический подход, 2-е издание. : Пер. с англ. - М.: Издательский дом «Вильямс», 2004. - 992 с.
. Гольденберг
Л. М. и др. Цифровая обработка сигналов: Учеб. пособие для вузов/ Л. М.
Гольденберг, Б. Д. Матюшкин, М. Н. Поляк. - 2-е изд. перераб. и доп.─ М.:
Радио и связь, 1990.─ 256 с.: ил.
Приложение.Текст программы
clc;;
%ПАРАМЕТРЫ ПО ЗАДАНИЮ=2*10^6;=2.1*10^6;=2.7*10^6;=2.8*10^6;=4*10^6;
%%%%%%%%%%%%%%%%%%%%%=178;%порядок фильтра -1=-R/2:1:R/2;=0:1:R;
%НОРМИРОВАНИЕ=2*fmax;n=f1/fd;n=f2/fd;n=f3/fd;n=f4/fd;
%ЧАСТОТЫ СРЕЗА=(f1n+f2n)/2;=(f3n+f4n)/2;
%ФОРМУЛА ИХ Рф=(sin(2*pi*fc1.*n)-sin(2*pi*fc2.*n))./(pi.*n);
h(R/2+1)=1+2*(fc1-fc2);
% ВЫВОД ИХ
figure; (x,h,'LineWidth',2);grid on;
title('ИХ Фильтра без окна');
xlabel('N');('A');;
figure; % Вывод АЧХ ФЧХ ИХ БЕЗ ОКНА
[H,w]=freqz(h,1,'whole');(3,1,1);=w/(2*pi);(fn,abs(H)); %ГрафикАЧХ
xlim([0 0.5]) %ограничение частоты для построения графика('АЧХ');('Нормированная частота'); ylabel('Амплитуда');(gca,'minor') %отображение сетки(3,1,2);(fn,20*log10(abs(H))); %График АЧХ в дБ([0 0.5]); %ограничение частоты для построения графика('АЧХ в дБ');('Нормированная частота'); ylabel('Амплитуда, дБ');(gca,'minor'); %отображение сетки(3,1,3);(fn,angle(H)); %График АЧХ в дБ([0 0.5]) %ограничение частоты для построения графика('АЧХ в дБ');('Нормированная частота'); ylabel('Амплитуда, дБ');(gca,'minor') %отображение сетки=3.384; %Параметр окна Кайзера=kaiser(R+1,b); %Задаем окно Кайзера=Wkz.';
%Построим окно Кайзера.
figure; (x,Wkz); % окноКайзерадляРФ;
title('Окно Кайзера');
%ПЕРЕМНОЖАЕМ ИХ И ОКНО=Wkz.*h;
% ВЫВОД ИХ C ПРИМЕНЕНИЕМ ОКНА
figure; (x,h1,'LineWidth',2);grid on;
title('ИХ с применением окна');
xlabel('N');('A');;
figure; % Вывод АЧХ ФЧХ ИХ С ОКНОМ
[H,w]=freqz(h1,1,'whole');(3,1,1)=w/(2*pi);(fn,abs(H)); %ГрафикАЧХ
xlim([0 0.5]) %ограничение частоты для построения графика('АЧХ');('Нормированная частота'); ylabel('Амплитуда');(gca,'minor') %отображение сетки(3,1,2)(fn,20*log10(abs(H))); %График АЧХ в дБ([0 0.5]) %ограничение частоты для построения графика('АЧХ в дБ');('Нормированная частота'); ylabel('Амплитуда, дБ');(gca,'minor') %отображение сетки(3,1,3)(fn,angle(H)); %График АЧХ в дБ([0 0.5]) %ограничение частоты для построения графика('АЧХ в дБ');('Нормированная частота'); ylabel('Амплитуда, дБ');(gca,'minor') %отображение сетки
% Реакция на типовые воздействия=100;i=0:1:Nt-1; % единичный импульс
ifi<=0;(i+1)=1;(i+1)=0;=0:1:length(U0)-1;;(2,1,1)(n,U0);('Единичный импульс');('N отчета'); ylabel('Значение');;=conv(h1,U0);=0:1:length(y)-1;(2,1,2)(n,y);('Реакция');('N отчета'); ylabel('Значение');;=0:1:Nt-1; % Прямоугольныйимпульсi<=50;(i+1)=1;(i+1)=0;=0:1:length(U1)-1;;(2,1,1)(n,U1);('Прямоугольный импульс');('N отчета'); ylabel('Значение');;=conv(h1,U1);=0:1:length(y)-1;(2,1,2)(n,y);('Реакция');('N отчета'); ylabel('Значение');;=0:1:100;=sin(2*pi*0.2*n)+sin(2*pi*0.3*n)+sin(2*pi*0.4*n);=conv(h1,U);;
[H,w]=freqz(U,1,'whole');(2,1,1)=w/(2*pi);(fn,abs(H)); %ГрафикАЧХ
xlim([0 0.5]) %ограничение частоты для построения графика('АЧХ');('Нормированная частота'); ylabel('Амплитуда');
gridon;(2,1,2)
[H,w]=freqz(y,1,'whole');=w/(2*pi);(fn,abs(H));
xlim([0 0.5]) %ограничение частоты для построения графика('АЧХ');('Нормированная частота'); ylabel('Амплитуда,');;