清空筛选条件
清空筛选条件

How to correctly apply dsp.FIRFilter

5 次查看（过去 30 天）
Kera Cai
Kera Cai 2024-4-3
评论： Kera Cai 2024-4-15
The code is simple:
N = 40;
fs = 20 * 10 ^ 9;
freq = [0, 1 * 10 ^ 9, 2 * 10 ^ 9, 3 * 10 ^ 9, 5 * 10 ^ 9, 7.5 * 10 ^ 9, 10 * 10 ^ 9]
mag = [1, 1, 1, 0.3, 0.1, 0.005, 0.001]
d = fdesign.arbmag('N,F,A', N, freq, mag, fs);
W = [1, 1, 1, 1, 1, 1, 1];
fir = design(d, 'equiripple', 'weights', W, 'SystemObject', true);
sents = ones(1, 300) * amplitude;
rcvds = fir(sents);
% plot
figure(1)
stem(sents)
hold on
stem(rcvds)
legend('sents', 'received')
By viewing freqz(fir), seems the filter is correctly designed, it's low pass:
But in the code I input a DC signal, the result does not seem to be correct, the output is almost 0.
Is the way I perform the filter wrong?
Seems I got correct result when apply the filter by function filter:
fir = fir.Numerator
rcvds = filter(fir, 1, sents)
figure(2)
stem(sents)
hold on
stem(rcvds)
legend('sents', 'received')

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

采纳的回答

Paul
Paul 2024-4-7
Ran in:
The input to fir() should be a column vector.
N = 40;
fs = 20 * 10 ^ 9;
freq = [0, 1 * 10 ^ 9, 2 * 10 ^ 9, 3 * 10 ^ 9, 5 * 10 ^ 9, 7.5 * 10 ^ 9, 10 * 10 ^ 9];
mag = [1, 1, 1, 0.3, 0.1, 0.005, 0.001];
d = fdesign.arbmag('N,F,A', N, freq, mag, fs);
W = [1, 1, 1, 1, 1, 1, 1];
fir = design(d, 'equiripple', 'weights', W, 'SystemObject', true);
sents = ones(1, 300);% * amplitude;
rcvds = fir(sents.');
% plot
figure(1)
stem(sents)
hold on
stem(rcvds)
legend('sents', 'received')

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Single-Rate Filters 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by