mywe.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 12 行

function [a,b,sig2,iflag] = mywe(x,n,ip,iq,epsilon) 
      
lag=ip+iq+1; 
r = correlat(n,lag,0,x,x); 
[a iflag] = toepeqn(r,ip,iq,epsilon);
if iflag == -1
	return 
end
y = czeroflt(x,n,ip,a);
l=round(n/5); 
[b sig2] = durbin(y,n,iq,l);