Спроектируйте конечную импульсную характеристику lowpass фильтр порядка 255 с переходной областью между $$0.25\pi$$ и $$0.3\pi$$. Использование fvtool для отображения величины и фазовых откликов фильтра.

b = firls(255,[0 0.25 0.3 1],[1 1 0 0]);
fvtool(b,1,'OverlayedAnalysis','phase')

Идеальный дифференциатор имеет частотную характеристику, заданную как $$D(\omega )=j\omega$$. Создайте дифференциатор порядка 30, который ослабляет частоты выше $$0.9\pi$$. Включите коэффициент $$\pi$$ в амплитуде, потому что частоты нормированы $$\pi$$. Отобразите нулевую фазовую характеристику фильтра.

b = firls(30,[0 0.9],[0 0.9*pi],'differentiator');

fvtool(b,1,'MagnitudeDisplay','zero-phase')

Создайте антисимметричный фильтр 24-го порядка с кусочно-линейными полосами пропускания.

F = [0 0.3 0.4 0.6 0.7 0.9]; 
A = [0 1.0 0.0 0.0 0.5 0.5];
b = firls(24,F,A,'hilbert');

Постройте график желаемой и фактической частотных характеристик.

[H,f] = freqz(b,1,512,2);
plot(f,abs(H))
hold on
for i = 1:2:6, 
   plot([F(i) F(i+1)],[A(i) A(i+1)],'r--')
end
legend('firls design','Ideal')
grid on
xlabel('Normalized Frequency (\times\pi rad/sample)')
ylabel('Magnitude')

Проектируйте конечную импульсную характеристику lowpass фильтр. Диапазон пропускания областей значений от постоянного тока до $$0.45\pi$$ рад/образец. Диапазон пробок колеблется от $$0.55\pi$$ рад/выборка на частоту Найквиста. Создайте три различных проекта, изменяя веса полос в метод наименьших квадратов.

В первом проекте сделать вес стопорной полосы выше веса полосы пропускания в 100 раз. Используйте эту спецификацию, когда очень важно, чтобы величина реакция в области значений была плоской и близкой к 0. Эта неравномерность в полосе пропускания примерно в ыше, чем пульсация стопорной полосы.

bhi = firls(18,[0 0.45 0.55 1],[1 1 0 0],[1 100]);

Во втором проекте реверсировать веса так, чтобы вес полосы пропускания в 100 раз превышал вес полосы упора. Используйте эту спецификацию, когда очень важно, чтобы величина реакция в полосе пропускания была плоской и близкой к 1. Пульсация стопорной полосы примерно в 100 раз выше, чем неравномерность в полосе пропускания.

blo = firls(18,[0 0.45 0.55 1],[1 1 0 0],[100 1]);

В третьем проекте придайте одинаковый вес обеим полосам. Результатом является фильтр с подобной пульсацией в полосе пропускания и полосе упора.

b = firls(18,[0 0.45 0.55 1],[1 1 0 0],[1 1]);

Визуализируйте величину ответы трех фильтров.

hfvt = fvtool(bhi,1,blo,1,b,1,'MagnitudeDisplay','Zero-phase');
legend(hfvt,'bhi: w = [1 100]','blo: w = [100 1]','b: w = [1 1]')