pi.m

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

%
% Function to calculate the objective function where
% input k = temporary vector
%	b = unstable vector

function [ value ] = PI(k, b)

Q=eye(length(k));
a=Levinson(k);
value=0.5*(a(2:length(a))-b)*Q*(a(2:length(a))-b)';