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