j=90*rand(1600,1);v=18*rand(1600,1)+2;A=[sqrt(sort(j)),sort(v)];%以上是模拟数据，将A换为自己的数据wv=2:0.5:20;ww=length(wv);R=[zeros(ww,1)wv'];j=A(:,1);v=A(:,2);forp=1:wwL=v1L=L&v>=wv(p-1);endR(p,1)=mean(j(L));endplot(R(:,2)
