water_salt.m

来自「matlab编写的计算机模拟程序 很方便使用 注释很清楚」· M 代码 · 共 37 行

h=1;
s0=2;
w0=2000;
r0=s0/w0;
s(1)=s0+0.5*6*h-4*h*r0;
w(1)=w0+2*h;
r(1)=s(1)/w(1);
t(1)=h;
y(1)=(2000000+3000000*h+3000*h^2+h^3)/(1000+h)^2;
for i=2:200
    t(i)=i*h;
    s(i)=s(i-1)+0.5*6*h-4*h*r(i-1);
w(i)=w(i-1)+2*h;
r(i)=s(i)/w(i);
    y(i)=(2000000+3000000*t(i)+3000*t(i)^2+t(i)^3)/(1000+t(i))^2;

m=floor(i/10);
if i/10-m<0.001
    tm(m)=m;
    sm(m)=s(i);
    rm(m)=r(i);
end

if r(i)>0.2
    if r(i)<0.2005
        t02=i*h;
        r02=r(i);
    end
end
end
[t02,r02]
x=[10*tm',sm',rm'];
plot(t,s,'b',t,y,'r')


XLSWRITE('C:\Documents and Settings\方如节\桌面\a.xls',x);