差分方程及matlab求解综述ppt课件.ppt
《差分方程及matlab求解综述ppt课件.ppt》由会员分享,可在线阅读,更多相关《差分方程及matlab求解综述ppt课件.ppt(102页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益1一个半球状雪堆,其体积融化的速率与半一个半球状雪堆,其体积融化的速率与半球面面积球面面积S成正比,比例系数成正比,比例系数k0。设融化中。设融化中雪堆始终保持半球状,初始半径为雪堆始终保持半球状,初始半径为R且且3小时小时中融化了总体积的中融化了总体积的7/8,问雪堆全部融化还需,问雪堆全部融化还需要多长时间?要多长时间?为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益为了规范事业单位聘用关系,建立和完善适应社
2、会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益差分方程初步差分方程初步为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益第一节第一节差分方程的基本概念差分方程的基本概念一、一、差分的概念差分的概念定义1 设函数yt=f(t)在t=,-2,-1,0,1,2,处有定义,对应的函数值为,y-2,y-1,y0,y1,y2,则函数yt=f(t)在时间t的一阶差分定义为yt=yt+1-yt=f(t+
3、1)-f(t)依此定义类推,有yt+1=yt+2-yt+1=f(t+2)-f(t+1),yt+2=yt+3-yt+2=f(t+3)-f(t+2),为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益一阶差分的性质(1)若yt=C(C为常数),则yt=0;(2)对于任意常数k,(kyt)=kyt;(3)(yt+zt)=yt+zt为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益定义2 函数yt=f(t)在时刻t的二阶差分定义为一阶差分的差分,即 2yt=(yt)=yt
4、+1-yt =(yt+2-yt+1)-(yt+1-yt)=yt+2-2yt+1+yt依此定义类推,有2yt+1=yt+2-yt+1=yt+3-2yt+2+yt+1,2yt+2=yt+3-yt+2=yt+4-2yt+3+yt+2,类推,计算两个相继的二阶差分之差,便得到三阶差分3yt=2yt+1-2yt=yt+3-3yt+2+3yt+1-yt,3yt+1=2yt+2-2yt+1=yt+4-3yt+3+3yt+2-yt+1,为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益一般地,k阶差分(k为正整数)定义为这里 为了规范事业单位聘
5、用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益二、二、差分方程差分方程定义3 含有未知函数yt=f(t)以及yt的差分yt,2yt,的函数方程,称为常差分方程(简称差分方程);出现在差分方程中的差分的最高阶数,称为差分方程的阶.n阶差分方程的一般形式为F(t,yt,yt,nyt)=0,其中F是t,yt,yt,nyt的已知函数,且nyt一定要在方程中出现 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益定义3 含有两个或两个以上函数值yt,yt+1,的函数方程,称为(常)差分方
6、程,出现在差分方程中未知函数下标的最大差,称为差分方程的阶 n阶差分方程的一般形式为F(t,yt,yt+1,,yt+n)=0,其中F为t,yt,yt+1,,yt+n的已知函数,且yt和yt+n一定要在差分方程中出现.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益三、三、差分方程的解差分方程的解定义4 如果将已知函数yt=(t)代入方程F(t,yt,yt+1,,yt+n)=0,使 其 对 t=,-2,-1,0,1,2,成 为 恒 等 式,则 称yt=(t)为方程的解.含有n个任意(独立)常数C1,C2,Cn的解yt=(t,C1
7、,C2,Cn)称为n阶差分方程的通解.在通解中给任意常数C1,C2,Cn以确定的值所得的解,称为n阶差分方程的特解.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益例如,函数yt=at+C(a为已知常数,C为任意常数)是差分方程yt+1-yt=a的通解.而函数yt=at,yt=at-1,均是这个差分方程的特解.由差分方程的通解来确定它的特解,需要给出确定特解的定解条件.n阶差分方程F(t,yt,yt+1,,yt+n)=0常见的定解条件为初始条件.y0=a0,y1=a1,,yn-1=an-1,这里a0,a1,a2,,an-1均为
8、已知常数为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益只要保持差分方程中的时间滞后结构不变,无论对t提前或推后一个相同的等间隔值,所得新方程与原方程是等价的,即二者有相同的解.例如,方程ayt+1-byt=0与方程ayt+2-byt+1=0都是相互等价的 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益四、四、线性差分方程及其基本定理线性差分方程及其基本定理形如yt+n+a1(t)yt+n-1+a2(t)yt+n-2+an-1(t)yt1+an(t)yt=f
9、(t)的差分方程,称为n阶非齐次线性差分方程.其中a1(t),a2(t),an-1(t),an(t)和 f(t)都 是 t的 已 知 函 数,且an(t)0,f(t)0.而形如yt+n+a1(t)yt+n-1+an-1(t)yt+1+an(t)yt=0 的 差 分 方 程,称 为 n阶 齐 次 线 性 差 分 方 程.其 中ai(t)(i=1,2,,n)为t的已知函数,且an(t)0.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益如果ai(t)=ai(i=1,2,n)均为常数(an0),则有yt+n+a1yt+n-1+a2y
10、t+n-2+an-1yt+1+anyt=f(t),yt+n+a1yt+n-1+a2yt+n-2+an-1yt+1+anyt=0 分别称为n阶常系数非齐次线性差分方程和n阶常系数齐次线性差分方程.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益定理定理1(齐齐次次线线性差分方程解的叠加原理性差分方程解的叠加原理)若y1(t),y2(t),ym(t)是齐次线性差分方程yt+n+a1yt+n-1+a2yt+n-2+an-1yt+1+anyt=0的m个特解(m2),则其线性组合y(t)=A1y1(t)+A2y2(t)+Amym(t)也
11、是方程的解,其中A1,A2,Am为任意常数定理定理2n阶齐次线性差分方程yt+n+a1yt+n-1+a2yt+n-2+an-1yt+1+anyt=0一定存在n个线性无关的特解为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益定理定理3(齐齐次次线线性差分方程通解性差分方程通解结结构定理构定理)如 果 y1(t),y2(t),yn(t)是 齐 次 线 性 差 分 方 程yt+n+a1yt+n-1+a2yt+n-2+an-1yt+1+anyt=0的n个线性无关的特解,则方程的通解为:yA(t)A1y1(t)+A2y2(t)+Anyn
12、(t),其中A1,A2,An为n个任意(独立)常数 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益定理定理4(非非齐齐次次线线性差分方程通解性差分方程通解结结构定理构定理)如果 (t)是非齐次线性方程yt+n+a1(t)yt+n-1+a2(t)yt+n-2+an-1(t)yt1+an(t)yt=f(t)的一个特解,yA(t)是其对应的齐次线性方程yt+n+a1yt+n-1+a2yt+n-2+an-1yt+1+anyt=0的通解,那么,非齐次线性差分方程的通解为:y(t)=yA(t)+(t)即y(t)=A1y1(t)+A2y2
13、(t)+Anyn(t)+(t),这里A1,A2,An为n个任意(独立)常数为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益第二节第二节一阶常系数线性差分方程一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为yt+1+ayt=f(t)和yt+1+ayt=0,其中f(t)为t的已知函数,a0为常数.分别称为一阶常系数非齐次线性差分方程和其对应的齐次差分方程.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益一、一、齐次差分方程的通解齐次差分方程的通解将方程yt
14、+1+ayt=0改写为:yt+1=-ayt,t=0,1,2,假定在初始时刻(即t=0)时,函数yt取任意值A,那么由上式逐次迭代,算得y1=-ay0=-aA,y2=-ay1=(-a)2A,方程的通解为yt=A(-a)t,t=0,1,2,如果给定初始条件t=0时yt=y0,则A=y0,此时特解为:yt=y0(-a)t 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益二、二、非齐次方程的通解与特解非齐次方程的通解与特解1.迭代法求通解迭代法求通解将方程改写为 yt+1=(-a)yt+f(t),t=0,1,2,逐步迭代,则有y1=(
15、-a)y0+f(0),y2=(-a)2y0+(-a)f(0)+f(1),y3=(-a)3y0+(-a)2f(0)+(-a)f(1)+f(2),为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益由数学归纳法,可得 yt=(-a)ty0+(-a)t-1f(0)+(-a)t-2f(1)+f(t-1)=(-a)ty0+,(t=0,1,2,),yA(t)=(-a)ty0为对应的齐次方程的通解.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益解解例例为了规范事业单位聘用关系
16、,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益方程的通解方程的通解为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益2.待定系数法求特解待定系数法求特解情形 f(t)为常数方程变为yt+1+ayt=b,a,b均为非零常数试以 (为待定常数)形式的特解代入方程得+a=(1+a)=b当a-1时,可求得特解当a=-1时,改设特解 (为待定系数),将其代入方程得 (t+1)+a t=(1+a)t+=b 求得特解为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用
17、制度,保障用人单位和职工的合法权益方程的通解为 解解例例为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益情形情形f(t)为为t的多项式的多项式不妨设f(t)=b0+b1t(t的一次多项式),即 yt+1+ayt=b0+b1t,t=1,2,,其中a,b0,b1均为常数,且a0,b10试以特解 =a+bt,(a,b为待定系数)代入方程得a+b(t+1)+a(a+bt)=b0+b1t,上式对一切t值均成立,其充分必要条件是:为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合
18、法权益当1+a0时,即a-1时,方程的特解为当a=-1时,改设特解 =(a+bt)t=at+bt2将其代入方程可求得特解为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益方程的通解为 解解例例为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益情形情形f(t)为指数函数为指数函数不妨设f(t)=bdt,b,d均为非零常数,方程变为 yt+1+ayt=bdt,t=0,1,2,求得特解当a+d0时,设方程有特解 =dt,为待定系数.将其代入方程得 dt+1+adt=bd
19、t,当a+d=0时,改设方程的特解 =tdt,为待定系数,将其代入方程可求得特解=btdt当a+d=0时,改设方程的特解 =tdt,为待定系数,将其代入方程可求得特解=btdt求得特解当a+d=0时,改设方程的特解 =tdt,为待定系数,将其代入方程可求得特解=btdt为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益方程的通解为解解例例为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益情形 f(t)为正弦、余弦型三角函数 设f(t)=b1cost+b2sint,
20、其中b1,b2,均为常数,且 0,b1与b2不同时为零.于是非齐次方程变为yt+1+ayt=b1cost+b2sint,a0,t=0,1,2,设方程有特解 =acost+bsint,a,b均为待定系数.将其代入方程得acos(t+1)+bsin(t+1)+aacost+absint =b1cost+b2sint,(acos+bsin+aa)cost+(-asin+bcos+ab)sint=b1cost+b2sint 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益(acos+bsin+aa)cost+(-asin+bcos+a
21、b)sint=b1cost+b2sint 上式对t=0,1,2,恒成立的充分必要条件是 其系数行列式为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益当D0时,则可求得其解当D=(a+cos)2sin2=0时,则有改设特解 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益代入方程并整理可得 方程的通解为 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益例例求差分方程求差分方程yt+1-2yt=c
22、ost的通解的通解解解对应齐次方程的通解为对应齐次方程的通解为yA(t)=A2t设设非非齐齐次方程的特解次方程的特解为为=cost+sint,其中其中,为为待定系数待定系数将其代入原方程将其代入原方程,并利用三角函数的和角公式并利用三角函数的和角公式,得得为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益所给方程的通解为所给方程的通解为为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益第三节第三节二阶常系数线性差分方程二阶常系数线性差分方程二阶常系数线性差分方程的
23、一般形式为yt+2+a1yt+1+a2yt=f(t),t=0,1,2,,其中f(t)为t的已知函数,a1,a2为已知常数,且a20,称为二阶常系数非齐次线性差分方程 特别地,当f(t)0时,方程变为yt+2+a1yt+1+a2yt=0 称为对应的齐次差分方程为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益一、齐次差分方程的通解称2a1+a2=0为二阶常系数非齐次线性差分方程或其对应的齐次差分方程的特特征征方方程程它的解(或根)称为方程的特征根特征根(值值)特征方程的两个根为(1)特征根为相异的两实根特征根为相异的两实根当0时,
24、1,2为两相异的实根.y1(t)=1t与y2(t)=2t是齐次差分方程的两个线性无关的特解.为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益齐次差分方程的通解 1,2由特征方程确定,A1,A2为两任意(独立)常数 例例求差分方程yt+2-7yt+1+12yt=0的通解解解特征方程为2-7+12=(-3)(-4)=0,有两相异实特征根 1=3,2=4原方程的通解为 为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益(2)特征根为两相等的实根特征根为两相等的实根当当
25、=0时时,=1=2=为为两相等的两相等的实实根根.方程的一个特解:方程的一个特解:yt(t)=t方程的另一个特解方程的另一个特解为为y(t)=t t,且与且与 t线线性无关性无关.方程的通解为方程的通解为为了规范事业单位聘用关系,建立和完善适应社会主义市场经济体制的事业单位工作人员聘用制度,保障用人单位和职工的合法权益例例求差分方程求差分方程yt+2-4yt+1+4yt=0的通解的通解.解解特征方程为特征方程为 2-4+4=(-2)2=0,方程有重特征根方程有重特征根=1=2=2原方程的通解原方程的通解为为yA(t)=(A1+A2t)2t,A1,A2为为任意常数任意常数为了规范事业单位聘用关系
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 方程 matlab 求解 综述 ppt 课件
限制150内