x6=[1.09 1.33 1.45 1.82 1.67]; x7=[0.12 0.17 0.39 0.70 1.18]; y=[2007 2008 2009 2010 2011];
plot(y,x1,y,x2,y,x3,y,x4,y,x5,y,x6,y,x7); gtext('一岗位'); gtext('二岗位'); gtext('三岗位'); gtext('四岗位'); gtext('五岗位'); gtext('六岗位'); gtext('七岗位');
title('各岗位级别与每年平均招生人数统计图')
2图
x1=[0.89 1.5 1.39 1.06 1.22]; x2=[0.769230769 1.538461538
0.923076923];
0.9375 1.125 1 0.9375];
0.569230769
0.953846154];
0.875 0.8125 0.5625 0.4375];
0.848484848 0.546448087
0.6060606060.469945355
x3=[0.4375 x5=[0.8125
1.153846154
1.076923077
x4=[0.353846154 0.4 0.615384615 x6=[0.636363636 0.696969697
0.242424242]; 0.459016393];
x7=[0.551912568 0.530054645 y=[2007 2008 2009 2010 2011];
plot(y,x1,y,x2,y,x3,y,x4,y,x5,y,x6,y,x7); gtext('一岗位'); gtext('二岗位'); gtext('三岗位'); gtext('四岗位'); gtext('五岗位'); gtext('六岗位'); gtext('七岗位');
title('各岗位级别与每年平均发表论文数统计图')
3图
x1=[0.944444444 0.555555556
0.222222222];
0.461538462
0.307692308
0.153846154];
x2=[0 0.153846154
0.611111111
0.222222222
x3=[0.125 0.125 0.375 0.4375 0.375]; x4=[0.046153846 0.123076923
0.092307692];
0.454545455
24
0.123076923 0.276923077
x5=[0.75 0.9375 0.6875 0.5625 0.0625]; x6=[0.424242424 0.484848485
0.575757576
0.121212121];
x7=[0.278688525 0.387978142 0.404371585
0.267759563
0.163934426];
y=[2007 2008 2009 2010 2011];
plot(y,x1,y,x2,y,x3,y,x4,y,x5,y,x6,y,x7); gtext('一岗位'); gtext('二岗位'); gtext('三岗位'); gtext('四岗位'); gtext('五岗位'); gtext('六岗位'); gtext('七岗位');
title('各岗位级别与每年平均申请专利数统计图')
4图
x1=[0.055555556 0.055555556 0.055555556 0.055555556 x2=[0 0 0 0 0]; x3=[0 0 0 0.125 0]; x4=[0.015384615 0.015384615 0.030769231 0.046153846 x5=[0 0 0 0.0625 0]; x6=[0 0 0 0 0]; x7=[0 0 0 0.049180328
0.016393443];
y=[2007 2008 2009 2010 2011];
plot(y,x1,y,x2,y,x3,y,x4,y,x5,y,x6,y,x7); gtext('一岗位'); gtext('二岗位'); gtext('三岗位'); gtext('四岗位'); gtext('五岗位'); gtext('六岗位'); gtext('七岗位');
title('各岗位级别与每年平均获奖数统计图')
5图
x=[30.80 41.98 71.67 25.42 32.31 97.06 33.01]; y=[1 2 3 4 5 6 7]; bar(y,x)
title('各岗位级别平均每位导师5年纵向经费统计图')
6图
x=[46.62 26.61 18.44 4.94 2.83 7.81 3.95]; y=[1 2 3 4 5 6 7]; bar(y,x)
title('各岗位级别平均每位导师5年横向经费统计图')
25
0];
0];
7图
x=[10/18 6/13 15/16 25/65 3/16 8/33 15/183]; y=[1 2 3 4 5 6 7]; bar(y,x)
title('各岗位级别平均每位导师5年获奖论文数统计图')
附录3：
function GM11(x) x=[313 392 453 598 608]; n=length(x); x0=x; for i=2:n
x0(i)=x(i)+x0(i-1); end x1=x; for i=2:n
x1(i)=0.5*x0(i)+0.5*x0(i-1); end
y=ones(n-1,1); for i=2:n
y(i-1)=x(i); end x2=y; for i=2:n
x2(i-1)=x1(i); end
X=[-x2 ones(n-1,1)]; A=X\\y;
a=A(1);b=A(2); x0_L=x; for i=2:n
x0_L(i)=(x(1)-b/a)*exp(-a*(i-1))*(1-exp(a)); end
x0e=abs(x0_L-x); mx=mean(x); sx=0; for i=1:n
sx=sx+(x(i)-mx)^2; end sx=sx/n; me=mean(x0e); se=0;
26
for i=1:n
se=se+(x0e(i)-me)^2; end se=se/n; c=se/sx; k=0; for i=1:n
if abs(x0e(i)-me)<=0.6745*sx k=k+1; end end p=k/n;
fprintf('*****灰色预测GM11模型******\\n'); fprintf('序号 原始数列 预测值 残差 \\n'); for i=1:n
fprintf('- %6.2f %6.2f %4.3f\\n',i,x(i),x0_L(i),x0e(i)); end
fprintf('优度检验值 \\n'); fprintf(' p c \\n');
fprintf('%4.3f %4.3f\\n',p,c); fprintf('优度等级为:'); if p>0.95&&c<0.35 fprintf('好');
elseif p>0.80&&c<0.45 fprintf('合格');
elseif p>0.70&&c<0.50 fprintf('勉强合格'); else
fprintf('合格'); end
if nargin==1 return end
27