2013年东三省数模A题第一问（Matlab绘制折线图、饼状图）

　　最近参加了东三省的数模竞赛，怎么说呢==怎一个坑爹了得。不过最后还是挣扎的把论文交上去了。因为我在数模中主要负责写代码的部分，正赶上最近有点时间，所以整理一下数模中学到的东西，也算对自己知识的一个梳理。PS：只是一些基本的知识和一些以前大家就碰到过的问题，希望大家不要吐槽。

　　首先是背景：

　　　　比赛：2013年东三省数模

　　　　所选题：A题（http://www.madio.net/thread-183416-1-1.html）

　　　　此题分析一下，大体可以归为大数据分析和数理统计的范畴。

下面是这道题的主要部分：

请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取），并根据这些资料来讨论：

1.      如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；

2.      从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；

3.      能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？

（下载数据后一共有一百九十多个excel和word表格==坑爹啊）

　　第一题就很坑爹，一个看似很简单的题目，确是最坑爹的题目，思路很简单，就是按年份统计分类，问题是数据量大啊==怎么办。。发扬吃苦耐劳精神，干呗。

 

　　统计完数据后，根据总数据绘制折线图，根据污染物绘制饼状图。不多说，直接上代码，代码中有注释：

%这是折线图的代码
clear all; 

%三年的数据
h_metal = [0.004160888  0           0.002053388 0.009210526 0.001730104   0.029742234 0.001589825 0.005432099 0.003926702 0.001432254   0           0           0.000789266 0.005124451 0.000937647 0.00426916  0.00426916];
microbe = [0.016643551  0.008429119 0.016427105 0.039473684 0.029411765   0.043621943 0.018282989 0.007901235 0.028795812 0.017187052   0.002024291 0.002922078 0.016574586 0.005124451 0.005625879 0.018296402 0.018296402];
additive = [0.030513176 0.004597701 0.008213552 0.022368421 0.012687428   0.04494382  0.032591415 0.007901235 0.016361257 0.008307075   0.000674764 0.006168831 0.009471192 0.00329429  0.003750586 0.003252694 0.003252694];
others = [0.06518724    0.002298851 0           0.040789474 0.009803922   0.009473452 0.016693164 0.019259259 0.018979058 0.010025781   0 0 0 0 0 0 0];

%转化成千分比
h_metal = h_metal * 100;
microbe = microbe * 100;
additive = additive * 100;
others = others * 100;

%三年，一共16个数据，2010年6个，2011年5个， 2012年5个
x = 1 : 17;
x=x';
%plot(x,a,
%'k-o',   黑色 实线 空心圆
%'linewidth',2,' 线宽  
%markersize',4);  标记点大小
plot(x,h_metal,'r-s','linewidth',1.5,'markersize',2.5);
%在同一副图上绘制四条曲线
hold on 
plot(x,microbe,'k-s','linewidth',1.5,'markersize',2.5);
plot(x,additive,'b-s','linewidth',1.5,'markersize',2.5);
plot(x,others,'g-s','linewidth',1.5,'markersize',2.5);
hold off

axis ([1 17 -1 10]); %对坐标轴进行设置 ，这点是需要注意的，点的数量要特别注意
set (gca,'xtick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17]);
set (gca,'xticklabel',{'2010' '' '' '' '' '2011' '' '' '' '' '2012' ''  '' '' '' '' '2013 (年)'});
set (gca,'ytick',[-1 0 1 2 3 4 5 6 7 8 9 10 ]);
set (gca, 'yticklabel', {'' '0' '1' '2' '3' '4' '5' '6' '7' '8' '9' '(%)'  });

set(gca,'FontSize',12, 'FontName','标楷体');

ylabel('千分比','fontname','标楷体','fontweight','bold','fontsize',12); %坐标轴标题
xlabel('时间','fontname','标楷体','fontweight','bold','fontsize',12);

%text(1,58,'(a)','Fontweight','bold','fontsize',12)  ; %文本框标注 添加三年的文本框标注
box off;
% plot(x,ma,'k--','linewidth',2);
legend 重金属 微生物 添加剂 其他 ; %图例
%legend('boxoff');
 
set(legend,'fontname','标楷体');
set(legend,'fontweight','bold');

然后是绘制饼状图的代码：

%2010添加剂

Names={'米面制品:';'休闲零食:';'糕点点心:';'熟肉制品:';... 4 
    '动物性水产品:';'干货:';'调味料:';'饮料:';'酒:'; '餐饮用具:'};

data=[4 3 13 4 1 1 1 1 3 4];


%食品种类总数
sum = 35;

data = data / sum;
%画饼状图
h =pie(data);
%下面部分是从网上抄的，具体功能是在饼状图上在各百分比处加上文字
colormap jet

textObjs = findobj(h,'Type','text');

oldStr = get(textObjs,{'String'});

val = get(textObjs,{'Extent'});

oldExt = cat(1,val{:});

newStr = strcat(Names,oldStr);%在名称后面加上x%

set(textObjs,{'String'},newStr)

val1 = get(textObjs, {'Extent'});

newExt=cat(1, val1{:});

offset = sign(oldExt(:,1)).*(newExt(:,3)-oldExt(:,3))/2;%sign返回一个同型矩阵，如果值大于0，对应位置则为1；等于0，对应位置为0；小于0，对应位置为-1

pos = get(textObjs, {'Position'});%标识的坐标应该是(x,y,z)这种形式，pos是一个cell，其中每个cell是一个1*3的向量

textPos =  cat(1, pos{:});%将cell转换成一个n*3的矩阵，其中n为数据点数

textPos(:,1) = textPos(:,1)+offset;%取第一列，即每个标识的横坐标

textPos(3,2)=textPos(3,2)-0.1;

textPos(7,2)=textPos(7,2)-0.1;

set(textObjs,{'Position'},num2cell(textPos,[3,2]))%将textPos转成每个元素都是1*3向量的cell---%num2cell将数转为cell,取cell的值，应该用{}，如Names{1}

ok，今天就先整理这么多吧，滚床单去撒~~