| 查看: 196 | 回复: 0 | |||
[求助]
matlab求解常微分问题
|
 问题如图： 以下是我的malab源码 function exm21 clear all,clc y0=[625*ones(1,11),0.0017*ones(1,11)] [t,y]=ode23(@Eqs,[0 500],y0); plot(t,y(:,[1：11])) %T随时间的关系 function y=Eqs(t,y0) N=10; dz=0.1; T=y0(1:N+1); PA=y0(N+2:end); dT(1)=0; dPA(1)=0; for i=2:N+1 dT(i)=-9.872e-4/dz*(T(i)-T(i-1))+3.0232e11*PA(i)*exp(-13636/T(i))-1.164e-2*(T(i)-625); dPA(i)=-3.548e-3/dz*T(i)*(PA(i)-PA(i-1))-PA(i)/T(i)*(T(i)-T(i-1))-2.543e6*T(i)*PA(i)*exp(-16363/T(i))+PA(i)/T(i)*dT(i); end y=[dT,dPA]'; |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab解微分方程
已经有10人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】如何求解下列偏微分方程
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人