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

% n=60; % меняем порядок n фильтра, найденный функцией ellipord
[b,a] = ellip(n,Rp,Rs,Wn,'high');    % расчет коэфф-ов
[h,f]=freqz(b,a);                % расчет АЧХ фильтра
subplot(1,2,1) % 1-ый график АЧХ с автоматическим выбором масштаба
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-60 10])      % диапазон
grid on          % включение сетки на графике
title ('АЧХ в децибелах');
subplot(1,2,2) %               2-ой график АЧХ для определения пульсаций
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-Rp Rp])      % диапазон Rp в децибелах
grid on          % включение сетки на графике
title ('АЧХ в децибелах');
disp порядок_фильтра ;
n
    % pause
hh=impz(b,a);       % расчет импульсной характеристики
figure
plot(hh)         % график импульсной характеристики
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)     % график входного сигнала от времени



                                            107