Как реализовать фильтр нижних частот MATLAB в C

хорошо, вы сказали, что вы получаете свои коэффициенты из фильтра [B, A] = oil (..) аналогично (попробуйте получить их в цифровом фильтре Z домена AKA), эти коэффициенты A, B соответствуют простой передаточной функции, которую вы знаете

H(z) = B(z)/A(z) = (b(1)+b(2) z^−1+⋯+b(n+1) z^−n)/(a(1)+a(2) z^−1+⋯+a(n+1) z^−n)

правильно?

вам просто нужно помнить, что выход y = H(z)*x или другими словами

y = B(z)/A(z) * x and finally A(z)*y = b(z)*x

и чему было равно x(t) * z^-1? да х(т-1)

это означает, что вы закончите с ecuation, подобным:

y(t)*a(1)+y(t-1)*a(2)+⋯+y(t-n)a(n+1) = x(t)*b(1)+x(t-1)*b(2)+⋯+x(t-n)b(n+1)

и нам нужно фактическое значение y (t) с известными значениями фактического x (t) и прошлого x (t-1) и т. д., а также с известными и сохраненными значениями прошлого y (t-1) и т. д. ..

y(t) = 1/a(1) * (x(t)*b(1)+x(t-1)*b(2)+⋯+x(t-n)b(n+1) - y(t-1)*a(2)-⋯-y(t-n)a(n+1))

это означает, что вам нужны два массива для x и y, и примените уравнение с массивами B и A, которые вы получили от Matlab...

к сожалению, это предполагает, что вы УЖЕ приняли во внимание время выборки в масле() (следовательно, Wn должно быть нормализовано) и убедитесь, что вы берете свои выборки в это точное время выборки (и в идеале также вычисляете свой результат в точное время)