shui2.m

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

clear all;clc
a=xlsread('2.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:476;
    a1=[a1;a(i,:)];
end
for i=2:17:476;
    a2=[a2;a(i,:)];
end
for i=3:17:476;
    a3=[a3;a(i,:)];
end
for i=4:17:476;
    a4=[a4;a(i,:)];
end
for i=5:17:476;
    a5=[a5;a(i,:)];
end
for i=6:17:476;
    a6=[a6;a(i,:)];
end
for i=7:17:476;
    a7=[a7;a(i,:)];
end
for i=8:17:476;
    a8=[a8;a(i,:)];
end
for i=9:17:476;
    a9=[a9;a(i,:)];
end
for i=10:17:476;
    a10=[a10;a(i,:)];
end
for i=11:17:476;
    a11=[a11;a(i,:)];
end
for i=12:17:476;
    a12=[a12;a(i,:)];
end
for i=13:17:476;
    a13=[a13;a(i,:)];
end
for i=14:17:476;
    a14=[a14;a(i,:)];
end
for i=15:17:476;
    a15=[a15;a(i,:)];
end
for i=16:17:476;
    a16=[a16;a(i,:)];
end
for i=17:17:476;
    a17=[a17;a(i,:)];
end
ku1=a1([1:5,13:17,25:28],:);
ping1=a1([6:7,18,19],:);
feng1=a1([8:12,20:24],:);
ku2=a2([1:5,13:17,25:28],:);
ping2=a2([6:7,18,19],:);
feng2=a2([8:12,20:24],:);
ku3=a3([1:5,13:17,25:28],:);
ping3=a3([6:7,18,19],:);
feng3=a3([8:12,20:24],:);
ku4=a4([1:5,13:17,25:28],:);
ping4=a4([6:7,18,19],:);
feng4=a4([8:12,20:24],:);
ku5=a5([1:5,13:17,25:28],:);
ping5=a5([6:7,18,19],:);
feng5=a5([8:12,20:24],:);
ku6=a6([1:5,13:17,25:28],:);
ping6=a6([6:7,18,19],:);
feng6=a6([8:12,20:24],:);
ku7=a7([1:5,13:17,25:28],:);
ping7=a7([6:7,18,19],:);
feng7=a7([8:12,20:24],:);
ku8=a8([1:5,13:17,25:28],:);
ping8=a8([6:7,18,19],:);
feng8=a8([8:12,20:24],:);
ku9=a9([1:5,13:17,25:28],:);
ping9=a9([6:7,18,19],:);
feng9=a9([8:12,20:24],:);
ku10=a10([1:5,13:17,25:28],:);
ping10=a10([6:7,18,19],:);
feng10=a10([8:12,20:24],:);
ku11=a11([1:5,13:17,25:28],:);
ping11=a11([6:7,18,19],:);
feng11=a11([8:12,20:24],:);
ku12=a12([1:5,13:17,25:28],:);
ping12=a12([6:7,18,19],:);
feng12=a12([8:12,20:24],:);
ku13=a13([1:5,13:17,25:28],:);
ping13=a13([6:7,18,19],:);
feng13=a13([8:12,20:24],:);
ku14=a14([1:5,13:17,25:28],:);
ping14=a14([6:7,18,19],:);
feng14=a14([8:12,20:24],:);
ku15=a15([1:5,13:17,25:28],:);
ping15=a15([6:7,18,19],:);
feng15=a15([8:12,20:24],:);
ku16=a16([1:5,13:17,25:28],:);
ping16=a16([6:7,18,19],:);
feng16=a16([8:12,20:24],:);
ku17=a17([1:5,13:17,25:28],:);
ping17=a17([6:7,18,19],:);
feng17=a17([8:12,20:24],:);

b1=mean(ku1);
b2=mean(ku2);
b3=mean(ku3);
b4=mean(ku4);
b5=mean(ku5);
b6=mean(ku6);
b7=mean(ku7);
b8=mean(ku8);
b9=mean(ku9);
b10=mean(ku10);
b11=mean(ku11);
b12=mean(ku12);
b13=mean(ku13);
b14=mean(ku14);
b15=mean(ku15);
b16=mean(ku16);
b17=mean(ku17);
c1=mean(ping1);
c2=mean(ping2);
c3=mean(ping3);
c4=mean(ping4);
c5=mean(ping5);
c6=mean(ping6);
c7=mean(ping7);
c8=mean(ping8);
c9=mean(ping9);
c10=mean(ping10);
c11=mean(ping11);
c12=mean(ping12);
c13=mean(ping13);
c14=mean(ping14);
c15=mean(ping15);
c16=mean(ping16);
c17=mean(ping17);
d1=mean(feng1);
d2=mean(feng2);
d3=mean(feng3);
d4=mean(feng4);
d5=mean(feng5);
d6=mean(feng6);
d7=mean(feng7);
d8=mean(feng8);
d9=mean(feng9);
d10=mean(feng10);
d11=mean(feng11);
d12=mean(feng12);
d13=mean(feng13);
d14=mean(feng14);
d15=mean(feng15);
d16=mean(feng16);
d17=mean(feng17);
b=[b1;b2;b3;b4;b5;b6;b7;b8;b9;b10;b11;b12;b13;b14;b15;b16;b17;c1;c2;c3;c4;c5;c6;c7;c8;c9;c10;c11;c12;c13;c14;c15;c16;c17;d1;d2;d3;d4;d5;d6;d7;d8;d9;d10;d11;d12;d13;d14;d15;d16;d17];
b(:,1)=[];
z=[];
biao=[];
for i=1:51;
    p=b(i,:);
    if p(1)>=7.5
       biao(1)=1;
    elseif p(1)>=6&p(2)<7.5
        biao(1)=2;
    elseif p(1)>=5&p(1)<6
        biao(1)=3;
    elseif p(1)>=3&p(1)<5
        biao(1)=4;
    elseif p(1)>=2&p(1)<3
        biao(1)=5;
    elseif p(1)>=0&p(1)<=2
        biao(1)=6;
    end
    if p(2)<=2
        biao(2)=1;
    elseif p(2)<=4&p(2)>2
        biao(2)=1;
    elseif p(2)<=6&p(2)>4
        biao(2)=3;
    elseif p(2)<=10&p(2)>6
        biao(2)=4;
    elseif p(2)<=15&p(2)>10
        biao(2)=5;
    elseif p(2)>15
        biao(2)=6;
    end
    if p(3)<=0.15
        biao(3)=1;
    elseif p(3)<=0.5&p(3)>0.15
        biao(3)=2;
    elseif p(3)<=1.0&p(3)>0.5
        biao(3)=3;
    elseif p(3)<=1.5&p(2)>1
        biao(3)=4;
    elseif p(3)<=2&p(2)>1
        biao(3)=5;
    elseif p(3)>2
        biao(3)=6;
    end
    biao=max(biao);
     z=[z;biao];
end
z1=z(1:17,:);
z2=z(18:34,:);
z3=z(35:51,:);
z1=z1'
z2=z2'
z3=z3'
xlswrite('z1',z1)
xlswrite('z2',z2)
xlswrite('z3',z3)