海
我需要在应用程序中使用这个 Kolmogorov 过滤器。您将一些测量数据放入其中,然后使用过滤器对其进行一些平滑处理。我尝试用“nchoosek”来做,但是当我尝试为 50 或更多的 I 做这个时,它需要的时间太长了。
有人知道如何以更快的方式做到这一点吗?
function [ filterd ] = kolmo(data, inter)
temp = 0;
temp1 = 0;
filterd(1:10, 1) = NaN;
for t=inter+1:(length(data)-inter)
for o=-inter:inter
temp = temp + (nchoosek(2*inter, (inter+o))*data(t+o));
temp1 = temp1 + nchoosek(2*inter, (inter+o));
end
filterd(t, 1) = temp/temp1;
temp = 0;
temp1 = 0;
end
end
谢谢安迪