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

s13 = A*sin(2*pi*f1*t).*(t>2/V).*(t<3/V);             % 3-ий битовый интервал
s14 = -A*sin(2*pi*f2*t).*(t>3/V).*(t<4/V);             % 4-ый битовый интервал
s15 = A*sin(2*pi*f1*t).*(t>4/V).*(t<5/V);             % 5-ый битовый интервал
s16 = -A*sin(2*pi*f2*t).*(t>5/V).*(t<6/V);            % 6-ой битовый интервал
s17 = A*sin(2*pi*f1*t).*(t>6/V).*(t<7/V);             % 7-ой битовый интервал
s18 = -A*sin(2*pi*f2*t).*(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) = 0.5; za1(1) =0; phi1(1) = 0;
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));
end
z21(m) = Fs*acos(z1)/(2*pi);    % частота входного сигнала
if abs(z21(m)) > 0.01 za1(m) = Fs*sqrt(s(m+1)^2-s(m)*s(m+2))/(2*pi*z21(m));
else za1(m) = 0.005;    % мгновенная амплитуда входного сигнала
end
 end



                                                119