6

我试图通过实现一阶 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.

结果： 在此处输入图像描述

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

提前致谢。

4

2 回答 2

3

您希望您的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

通过此校正，响应变为

在此处输入图像描述

于 2012-08-04T22:44:21.817 回答
2

Ben Voigt 给出了正确的答案。他只是没有展示如何更改 bitcyber 的代码以使代码生成他的图表。我需要这个StackExchange 答案来理解 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.

谢谢你的回答本。只是我比大多数人慢一点。

于 2015-11-28T00:30:59.980 回答