clear;clc;
fid=fopen（‘D:\Ship\Clutter\2‘‘rb‘）;
width=512;
height=512;
imgdata=fread（fid[widthheight]‘uint16‘）;
fclose（fid）;

L=2;?%視數
imgdata=double（imgdata）;
imgdata2=imgdata.^2;
mean1=mean（mean（imgdata））;
mean2=mean（mean（imgdata2））;
vmuc=1/（mean2/mean1/mean1/（1+1/L）-1）;?%形狀參數

y=imgdata（:）;
[hw]=hist（ymax（y））;
bar（wh/sum（h））;?%直方圖
hold?on;

%高斯分布
[fs]=normfit（y）;
y1=normpdf（wfs）;
[h1p1k1cv1]=kstest（y[ynormcdf（yfs）]）;?%ks檢驗
hold?on;
%瑞利分布
[pn]=raylfit（y0.05）;
y2=raylpdf（wp）;
[h2p2k2cv2]=kstest（y[yraylcdf（yp）]）;??%ks檢驗
hold?on;
%威布爾分布
a=wblfit（y）;
y3=wblpdf（wa（1）a（2））;
[h3p3k3cv3]=kstest（y[ywblcdf（ya（1）a（2））]）;??%ks檢驗
hold?on;
%對數正態分布
par=lognfit（y）;
y4=lognpdf（wpar（1）par（2））;
[h4p4k4cv4]=kstest（y[ylogncdf（ypar（1）par（2））]）;??%