How to write my own filter function with initial condition similar to Filter.m?

2013 2 5
2 个回答
回答已采纳
更新时间：2019 10 7
5 次查看（30 天）

0 个投票

I have read the documents on how MATLAB implements Filter(b,a,y,z0). To fully understand the concept and learn more about Signal Processing, I tried to implement it by myself, but the answer is different from MATLAB filter.m
Could anyone let me know what I am missingg?! :-(
function y =myfilter(num, den, x, zi)
N = length(x);
% assuming order of Num and Den are the same.
ord = (length(num)-1)/2;
y = zeros(N,1);
for i=1:N
c = 0;
for j=1:2*ord
if i-j < 1
c = c+1;
xo(j) = zi(c);
yo(j) = 0;
else
xo(j) = x(i-j);
yo(j) = y(i-j);
end
end
y(i) = num(1)*x(i) + sum(num(2:end).*xo') - sum(den(2:end).*yo');
end

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

 采纳的回答

Hossein
Hossein 2013-2-6

0 个投票

Ok, I searched online and found out that I was wrong on the concept :) Here' the correct formulation of 1D filtering with delay initial values: SciPY Filtering Documentation, and it outputs the same result as MATLAB filter() does.

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Mariasole Laureti
Mariasole Laureti 2019-10-7
Hi, I have the same problem, can I ask you if you can show me your code? I would be very grateful. Thank you very much

请先登录，再进行评论。

更多回答（1 个）

ssp
ssp 2017-2-12

0 个投票

Hi! Greetings!
I was also trying to write my filter function. However y=filter(b,a,x) seems easy to implement, but I am stuck at implementing y=filter(b,a,x,zi). What would be the new sequence of the following sequence after applying initial condition zi? 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)
Thanks! SSP

类别

在 帮助中心 和 File Exchange 中查找有关 Signal Generation, Manipulation, and Analysis 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by