《2011全国数学建模大赛A题相关matlab程序.pdf》由会员分享，可在线阅读，更多相关《2011全国数学建模大赛A题相关matlab程序.pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、20112011 全国数学建模大赛全国数学建模大赛 A A 题相关题相关 matlabmatlab 程序程序(2011-09-16 16:42:56)这是我在今年全国大学生数学建模大赛中解答 A 题的相关 MATLAB 程序（供大家参考）,同时晒两张相关图片：%绘单个元素空间分布图程序clear all;load shuju.txta=shuju(1:end,8);%b=a/69;%load data.txtx=data(1:end,1);y=data(1:end,2);z=data(1:end,3);scatter3(x,y,z,50,b,filled)xlabel(X)ylabel(Y)z
2、label(海拔)%算各区污染平均指数程序clear all;load shuju.txtfor i=1:8a=shuju(1:end,i); %i 随数据文本列数决定p=3.6 130 31 13.2 35 12.3 31 69b(:,i)=a/p(i); %i 随数据文本列数决定endload data.txtx=data(1:end,4);m=1;sum=0;for i=1:319if x(i)=5cb(m)=b(i,8);sum=sum+cb(m);m=m+1;endendav=sum/(m-1)%层次分析法一致性检验程序clearclose allclcticdisp(please
3、input Matrix A)A=input(A=n);%A 矩阵为需要分析的 n 阶方阵！n,n=size(A);x=ones(n,100);y=ones(n,100);m=zeros(1,100);m(1)=max(x(:,1);y(:,1)=x(:,1);x(:,2)=A*y(:,1);m(2)=max(x(:,2);y(:,2)=x(:,2)/m(1);p=0.0001;i=2;k=abs(m(2)-m(1);while kp;i=i+1;x(:,i)=A*y(:,i-1);m(i)=max(x(:,i);y(:,i)=x(:,i)/m(i);k=abs(m(i)-m(i-1);end
4、a=sum(y(:,i);w=y(:,i)/a;t=m(i);disp(权重向量 w=)disp(w)%fprintf(权重向量 w=%fn,w);fprintf(最大特征根 t=n %fn,t);%以下为该方阵的一致性检验！CI=(t-n)/(n-1);RI=0 0 .58 .90 1.12 1.24 1.32 1.41 1.45 1.49 1.51;CR=CI/RI(n);if CR0.1;disp(CI=)disp(CI)disp(RI=)disp(RI(n)disp(CR=)disp(CR)fprintf(所以n CR0.1nn);disp(因此该方阵可以通过一致性检验！)else d
5、isp(因此该方阵不可以通过一致性检验！)end%总权重计算程序clear all;clcload chu.txta=chu(1,:);%b=chu(2,:);%c=chu(3,:);%d=chu(4,:);%e=chu(5,:);%k=chu(5,:);sum=0;for i=1:4 %根据总矩阵列数决定 4 所在位置的数据l(i)=a(i)*b(i);sum=sum+l(i);endsum1=0;for i=1:4l(i)=a(i)*c(i);sum1=sum1+l(i);endsum2=0;for i=1:4l(i)=a(i)*d(i);sum2=sum+l(i);endsum3=0;for i=1:4l(i)=a(i)*e(i);sum3=sum+l(i);endsum4=0;for i=1:4%3 这一位的数由阵列列数决定l(i)=a(i)*k(i);sum4=sum+l(i);endsumsum1sum2sum3sum4