|
|
[求助]
怎么用matlab求解这个偏微分方程组,谢谢大家指点?
怎么用matlab求解这个偏微分方程组,谢谢大家指点?
这是我写得两个程序,但是运行的时候总是出错,大家也帮忙看一下
程序一:
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 |
» 猜你喜欢
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
各位道友,我要去昆明玩几天,回来见。
已经有8人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有7人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
是这周出结果还是下周出结果?
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
|