www.gusucode.com > 步态识别代码视频序列ppt展示及相关文献matlab源码程序 > lslp/forback.m
function a = forback( x, m ) % % FORBACK(x, m) - compute m linear prediction coefficients for signal x % n = size(x,2); Xfor = []; for i = m:n-1 Xfor = [ Xfor x(i:-1:i-m+1 ) ]; end Xbac = []; for i = 2:n-m+1 Xbac = [ Xbac x(i:i+m-1 ) ]; end Xforback = [ Xfor Xbac ]; yforback = [ x(m+1:n)' x(1:n-m)' ]; % least squares solution a = pinv( Xforback ) * yforback; a = a'; end