matlab中ode15s函数的用法,求解n个常微分方程构成的方程组，如何使用ode15s

深藏阁楼爱情的钟 2023-01-14 10:54 1阅读 0赞

我的方程组涉及的方程重复性较高，而且数量较多，故采用i和for循环来编写方程，不知道哪里有问题，下面是我的代码

function dy = ad\_mass (t,y)

N=100;%间隔数

e=0.44;%床孔隙率

Dl=11;%轴向弥散系数，cm^2/s

us=40.7;%表观速度,cm/s

u=us/e;%间隙速度,cm/s

leng=25.5;%吸附剂填充长度,cm

z=leng/N;%间隔长度,cm

R=8.314;%气体常数，kPa\*cm^3/mmol/K

P=101.325;%气体压力，kPa

T=303;%温度，K

mdp=0.72;%颗粒物密度，g/cm^3

c0=0.003146;%初始浓度，mol/mol

mT=2.889\*exp(79.11/T);%langmuir常数，mmol/g

bT=5.429E-7\*exp(5292/T);%langmuir常数，1/kPa

ye=zeros(N,1);%平衡吸附量定义，mmol/g

for i=1:N+1

ye(i)=mT\*bT\*P\*y(i)/(1+bT\*P\*y(i));

if i==1

dy(i)=c0;

dy(i+N+1)=k\*(ye(i)-y(i+N+1));

else

if i==N+1

dy(i)=dy(i-1);

dy(i+N+1)=k\*(ye(i)-y(i+N+1));

else

dy(i)=Dl\*(y(i+1)-2\*y(i)+y(i-1))/z^2-u\*y((i+1)-y(i))/z-R\*T/P\*(1-e)/e\*mdp\*Dx(i+N+1);

dy(i+N+1)=k\*(ye(i)-y(i+N+1));

end

end

end

end

调用函数为

function sol= sy

y0=zeros;

sol=ode15s(@ad\_mass,\[0,3600\],y0);

出现以下错误

Error using feval

Error: File: ad\_mass.m Line: 22 Column: 5

"i" previously appeared to be used as a function or command, conflicting with its use here as

the name of a variable.

A possible cause of this error is that you forgot to initialize the variable, or you have

initialized it implicitly using load or eval.

Error in odearguments (line 87)

f0 = feval(ode,t0,y0,args\{:\});   % ODE15I sets args\{1\} to yp0.

Error in ode15s (line 148)

\[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in sy (line 3)

sol=ode15s(@ad\_mass,\[0,3600\],y0);

发表评论取消回复

表情：

评论列表（有 0 条评论，1人围观）

还没有评论，来说两句吧...

相关阅读

相关求解非刚性常微分方程（ODE）的方法

求解非刚性常微分方程（ODE）的方法非刚性常微分方程是指其中的参数或变量可以随时间变化的常微分方程。在Matlab中，我们可以使用ode45函数来求解非刚性ODE。ode4

秒速五厘米/ 2024年03月04日 03:50/ 0 赞/ 96 阅读

相关 matlab 多元多次函数,fslove - Matlab求解多元多次方程组

fslove - Matlab求解多元多次方程组简介：之前看到网上的一些资料良莠不齐，各种转载之类的，根本无法解决实际问题，所以我打算把自己的学到的总结一下，以实例出发讲

男娘i/ 2023年01月19日 01:56/ 0 赞/ 36 阅读

相关 matlab求解常微分方程,matlab 求解常微分方程式

MATLAB解常微分方程式的语法是dsolve('equation','condition')，其中equation代表常微分方程式即y'=g(x,y)，且须以Dy代表一阶微分

Bertha 。/ 2023年01月16日 06:20/ 0 赞/ 182 阅读

相关 matlab中如何画隐函数方程组,matlab隐函数方程组

用matlab 求解常微分方程在 MATLAB 中,由函数 dsolve()解决常微分方程(组)的求解问题,其具体格式如下: r = dsolve('eq1,eq2,...'

Dear 丶/ 2023年01月13日 10:29/ 0 赞/ 161 阅读

相关 matlab求解极限问题（limit函数的用法）

本文介绍利用MATLAB求解函数或序列的极限问题，顺便介绍`limit`函数的用法。内容主要包括单变量函数的极限和多变量函数的极限。目录单变

迈不过友情╰/ 2022年11月25日 04:12/ 0 赞/ 267 阅读

相关 java常微分方程数值解_如何求解R中具有时间相关参数的常微分方程？

我试图复制一个现有的数学模型，以获得一个独立研究项目的"control"数据 . Link到文章 . 我想用deSolve包解决R中五个常微分方程的系统 . 但是，大多数参数

偏执的太偏执、/ 2022年11月06日 12:51/ 0 赞/ 219 阅读

相关 MATLAB求解微分方程

求微分方程的解自牛顿发明微积分以来，微分方程在描述事物运动规律上已发挥了重要的作用。实际应用问题通过数学建模所得到的方程，绝大多数是微分方程。由于实际应用的需要，人们必

女爷i/ 2022年10月15日 04:58/ 0 赞/ 263 阅读

相关 matlab中微分方程的求解

1、求解一阶微分方程 ![70][] > clc,clear > > dy=@(t,y)y\sin(t)+cos(t);%定义微分方程右端项的匿名函数 > > s=od

╰半夏微凉°/ 2022年05月14日 05:39/ 0 赞/ 339 阅读

相关常微分方程数值求解【python】

简述这里只考虑最为简单的一种常微分方程 d y d x = f ( x , y ) \\frac\{dy\}\{dx\} = f(x,y) dxdy=f(x,y)

短命女/ 2022年04月13日 06:13/ 0 赞/ 531 阅读

相关利用matlab求解方程组的解

定解方程组 ![221250451684152.png][] 计算过程： ![221250454804122.png][] 注意：矩阵系数中，逗号和分号的使

- 日理万妓/ 2021年12月15日 04:09/ 0 赞/ 292 阅读