非线性电路--微分方程数值解法

非线性电路中的微分方程解法

非线性电路理论及应用

周波 电路研-11 2011307080114

非线性电路中的微分方程解法

微分方程数值解法初值： 初值： 所谓初值问题， 所谓初值问题，是函数及其必要的导数在积分的起始 点为已知的一类问题，一般形式为： 点为已知的一类问题，一般形式为：

y

(n)

= f ( x)或y

(n)

= f ( x, y′, y′′,L, y

( n 1)

)

我们先介绍 y′(x) = f (x, y(x)) 简单的一阶问题： 简单的一阶问题：

a≤x≤b

y(a) =y0

(8 1 )

只要f ( x, y )满足(李卜希兹)( Lipschitz条件) : f ( x, y ) f ( x, y ) ≤ L y y

非线性电路中的微分方程解法

由常微分方程的理论可知：上述问题的解唯一存在。 的理论可知：上述问题的解唯一存在。 常微分方程求解求什么？应求一满足初值问题(8—1) 求解求什么？应求一满足初值问题( 的解函数y 如对下列微分方程： 的解函数y = y(x) ,如对下列微分方程：

第八章 序

y ′ = xy 2 dy dy 2 x2 = xy 2 2 = xdx = +c dx y 2 y (0) = 1 y yx =0

2 x2 1 4 x2 4 = 1 c = 2 = 2 = y= y 2 y 4 4 x2

《高等数学》中，微分方程求解，如对一阶微分方程： 高等数学》 微分方程求解，如对一阶微分方程： y′ =f(x,y)是求解解函数y = y(x) ,使满足上述方程。但能够 =f(x,y)是求解解函数 是求解解函数y 使满足上述方程。 求出准确的解析函数y(x)的微分方程是很少的 的微分方程是很少的， 高数》 求出准确的解析函数y(x)的微分方程是很少的，《高数》 中研究微分方程的求解， 分门别类讨论， 中研究微分方程的求解，是分门别类讨论，对不同类型的 微分方程，求解方法不一样，因此，要求解微分方程， 微分方程，求解方法不一样，因此，要求解微分方程，首 先必须认清类型。 先必须认清类型。

非线性电路中的微分方程解法

微分方程 数值解法而常微分方程 初值问题的数值解法，是要寻求解函数 初值问题的数值解法， y(x) 在一系列点y(x 离散点) 在一系列点y(xi ) (离散点):

x0 < x1 < L < xn < L

上 y(xi )的近似值yi( i=1,2,…,n),并且还可由这些(xi,yi) 近似值y i=1,2,…,n),并且还可由这些( ),并且还可由这些(i=1,2,…,n)构造插值函数作为近似函数。上述离散点 i=1,2,…,n)构造插值函数作为近似函数。 相邻两点间的距离h 称为步长， 相邻两点间的距离hi=xi-1-xi 称为步长，若hi 都相等为一定 定步长， 变步长。 则称为定步长 否则为变步长 数h, 则称为定步长，否则为变步长。 a≤ x≤b y ′( x) = f ( x, y ( x)) 本章重点讨论如下 y (a ) = y0 一阶微分方程： 一阶微分方程： 在此基础上介绍一阶微分方程组与 高阶微分方程的数值解法。 高阶微分方程的数值解法。

非线性电路中的微分方程解法

§

欧拉(Euler)法

以Euler法及其改进方法为例,说明 Eul

er法及其改进方法为例 法及其改进方法为例, 常微分方程初值问题数值解法的一般概 Euler法很简单 准确度也不高， 法很简单， 念，Euler法很简单，准确度也不高， 介绍此方法的目的， 介绍此方法的目的，是由于对它的分析 讨论能够比较清楚地显示出方法的一些 特点， 特点，而这些特点及基本方法反映了其 它方法的特点。 它方法的特点。Euler法用于求 Euler法用于求 解一阶微分方 程初值问题： 程初值问题： y ′( x) = f ( x, y ( x)) y (a) = y0 a≤ x≤b

非线性电路中的微分方程解法

