Формирование импульсной характеристики фильтра

function h = build_SCLF(Fs, flag_draw)
 
% Fs - крайняя частота в массиве АЧХ
% flag_draw - флаг прорисовки кривых (АЧХ и ИХ)
% flag_draw = 1 - рисует графику;
% flag_draw = 0 - не рисует графику;
 
%% Проектирование ФНЧ SCLF-10+ для ЦАП
% Частота среза
% Fc = 12.2e6;
 
% Шаг по времени
dt = 1/Fs;
 
% Сетка частот АЧХ
f = [0.5 1 5 8.5 10 11 12.2 12.4 12.6 13.8 15 20, ...
     50 100 150 200 250 300 350 370 400]*1e6;
 
i = find(f == Fs);
 
f = f(1:i); 
 
% f = [0.5 1 5 8.5 10 11 12.2 12.4 12.6 13.8 15 20, ...
%      50 100 150 200]*1e6;
 
% Ослабление на сетке частот, дБ 
attention_dB = [0.06 0.07 0.17 0.35 0.56 2.7 3 14.39 16.31 27.59 39.05, ...
             62.85 72.81 63.39 59.75 72.82 38.33 47.06 59.23 42.02 14.91];
 
% attention_dB = [0.06 0.07 0.17 0.35 0.56 2.7 3 14.39 16.31 27.59 39.05, ...
%              62.85 72.81 63.39 59.75 72.82];
 
attention_dB = attention_dB(1:i);
 
% Перевод ослабления в разы
a = 10.^(attention_dB/20);
 
% АЧХ
H = 1./a;
 
% Значения АЧХ, дБ       
H_dB = 20*log10(H);
 
%% Нахождение импульсной характеристики (ИХ) ФНЧ
% Порядок фильтра 234
Norder = 234;
% Окно 
window = ones(1, Norder+1);
% Нормирование частоты относительно Fs
f_norm = [0 f/Fs];
H = [1 H];
% ИХ
h = fir2(Norder, f_norm, H, window);
 
%% Графический вывод АЧХ
if flag_draw == 1
 
    fig1 = figure;
    scrsz = get(0, 'ScreenSize');
    set(fig1 , 'Position', [400 150 scrsz(3)/1.6 0.7*scrsz(4)]);
    datacursormode on;
 
    subplot(3,1,1);
    plot(f, attention_dB, '-om', 'LineWidth', 2);
    grid on;
    xlabel('{\itf}, МГц'); ylabel('Ослабление, дБ', 'Rotation', 0, 'HorizontalAlignment', 'right');
    title('Ослабление ФНЧ SCLF-10+');
    set(gca, 'Xlim', [0 100]*1e6);
 
    subplot(3,1,2);
    plot(f, H_dB, '-ob', 'LineWidth', 2);
    grid on;
    xlabel('{\itf}, МГц'); ylabel('|{\itH_{ФНЧ}}({\itf})|, дБ', 'Rotation', 0, 'HorizontalAlignment', 'right');
    title('АЧХ ФНЧ SCLF-10+');
    set(gca, 'Xlim', [0 Fs]);
 
    subplot(3,1,3);
    plot(dt*(0:length(h)-1), h, '-r', 'LineWidth', 2);
    grid on;
    xlabel('{\itt}, c'); ylabel('{\ith_{ФНЧ}}({\itt})', 'Rotation', 0, 'HorizontalAlignment', 'right');
    title('ИХ ФНЧ SCLF-10+');
    legend(['Порядок фильтра ', num2str(Norder+1)]);
    
end;