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

subplot(1,2,2) %                2-ой график АЧХ для определения пульсаций
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-Rp Rp])     % диапазон Rp в децибелах
grid on          % включение сетки на графике
title ('АЧХ в децибелах');
disp порядок_фильтра ;
n
hh=impz(b,a);      % расчет импульсной характеристики
figure
plot(hh)         % график импульсной характеристики
xlim([0 150])
title ('Импульсная характеристика');
interv=8; % количество битовых интервалов
t = 0:1/Fs:interv/V;   % диапазон времени равен interv битовым интервалам
s11 = A*sin(2*pi*f1*t).*(t<1/V);    % сигнал ЧМ на 1-ом битовом интервале
s12 = 0.*(t>1/V).*(t<2/V);              % 2-ой битовый интервал
s13 = A*sin(2*pi*f1*t).*(t>2/V).*(t<3/V);      % 3-ий битовый интервал
s14 = 0.*(t>3/V).*(t<4/V);              % 4-ый битовый интервал
s15 = A*sin(2*pi*f1*t).*(t>4/V).*(t<5/V);      % 5-ый битовый интервал
s16 = 0.*(t>5/V).*(t<6/V);              % 6-ой битовый интервал
s17 = A*sin(2*pi*f1*t).*(t>6/V).*(t<7/V);      % 7-ой битовый интервал
s18 = 0.*(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-6)/V interv/V])   % вых-ой сигнал на 6 последних интервалах
grid on
title ('выходной сигнал на 6 интервалах');
     % расчет мгновенной частоты по трем точкам входного сигнала s(t)



                                             124