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

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


                                                                          ПРИЛОЖЕНИЕ 3

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


%         Полосовой КИХ-фильтр, фазовая манипуляция
%         окно Кайзера
clear all; % уничтожаем все переменные
A = 5; Fs =32000; V = 600; % амплитуда, частота дискрет - ии, скорость
Rp = 1.5; Rs = 28; % пульсации и минимальное ослабление в децибелах
Kp =(10^(Rp/10)-1)^0.5; Km = 10^(-Rs/20); % пересчет Rp, Rs в абсолютные значения
f1 = 4.0*V; f2 = f1; % f1 - частота входного сигнала
f3 = 0.40*f1; f4=0.85*f1; f5 = 1.35*f2; f6 = 1.75*f2;       % точки АЧХ ПФ
[n,Wn,beta,ftype] = kaiserord ([f3,f4,f5,f6], [0 1 0], [Km,Kp,Km], Fs);
% определен порядок n фильтра с окном Кайзера


                                             109