时间离散实现一阶RC滤波器

Question

我试图通过实现一阶RC滤波器的时间离散版本来确保理解我的（数字）信号处理知识。 （背景是我正在尝试在软件中实现用于SDR目的的PLL，但这是一个不同的故事...）

我的问题是我认为我理解如何为这样的滤波器创建差分方程，因此得出它的系数。然而，当我使用freqz函数在MATLAB中绘制响应时 - 使用计算的a和b系数 - 我看起来不像RC滤波器响应。

我引用了关于这个主题的维基百科页面（http://en.wikipedia.org/wiki/Low-pass_filter#Discrete-time_realization），只是为了确保我没有完全脱离杂草，但它仍然无济于事。这将差异方程详述为：

yi = alpha * xi + ( 1 - alpha ) * yi-1
where: alpha = sample period / ( RC + sample period )

一个例子：

fs = 96000.0;                         % Sample rate.
delta_t = 1.0 / fs;                   % Sample period.
fc = 5000.0;                          % Filter cut off frequency.
tau = 1 / ( 2 * pi * fc );            % Time constant of filter.
alpha = delta_t / ( tau + delta_t );  % Smoothing factor per Wikipedia page.
b = [ alpha ];                        % 'b' coefficients
a = [ 1.0, ( 1 - alpha ) ];           % 'a' coefficents
freqz( b, a, 1024, fs );              % 1024 point FFT used.

结果：

关于我哪里出错的任何想法？我完全误解了什么吗？

提前致谢。

Answer 1

您希望a(2)系数为负数，因为a表示出现在等式左侧的系数。

a(1)*y(n) + a(2)*y(n-1) - ... + a(na+1)*y(n-na) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb)

或等效，

a = a ./ a(1)
y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb)
                 - a(2)*y(n-1) - ... - a(na+1)*y(n-na)

请参阅filter

的文档

---

通过此更正，响应变为

Answer 2

Ben Voigt给出了正确的答案。他只是没有展示如何改变bitcyber的代码以使代码生成他的图形。我需要这个StackExchange answer来理解Ben所说的，然后我将第三行改为“a = [1.0，（alpha -1）];”。

fs = 96000.0;                         % Sample rate.
delta_t = 1.0 / fs;                   % Sample period.
fc = 5000.0;                          % Filter cut off frequency.
tau = 1 / ( 2 * pi * fc );            % Time constant of filter.
alpha = delta_t / ( tau + delta_t );  % Smoothing factor per Wikipedia page.
b = [ alpha ];                        % 'b' coefficients
a = [ 1.0, ( alpha -1 ) ];           % 'a' coefficents
a = a ./ a(1);
freqz( b, a, 1024, fs );              % 1024 point FFT used.

感谢您的回答。只是我比大多数人都要慢一些。