shui3.m

05年全国大学生数学建模大赛a题长江水质预测代码

clear all;clc
a=xlsread('shui.xls');
m=mean(a);
a1=[];a2=[];a3=[];a4=[];a5=[];a6=[];a7=[];a8=[];a9=[];a10=[];a11=[];a12=[];a13=[];a14=[];a15=[];a16=[];a17=[];
for i=1:17:221;
    a1=[a1;a(i,:)];
end
for i=2:17:221;
    a2=[a2;a(i,:)];
end
for i=3:17:221;
    a3=[a3;a(i,:)];
end
for i=4:17:221;
    a4=[a4;a(i,:)];
end
for i=5:17:221;
    a5=[a5;a(i,:)];
end
for i=6:17:221;
    a6=[a6;a(i,:)];
end
for i=7:17:221;
    a7=[a7;a(i,:)];
end
for i=8:17:221;
    a8=[a8;a(i,:)];
end
for i=9:17:221;
    a9=[a9;a(i,:)];
end
for i=10:17:221;
    a10=[a10;a(i,:)];
end
for i=11:17:221;
    a11=[a11;a(i,:)];
end
for i=12:17:221;
    a12=[a12;a(i,:)];
end
for i=13:17:221;
    a13=[a13;a(i,:)];
end
for i=14:17:221;
    a14=[a14;a(i,:)];
end
for i=15:17:221;
    a15=[a15;a(i,:)];
end
for i=16:17:221;
    a16=[a16;a(i,:)];
end
for i=17:17:221;
    a17=[a17;a(i,:)];
end
b1=mean(a1);
b2=mean(a2);
b3=mean(a3);
b4=mean(a4);
b5=mean(a5);
b6=mean(a6);
b7=mean(a7);
b8=mean(a8);
b9=mean(a9);
b10=mean(a10);
b11=mean(a11);
b12=mean(a12);
b13=mean(a13);
b14=mean(a14);
b15=mean(a15);
b16=mean(a16);
b17=mean(a17);
b=[b1;b2;b3;b4;b5;b6;b7;b8;b9;b10;b11;b12;b13;b14;b15;b16;b17];
gao=b(:,3);
dan=b(:,4);
gao=gao';
dan=dan';
for i=1:17;
    for j=1:17;
       ju1(i,j)=gao(i)/gao(j);
       ju2(i,j)=dan(i)/dan(j);
    end
end
he1=sum(ju1(:,1));
guiyi1=ju1(:,1)/he1
he2=sum(ju2(:,1));
guiyi2=ju2(:,1)/he2
guiyi=[guiyi1,guiyi2];
hang=[0.5;0.5];
paixu=guiyi*hang
xlswrite('ju1.xls',ju1)
xlswrite('ju2.xls',ju2)