// t y
2000 1248 2001 1223 2002 1134 2003 1465 2004 1430 2005 1452 2006 1466 2007 1578
对数据进行灰色预测的源程序 1.灰色预测的主程序：
function [X0 detaeq gama C p]=di(x0) n=size(x0',1); x1=zeros(1,n); z1=zeros(1,n); X0=zeros(1,n); deta=zeros(1,n); if any(x0<0)
min1=min(x0); else min1=0; end
for i=1:n
x0(i)=x0(i)+2*abs(min1); end
for i=1:n
x1(i)=sum(x0(1:i)); end
for i=2:n
z1(i)=(x1(i)+x1(i-1))/2; end
B=[-z1(2:n)' ones(n-1,1)]; Y=x0(2:n)';
a=inv(B'*B)*B'*Y;
X0(1)=(x0(1)-a(2)/a(1))*exp(-a(1)*0)+a(2)/a(1)-2*abs(min1); for t=1:n+9
X1=(x0(1)-a(2)/a(1))*exp(-a(1)*t)+a(2)/a(1); X12=(x0(1)-a(2)/a(1))*exp(-a(1)*(t-1))+a(2)/a(1); X0(t+1)=X1-X12-2*abs(min1); end
%求deta平均值。 for i=1:n
deta(i)=abs(X0(i)-x0(i))/x0(i); end
detaeq=mean(deta);
21
p=(1-detaeq)*100/100; [gama C p]=f3(x0,X0); 求关联度，C，p的程序
function [gama C p]=f3(x0,X0) %以下求gama n=size(x0',1); gama=0; k=1;
while k<=n-1
gama=gama+1/(n-1)*1/(1+abs((x0(k+1)-x0(k))/x0(k+1)-(X0(k+1)-X0(k))/X0(k+1))); k=k+1; end
%以下求C temp=0; j=1;
while j<=n
temp=temp+1/n*x0(j); j=j+1; end %求s1 s1=0; k=1;
while k<=n
s1=s1+1/n*(x0(k)-temp)^2; k=k+1; end %求s2
ebq=0;?q相对误差的均值。 k=1;
while k<=n
ebq=ebq+1/n*(x0(k)-X0(k));?q相对误差的均值。 k=k+1; end s2=0; k=1;
while k<=n
s2=s2+1/n*(x0(k)-X0(k)-ebq)^2; k=k+1; end
C=(s2/s1)^(1/2); %以下求p r=0;w=0; for k=1:n
if abs((x0(k)-X0(k)-ebq))<0.6745*s1
22
r=r+1; else
w=w+1; end end
p=r/(r+w);
第二问源程序
模糊评判模型源程序
function [r,w]=pingjia(c,s) for i=1:1:4
if c(i) if c(i) r(i,1)=(-c(i)+s(i,2))/(-s(i,1)+s(i,2)); r(i,2)=1-r(i,1); r(i,3:5)=0; end if c(i) r(i,2)=(-c(i)+s(i,3))/(-s(i,2)+s(i,3)); r(i,3)=1-r(i,2); r(i,4:5)=0;r(i,1)=0; end if c(i) r(i,3)=(-c(i)+s(i,4))/(-s(i,3)+s(i,4)); r(i,4)=1-r(i,3); r(i,5)=0;r(i,1:2)=0; end if c(i)>=s(i,4) r(i,5)=1; r(i,1:4)=0; end se=mean(s(i,1:4)); I(i)=c(i)/se; end for i=1:1:4 w(i)=I(i)/sum(I); end end 求解程序：其中p为我们要求的结果。 s=[50 140 395 696.5;50 97.5 146 206;3400 4750 5710 6721;171.5 200 299 314]; load g.txt %g为评判指标值所构成的矩阵。 for i=1:1:size(g,1) 23 c=g(i,:); [r,w]=pingjia(c,s); b=w*r; p(i)=b*[1 2 3 4 5]'; end 第三问源程序 max=(0.2727*x-1071)/(1994.7935*A) ; 1.227*x-7282.3<=6556.12*(1+u); A<=2.2*10000*b; A=13750+0.2*x; b=1; !u可调，取0.1—0.6; u=0.5; 表23 各界世博会的规模数据 年份 举办国城市 名称 类型 主题 场馆面积（公顷） 10.4 参展国家和组织 25 参观人数（万人） 604 天 数 1851 英国伦敦 美国费城 法国巴黎 澳大利亚墨尔本 法国巴黎 伦敦万国工业产品大博览会 美国独立百年博览会 第3届巴黎世界博览会 万国工农业，制造与艺术博览会 世界博览会（1889） 综合 万国工业 190 1876 综合 庆祝美国百年独立 农业 115 35 800 180 1878 综合 75 36 1616 170 1880 综合 万国工农业 25 33 1200 210 1889 综合 法国大革命百年，艾菲尔铁塔落成 哥伦布发现新大陆四百年 世纪回顾 96 35 2512 182 1893 美国芝加哥 法国巴黎 芝加哥哥伦布纪念博览会 第5届巴黎世界博览会 综合 290 19 2700 183 1900 综合 120 58 5000 210 1904 美国圣圣路易斯百路易斯 周年纪念博览会 综合 该市成立百年 500 60 1969 185 24 1915 美国旧金山 美国芝加哥 美国纽约 比利时布鲁塞尔 加拿大蒙特利尔 日本大坂 旧金山巴拿马太平洋博览会 芝加哥万国博览会 纽约世界博览会 布鲁塞尔世界博览会 加拿大世界博览会 日本万国博览会 综合 庆祝巴拿马运河通航 进步的世纪 建设明天的世界 科学 254 32 1883 288 1933 1939 1958 综合 综合 综合 170 500 200 21 25 42 2257 4500 4150 170 340 186 1967 综合 人类与世界 287 62 5031 185 1970 1992 2000 2005 2010 综合 综合 综合 综合 人类的进步与和谐 发现的时代 人类 自然的睿智 城市，让生活更美好 865 215 160 173 528 81 114 172 120 240 6422 4100 1800 2200 7500（预计） 183 176 153 185 184 西班牙塞维亚世界塞维亚 博览会 德国汉诺威 日本爱知 中国上海 汉诺威世界博览会 爱知地球博览会 上海世博会 综合 25