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

title ('АЧХ в децибелах');
disp порядок_фильтра ;
n
hh=impz(b,a);     % расчет импульсной характеристики
figure
plot(hh)         % график импульсной характеристики
xlim([0 100])     % диапазон графика
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-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;



                                             113