How can I write an FIR filter function

9 次查看（过去 30 天）
Aliyu Suleiman Liman
评论： Aliyu Suleiman Liman 2019-2-12
Hi! I'm trying to write a function that performs FIR filtering without having to use the built-in function.
so I have this code at the moment.
function y = myFIR(x,h)
M = length(x);
L = length(h);
y = zeros(0,((M+L)-1));
for n = 1:M+1
for k = 1:L+1
Y1 = (h(k).*x(n-k));
y(n)= sum(Y1);
end
end
where (h) is an array consisting of the impulse response, (x) is an array consisting of the input sequence and (y) is a vector containing the filtered samples.
but I have been getting the error below;
Array indices must be positive integers or logical values.
Error in myFIR (line 7)
Y1 = (h(k).*x(n-k));
I have really tried everything I could but it just doesnt seem to go. I will really be glad if anyone can help figure out the problem.
Thank you.
  2 个评论
KALYAN ACHARJYA
KALYAN ACHARJYA 2019-2-12
编辑：KALYAN ACHARJYA 2019-2-12
Can you define x and h inputs, examples?
Aliyu Suleiman Liman
Aliyu Suleiman Liman 2019-2-12
Yh sure!
h = { 0.2 0.3 0.2}
x = {1 2 3 4 5 6 7 8 9 10}

请先登录，再进行评论。

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

回答（1 个）

KALYAN ACHARJYA
KALYAN ACHARJYA 2019-2-12
编辑：KALYAN ACHARJYA 2019-2-12
My suggestion to avoid the error, In the loop, during first iterations when when n=1 and k=1, then following statement
Y1 = (h(k).*x(n-k));
x(n-k) becomes x(0), which is not allowed in indexing vectors. x(Positives only)

类别

在 Help Center 和 File Exchange 中查找有关 Digital Filtering 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by