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

title ('входной сигнал от времени');
s1 = filter (b,a,s); % s1 - выходной сигнал фильтра
subplot(1,2,2)
plot(t,s1)


xlim([(interv-4)/V interv/V])     % вых-ой сигнал на 4 последних интервалах
grid on
title ('выходной сигнал на 4 интервалах');
      % расчет мгновенной частоты по трем точкам входного сигнала s(t)
maxt =length(t)-5;   % последние 5 точек вектора t не используются
m = 1:maxt; % m - количество точек для расчета частоты
con1=length(m);
z1=0.5;z21(1) = 0;z24(1)= 0;
for m = 2:maxt
if abs(s(m+1)) >0.02 z1 =0.5*(s(m)+s(m+2))/s(m+1);else z1 =0.5;
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);      %    частота входного сигнала
% расчет амплитуды входного сигнала
z3=sqrt(abs(s(m+1)^2-s(m)*s(m+2)));
 if abs(z21(m)) > 0.02 z24(m)=z3*Fs/(2*pi*z21(m));
 else z24(m)=z24(m-1) ;%         амплитуда входного сигнала
end
end
m=1:maxt;
con2=length(z21);
figure
subplot(1,2,1)
plot(m,z24)                     % график мгновенной амплитуды A_ВХ (m)
grid on                         % время t = m/FS
title ('амплитуда входного сигнала от m');



                                             108