参赛队号#1893
参考文献
[1]司守奎,孙玉箐，数学建模算法与应用,北京:国防工业出版社[M],193-261,2011.08 [2]张福顺，不同干扰下荒漠啮齿动物群落格局动态及敏感研究[D],内蒙古农业大学,26-30, 2008.05
[3]袁帅,荒漠区啮齿动物群落及其优势种对不同干扰的响应机制[D],内蒙古农业大学, 70-76,2013.06
[4]王利清,不同干扰方式下荒漠啮齿动物群落生物量与生境因子的关系[D],内蒙古农业大学,39-49,2006.05
[5]王利清,不同干扰下荒漠啮齿动物群落组成种数量和生物量格局[D],内蒙古农业大学,53-60,2009.04
[6]张福顺,不同干扰条件下荒漠啮齿动物种群和群落动态研究[D],内蒙古农业大学,58-65, 2011.04
[7]王利清,武晓东,不同干扰方式对荒漠植物与啮齿动物群落的作用[J],内蒙古大学学报,27(3), 47-52,2006.09
[8]刘保柱等,《MATLAB 7.0 从入门到精通》,北京:人民邮电出版社[M],19-46,2010.05 [9]乌云噶,查木哈,荒漠区啮齿动物群落与植物因子的冗余分析[J],草业科学,31(12), 2323-2332,2014.12
22
参赛队号#1893
附录1：
建议信
尊敬的政府部门你们好：
我们都知道，经济的发展离不开环境的良性循环。干旱区作为全球生态系统中的重要类型之一，也是目前全球开发较晚的区域之一，因此，积极开展和深化干旱区的生态学研究，对于该区域的经济发展与生态环境保护具有深远的理论意义与时间价值。
随着人为干扰不断加重，干旱区受干扰比较严重。对干旱区的健康状态的评估与该地区的植物生物量和动物生物量都有关系，干旱区研究物种主要植物有霸王、白刺、猫头刺、藏锦鸡儿等，优势动物为啮齿动物三趾跳鼠、子午沙鼠、小毛足鼠，对干旱区不同地方做取样处理，统计数据得到的干旱区植物动物数据统计，对该数据表进行数据处理分析，我们得到一下结论：
（1）过牧降低了草本植物生物量，有利于灌木植物生物量的提高，轮牧降低了啮齿动物的总丰富度；
（2）过度放牧减少了子午沙鼠和三趾跳鼠的丰富度，轮牧降低了小毛足鼠的数量； （3）影响过牧区、轮牧区啮齿动物群落多样性和动物个体数量的最大直接效应均来自草本；
（4）气候会影响草本植物、灌木植物的生物量，7月份左右适合过牧，10月份前后应该以轮牧的方式放牧。
因此针对以上得出的信息，对干旱区生态环境的管理给出以下相应的建议为： ?在过牧区适当减少放牧，减少人为干扰，使过牧区进入生态的恢复调养；
?在植物生长代谢旺盛的季节（6月、7月、8月）可以适当的增加放牧，提高干旱区植物生态资源的利用率；
?过牧区，可以适量增加草本植物的种植量，移植过多的灌木植物到轮牧区域内，可以缓解过牧区植物之间的竞争，提高草本植物对于空间、食物的利用率，促进过牧区生态环境的恢复；
?在植物生长代谢缓慢的季节（10月、11月、12月），应当对放牧区进行轮牧规划，依照湿润区、干旱区、半干旱区将放牧区进行合理规划，在使牧场资源运用效率最大化时不至于破坏牧区；
?轮牧区是一种保险的放牧方式，但是不能合理分配，往往不能满足广大的牧场资源需求，因此应当合理规划。
以上为我们团队经研究给出的观点和主张，由于我们团队条件、能力的限制，使得观点存在很大的局限，希望这次研究的微小成果，能给政府解决干旱区问题尽微薄之力。
23
参赛队号#1893
附录2：
1.典型相关分析MATLAB代码 clc
a=[1,2,5/2; 1/2,1,1; 2/5,1,1];
[x,y]=eig(a);eigenvalue=diag(y);lamda=eigenvalue(1); ci1=(lamda-3)/2;cr1=ci1/0.58 w1=x(:,1)/sum(x(:,1)) b1=[1,5,4,5,7; 1/5,1,1/2,3,3; 1/4,2,1,5,5; 1/5,1/3,1/5,1,3; 1/7,1/3,1/5,1/3,1];
[x,y]=eig(b1);eigenvalue=diag(y);lamda=eigenvalue(1); ci21=(lamda-5)/4;cr21=ci21/1.12 w21=x(:,1)/sum(x(:,1))
b2=[1,3,1/4,1/5,3; 1/3,1,1/4,1/6,1/4; 3,4,1,1/6,4; 5,6,6,1,4;
1/3,4,1/4,1/4,1];
[x,y]=eig(b2);eigenvalue=diag(y);lamda=eigenvalue(1); ci22=(lamda-5)/4;cr22=ci22/1.12 w22=x(:,1)/sum(x(:,1))
b3=[1,2,1/4,4,2; 1/2,1,1/5,4,1/2; 4,5,1,6,5;
1/4,1/4,1/6,1,1/4; 1/2,2,1/5,4,1];
[x,y]=eig(b3);eigenvalue=diag(y);lamda=eigenvalue(1); ci23=(lamda-5)/4;cr23=ci23/1.12 w23=x(:,1)/sum(x(:,1))
w_sum=[w21,w22,w23]*w1 ci=[ci21,ci22,ci23];
cr=ci*w1/sum(0.58*w1)
24
参赛队号#1893
2.时间序列MATLAB代码
function [yucezhi]=shijianxulie(data,cishu); [hang lie]=size(data); yucezhi=zeros(hang,lie); z=zeros(hang,cishu); for nt=1:hang x=data(nt,:); xx=x;
x=[0 diff(x)];
[NnN Nn]=size(x); %%%%%%用AR模型 for n=1:sqrt(Nn) x2=x((n+1):Nn); for j=1:(Nn-n) for ji=1:n
x1(j,ji)=x(n-ji+j); end end
dd=pinv((x1'*x1))*x1'*x2';
p=dd;%%%%%%%%p为行向量 for j=n+1:Nn for ii=1:n
ax(ii)=p(ii)*x(j-ii);%%%%拟合出的xt end
ct(j)=x(j)-sum(ax);%%%%%%%% end t=[];
t=ct(n+1:Nn).^2; b(n)=sum(t)/(Nn-n);
AIC(n)=log(b(n))+2*n/Nn;%%%%%%%% end
%%%%%%%%%%%%%%%%% pp=AIC(1); nn=1;
for i=1:sqrt(Nn) if pp>AIC(i) pp=AIC(i); nn=i; end end
for j=nn+1:Nn for ii=1:nn
ax(ii)=p(ii)*x(j-ii);
25
参赛队号#1893
end
ct(j)=x(j)-sum(ax); at0(nn+1:Nn)=ct(nn+1:Nn);
%%%%%%%%%%%%%估计参数aa bb jies=12; if jies>=Nn jies=Nn-1; end
for p=1:jies
for q=1:jies
hh=[p,q,nn]; L=max(hh);
X=[];e1=[];YY=[];Y=[];m=[];n=[];aabb=[];pa=[];pb=[]; for j=1:Nn-L for km=1:p
X(j,km)=x(L+j-km); end end
26
参赛队号#1893
27