一次阶跃方程 matlab,使用初始条件阶跃函数求解 PDE 方程组

编写方程代码

在编写方程代码之前，您需要确保它的形式符合 pdepe 求解器的要求：

c(x,t,u,∂u∂x)∂u∂t=x-m∂∂x(xmf(x,t,u,∂u∂x))+s(x,t,u,∂u∂x).

由于该 PDE 方程组中有两个方程，PDE 方程组可以重写为

[1001]∂∂t[nc]=∂∂x[d∂n∂x-an∂c∂x∂c∂x]+[Srn(N-n)S(nn+1-c)].

则方程中系数的值为

m=0

c(x,t,u,∂u∂x)=[11](仅对角线值)

f(x,t,u,∂u∂x)=[d∂n∂x-an∂c∂x∂c∂x]

s(x,t,u,∂u∂x)=[Srn(N-n)S(nn+1-c)]

现在，您可以创建一个函数以编写方程代码。该函数应具有签名 [c,f,s] = angiopde(x,t,u,dudx)：

x 是独立的空间变量。

t 是独立的时间变量。

u 是关于 x 和 t 微分的因变量。它是二元素向量，其中 u(1) 是 n(x,t)，u(2) 是 c(x,t)。

dudx 是偏空间导数 ∂u/∂x。它是二元素向量，其中 dudx(1) 是 ∂n/∂x，dudx(2) 是 ∂c/∂x。

输出 c、f 和 s 对应于 pdepe 所需的标准 PDE 形式中的系数。

因此，此示例中的方程可由以下函数表示：

function [c,f,s] = angiopde(x,t,u,dudx)

d = 1e-3;

a = 3.8;

S = 3;

r = 0.88;