filter函数看名字以为是滤波函数,但又对滤波的概念不是太清晰,所以是一头雾水。。。
看了百度知道里面的一个最佳答案,自己看了半天,没弄明白,不知道怎么那么多人顶,自己试了好多数据,都一次次的推翻得到的结果
最后还是灵机一现,看了出来,哈哈...
先看别人怎么解释的:
看楼主的意思主要是想实现差分方程
先从简单的说起:
filter([1,2],1,[1,2,3,4,5])
实现 y[k]=x[k]+2*x[k-1]
y[1]=x[1]+2*0 %(x[1]之前状态都用0)
=1
y[2]=x[2]+2*x[1]=2+2*1=4
...
然后说你感兴趣的filter(b,a,x)
这里a,b都是矢量
你提到输入[zeros(1,100),ones(1,101)]
差分方程 2y[k]-y[k-1]-3y[k-2]=2x[k]-x[k-1]
往里面代
[1]之前的状态都为0
y[1]=1/2*(2*x[1])=0
y[2]=1/2*(2*x[2]-x[1]+y[1])=0
........
显示的数为差分方程的输出y
这个函数和还有其它的用法,再次我就不说了
matlab的帮助中说的比较详细,你可以自己看一下
我个人得到的结论是:
f=(1,1,1,1,1);
x=(1,2,3,4,5);
y=filter(f,1,x);
y = 1 3 6 10 15
用公式表示:
y1 = x1
y(i) = x(i-1)+f(i)+f(i-1)+...+f(1);
如果对于y=filter(f,a,x)
y1 = x1
y(i)=(x(i-1)+f(i)+f(i-1)+...+f(1))/a;