Оценка частотной характеристики с помощью фильтрации
[ вычисляет оценку частотной характеристики входного фильтра System object™. Функция оценивает, выполняя входные данные, составленные из синусоид с равномерно распределенными случайными частотами через фильтр и формируя отношение между выходными данными и входными данными.h,w] = freqrespest(sysobj)
freqrespest функция вычисляет оценку частотной характеристики объекта входного фильтра, которая близка к частотной характеристике, полученной с использованием только квантованных коэффициентов, при условии, что:
При выполнении фильтрации с фиксированной точкой переполнение не происходит
В сложениях и умножениях с фиксированной точкой не происходит значительного квантования.
freqrespest может использоваться в качестве диагностического инструмента для реализации с фиксированной точкой. Если freqrespest значительно отличается от freqz (что учитывает только квантование коэффициентов), это указывает на то, что, по меньшей мере, одно из двух условий может быть верным.
[ использует h,w] = freqrespest(sysobj,L,Name,Value)Name,Value пар в качестве входных аргументов для указания необязательных параметров теста. Введите имя параметра в качестве входного аргумента в одинарных кавычках, за которым следует возможное значение параметра.
freqrespest( без выходного аргумента запускает FVTool и показывает оценку отклика величины объекта фильтра.sysobj,___)
Оцените частотную характеристику FIR-фильтра с фиксированной точкой.
firFilt = design(fdesign.lowpass(.4,.5,1,60),'equiripple','Systemobject',true); dataIn = fi(randn(16,15),1,16,15); dataOut = firFilt(dataIn); %#ok [h,w] = freqrespest(firFilt); %#ok % This should be about the same as freqz. release(firFilt);
Продолжение работы с объектом фильтра firFilt, измените FullPrecisionOverride свойство для false и затем укажите длины слов и точность (длины дробей), применяемые к выходам из операций внутреннего сложения и умножения. После установки длин слов и дробей используйте freqrespest функцию для вычисления оценки частотной характеристики для фильтра с фиксированной точкой.
firFilt.FullPrecisionOverride = false; firFilt.ProductDataType = 'Custom'; firFilt.CustomProductDataType = numerictype(1,16,15); firFilt.AccumulatorDataType = 'Custom'; firFilt.CustomAccumulatorDataType = numerictype(1,16,15); firFilt.OutputDataType = 'Same as accumulator'; dataOut = firFilt(dataIn); [h,w] = freqrespest(firFilt,2); [h2,w2] = freqz(firFilt,512); plot(w/pi,20*log10(abs([h,h2]))) legend('Frequency response estimated by filtering',... 'Freq. response computed by quantizing coefficients only'); xlabel('Normalized Frequency (\times\pi rad/sample)') ylabel('Magnitude (dB)')
