...kiyb注:一般有 其中Bk 与yi1 计算公式中 fi … fi?k 各项的系数均可查表得到 fi…1fi…Step 2: 用Adams 显式计算预估值Modified final value yi1...0kyb(1?=yb7注:上式也可通过数值积分导出即将 在区间
单击此处编辑母版标题样式单击此处编辑母版文本样式第二级第三级第四级第五级偏微分方程数值解第十章典型的偏微分方程介绍一. 椭圆型方程2Laplace 方程1Poisson方程 二. 抛物型方程1热传导方程 其中a是常数它表示长度为L的细杆内物体温度分布的规律. 2土壤水运动方程: 3溶质运移方程: (水流稳态)(瞬态)三.双曲型方程1波动方程它表示长度为L的弦振动规律 定解问题定义:决定方程唯一解所
本章讨论常微分方程的数值解法设解函数在节点的近似为② 误差估计做等距分割利用数值微分代替导数项梯形公式隐式的显隐尤拉隐式公式2阶二步注:? 龙格-库塔法的主要运算在于计算 Ki 的值即计算 f 的值Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系:方程组和高阶方程的数值解法精确解
单击此处编辑母版标题样式单击此处编辑母版文本样式第二级第三级第四级第五级第五章 常微分方程数值解 Numerical Methods for Ordinary Differential Equations ? 待求解的问题:一阶常微分方程的初值问题 Initial-Value Problem : 解的存在唯一性(常微分方程理论):只要 f (x y) 在[a b] ? R1
() 假设初值问题()的解y=y(x)唯一存在且足够光滑.对求解区域[ab]做剖分 o称为Euler中点公式或称双步Euler公式.hh= 实际上常将Euler公式与梯形公式结合使用: 称之为改进的Euler方法. 这是一种单步显式方法.012345678910 可见公式的局部截断误差为: y(xn1)-yn1=O(hp1). 三阶R-K
数值计算常微分方程初值问题常微分方程边值问题边值问题:在自变量两端均给定附加条件一般形式:边值问题可能有解也可能无解可能有唯一解也可能有无数解边值问题有3种基本解法迭加法打靶法松弛法指令odefile普通变阶法解ODEodephas3odefile求解初值问题: function f=fun(xy)dy1dx = (1-y(1))-(1-y(2)).y(1)(1-y(2)).2 dy2dx = -
Numerical Analysis J. G. Liu School of Math. Phys. North China Elec. P.U.第八
第一章 基本概念 常微分方程基本概念 常微分方程初值问题 Euler法及其基本问题 线性多步方法 数值稳定性 Runge-Kutta方法:未知函数是一元函数例如:— 不含任意常数的解. 2) n 阶方程的边界条件(或边值条件): 问题 上上世纪初英国物理学家Rutherford发现放射性元素的原子是不稳定的在每一段时间内总有一定比例的原子自然衰变而形成新元素的原子. 记t时刻放射性
2.欧拉法与改进欧拉法1424
() 假设初值问题()的解y=y(x)唯一存在且足够光滑.对求解区域[ab]做剖分 o称为Euler中点公式或称双步Euler公式.hh= 实际上常将Euler公式与梯形公式结合使用: 称之为改进的Euler方法. 这是一种单步显式方法.012345678910 如果单步差分公式的局部截断误差为O(hp1)则称该公式为p阶方法.这里p为非负整数
违法有害信息,请在下方选择原因提交举报