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

title ('Импульсная характеристика');
  % pause
interv=6; % количество битовых интервалов
t = 0:1/Fs:interv/V;    % диапазон времени равен interv битовым интервалам
s11=A*sin(2*pi*f1*t).*(t<1/V); phi1=2*pi*f1/V; % сигнал ЧМ на 1-ом интервале
s12=A*sin(2*pi*f2*(t-1/V)+phi1).*(t>1/V).*(t<2/V);phi2=2*pi*f2/V;% 2-ой интервал
s13=A*sin(2*pi*f1*(t-2/V)+phi1+phi2).*(t>2/V).*(t<3/V);       % 3-ий интервал
s14=A*sin(2*pi*f2*(t-3/V)+2*phi1+phi2).*(t>3/V).*(t<4/V);       % 4-ый интервал
s15=A*sin(2*pi*f1*(t-4/V)+2*phi1+2*phi2).*(t>4/V).*(t<5/V);      % 5-ый интервал
s16=A*sin(2*pi*f2*(t-5/V)+3*phi1+2*phi2).*(t>5/V).*(t<6/V);      % 6-ой интервал
s17=A*sin(2*pi*f1*(t-6/V)+3*phi1+3*phi2).*(t>6/V).*(t<7/V);      % 7-ой интервал
s18=A*sin(2*pi*f2*(t-7/V)+4*phi1+3*phi2).*(t>7/V).*(t<8/V);      % 8-ой интервал
s = s11+s12+s13+s14+s15+s16+s17+s18;          %    расчет входного сигнала
figure
subplot(1,2,1)
plot(t,s)    % график входного сигнала от времени
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) =(f1+f2)/2;
for m = 2: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));



                                             130