海杂波统计建模 K分布

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））]）;  %