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

grid on          % включение сетки на графике
title ('АЧХ в децибелах');
disp порядок_и_параметр_фильтра ;
n,beta
  % pause
hh=impz(b);      % расчет импульсной характеристики
figure
plot(hh)         % график импульсной характеристики
title ('Импульсная характеристика');
  % pause
interv=8; % количество битовых интервалов
t = 0:1/Fs:interv/V;    % диапазон времени равен interv битовым интервалам
s11 = A*sin(2*pi*f1*t).*(t<1/V);    % сигнал ЧМ на 1-ом битовом интервале
s12 = -A*sin(2*pi*f2*t).*(t>1/V).*(t<2/V);         % 2-ой битовый интервал
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,1,s); % s1 - выходной сигнал фильтра
subplot(1,2,2)
plot(t,s1)
xlim([(interv-4)/V interv/V])   % вых-ой сигнал на 4 последних интервалах
grid on
title ('выходной сигнал на 4 интервалах');
 % pause
    % расчет мгновенной частоты по трем точкам входного сигнала s(t)
maxt =length(t)-5;     % последние 5 точек вектора t не используются



                                             127