var_my.m

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

function [sig2] = var_my(x,n,ip,iq,a,b)
      
y = polefilt(x,n,iq,b);
e = zerofilt(y,n,ip,a);
sum=0.;
for i=1:n 
	sum=sum+e(i)*e(i);
end
sig2=sum/n;