How to implement a generic 'filter' function using direct form I?

Question

vinicius lanziotti 2019-6-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/467020-how-to-implement-a-generic-filter-function-using-direct-form-i

回答： Walter Roberson 2019-6-14

I need a code for the filter function in direct form I, like as I implemented a filter function in direct form II transposed below:

function[y] = filtro(a,b,x)
% Generic filter function in direct form II transposed
    N = size(b,2) - 1;% filter order 
    y = zeros(size(x)); % Initializes y with zeros
    %% FOR THE CASE FIR
    if not(size(a,2) == size(b,2))
        a(1) = 1;
        a = zeros(size(b));
    end
    %% Recursive implementation of equations to differences
    aux = zeros(1,N); 
    for r = 1 :size(x,2)      
        y(r) = b(1)*x(r) + aux(1);
        for n = 2 : N   
            aux(n-1) = b(n)*x(r) + aux(n);
            aux(n-1) = aux(n-1)- a(n)*y(r);
        end
        aux(N) = b(N+1)*x(r) - a(N+1)*y(r);
    end
end % END OF FUNCTION

Thanks!

2 个评论
显示无隐藏无

Walter Roberson 2019-6-14

What is the difference between this question and the other one you asked about these filters?

vinicius lanziotti 2019-6-14

在 MATLAB Online 中打开

The difference is that if I use the command "edit(which('dfilt.df1'))" the function has only 2 inputs (num, den), and use others parameters already defined in Matlab, like (dfilt).

% Using edit(which('dfilt.df1'))
function Hd = df1(num,den)
Hd = dfilt.df1;
Hd.ncoeffs = [1 1];
Hd.FilterStructure = 'Direct-Form I'; 
Hd.tapIndex = [0 0];
if nargin>=1
  Hd.Numerator = num;
end
if nargin>=2
  Hd.Denominator = den;
end

But I need of 3 inputs, like the filter function. But this function was made using direct-form II transposed, and I want the direct-form I from function "filter". And no one parameters already defined in Matlab.

请先登录，再进行评论。

请先登录，再回答此问题。