ÔõôÓÃmatlabÇó½âÕâ¸öƫ΢·Ö·½³Ì×飬лл´ó¼ÒÖ¸µã£¿
ÕâÊÇÎÒдµÃÁ½¸ö³ÌÐò£¬µ«ÊÇÔËÐеÄʱºò×ÜÊdzö´í£¬´ó¼ÒÒ²°ïæ¿´Ò»ÏÂ
³ÌÐòÒ»£º
function dx=weifen_fun(t,x)
Dx-Dy*(0.6/(1+0.4*t))=t^2*69.888/(1+0.4*t)+t*2096.64/(1+0.4*t)+4914/(1+0.4*t)-x*(4.4+0.16*t)/(1+0.4*t);
Dy-y^2*1.8*exp(38)*exp(1-1/x)+y*1.8*exp(38)*exp(1-1/x)*(t+1)=1.8*exp(38)*exp(1-1/x)*t;
clear,clc;
t12=linspace[0 1];
x0=0
[t,x]=ode45(@weifen_fun,tspan,0);
plot(x,t);
grid on
³ÌÐò¶þ£º
function dz=weifeng_fun(t,z)
dz(2,1)=dz(1,1)*(0.6/(1+0.4*t))+t^2*69.888/(1+0.4*t)+t*2096.64/(1+0.4*t)+4914/(1+0.4*t)-z(2)*(4.4+0.16*t)/(1+0.4*t);
dz(1,1)=z(1)^2*1.8*exp(38)*exp(1-1/z(2))-z(1)*1.8*exp(38)*exp(1-1/z(2))*(t+1)+1.8*exp(38)*exp(1-1/z(2))*t;
clear,clc;
t12=linspace(0,0.2,1);
[t,z]=ode45('weifeng_fun',t12,[2,0]);
plot(t,z(:,1))
grid on
![ÔõôÓÃmatlabÇó½âÕâ¸öƫ΢·Ö·½³Ì×飬лл´ó¼ÒÖ¸µã£¿]()
΢·Ö·½³Ì×é.png |