当前位置: 首页 > news >正文

差分及其性质不变算子移位算子

目录

  • 差分及其性质
  • 不变算子&移位算子

差分及其性质

设函数 y = f ( x ) =f(x) =f(x)在等距节点 x k = x 0 + k h x_k=x_0+kh xk=x0+kh ( k = 0 , 1 , ⋅ ⋅ ⋅ , n k=0,1,\cdotp\cdotp\cdotp,n k=0,1,⋅⋅⋅,n)上的值 f k = f ( x k ) f_k=f(x_k) fk=f(xk)为已知,这里 h h h 为常数,称为步长。

Δ f k = f k + 1 − f k \Delta f_{k}=f_{k+1}-f_{k} Δfk=fk+1fk

∇ f k = f k − f k − 1 \nabla f_k=f_k-f_{k-1} fk=fkfk1
δ f k = f ( x k + h / 2 ) − f ( x k − h / 2 ) = f k + 1 2 − f k − 1 2 \delta f_k=f(x_k+h/2)-f(x_k-h/2)=f_{k+\frac12}-f_{k-\frac12} δfk=f(xk+h/2)f(xkh/2)=fk+21fk21

分别称为 f ( x ) f(x) f(x) x k x_k xk 处以 h h h 为步长的向前差分、向后差分及中心差分。符号 Δ , ∇ , δ Δ,∇,δ Δ,,δ 分别称为向前差分算子、向后差分算子及中心差分算子。

利用一阶差分可定义二阶差分为
Δ 2 f k = Δ f k + 1 − Δ f k = f k + 2 − 2 f k + 1 + f k . \Delta^2f_k=\Delta f_{k+1}-\Delta f_k=f_{k+2}-2f_{k+1}+f_k. Δ2fk=Δfk+1Δfk=fk+22fk+1+fk.
一般地,可定义 m m m 阶差分为

Δ m f k = Δ m − 1 f k + 1 − Δ m − 1 f k \Delta ^mf_k= \Delta ^{m- 1}f_{k+ 1}- \Delta ^{m- 1}f_k Δmfk=Δm1fk+1Δm1fk
∇ m f k = ∇ m − 1 f k − ∇ m − 1 f k − 1 \nabla ^mf_k= \nabla ^{m- 1}f_k- \nabla ^{m- 1}f_{k- 1} mfk=m1fkm1fk1

因中心差分 δ f k \delta f_k δfk 用到 f k + 1 2 f_{k+\frac12} fk+21 f k − 1 2 f_{k-\frac12} fk21这两个值,实际上不是函数表上的值。

如果用函数表的值,一阶中心差分应写成

δ f k + 1 2 = f k + 1 − f k \delta f_{k+ \frac 12}= f_{k+ 1}- f_k δfk+21=fk+1fk
δ f k − 1 2 = f k − f k − 1 \delta f_{k- \frac 12}= f_k- f_{k- 1} δfk21=fkfk1

二阶中心差分为
δ 2 f k = δ f k + 1 2 − δ f k − 1 2 \delta^2f_k=\delta f_{k+\frac12}-\delta f_{k-\frac12} δ2fk=δfk+21δfk21



不变算子&移位算子

除了已引人的差分算子外,常用算子符号还有不变算子 I I I 及移位算子 E E E
I f k = f k {I} f_k= f_k Ifk=fk
E f k = f k + 1 {E} f_k= f_{k+ 1} Efk=fk+1

于是,由
Δ f k = f k + 1 − f k = E f k − I f k = ( E − I ) f k \Delta f_k=f_{k+1}-f_k=Ef_k-If_k=(E-I)f_k Δfk=fk+1fk=EfkIfk=(EI)fk
可得
Δ = E − I \Delta=E-I Δ=EI

同理,可得
∇ = I − E − 1 \nabla=I- \mathrm{E} ^{- 1} =IE1
δ = E 1 2 − E − 1 2 \delta = \mathrm{E} ^{\frac 12}- \mathrm{E} ^{- \frac 12} δ=E21E21



相关文章:

  • 深入探索C语言中的字符串处理函数:strstr与strtok
  • 磁盘与mysql的故事
  • 聊一聊FutureTask源码中体现的“自旋锁”思想
  • 【无线感知会议系列-22 】Vivisecting Mobility Management in 5G Cellular Networks
  • 类型通配符上限
  • 《pyqt+open3d》第三章——icp配准点对点
  • 交叉编译foxy版ros2部署到ARM上运行
  • Hive的动态分区的原理
  • web第三次作业
  • 初学 mybatis
  • 第J7周:对于ResNeXt-50算法的思考
  • SpringBoot的日志框架
  • 矩阵碰一碰发视频的源码技术开发,支持OEM
  • 华为云之CodeArts IDE的使用体验
  • Git 分布式版本控制
  • VSCode配合cline实现自动编程
  • 完善sql盲注中的其他函数 dnslog+sqlmap外带数据
  • minio在上传pdf文件时设置Content-Type: application/pdf有什么作用
  • 论面向服务的架构设计
  • Qt QSpinBox 总结
  • 建设智能网站/网店营销
  • 政府网站模板免费下载/成都业务网络推广平台
  • 网站建设背景/微博指数查询入口
  • 新疆网站开发报价/黄页88网络营销宝典
  • 做直播网站需要那些技术/品牌策略
  • 一个网站占空间有多少g/最新国际要闻