我想要：四组图在一起的那种
x=[10 11 12 13 14 15 16];
y1=[6.51 6.32 6.49 7.09 6.82 7.15 7.06];
y2=[7.39 6.59 2.00 2.26 1.85 3.09 2.95];
y3=[7.1 7.1 7.5 7.7 8.3 7.4 7.3];
y4=[1.52 1.35 2.03 1.71 1.56 1.93 1.56];
scatter (x,y,'b.') %最大流量与年份散点图

w=polyfit(x,y,3);
plot(w)


x=[10 11 12 13 14 15 16];
y=[6.51 6.32 6.49 7.09 6.82 7.15 7.06];
y2=[7.39 6.59 2.00 2.26 1.85 3.09 2.95];
y3=[7.1 7.1 7.5 7.7 8.3 7.4 7.3];
y4=[1.52 1.35 2.03 1.71 1.56 1.93 1.56];
scatter (x,y,'b.') %最大流量与年份散点图

plot(x,y,'b-*') %最大流量与年份折线图

p=polyfit(x,y,6);  % 拟合出的六次函数的系数
ye=y-polyval(p,x);  % 计算误差
ye2s=sum(ye.^2); % 误差的平方和
disp(sprintf('误差的平方和=%d',ye2s));
xx=linspace(min(x),max(x));  % 绘图用到的点的横坐标
yy=polyval(p,xx);   % 拟合曲线的纵坐标
plot(x,y,'*-',xx,yy);  % 绘图，原始数据+拟合曲线
legend('原始数据','拟合曲线');  % 图示
s=char(vpa(poly2sym(p,'x'),5));  % 二次函数式转换为字符串，vpa转换小数，保留5位有效数字
title('赣江南昌站最小流量图');
        xlabel('年份');
        ylabel('流量');
        text(['y=' s]);



randn('state',1),rand('state',31) %初始化
x=randn(100,1);y=rand(100,1); %生成正态和均匀分布实验样本
%观察正态数据组的频数直方图5.9.2-1 在不同区间分段数时的变化
subplot(1,2,1),hist(x,7) %7区间情况
subplot(1,2,2),histfit(x,20) %15 区间情况（带正态拟合线）