sisodecon.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function y=sisodecon(s)

n=32;
h=8;
r=2^(-1/h);
for i=1:n
    a(i)=r^(i-1);
end
x=filter(a,1,s);