bispeci 直接法用于双谱估计

直接法用于双谱估计，比较简单，适合初学者使用

function [Bspecwaxis]=bispeci（ynlagnsampoverlapflagnfftwind）
%Bspec:由输入数据y估计出的双谱
%nlag：累积量计算的最大延迟
%nfft：FFT计算的长度，缺省值为128
%wind：窗口选择，wind=0，采用Parzen窗（缺省时），否则为六角形窗
%nsamp：每个数据分段的长度，缺省值为信号长度的1/8
%overlap：每个数据段间重叠的点段（0~99），缺省值为50
%flag：用于指定估计是有偏的还是无偏的
%——————————————Parameter checks—————————————————
[lynrecs]=size（y）;
if（ly==1）                  %如果输入信号y为1行，则将数据矩阵由行矩阵变为列矩阵原来的列数变为行数，现在的列数为1
    y=y（:）;  
    ly=nrecs;nrecs=1; 
end
if （exist（‘overlap‘）~=1）         %如果overlap是一个数，则不进入条件句
    overlap=0; 
end
overlap=min（99max（overlap0））;  %不进入条件语句的情况下，overlap等于设定值
if（nrecs>1）              
    overlap=0;          
end
if（exist（‘nsamp‘）~=1）    
    nsamp=ly;           
end
if（nsamp>ly|nsamp<=0）           %数据分段的长度不能大于输入数据本身，也不能小于0，否则取长度为输入数据长度本身
    nsamp=ly;      

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4110  2012-05-03 19:41  bispeci.m

----------- ---------  ---------- -----  ----

                 4110                    1