恒温箱数学建模(4)

ir(i))]';

K3(:,1)=[h*func1(y(1,i)+K2(1,1)/2,y(2,i)+K2(2,1)/2,Tair(i)),h*func3(y(1,i)+K2(1,1)/2,y(2,i)+K2(2,1)/2,Tair(i))]';

K4(:,1)=[h*func1(y(1,i)+K3(1,1),y(2,i)+K3(2,1),Tair(i)),h*func3(y(1,i)+K3(1,1),y(2,i)+K3(2,1),Tair(i))]'; y(:,i+1)=y(:,i)+1/6*(K1+2*K2+2*K3+K4);%加热循环 i=i+1; end

while (i<=n)&&(y(1,i)>67.8)

K1(:,1)=[h*func2(y(1,i),y(2,i),Tair(i)),h*func3(y(1,i),y(2,i),Tair(i))]';

K2(:,1)=[h*func2(y(1,i)+K1(1,1)/2,y(2,i)+K1(2,1)/2,Tair(i)),h*func3(y(1,i)+K1(1,1)/2,y(2,i)+K1(2,1)/2,Tair(i))]';

K3(:,1)=[h*func2(y(1,i)+K2(1,1)/2,y(2,i)+K2(2,1)/2,Tair(i)),h*func3(y(1,i)+K2(1,1)/2,y(2,i)+K2(2,1)/2,Tair(i))]';

K4(:,1)=[h*func2(y(1,i)+K3(1,1),y(2,i)+K3(2,1),Tair(i)),h*func3(y(1,i)+K3(1,1),y(2,i)+K3(2,1),Tair(i))]'; y(:,i+1)=y(:,i)+1/6*(K1+2*K2+2*K3+K4);%不加热循环 i=i+1; end end

hwxwd=y(1,:); grcwd=y(2,:); end

微分方程1：

functiondelta_hwx= func1( Tbox,Tgere,Tair ) %加热方程时的恒温箱温度微分方程

delta_hwx=1/63.9245*(22.4495-1.0366*Tbox+Tgere+0.0366*Tair); end

微分方程2：

function delta_hwx1= func2( Tbox,Tgere,Tair ) %不加热方程时的恒温箱温度微分方程

delta_hwx1=1/63.9245*(-1.0366*Tbox+Tgere+0.0366*Tair); end

微分方程3：

functiondelta_grc= func3( Tbox,Tgere,Tair ) %隔热层温度微分方程

delta_grc=1/122.5048*(Tbox-1.4669*Tgere+0.4669*Tair); end