[an error occurred while processing this directive]
|
радиотехнической практикой в смысле радиосигналов я не занимался вообще. Я больше по проводам. fir2 кстати это оконный метод, и мало чем отличается от того, что Вы ранее гоняли вручную через fft. А интересующие Вас функции - это вероятно firlpnorm - где можно задавать весовые коэффициенты для ошибки для разных полос, задавать переходные полосы, не подлежащие оптимизации и нормой ошибки играться с режимом оптимизации от минимума квадрата почти до минимаксной. Но там все кусочно-линейно.
Ну или ремез по имени gremez, которому можно скормить как параметр произвольную функцию АЧХ, определить особые точки, задавать разные весовые коэффициенты в разных полосах, и вообще чего только не делать.
b = gremez(n,f,fresp,w) returns a length N+1 FIR filter which has the best approximation to the desired frequency response as returned by the user-defined function fresp. gremez uses the following syntax to call fresp
тута n - порядок, f - вектор с определением полос, w весовые коэффициенты в этих полосах. Ну а [dh,dw] = fresp(n,f,gf,w) - это функция, Вами написанная, которая рассчитывает ЧХ по входному вектору частот gf. При этом выходными параметрами являются dh - это собственно вектор коэффициентов передачи. и dw - вектор весовых коэффициентов.
E-mail: info@telesys.ru