Сигнальные процессоры. Аксенов В.П. - 129 стр.

title ('частота выходного сигнала от m');


                                                               ПРИЛОЖЕНИЕ 10

         Синтез полосового БИХ-фильтра при частной манипуляции


%          Полосовой БИХ-фильтр, частотная манипуляция
%          функция - ellip (эллиптический фильтр)
clear all; % уничтожаем все переменные
A = 5; Fs =32000; V = 600; % амплитуда, частота дискрет - ии, скорость
Rp = 1.5; Rs = 25; % пульсации и минимальное ослабление в децибелах
f1 = 2.0*V; f2 = 3.0*V; % f1,f2 - частоты входного сигнала
f3 = 0.40*f1; f4=0.85*f1; f5=1.35*f2; f6=1.75*f2;   % точки АЧХ фильтра
w3 = 2*f3/Fs; w4 = 2*f4/Fs; w5 = 2*f5/Fs; w6 = 2*f6/Fs; % нормированная частота
[n, Wn] = ellipord([w4,w5],[w3,w6], Rp,Rs);
% определен порядок n фильтра
[b,a] = ellip (n,Rp,Rs,Wn);    % расчет коэфф-ов
[h,f]=freqz(b,a);              % расчет АЧХ фильтра
subplot(1,2,1) % 1-ый график АЧХ с автоматическим выбором масштаба
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-60 10])      % диапазон АЧХ
grid on          % включение сетки на графике
title ('АЧХ в децибелах');
subplot(1,2,2) %               2-ой график АЧХ для определения пульсаций
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-Rp Rp])      % диапазон Rp в децибелах
grid on          % включение сетки на графике
title ('АЧХ в децибелах');
disp порядок_фильтра ;
n
    % pause
hh=impz(b,a);       % расчет импульсной характеристики
figure
plot(hh)         % график импульсной характеристики



                                            129