1、Anscombe’s Quartet的数据
数据的来源 http://www.matrix67.com/blog/archives/2308
文本格式的数据直接去上面的网站复制就好啦,或者下面是我在txt里面粘贴过来的
10.0 8.04 10.0 9.14 10.0 7.46 8.0 6.58
8.0 6.95 8.0 8.14 8.0 6.77 8.0 5.76
13.0 7.58 13.0 8.74 13.0 12.74 8.0 7.71
9.0 8.81 9.0 8.77 9.0 7.11 8.0 8.84
11.0 8.33 11.0 9.26 11.0 7.81 8.0 8.47
14.0 9.96 14.0 8.10 14.0 8.84 8.0 7.04
6.0 7.24 6.0 6.13 6.0 6.08 8.0 5.25
4.0 4.26 4.0 3.10 4.0 5.39 19.0 12.50
12.0 10.84 12.0 9.13 12.0 8.15 8.0 5.56
7.0 4.82 7.0 7.26 7.0 6.42 8.0 7.91
5.0 5.68 5.0 4.74 5.0 5.73 8.0 6.89
2、实验环境
win10 + MATLABR2014b
3、实验代码
figure()
data= load('Anscombe.txt');%读入文件
x=0:2:20; % x 坐标从0 到20间隔为2
for i =1:2:8%从1开始到8,间隔为2[1,3,5,7]
subplot(2,2,(i+1)/2);%子图的编号
plot(data(:,i),data(:,(i+1)),'r.', 'markersize',20);%画图
hold on
p = polyfit(data(:,i),data(:,(i+1)), 1);%拟合
plot(x, polyval(p, x));%画拟合的图
axis([0 20 0 20])%坐标轴的范围
set(gca,'linewidth',2,'fontsize',10,'fontname','Times');
set(gca,'XTick',[0:2:20]) %改变x轴坐标间隔显示 这里间隔为2
set(gca,'YTick',[2:2:20]) %改变x轴坐标间隔显示 这里间隔为2
xlabel(['x',num2str((i+1)/2)],'Fontname', 'Times New Roman','FontSize',12);
ylabel(['y',num2str((i+1)/2)],'Fontname', 'Times New Roman','FontSize',12);
end