x =[2 3 5 6 6 8 10 11;3 6 9 9 13 9 9 3];
x1=x(1,:);
x11=x1;
x2=x(2,:);
d=x(1,1);
for i=1:6
f1=find(x1<=d+3);
f2=find(x1
break
end
x1(f2)=1000;
jz(i)=mean(x2(f1));
d=d+3;
disp(num2str([x11(f1);x2(f1)]'));
disp(['对应的均值是: ' num2str(jz(i)) char(13)]);
end