计算方法-95线性多步法资料
- 格式:ppt
- 大小:838.00 KB
- 文档页数:8
多步法应用于常微分方程的数值解。
从概念上讲,数值方法从初始点开始,然后在时间上向前迈出一小步,以找到下一个求解点。
该过程的下一步是绘制解决方案。
一步法(例如Euler方法)仅引用前一点及其导数来确定当前值。
诸如Runge Kutta之类的方法采取一些中间步骤(例如,半个步骤)来获得高阶方法,但是在进行第二步之前会丢弃所有先前的信息。
多步尝试通过保留和使用先前步骤中的信息而不是丢弃信息来提高效率。
因此,多步法涉及前几个要点和导数。
在多步的情况下,使用先前点和导数的线性组合。
简单的介绍多步法应用于常微分方程的数值解。
从概念上讲,数值方法从初始点开始,然后在时间上向前迈出一小步,以找到下一个求解点。
该过程的下一步是绘制解决方案。
一步法(例如Euler方法)仅引用前一点及其导数来确定当前值。
诸如Runge Kutta之类的方法采取一些中间步骤(例如,半个步骤)来获得高阶方法,但是在进行第二步之前会丢弃所有先前的信息。
多步尝试通过保留和使用先前步骤中的信息而不是丢弃信息来提高效率。
因此,多步法涉及前几个要点和导数。
在多步的情况下,使用先前点和导数的线性组合。
[1-3]具体定义常微分方程的数值方法近似地解决了形式初值问题结果是离散时间的Ti的Y(T)的近似值其中h是时间步长,而I是整数。
Multistep使用上一步中的信息来计算下一个值。
特别地,多步法使用Yi和f(Ti,Yi)来计算所需当前步长的Y值。
因此,多步方法是以下形式的方法:确定系数AI和Bi。
该方法的设计者选择系数平衡了对实际解决方案的需求,以便获得一种易于使用的方法。
通常,许多系数为零以简化该方法。
显式和隐式方法可以区分。
如果Bi = 0,则该方法称为“显式”,因为它可以直接计算yn + s。
如果Bi≠0,则该方法称为“隐式”,因为YN + s的值取决于f(TN + s,yn + s),并且必须为yn + s。
迭代方法(例如牛顿法)通常用于求解隐式公式。
摘要本文主要研究线性三步法的性质及其应用问题,在已有线性多步法基本公式的及线性二步法的基础上,本文又推导出了一个线性三步法公式,并对其进行性质分析验证。
对构造出的线性三步法公式进行相容性、稳定性、收敛性的判断。
对于一些简单而典型的微分方程模型,是可以设法求出其解析解的,并有理论上的结果可利用。
但在数学模型中遇到的常微分方程初值问题模型,通常很难直接求出结果,甚至根本无法求出其解析解,而只能求其近似解。
因此,研究其数值方法,以便快速求得数值解有其重大意义。
对此,本文对常微分方程初值问题模型用线性三步法进行了计算机实现。
本文工作如下:首先,介绍线性多步法公式的基本概念、构造方法、误差分析。
然后,在已有的线性多步法公式特别是线性二步法的基础上,推导出线性三步法的公式,并对其性质进行分析判断。
最后,对构造的线性三步法公式进行应用,主要分析出口服药物在体内吸收变化的情况,用Matlab程序对饮食和非饮食两种情况进行作图比较。
关键词:线性三步法,常微分方程数值解,初值问题,口服药AbstractThis paper studies the nature of the linear three-step method and its application,the existing basic formula of linear multi-step and linear two-step method,the paper has derived a linear three-step formula,and verify the nature of their conduct.Of the constructed linear three-step formula for compatibility,stability,convergence of the judge.For some simple and typical differential equation model,is to derive its analytical solution,and the results are theoretically available.However,mathematical models encountered in the ODEs model,the results are usually difficult to acquire,or even impossible to derive its analytical solution,but can only seek its approximate solution. Therefore,to study the numerical method to quickly obtain the numerical solution to be of significance.In this regard,this paper model of Ordinary Differential Equation of linear three-step method using a computer to achieve.This works as follows:First,the introduction of linear multi-step formula the basic concepts,construction methods,error analysis.Then,in the existing formula,especially linear multi-step linear two-step method based on the derived formula of linear three-step method,and the nature of its judgments.Finally,structural formula of linear three-step application,the main export services of the in vivo absorption of changing circumstances,using Matlab program on food and non food plot comparison of two situations.Key words:Linear three-step method,Numerical Solution of Ordinary Differential Equations,Initial Value Problem,Oral目录第一章绪论 (1)第二章线性多步法的基本理论 (3)2.1常微分方程的数值解法 (3)2.2线性三步法的构造 (4)第三章线性三步法相容性、稳定性、收敛性的研究 (7)3.1相容性 (7)3.2稳定性 (7)3.3收敛性 (8)第四章口服药物在体内的变化 (10)4.1问题的基本概述 (10)4.2建立口服药物的吸收模型 (11)4.2.1问题的提出 (11)4.2.2模型的假设 (11)4.2.3模型的符号及意义 (12)4.2.4应用线性三步法求解 (12)第五章结论与展望 (16)5.1结论 (16)5.2进一步展望 (16)参考文献 (17)致谢 (18)附录 (19)声明 (22)第一章绪论自然界和工程技术中的很多现象,例如自动控制系统的运行、电力系统的运行、飞行器的运动、化学反应的过程、生态平衡的某些问题等,都可以抽象成为一个常微分方程初值问题。
线性多步法:线性多步法(linear multistep method)是1993年发布的数学名词。
线性:线性特性是卷积运算的性质之一,即设a,b为任意常数,则对于函数f(z,y),h(x,y)和g(x,y),{af(x,Y)+bh(z,y)}*g(z,y)=af(x,y)*g(x,y)+bh(x,y)*g(z,y)。
同样有:f(x,y)*{ah(x,y)+bg(x,y)=af(x,y)*h(x,y)+bf(x,y)*g(x,y)。
定义:卷积(Convolution)既是一个由含参变量的无穷积分定义的函数,又代表一种运算。
其运算性质在线性系统理论、光学成像理论和傅里叶变换及其应用中经常用到。
卷积的运算性质有线性特性,复函数的卷积,可分离变量,卷积符合交换律,卷积符合结合律,坐标缩放性质,卷积位移不变性,函数f(x,y)与函数的卷积。
其中线性特性可描述为:设a,b为任意常数,则对于函数f(z,y),h(x,y)和g(x,y),{af(x,Y)+bh(z,y)}*g(z,y)=af(x,y)*g(x,y)+bh(x,y)*g(z,y)。
同样有:f(x,y)*{ah(x,y)+bg(x,y)=af(x,y)*h(x,y)+bf(x,y)*g(x,y)。
多步法:多步法用于普通微分方程的数值解。
从概念上讲,一个数值方法从一个初始点开始,然后在时间上向前迈出一小步,找到下一个解点。
该过程以后的步骤来绘制解决方案。
单步方法(如欧拉方法)只指一个前一点及其导数来确定当前值。
诸如Runge-Kutta的方法采取一些中间步骤(例如,半步)来获得更高阶的方法,但是在进行第二步之前丢弃所有先前的信息。
多步法尝试通过保留和使用先前步骤的信息而不是丢弃它来提高效率。
因此,多步法是指前几个点和导数值。
在多步法的情况下,使用先前点和导数值的线性组合。
常微分方程数值解的多步法。
从概念上讲,一种数值方法是从一个初始点开始的,然后在时间上向前迈出一小步,以找到下一个求解点。
以下过程绘制解决方案。
单步方法(例如欧拉方法)仅参考前一点及其导数来确定当前值。
诸如Runge-Kutta之类的方法采取了一些中间步骤(例如,半步骤)来获得高阶方法,但是在进行第二步之前会丢弃所有先前的信息。
多步方法试图通过保留和使用先前步骤的信息而不是丢弃信息来提高效率。
因此,多步法是指前几个点和导数值。
在多步法的情况下,使用先前点和导数值的线性组合。
常微分方程数值解的多步法。
从概念上讲,一种数值方法是从一个初始点开始的,然后在时间上向前迈出一小步,以找到下一个求解点。
以下过程绘制解决方案。
单步方法(例如欧拉方法)仅参考前一点及其导数来确定当前值。
诸如Runge-Kutta之类的方法采取了一些中间步骤(例如,半步骤)来获得高阶方法,但是在进行第二步之前会丢弃所有先前的信息。
多步方法试图通过保留和使用先前步骤的信息而不是丢弃信息来提高效率。
因此,多步法是指前几个点和导数值。
在多步法的情况下,使用先前点和导数值的线性组合。
具体定义常微分方程的数值方法近似地解决了形式初值的问题结果是离散时间ti处y(t)的近似值:其中h是时间步长,而i是整数。
多步方法使用上一个S步骤的信息来计算下一个值。
特别地,多步方法使用yi和f(ti,yi)来计算当前步骤所需的y值。
因此,多步方法是一种具有以下形式的方法:确定系数ai和bi的方法。
该方法的设计者选择系数来平衡对实际解决方案的需求,从而获得一种易于使用的方法。
通常,许多系数为零以简化方法。
可以区分显式和隐式方法。
如果bi = 0,则此方法称为“显式”,因为此公式可以直接计算yn + s。
如果bi≠0,则此方法称为“隐式”,因为yn + s的值取决于f(tn + s,yn + s),并且必须为yn + s。
迭代方法(例如牛顿法)通常用于求解隐式公式。
95置信区间计算公式
1、样本数量少的话可以直接算:可信区间为阳性样本平均值±标准差(X±SD) 。
2、可信区间介绍:按一定的概率或可信度(1-α)用一个区间来估计总体参数所在的范围,该范围通常称为参数的可信区间或者置信区间(confidenceinterval,CI),预先给定的概率(1-α)称为可信度或者置信度(confidencelevel),常取95%或99%。
3、总体参数的估计,是统计学一大重要的应用。
主要为均数和率的估计,本期做了一个简单的小结,实现该项功能,希望对大家有用。
SPSS对总体均数在探索里是默认实现的,然而对于率却不可以,本例采用比率方法实现。
扩展资料
例:估计该县成年人HBsAg阳性率的95%置信区间。
本例n=100,p=0.12,可采用正态近似法估计总体率的置信区间。
阳性率的95%的置信区间按式(p -Zα/2Sp,p+Zα/2Sp)计算:
下限:p-1.96Sp=0.12-1.96×0.0325=0.0563
上限:p+1.96Sp=0.12+1.96×0.0325=0.1837
所以该县成年人HBsAg阳性率的95%置信区间为(5.63%,18.37%)。