huidu.m

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

clear all;clc
x0=[3
1.9
2.6
1.7
6.2
4.4
5.5
3.2
5.8
5.9];
x0=x0'+100;
for j=10:19;
for i=2:length(x0);
    jian(i-1)=x0(i-1)/x0(i);
end
jian;
jian1=exp(-2/(j+1));
jian2=exp(2/(j+1));
x1=cumsum(x0);
for i=2:length(x0);
    z(i-1)=0.5*(x1(i)+x1(i-1));
end
z;
duizhao=x0;
duizhao(1)=[];
aa=polyfit(z,duizhao,1);
a=-aa(1);
b=aa(2);
x=(x0(1)-b/a)*(exp(-a*i)-exp(-a*(i-1)))
x0=[x0,x];
end
x0=x0-100