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

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



                                            128