function?[hurst]=estimate_hurst_pic（nopolinput_data）
input_dir=‘.......‘;%%%%%%設置數據讀取路徑
data_name=[input_dirpol‘_‘input_data];
data0=importdata（data_name）;
for?i=1:14
????data=data0（:i）;
%產生hurst指數時間尺度圖并計算Hurst指數
[MN]=size（data‘）;
data=data‘-ones（MN）*mean（data‘）;%減去均值
[da]=wavedec（data16‘db4‘）;%16級小波分解
[cd1cd2cd3cd4cd5cd6cd7cd8cd9cd10cd11cd12cd13cd14cd15cd16]=detcoef（da[12345678910111213141516]）;%提取小波分解細節分量

%%%%%%%%%%%經驗證4—12尺度的分形維數效果較好

yvals=[var（cd4）var（cd5）var（cd6）var（cd7）var（cd8）var（cd9）var（cd10）var（cd11）var（cd12）];%取方差
xvals=[456789101112];
x=xvals;
y=log2（yvals）;

p2=polyfit（xy1）;%線性擬合
hurst（i）=abs（（p2（1）-1）/2）;?%?Hurst?exponent?is?the?slope?of?the?linear??r（r=2H+1）
end
?return;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????845??2012-02-20?16:19??基于空間分形特性差異的海雜波目標檢測\estimate_hurst_pic.m

?????文件???????4808??2012-02-20?16:18??基于空間分形特性差異的海雜波目標檢測\time_fractal.m

?????目錄??????????0??2012-02-20?16:20??基于空間分形特性差異的海雜波目標檢測

-----------?---------??----------?-----??----

?????????????????5653????????????????????3