微分动力学问题参数确定.pdf

微分动力学问题参数确定基本方法 微分动力学问题参数确定基本方法 吴树林 吴树林 单位：四川理工学院/理学院（讲师） 单位：四川理工学院/理学院（讲师） 电子科技大学/数学科学学院（博士后） 电子科技大学/数学科学学院（博士后） Email: wushulin_sh@163.com Email: wushulin_sh@163.com 主要内容 Part 1: 数值积分基本知识 Part 2: 微分方程数值解初步 Part 3: 必要的优化技术 Part 1 数值积分初步 §1.1 构造数值积分公式的基本方法 §1.2 Newton-Cotes公式 §1.3 复化求积公式 §1.4 龙贝格求积算 Part 1 引言 F (x ) f (x ) 如果 是 的一个原函数则可用牛顿-莱布尼兹公式 b 计算定积分 I ∫a f (x )dx F (b ) −F (a ) . 但在实际计算中常常会碰到一些困难： ①有些函数的原函数均不能用初等函数表达； f (x ) 的原函数表达式太复杂，计算量太大； ② f (x ) 没有解析表达式，仅知道它在某些离散点处的值。 ③ Part 1 引言 数值积分基本形 b n f ( x )dx ≈ A f ( x ) ∫a ∑ k k k 0 其中x 为求积节点, A 为求积系数. k k 为了便于讨论，我们假定被积函数在闭区间上连续。 Part 1 §1.1 构造数值积分公式的基本方法 [ ,a]b a x x x b ≤ ... ≤ 在积分区间 上取有限个点 0 1 n f (x ) 的n 次插值多项式 作 n L x( ) f x( l) x(n ) ∑ k k k 0 l (x ) (k 0,1,...,n ) 其中 k 为 n 次插值基函数。 L x( ) f (x ) 用 n 近似代替被积函数 ，则得 b b n b f( x) dx L ( x) dx f( x