b a Euler公 Euler公 对等距节点 x n = x 0 + nh ( n = 0 ,1, L , N ), h = N 式为： 式为： y = y + hf ( x , y ) n +1 n n n n = 1, 2 , L , N (8 2 ) y ( x0 ) = y 0 出发 ,…,x 由x0出发 x1,x2,…,xN,而利用此式可算出对应的 y1,y2,…,yN,式(8-2)称为差分方程(序列{yn}满足的方 ,…,y 称为差分方程 序列{ 差分方程( 程)。下面是Euler公式的推导 ： 下面是Euler公式的推导 一、从几何意义出发：y′=f (x,y)的解函数y=y (x) 在xoy平 从几何意义出发： 的解函数y=y y= xoy平 面上是一族解曲线， 而初值问题则是其中一条积分曲线， 面上是一族解曲线， 而初值问题则是其中一条积分曲线， 假定y 的曲线如图8 从给定的点P 出发， 假定y = y(x)的曲线如图8-1从给定的点P0(x0,y0) 出发，以 为切点，作切线，切线斜率为曲线y P0为切点，作切线，切线斜率为曲线y(x)的切线斜率 y′ =f (x0,y0),因此可 P1 P2 得切线：(点斜式) :(点斜式 得切线：(点斜式)

1.1 Euler法及其简单改进

y(x)

y = y0 + f ( x0 , y0)x x0 ) (

P0

x0 x1 x2

非线性电路中的微分方程解法

切线与x = x1交于P1(x1,y1),在[x0,x1]上以切线 切线与x Euler公式的推导(续1) P 0 P 1 交于P 近似曲线， 近似曲线，y1 = y0 + f ( x0 , y 0 )( x1 x0 ) = y0 + hf ( x0 , y0 )

以y1近似y ( x1 ),当x = x1时，可求y1:

再以P1 ( x1 , y1 )为切点，以 f ( x1 , y1 )为斜率，可得切线： y = y1 + f ( x1 , y1 )( x x1 ) 与x = x2 相交于P2 ( x2 , y 2 )可得：y 2 = y1 + f ( x1 , y1 )( x2 x1 ) = y1 + hf ( x1 , y1 ) 而用P1 P2 近似曲线(在[ x1 , y1 ]上)

重复上述过程，即可求出在x0 , x1 ,L, x N ,L与所对应的微分方程 初值问题的y = y ( x)的y ( xi )(i = 1,2,..., N ,...)的近似值y0 , y1 , L, y N ,L一般地假定已求出 n , yn ), yn ≈ y( xn )可作直线: y = yn + f ( xn , yn )(x xn ) (x 与x = xn+1相交得: yn+1 = yn + f ( xn , yn )(xn+1 xn ) = yn + hf ( xn , yn )

几何意义：用折线近似解曲线y 几何意义：用折线近似解曲线y = y(x),折线不会偏离太 因为每项以f 斜率)修正。 远 ，因为每项以f (x, y)(斜率)修正。

非线性电路中的微分方程解法

微分方程数值解法1欧

欧

:

非线性电路中的微分方程解法

微分方程数值解法

非线性电路中的微分方程解法

微分方程数值解法果 时间区间 为 均 则称 欧

其中K=1,2,3

………..

非线性电路中的微分方程解法

微分方程数值解法

非线性电路中的微分方程解法

微分方程数值解法果 时间区间 为 均 则称 欧其中K=1,2,3………..

非线性电路中的微分方程解法

微分方程数值解法果 时间区间 为 均 则称 欧其中K=1,2,3………..

非线性电路中的微分方程解法

微分方程数值解法果 时间区间 为 均 则称 欧其中K=1,2,3………..

非线性电路中的微分方程解法

§2 龙格-库塔(Runge-kutta)方法龙格龙格-库塔方法的基本思想： y ( xn +1 ) y ( xn ) 我们从研究差商 开始,由微分中值定理 h y ( x n +1 ) y ( xn ) = y ′( xn + θh) (0 < θ < 1) h 于是由微分方程 y ′ = f ( x, y )得到解 y ( xn +1 ) = y ( x n ) + hf ( x n + θh, y ( xn + θh) ) (8 - 11)这里f ( xn + θh, y ( xn + θh))称作区间( x n , xn +1 )上的平均斜率 , 记作k *,即 : k * = f ( x n + θh, y ( xn + θh))

因此只要对平均斜率k 提供一种算法， 因此只要对平均斜率k*提供一种算法，由(8-11)式 11) …… 此处隐藏：3336字，全部文档内容请下载后查看。喜欢就下载吧 ……