matlab建模文献,建模论文matlab附录-主要程序代码

主要程序代码

clear;

a=10;

b=0;

h=(a-b)/500;

w1=100;

w2=50;

w3=0.001;

for i=2:501

k11=h*((1-1*w3)*w1-(1-1*w3)*w1*w1/1000-0.05*w2);

k12=h*((1-1.5*w3)*w2-(1-1.5*w3)*w2*w2/(0.1*w1));

k13=-h*(0.00001*w1*w3+0.00002*w2*w3);

k21=h*((1-1*(w3+0.5*k13))*(w1+0.5*k11)-(1-1*(w3+0.5*k13))*(w1+0.5*k11)*(w1+0.5*k11)/1 000-0.05*(w2+0.5*k12));

k22=h*((1-1.5*(w3+0.5*k13))*(w2+0.5*k12)-(1-1.5*(w3+0.5*k13))*(w2+0.5*k12)*(w2+0.5*k1 2)/(0.1*(w1+0.5*k11)));

k23=-h*(0.00001*(w1+0.5*k11)*(w3+0.5*k13)+0.00002*(w2+0.5*k12)*(w3+0.5*k13));

k31=h*((1-1*(w3+0.5*k23))*(w1+0.5*k21)-(1-1*(w3+0.5*k23))*(w1+0.5*k21)*(w1+0.5*k21)/1 000-0.05*(w2+0.5*k22));

k32=h*((1-1.5*(w3+0.5*k23))*(w2+0.5*k22)-(1-1.5*(w3+0.5*k23))*(w2+0.5*k22)*(w2+0.5*k2 2)/(0.1*(w1+0.5*k21)));

k33=-h*(0.00001*(w1+0.5*k21)*(w3+0.5*k23)+0.00002*(w2+0.5*k22)*(w3+0.5*k23));

k41=h*((1-1*(w3+k33))*(w1+k31)-(1-1*(w3+k33))*(w1+k31)*(w1+k31)/1000-0.05*(w2+k32)); k42=h*((1-1.5*(w3+k33))*(w2+k32)-(1-1.5*(w3+k33))*(w2+k32)*(w2+k32)/(0.1*(w1+k31))); k43=-h*(0.00001*(w1+k31)*(w3+k33)+0.00002*(w2+k32)*(w3+k33));

w1=w1+(k11+2*k21+2*k31+k41)/6;

w2=w2+(k12+2*k22+2*k32+k42)/6;

w3=w3+(k13+2*k23+2*k33+k43)/6;

p1(i)=w1;

p2(i)=w2;

p3(i)=w3;

pp(i)=i;

end

>> plot(pp,p1)

>> plot(pp,p2)

>> plot(pp,p3)

>> plot(pp,p3)

>>