Тип линейного фильтра фазы FIR
t = firtype(b)t = firtype(d)Разработайте два КИХ-фильтра с помощью метода окна, одного из даже порядка и другого нечетного порядка. Определите их типы и постройте их импульсные ответы.
subplot(2,1,1) b = fir1(8,0.5); impz(b), title(['Type ' int2str(firtype(b))]) subplot(2,1,2) b = fir1(9,0.5); impz(b), title(['Type ' int2str(firtype(b))])
Разработайте два equiripple Гильбертовых преобразователя, один из даже порядка и другой нечетный порядок. Определите их типы и постройте их импульсные ответы.
subplot(2,1,1) b = firpm(8,[0.2 0.8],[1 1],'hilbert'); impz(b), title(['Type ' int2str(firtype(b))]) subplot(2,1,2) b = firpm(9,[0.2 0.8],[1 1],'hilbert'); impz(b), title(['Type ' int2str(firtype(b))])
Используйте designfilt, чтобы разработать фильтры от предыдущего примера. Отобразите их типы.
d1 = designfilt('lowpassfir','DesignMethod','window', ... 'FilterOrder',8,'CutoffFrequency',0.5); disp(['d1 is of type ' int2str(firtype(d1))])
d1 is of type 1
d2 = designfilt('lowpassfir','DesignMethod','window', ... 'FilterOrder',9,'CutoffFrequency',0.5); disp(['d2 is of type ' int2str(firtype(d2))])
d2 is of type 2
d3 = designfilt('hilbertfir','DesignMethod','equiripple', ... 'FilterOrder',8,'TransitionWidth',0.4); disp(['d3 is of type ' int2str(firtype(d3))])
d3 is of type 3
d4 = designfilt('hilbertfir','DesignMethod','equiripple', ... 'FilterOrder',9,'TransitionWidth',0.4); disp(['d4 is of type ' int2str(firtype(d4))])
d4 is of type 4