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

subplot(1,2,1)
plot(m,z21)                   % график мгновенной частоты FВХ (m)
grid on                       % время t = m/FS
title ('частота входного сигнала от m');
      % расчет мгновенной частоты по трем точкам выходного сигнала s1(t)
z1=0.5; z22(1) = (f1+f2)/2;
for m = 1:maxt
if abs(s1(m+1)) > 0.02 z1 =0.5*(s1(m)+s1(m+2))/s1(m+1); else z1 = 0.7;
end
if m > 1 z22(m) = z22(m-1);
end
if abs(z1) > 1 z1 = sign(0.5*(s1(m)+s1(m+2))/s1(m+1));
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');


                                                                    ПРИЛОЖЕНИЕ 2

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


%         БИХ-фильтр верхних частот, амплитудная манипуляция
%         функция - ellip
clear all; % уничтожаем все переменные
A = 5; Fs =32000; V = 600; % амплитуда, частота дискрет - ии, скорость
Rp = 1.5; Rs = 25; % пульсации и минимальное ослабление в децибелах
f1 = 6.0*V; f2 = f1; % f1 - частота входного сигнала
f3 = 0.50*f1; f4=0.85*f1;      % точки АЧХ ВЧ-фильтра
w3 = 2*f3/Fs; w4 = 2*f4/Fs; % нормированная частота
[n, Wn] = ellipord(w4,w3, Rp,Rs);
% определен порядок n фильтра


                                            106