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

end
z22(m) = Fs*acos(z1)/(2*pi); % расчет частоты выходного сигнала
end
m=1:maxt;
subplot(1,2,2)
plot(m,z22)                  % график мгновенной частоты FВЫХ (m)
grid on                      % время t = m/FS
title ('частота выходного сигнала от m');


                                                                 ПРИЛОЖЕНИЕ 4

 Синтез режекторного БИХ-фильтра при амплитудной манипуляции


%           Режекторный БИХ-фильтр, амплитудная манипуляия
%           функция ellip
clear all; % уничтожаем все переменные
A = 5; Fs =32000; V = 600; % амплитуда, частота дискрет - ии, скорость
Rp = 1.5; Rs = 25; % пульсации и минимальное ослабление в децибелах
f1 = 6.0*V; f2 = f1; % f1 - частота входного сигнала
f3 = 1.20*f2; f4=1.65*f2; f5=1.95*f2; f6=2.35*f2;   % точки АЧХ ВЧ-фильтра
w3 = 2*f3/Fs; w4 = 2*f4/Fs; w5 = 2*f5/Fs; w6 = 2*f6/Fs; % нормированная частота
[n, Wn] = ellipord([w3,w6],[w4,w5],Rp,Rs);
% определен порядок n фильтра
% n=3; % меняем порядок n фильтра, найденный функцией ellipord
[b,a] = ellip (n,Rp,Rs,Wn,'stop');   % расчет коэфф-ов
[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          % включение сетки на графике


                                             112