How do filter function work?
显示 更早的评论
I understand how the filter(b,a,X) works. What confuses me is the function with the initial conditions. The function filter(b,a,X,zi).
for example: y(n)=-0.5y(n-1)+x(n) x(n)=u(n) y(-1)=1
how will the filter function solve this? i tried to do it manually, but i get different output.
采纳的回答
Honglei Chen
2012-8-26
编辑：Honglei Chen
2012-8-26
MATLAB's filter() function an implementation of direction form II, the detailed algorithm can be found at
For your case, you need to first figure out what the initial condition is. You can use filtic function.
e.g.,
b = 1;
a = [1 0.5];
y = filter(b,a,ones(10,1),filtic(b,a,1));
3 个评论
Treant
2012-8-27
Honglei Chen
2012-8-27
It basically perform a reverse filtering operation. If you want the details, you can look at filtic code.
Sergio Cuadros
2020-6-18
编辑：Sergio Cuadros
2020-6-18
Hello Mr. Honglei, Could you check my question. I have a issue.
Best regards.
更多回答（1 个）
You find a Matlab implementation of the FILTER command at Answers: hard-code-filter. This should reveal the meaning of the initial and final conditions.
Splitting a signal into parts might be helpful also:
signal = rand(1,1000);
[b,a] = butter(3, 0.5); % Arbitrary filter parameters
inital_cond1 = zeros(1, length(a) - 1);
[filtered1, final_cond1] = filter(b, a, signal(1:500), initial_cond1);
[filtered2, final_cond2] = filter(b, a, signal(501:1000), final_cond1);
Joined = [filtered1, filtered2];
% Now filter in one step:
OnePiece = filter(b, a, signal, initial_cond1);
% Compare the results - should be tiny:
disp(mean(abs(Joined - OnePiece)));
This shows, that the inital conditions are the values of the internal status of the filter. When the signal is split into pieces, the final conditions of the former part are the inital conditions of the current part.
类别
在 帮助中心 和 File Exchange 中查找有关 Matched Filter and Ambiguity Function 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
