Лабораторные работы по курсу "Базы данных и знаний"

[b1,a1]=lp2lp(b,a,W₀)
где, b, а – коэффициенты исходного фильтра низких частот; W₀ – желаемая частота среза проектируемого фильтра; b1, a1 – коэффициенты спроектированного фильтра с новыми частотами среза.
В Программе 2 приведены примеры определения передаточных функций ФНЧ (фильтра Баттерворта) для трех частот среза. Командой buttap определяют нули, пояса и коэффициент усиления фильтра Баттерворта заданного порядка n. Команда zp2tf преобразует математическую модель ФНЧ: от информации заданной нулями, полосами и коэффициентом усиления, к информации представленной передаточными функциями.
%Программа 2
wr=0.1; %Частота пропускания в относительных единицах.
ws=0.6; %Частота задерживания в относительных единицах.
rp=6; %Допустимое подавление в полосе пропускания.
rs=40; %Допустимое подавление в полосе задерживания.
[n1,wc1]=buttord(wr,ws,rp,rs) %Определение параметров фильтра
%Баттерворта.
[z1,p1,k1]=buttap(n1); %Определение фильтра Баттерворта в форме ZPK
[b1,a1]=zp2tf(z1,p1,k1)%Определение полинома числителя и знаменателя
%фильтра Баттерворта.
h1=tf([b1],[a1]) %Определение фильтра Баттерворта в форме TF.
w1=10; %Частота среза проектируемого фильтра Баттерворта.
[b10,a10]=lp2lp(b1,a1,w1)Параметры фильтра Баттерворта на частоте w1.
w2=100; %Частота среза проектируемого фильтра Баттерворта.
[b100,a100]=lp2lp(b1,a1,w2) %Параметры фильтра Баттерворта
%с частотой среза w2.
q1=tf([b10],[a10]) %Передаточная функция фильтра Баттерворта
%с частотой среза w1
q2=tf([b100],[a100]) %Передаточная функция фильтра Баттерворта
%с частотой среза w2.
figure(1) %Логарифмические характеристики проектируемых фильтров
bode(h1,'g',q1,'k',q2,'r'),grid on %при разных частотах среза