隐式差分方程c语言,科学网—显式与隐式差分 - 张凌的博文

数值差分方法一般分为显式差分与隐式差分,对二者的区别我查阅了相关的资料。我后面发现关于二者区分不是有一些差异。如Wiki:

Explicit methodscalculate the state of a system at a later time from the state of the system at the current time, whileimplicit methodsfind a solution by solving an equation involving both the current state of the system and the later one. Mathematically, if{displaystyle Y(t)}

b160c901a32a90bdf180749f101db13c894cbed2is the current system state and{displaystyle Y(t+Delta t)}

4984508fe6544e1a596d92cdb4d3cbcbf7fd0230is the state at the later time ({displaystyle Delta t}

8c28867ecd34e2caed12cf38feadf6a81a7ee542is a small time step), then, for an explicit method

{displaystyle Y(t+Delta t)=F(Y(t)),}

4d72c127fbdfdfd3b11386bcc2e8b1672b59e8c4

while for an implicit method one solves an equation

{displaystyle G{Big (}Y(t),Y(t+Delta t){Big )}=0qquad (1),}

3e5330e3b71d7ad310f18f2eeb4412429e873228

to find{displaystyle Y(t+Delta t).}

57166f3b307546bc92c510165b79ac653df805a8

下面几个链接是其它的一些讨论:

这几个讨论,认为如果下一时刻状态值,能通过当前时刻值进行直接计算,则表示显式差分,如果下一时刻状态值,不能通过当前时刻值进行直接计算,而必须求解(1)式来计算,则表示隐式差分。

这个定义似乎很简单,也容易理解,但现在有一个情况是,式(1)如果是线性的,并且通过式(1),我们能反解出Y(t+dt),使其能够写出显示差分的格式。举一个例子吧,

du/dt = f(u)where u is a vector and f is a known function. Then we may approximate the time derivative over a time step h as(implicit) u(n+1) = u(n) + h.f(u(n+1))

if f(u) = Au where A is a constant, then u(n+1)=u(n)+Ahu(n+1),then, u(n+1)=u(n)/(1-Ah)

这种情况下,差分方法是显式还是隐式呢? 本质上我们也求解方程(1), 但是我们得到了显示格式

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值