pdf.rar_DAA_K海杂波_海尖峰_海杂波模型_海杂波的模型拟合分布

% xn=R(:,1); % xn1=double(xn); % y=xn1(:); % % [pdft,x]=ksdensity(y); % [pdft,x] = ksdensity(y); % plot(x,pdft); % axis([0 10 0 2]); xn=R(:,1); [y,z] = ksdensity(xn); % plot(z,y); hold on % % 正态分布 % [f,s]=normfit(xn); % y1=normpdf(xi,f,s); % [h1,p1,k1,cv1]=kstest(xn,[xn,normcdf(xn,f,s)]); %ks检验 % plot(xi,y1); % axis([-1 10 0 1]); % hold on; % %瑞利分布 % [p,n]=raylfit(xn,0.05); % y2=raylpdf(xi,p); % [h2,p2,k2,cv2]=kstest(xn,[xn,raylcdf(xn,p)]); %ks检验 % % axis([-1 10 0 1]); % plot(xi,y2); % hold on; z = 1e-9:0.1:(9.9+1e-9); seta = rms(R(:,1))/sqrt(2); fz0 = z/(seta.^2) .* exp(-(z.^2)/(2*(seta.^2))); axis([-1 10 0 1]); % plot(z,fz0); hold on % %威布尔分布 % a=wblfit(xn); % y3=wblpdf(xi,a(1),a(2)); % [h3,p3,k3,cv3]=kstest(xn,[xn,wblcdf(xn,a(1),a(2))]); %ks检验 % plot(xi,y3); % axis([0 10 0 1]); % z = 1e-9:0.1:(9.9+1e-9); A1 = mean(R(:,1)) * mean(R(:,1).^(-1)); % Weibull distribution parameter estimation A2 = mean(R(:,1).^(2)) * mean(R(:,1).^(-2)); B1 = mean(R(:,1)); alpha = pi*A2/(sqrt(A2^2-A1^2)*A1); beta = B1/gamma(1+1/alpha); fz1 = alpha/beta * ((z/beta).^(alpha-1)) .* exp(-(z/beta).^alpha); % plot(z,fz1); axis([0 10 0 1]); hold on; plot(z,y,z,fz0,z,fz1); % %Gamma分布 % pha=gamfit(xn); % y4=gampdf(xi,pha(1),pha(2)); % [h5,p5,k5,cv5]=kstest(xn,[xn,gamcdf(xn,pha(1),pha(2))]); %ks检验 % axis([0 10 0 1]); % hold on; %对数正态分布 par=lognfit(xn); y5=lognpdf(z,par(1),par(2)); [h4,p4,k4,cv4]=kstest(xn,[xn,logncdf(xn,par(1),par(2))]); %ks检验 axis([0 10 0 1]); hold on; % % K分布 % m2 = mean(R(:,1).^2);% K distribution parameter estimation % m4 = mean(R(:,1).^4); % v = (m4/(2*(m2^2))-1)^(-1); % b = 2 * sqrt(v/m2); % fz2 = 2*b/gamma(v) .* (b*z/2).^v .* besselk(v-1,b*z); %K distribution % axis([0 10 0 4]); % hold on plot(z,y,z,fz0,z,fz1,z,y5); % axis([0 10 0 4]); legend('实测数据概率密度分布','Rayleigh','Weibull','Lognormal') title('纯杂波拟合PDF'); xlabel('杂波幅度');ylabel('概率密度值'); % plot(xi,y,'r-',xi,y1,'k:',xi,y2,'m-.',xi,y3,'y--',xi,y4,'c',xi,y5,'k');