22-3.m

高阶谱分析包

load nl1
[h1,q1]=nlpow (x,y,128);