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

b = fir1(n, Wn, ftype, kaiser (n+1, beta), 'noscale');      % расчет коэфф-ов ПФ
[h,f]=freqz(b);                % расчет АЧХ фильтра
subplot(1,2,1) % 1-ый график АЧХ с автоматическим выбором масштаба
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
grid on           % включение сетки на графике
title ('АЧХ в децибелах');
subplot(1,2,2) %                 2-ой график АЧХ для определения пульсаций
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-Rp Rp])      % диапазон Rp в децибелах
grid on           % включение сетки на графике
title ('АЧХ в децибелах');
disp порядок_и_параметр_фильтра ;
n,beta
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 - выходной сигнал фильтра



                                                110