Re: Целочисленное проектирование фильтров
Добавлено: 03 дек 2018, 20:59
Здравствуйте, Как в OCTAVE Минимальная фаза ких-фильтр делать?
Теория и практика цифровой обработки сигналов
https://ru.dsplib.org/forum/
Код: Выделить всё
clear all; close all; clc;
pkg load signal;
ORD = 32; % lowpass filter order
%linear phase FIR
h = fir1(ORD, 0.3);
% calculate linear phase FIR roots
z = roots(h);
% all roots outside unit circle will be reflected inside unit circle
y = 1;
for k = 1:length(z)
if(abs(z(k))>1)
z(k) = 1/z(k);
end
y = conv([1, -z(k)], y);
end
% minimum phase FIR
y = y/sum(y)
%plot
figure; plot(0:ORD, h, 0:ORD, y)
хуссейн макассис писал(а): ↑14 фев 2019, 20:18Good day
The coefficients of the FIR filter of the 24th order in the cascade form of construction
have the following below
Звено .. B0 ...... В1 ...... В2
1) -0.461879 -0.639376 -0.528379
2) 0.037240 -0.818323 -0.089716
3) 0.692334 0.238436 0.693423
4) 0.364908 -0.270356 -0.510364
5) -0.385463 0.858949 -0.731630
6) 0.735379 0.426089 0.742754
7) 0.945038 0.141224 -0.146544
8) -0.016565 0.645663 0.771894
9) 0.698864 0.225619 -0.260442
10) 0.525961 0.935365 0.580976
11) 0.644696 0.582789 0.595606
12) -0.233478 -0.014993 0.890144
can you help me to calculate the frequency response, phase response and z-plane of this filter in MATLAB and make it minimally phase .
Несколько слов по важным вопросам «инерционности» цифрового фильтра, поднятым на форуме. Вопросы инерционности, быстродействия, скорости цифровой обработки (фильтрации сигнала) имеют смысл, естественно, только при физической реализации, при работе ЦФ в реальном времени и определяются общим временем расчёта отклика to при реализации фильтра на конкретной цифровой, как вы говорите, платформе - это MCU, DSP, FPGA, БMK, заказные или полузаказные СБИСы. А оно, естественно, является суммой tгз - времени группового запаздывания (ГВЗ) и времени tфр - времени «физической реализации» фильтра на конкретном «железе» - рис.1.
На вашей картинке вы никогда не уменьшите время задержки до величины меньшей tгз.
Здравствуйте!Vlad27 писал(а): ↑07 июн 2019, 12:50Несколько слов по важным вопросам «инерционности» цифрового фильтра, поднятым на форуме. Вопросы инерционности, быстродействия, скорости цифровой обработки (фильтрации сигнала) имеют смысл, естественно, только при физической реализации, при работе ЦФ в реальном времени и определяются общим временем расчёта отклика to при реализации фильтра на конкретной цифровой, как вы говорите, платформе - это MCU, DSP, FPGA, БMK, заказные или полузаказные СБИСы. А оно, естественно, является суммой tгз - времени группового запаздывания (ГВЗ) и времени tфр - времени «физической реализации» фильтра на конкретном «железе»