微分方程的多解与变号解
数学分析知识点总结(微分方程)
2.7．微分方程初步2.7.1 概说涉及到量的变化率满足的制约关系，通常是含有导数的方程——微分方程。
简单例子：（1）放射性物质，在每一时刻t ，衰变的速率dm dt -（由于是减少，因此0dm dt<，速率为标量，是正值）正比于该放射性物质尚存的质量，因此质量应满足一下微分方程。
dmkm dt-= （2）质量为m 的物体自由落体，取坐标轴沿竖直方向指向地心，下落距离()y y t =应该满足牛顿第二定律F ma =，即22d ymg m dt=（3）质量为m 的跳伞员下落，所受空气阻力正比下降的速度，取坐标轴沿竖直方向指向地心，则t 时刻下降距离()y y t =满足22dy d y mg k m dt dt-=（1）如下图所示，钢球在以水平光滑杆上，受到弹力而来回整栋，原点位置为O ，钢球在t 时刻的坐标()x x t =满足微分方程()22d x kx m dt-=如果钢球还受到一个与速度成正比，方向与速度相反的阻尼力的作用，那么它所满足的微分方程是22dx d xkx h m dt dt--=总结：最简单的一阶微分方程是()dxf t dt= 其中t 是自变量，上述方程的一般解应该是()x f t dt C =+⎰最简单的n 阶方程()n nd xf t dt = 它等价于说11n n d xdt--是()f t 的原函数，即11()n n d xf t dt C dt --=+⎰则再次积分，一直积分下去得到111()(1)!n n n t x f t dtdt C C t C n --=++++-⎰⎰2.7.2 一阶线性微分方程考察下面的方程()()dxa t xb t dt+= 方程中有未知函数的一阶导数，且其一阶导数的系数为常数，其余部分未知函数最高层次数为一次，称为线性，上述方程为一阶线性微分方程。
如果()0a t =，则称为一阶线性常微分方程。
试着求解上述方程，方程两端都乘以()a t dte ⎰，得到()()()()()a t dta t dt a t dt dxe a t e x b t e dt⎰⎰⎰+= 即为下面的形式()()()()a t dta t dta t dt d e dxe x b t e dt dt ⎛⎫⎰ ⎪⎝⎭⎰⎰+=即()()()a t dta t dt d xeb t e dt⎛⎫⎰ ⎪⎝⎭⎰=于是有()()()a t dta t dtxe b t e dt C ⎰⎰=+⎰那么有()()()a t dt a t dt x e b t e dt C -⎛⎫⎰⎰=+ ⎪⎝⎭⎰ 这就是一阶线性微分方程的一般解。
第X章 代数(微分)方程(组)
符号表达式的操作—expand
expand(EXPR)—对EXPR表达式进行多项式、
三角函数、指数对数函数等展开。
>> syms x y; >> expand((x-2)*(x-4)) ans =x^2-6*x+8 >> expand(cos(x+y))
ans =cos(x)*cos(y)-sin(x)*sin(y)
>> expand(exp((x+y)^2)) ans =exp(x^2)*exp(x*y)^2*exp(y^2) >> expand([sin(2*x),cos(2*x)]) ans =[ 2*sin(x)*cos(x),2*cos(x)^2-1]
23
符号表达式的操作—factor
factor(EXPR)—对EXPR表达式(或正整数)进
sol：输出的求解结果，为结构体类型数据，sol.x表示 时间变量采样值，sol.y表示函数值分量上的取值。 ddefun：表示待求时滞微分方程组。 lags：表示时间延迟常数。 history：表示求解时间区间上的变量值函数。 tspan：表示求解时间范围。
18
dde23函数函数举例
‘eq1’,’eq2’,…,’eqn’是字符串表达的方程或是
字符串表达式,’v1’,’v2’,…,’vn’是求解变量名。 eq1,eq2,…,eqn是不含‚等号‛的表达式，则是对 eq1,eq2=0,…,eqn=0求解。 exp1,…,expn只能是符号表达式，v1,…vn是求解的符号 变量。
>> syms x ;……………………………………..(1) >> f1=x^4-5*x^3+5*x^2+5*x-6; >> factor(f1) ans = (x-1)*(x-2)*(x-3)*(x+1)
信号与系统第二章第一讲
则相应于1的k阶重根，有k项：
( A1t k 1 A2t k 2 Ak 1t Ak )e1t ( Ai t k i )e1t
i 1
k
例2-3
信 号 与 系 统
求如下所示的微分方程的齐次解。
Hale Waihona Puke d3 d2 d r (t ) 7 2 r (t ) 16 r (t ) 12r (t ) e(t ) 3 dt dt dt
等式两端各对应幂次的系数应相等，于是有：
信 号 与 系 统
特解为： 联立解得：
3B1 1 4 B1 3B2 2 2 B 2 B 3 B 0 2 3 1
统
线性时不变系统
线性的常系数微分方程
按照元件的约束特性及 系统结构的约束特性
也即：
具体系统物理模型
常系数微分方程建立
（1）元件端口的电压与电流约束关系
iR (t ) R
信 号 与 系 统
vR (t )
C
vR (t ) iR (t ) R
dvC (t ) iC (t ) C dt
vR (t ) Ri R (t )
与
时域经典法就是直接求解系统微分方程的方法。这种方 系 法的优点是直观，物理概念清楚，缺点是求解过程冗繁，应 用上也有局限性。所以在20世纪50年代以前，人们普遍喜欢 统 采用变换域分析方法(例如拉普拉斯变换法)，而较少采用时 域经典法。20世纪50年代以后，由于δ(t)函数及计算机的普 遍应用，时域卷积法得到了迅速发展，且不断成熟和完善， 已成为系统分析的重要方法之一。时域分析法是各种变换域 分析法的基础。
信 号 与 系 统
is (t )
信号与系统3-3微分方程的拉氏变换解课件
U (s)
u(t)
U (s)
sL
I (s)
3.4 微分方程的拉氏变换解
返回
用拉氏变换求解线性常系数微分方程，主要用到拉氏变换 的微分性质：
• 对于一阶导数： L [ y(t)] sY (s) y(0 )
• 对于二阶导数：L [ y(t)] s[sY (s) y(0 )] y(0 )
• 对于三阶导数：
s2Y (s) s y(0 ) y(0 )
4
系统函数包含了两层含义
系统函数与冲激响应
yzs (t) h(t) f (t)
Yzs (s) H (s)F (s)
可见系统函数可h视(t)为系统H对(s复) 指数信号的加权
系系统数函，数它与与输复入指无数关信，号反映系统本身特性。只
不复过 频h域(yt的)zs是(描t)系述统h。(在t) 时f域(t)的描h(述t) ，eHst (s)是对系统在
(s2 5s 6)Y (s) 3 s 4 s 1
3Y (s)s2s 1 5s 6
s2
s4 5s
6
零状态响应
零输入响应
y(t) yzs (t) yzi (t)
8
例 3.21
解2
系统函数为
H (s)
s2
3 5s
6
零状态响应为
Yzs
(s)
H (s)F(s)
(s
3)( s
3
2)( s
定义
H (s)
零状态响应的拉氏变换 激励信号的拉氏变换
Yzs (s) F (s)
二阶系统零状态响应
Yzs (s)
b2 s 2 a2s2
b1s b0 a1s a0
F(s)
常微分方程的符号解
一阶常微分方程是指未知数的导数为一次方 的方程，如 dy/dx = f(x, y)。二阶及高阶常 微分方程则是导数的次数更高的方程，如 d^2y/dx^2 = f(x, y, dy/dx)。根据具体问 题的需求，可以选择不同类型和阶数的常微 分方程来描述。
常微分方程的应用
总结词
常微分方程在物理学、工程学、经济学等领域有广泛的应用。
符号计算库的完善
现有的符号计算库如SymPy等将继续完善，提 供更多实用的功能和算法。
符号解法与其他方法的结合
符号解法将与数值解法、计算机代数等其他方法相结合，形成更全面、高效的 求解方案。
符号解法与其他方法的结合
01
符号解法与数值解法的结合
通过混合符号计算和数值计算的方法，可以在保证析方法 来求解常微分方程的方法，它能够给 出方程的解析解，即以数学表达式形 式给出的解。
符号解法通常使用符号计算软件（如 SymPy）进行计算，能够处理包含任意 常数、函数和它们的导数的复杂方程。
符号解法的步骤
建立微分方程
导入符号计算库
使用Python的SymPy库或其他符 号计算软件进行符号计算。
详细描述
在物理学中，常微分方程被用来描述物体的运动规律、电磁波的传播等。在工程学中，常微分方程用于分析机械、 电路、控制系统等各种实际系统的动态行为。此外，经济学中用于研究经济变量的变化规律，如供需关系、人口 增长等。通过建立适当的常微分方程模型，可以深入了解各种实际问题的内在机制。
02
常微分方程的符号解法
05
常微分方程的数值解法的应用
在物理中的应用
描述物体运动规律
01
常微分方程可以用来描述物体的运动规律，例如牛顿第二定律、
信号与系统(第四版)陈生潭第四章课后答案
五、单边拉氏变换与傅里叶变换的关系
F (s) f (t) est d t 0
Re[s]>0
F (j) f (t) e j t d t
要讨论其关系，f(t)必须为因果信号。
根据收敛坐标0的值可分为以下三种情况：
（1）0<0，即F(s)的收敛域包含j轴，则f(t)的傅里叶
变换存在，并且
F(j)=F(s) s=j
如f(t)=e-2t(t) ←→F(s)=1/(s+2) , >-2；
则 F(j)=1/( j+2)
第5-13页
■
©西安电子科技大学电路与系统教研中心
信号与系统 电子教案
4.1 拉普拉斯变换
（2）0 =0，即F(s)的收敛边界为j轴，
F(j) lim F(s) 0
■
©西安电子科技大学电路与系统教研中心
信号与系统 电子教案
4.1 拉普拉斯变换
例3 双边信号求其拉普拉斯变换。
f3 (t)
f1 (t)
f
2
(t
)
e e
t t
, ,
t0 t 0
求其拉普拉斯变换。
解 其双边拉普拉斯变换 Fb(s)=Fb1(s)+Fb2(s)
jω
仅当>时，其收敛域为
例1:已知因果信号f(t)的象函数F(s)=
s s2 1
求e-tf(3t-2)的象函数。
解：e-tf(3t-2) ←→
(s
s 1 1)2
9
e
2 (s1) 3
例2：
e-2t cos 3t
s 2 (s 2)2 9
e 2t sin 3t
微分方程的应用
微分方程的应用嘿，朋友！想象一下，你正走在繁忙的街道上，周围的人们行色匆匆。
突然，你的脑海中蹦出了一个奇怪又神秘的词——微分方程。
是不是觉得这东西好像跟咱们的日常生活风马牛不相及？那你可就大错特错啦！其实，微分方程就像一个隐藏在幕后的超级英雄，默默为我们解决着各种看似无解的难题。
比如说，当医生在研究药物在人体内的代谢过程时，微分方程就派上了大用场。
你看，药物进入身体后，它的浓度会随着时间不断变化。
这时候，微分方程就像一个精准的导航仪，能够帮助医生预测在不同时间点药物的浓度分布，从而确定最佳的用药剂量和时间间隔。
这难道不神奇吗？再想想咱们每天都离不开的手机。
手机电池的充电和放电过程，其实也可以用微分方程来描述。
要是没有对这个过程的精确理解和控制，咱们的手机可能一会儿就没电啦，那得多让人崩溃！还有交通领域，你有没有想过为什么交通信号灯能合理地控制车流量？这背后也有微分方程的功劳。
它能够模拟车辆的流动情况，帮助交通规划师设计出更高效的道路系统和信号灯设置。
就像你在玩游戏时，需要掌握规则才能玩得好，微分方程就是这些现实问题中的“游戏规则”。
你可能会疑惑，这东西这么厉害，难道就没有搞不定的时候？嘿，还真有！但这并不妨碍它在大多数情况下大显身手。
想象一下，如果没有微分方程，我们的世界会变成什么样？可能医生治病就像盲人摸象，交通会乱成一锅粥，各种技术发展也会停滞不前。
所以说，微分方程可不仅仅是那些数学课本里枯燥的公式和符号，它是实实在在影响着我们生活的重要工具。
它在我们看不见的地方默默发挥着作用，让我们的生活变得更加便捷、高效和美好。
朋友，现在你还觉得微分方程离我们很遥远吗？相信你已经明白，它就在我们身边，无处不在，为我们的生活保驾护航！。
matlab求解微分方程解析解
matlab求解微分方程解析解在数学和工程学科中，微分方程是一种重要的数学工具，它涉及到很多实际问题的模型和解决方法。
而Matlab作为一款强大的数学软件，可以方便地求解微分方程的解析解。
Matlab中求解微分方程的一种常见方法是使用符号计算工具箱（Symbolic Math Toolbox），它可以处理符号表达式和符号函数，包括微积分、代数、矩阵、符号等数学操作。
首先，我们需要定义微分方程的符号变量和初值条件。
例如，我们假设要求解的微分方程为dy/dx = x^2，初值条件为y(0)=1，则可以使用如下代码：syms x yode = diff(y,x) == x^2;cond = y(0) == 1;然后，我们可以将微分方程和初值条件作为参数传递给dsolve函数来求解微分方程的解析解。
例如：sol = dsolve(ode, cond);其中，sol为求解得到的符号表达式，可以使用vpa函数将其转换为数值解。
例如：sol_num = vpa(sol, 5);这样，我们就得到了微分方程的解析解，并将其转换为5位有效数字的数值解。
除了使用符号计算工具箱，Matlab还提供了许多数值方法来求解微分方程的数值解。
例如，可以使用ode45函数来求解微分方程的数值解。
例如，求解dy/dx = x^2，y(0)=1的数值解可以使用如下代码：fun = @(x,y) x^2;[t,y] = ode45(fun, [0,1], 1);其中，fun为微分方程的函数句柄，[0,1]为求解区间，1为初值条件。
t和y分别为求解得到的时间序列和解向量。
总之，Matlab提供了多种方法来求解微分方程的解析解和数值解，可以根据实际问题的需要选择不同的方法来求解。
变微分方程
变微分方程
变微分方程是指对已知的微分方程进行变换，以便将其转化为更简单或更易于解析的形式。
这种变换可以通过代换、变量替换、参数化等方式进行。
下面是一些常见的变微分方程的方法：
1. 代换法：通过引入一个新的变量或函数，将原微分方程转化为一个新的微分方程。
这个新的微分方程可能更容易求解或转化为已知的标准形式。
常见的代换包括指数代换、三角函数代换等。
2. 变量替换法：通过引入新的变量，将原微分方程转化为关于新变量的微分方程。
这样可以改变微分方程的形式，使其更容易求解或分离变量。
常见的变量替换包括极坐标替换、球坐标替换等。
3. 参数化方法：将原微分方程的解表示为一个参数方程，通过引入参数，将微分方程转化为一个关于参数的方程。
这样可以将原微分方程的求解问题转化为参数方程的求解问题。
4. 齐次化方法：对于非齐次微分方程，可以通过适当的变量替换将其转化为齐次微分方程。
这样可以简化求解过程，因为齐次微分方程的解结构更简单。
5. 线性化方法：对于非线性微分方程，可以通过适当的变换将其转化为线性微分方程。
线性微分方程的求解通常更为直接和简单。
需要注意的是，变微分方程的方法取决于具体的微分方程形
式和求解目标。
不同的微分方程可能需要使用不同的方法进行变换。
因此，在变微分方程之前，需要对原微分方程的形式和性质进行分析，并选择合适的变换方法。
偏微分方程
关键词 ：无 限级型 函数 ；Ｂ ｒｌ 向；零 点收敛指数 ｏｅ 方
０００８ ８ ９０ ３ １ ０・ ４ 常微 分 方 程 １ ４
ｏ ｃａｓ ｆ ｉｇｌ ｃ ｎ —ｒｅ ｂ ｕ ｄｒ ａ ｅｐｏ ｌｍｓ［ ， ｆ ｓ ｎ ｕａ ｓ ｏ ｄｏｄｒ ｏ ｎ ａｙｖ ｕ ｒｂｅ ａｌ ｏｓ ｒ ｅ ｌ 刊 报 ． ２ ０ ，５ （） １ ５ ￣ １６ 一 ０ ７ ０６． ３ ７ ３ ４ 一
Ｉ （ ＝ｆ ｔ ｌ） ２ ） ３ ） ， ｔ ［， ， “ ｔ ａ， （， （， （， ｆ） “ ｆ“ ｆ “ ｆ ） ∈ ０１ 】
（ ） ｉ０ ＝Ｕ （ ：０ ０ ＝Ｕ （ ） ｉ ） ， １ ｉ ， ， ＝１２ ３
１ ０ ・４ １ ４
类带不 变号权函数 二阶超线性 方程 的周 期碰撞解 ＝Ｐｒ ｄｃ ｅｏｉ ｉ
［ 中］黄建科（ 刊， ／ 西安陆军学 院军事运筹教研 室， 西安 ７０ ０ ） １ １８， 吴筱宁 ∥纺织高校基础科 学学报 ． ２ ０ ，２ （） ３２ ９ 一 ０ ７ ０４． ９ ￣３ ５ － 研究 了一类具有 Ｈｏｌ ｇ Ｉ 类功 能反应且周期系统 的三维非 自 ｌｎ Ｉ ｉ Ｉ 治捕食系统 ，利用微分 不等式给 出了系统永久持续生存 的充分 条件 ；同时 ，通过构造 Ｌ ａｕ ｏ 函数 ，建立 了系统正周期解 ｙｐ ｎ ｖ 存在惟一及全局渐 近稳 定的充分判据．参 ６ 关键词：捕食系统 ；Ｈ ｌｎ Ｉ ｏｌ ｇＩ 类功能反应 ；永久持续生存；周 ｉ Ｉ
０ １Ｉ ５ ８ ０Ｄ ８
一
［ ，中］孙 忠 民（ 刊 ／ 山东省潍坊 教育学 院，青州 ２ ２ ０ ） ６ ５ ０ ，赵增 勤． ∥系 统 科 学 与 数 学 ． ２０ ，２ （） ８ １ ８ ９ 一 ０ ７ ７６． １ ～ １ 一 利用 Ｋ ａｎｓｌ ｉ 不动点 定理 ，结合 Ｌ ｒｙＳ ｈｕ ｅ 度 ，研究 ｒｓｏｅｓ ｉ ｋ ｅａ —ｃａ ｄｒ 下 列 三 阶 微 分 方 程 组 边 值 问题
向量微分方程x′=ax的一种求解方法
向量微分方程x′=ax的一种求解方法，在数学中，向量微分方程是一种表示变量及其变化随时间的场合的常用方法。
例如，x′=ax，表示变量a随时间x的变化，即x(t+delta t)=a*x(t)。
解决向量微分方程，可以使用变换方法，解析解，积分法等多种求解方法。
变换法求解向量微分方程x′=ax，其基本思想是把原方程变换成能够解决的更容易的方程，也就是通过变换，将偏微分方程变换成通常的积分形式。
变换法主要依靠选择恰当的变换函数，从而使得变换前后的微分方程之间存在简单的映射，本质上在解微分方程等号前和等号后分别进行积分操作，计算各个积分量的乘积，其结果就是需要求解的微分方程的解析形式。
同样，解析解法也可以求解向量微分方程x′=ax。
解析解法的思想与变换法相似，也是在条件的情况下，可以将微分方程变换成积分形式，这时就可以用积分来求得方程解析解。
可以各元素分别求解，将求得的结果相加，就可以给出方程x′=ax的解析解。
积分法也可以用来求解向量微分方程x′=ax。
它将通常微分方程转为积分形式，通过一定的技巧，把积分中的区间拆分为许多小区间，然后用泰勒展开，将非定积分拆分成定积分，计算各定积分的和，就可以给出方程x′=ax的解。
通过以上三种求解方法，我们可以容易地求解向量微分方程x′=ax。
只要选定合适的变换函数，或者把原微分方程求解成积分形式，分解小区间，就可以求得方程的解析形式。
这些求解方法也广泛用于数值计算，因为无论什么情况下都可以利用它们计算出准确的解。
此外，这些求解方法也可以被应用到其他的微分方程中，如热力学，磁现象等，使研究者及工程师们更好地了解问题的物理规律，也能更容易地提出合适的解决方案。
北邮信号与系统复习： 典型例题 微分方程3种解法——含冲激函数匹配法
k 1 n
特解 强迫响应
11
2.5节
总结
P58页，理解问题核心。 缺少如下补充，请大家抄写到58页公式2-41旁边 零输入响应:rzi(0-)=rzi(0+);-> 对应于起始储能 rzi’(0-)=rzi’(0+);
零状态响应:rzs(0-) = 0 ;->从定义上与储能无关 rzs’(0-) = 0;
i)写出完全解形式：
ii)冲激函数匹配法求跳变值：根据t=0时刻微分方程左 右两端的 及其各阶导数应该平衡相等
系统用微分方程表示时，系统地0-状态到0+状
态有无跳变决定于微分方程的右端自由项是否
包含
及其高阶导数。有则跳变。 (t )
5
例题2：经典法
带入
在 考虑换路时情况，即t=0时刻e(t)有2 u(t) 变化，得
2
例题2：经典法
，
，
求完全响应。
3
例题2：经典法
解：①由特征根写出齐次解形式 i)特征方程： 特征根： ii)齐次解形式： ②求特解
i) t>0时自由项=16[ຫໍສະໝຸດ 2+2u(t) ] * 4
ii)0不是特征根，设特解为：
iii)代入方程左边解得：B=8/5
4
例题2：经典法
③求完全解中的齐次解待定系数
信号与系统(郑君里)复习要点
信号与系统复习书中最重要的三大变换几乎都有。
第一章 信号与系统 1、信号的分类①连续信号和离散信号 ②周期信号和非周期信号 连续周期信号f (t )满足f (t ) = f (t + m T )， 离散周期信号f(k )满足f (k ) = f (k + m N )，m = 0,±1,±2,…两个周期信号x(t)，y(t)的周期分别为T 1和T 2，若其周期之比T 1/T 2为有理数，则其和信号x(t)+y(t)仍然是周期信号，其周期为T 1和T 2的最小公倍数。
③能量信号和功率信号 ④因果信号和反因果信号2、信号的基本运算（+ - × ÷） 2.1信号的（+ - × ÷）2.2信号的时间变换运算 （反转、平移和尺度变换） 3、奇异信号3.1 单位冲激函数的性质f (t ) δ(t ) = f (0) δ(t ) ， f (t ) δ(t –a) = f (a) δ(t –a)例： 3.2序列δ(k )和ε(k ) f (k )δ(k ) = f (0)δ(k ) f (k )δ(k –k 0) = f (k 0)δ(k –k 0) 4、系统的分类与性质4.1连续系统和离散系统4.2 动态系统与即时系统 4.3 线性系统与非线性系统 ①线性性质 T [a f (·)] = a T [ f (·)](齐次性) T [ f 1(·)+ f 2(·)] = T[ f 1(·)]+T[ f 2(·)] （可加性）②当动态系统满足下列三个条件时该系统为线性系统：)0(d )()(f t t t f =⎰∞∞-δ)(d )()(a f t a t t f =-⎰∞∞-δ?d )()4sin(91=-⎰-t t t δπ)0('d )()('f t t f t -=⎰∞∞-δ)0()1(d )()()()(n n n f t t f t -=⎰∞∞-δ4)2(2])2[(d d d )(')2(0022=--=--=-==∞∞-⎰t t t t tt t t δ)(1||1)()()(t a a at n n n δδ⋅=)(||1)(t a at δδ=)(||1)(00a t t a t at -=-δδ)0()()(f k k f k =∑∞-∞=δy (·) = y f (·) + y x (·) = T[{ f (·) }, {0}]+ T[ {0}，{x (0)}] (可分解性) T[{a f (·) }, {0}] = a T[{ f (·) }, {0}]T[{f 1(t ) + f 2(t ) }, {0}] = T[{ f 1 (·) }, {0}] + T[{ f 2 (·) }, {0}](零状态线性)T[{0},{a x 1(0) +b x 2(0)} ]= aT[{0},{x 1(0)}] +bT[{0},{x 2(0)}](零输入线性) 4.4时不变系统与时变系统T[{0}，f (t - t d )] = y f (t - t d)(时不变性质)直观判断方法：若f (·)前出现变系数，或有反转、展缩变换，则系统为时变系统。
微分方程符号解
ans =
C1*sinh(5^(1/2)*t)+C2*cosh(5^(1/2)*t)
例2 求解下列微分方程组的通解 (1)dy/dx=3y-2z; dz/dx=2y-z; >> [Y,Z]=dsolve('Dy=3*y-2*z','Dz=2*y-z','x')
Y= -exp(x)*(-C1-2*x*C1+2*x*C2) Z= -exp(x)*(-2*x*C1-C2+2*x*C2) (2)(d^2)x/dt^2+dy/dt-x=e^(t); (d^2)y/dt^2+dx/dt+y=0. >> [X,Y]=dsolve('D2y+Dy-x=exp(t)','D2y+Dx+y=0')
dx dt 3 x y z , dy x 5 y z, dt dz x y 3 z dt
( 6)
精品课件！
精品课件！
2．求下列常微分方程或方程组的特解：
(1)
()
y'' 3 y' 2 y 0; y(0) 2, y' (0) 3;
例3 求下列常微分方程或方程组的特解 (1) >> dsolve('(1+x^2)*D2y=2*x*Dy','y(0)=1','Dy(0)=3','x')
ans =
1+x^3+3*x
(2) >>[X,Y]=dsolve('Dx+Dy-x+3*y=exp(-t)1','Dx+2*x+Dy+y=exp(2*t)+t','x(0)=48/49','y(0)=95/98') X= 3/7*t-1/49+5/17*exp(2*t)+12/17*exp(-7/5*t)
信号与系统第5章
s a n 1 s
n 1
... a 1 s a 0
若m≥n （假分式）,可用多项式除法将象函数F(s)分 解为有理多项式P(s)与有理真分式之和。
F (s) P (s)
第5-9页
■
B0 (s) A(s)
©南京信息工程大学滨江学院
信号与系统
F (s) s 8 s 25 s 31 s 15
5.3
拉普拉斯逆变换
直接利用定义式求反变换---复变函数积分，比较困难。 通常的方法 （1）查表:直接利用拉普拉斯逆变换表 （2）利用性质 （3） 部分分式展开 -----结合 若象函数F(s)是s的有理分式，可写为
F (s) bm s
n m
b m 1 s
m 1
.... b1 s b 0
F (s) 1 e
sT
sT
e
2 sT
e
3 sT
＋)
特例：T(t) ←→ 1/(1 – e-sT)
第5-5页
■
©南京信息工程大学滨江学院
信号与系统
5.2
拉普拉斯变换性质
四、复频移（s域平移）特性
若f(t) ←→F(s) , Re[s]>0 , 且有复常数sa=a+ja, 则f(t)esat ←→ F(s-sa) , Re[s]>0+a 例1:已知因果信号f(t) 的象函数F(s)=
第5-1页
■
©南京信息工程大学滨江学院
信号与系统
5.1
拉普拉斯变换
四、常见函数的拉普拉斯变换
1、(t) ←→1，> -∞
’(t) ←→s，> -∞
2、(t)或1 ←→1/s ，> 0 3、指数函数e-s0t ←→
信号与系统 第二章
例2 描述某系统的微分方程为 y”(t)+5y’(t)+6y(t)= f”(t) + 2f’(t) + 3f(t) 求其冲激响应h(t)。 解： 选新变量y1(t) 它满足 y1”(t)+5y1’(t)+6y1 (t)=f(t) 设 h1(t)为冲激 响应，得到系统的冲激响应 h(t)= h1”(t) + 2 h1’(t) + 3 h1(t) 式1 求h1？同例1，得其冲激响应 h1 (t)=(e-2t -e-3t)ε(t) 再求h1”(t) 、 h1’(t) ，并带入 式1 ， 等到系统的冲激响应为 h(t)= δ(t)+（3e-2t -6e-3t)ε(t)
ˆ (t ) f
n
f (n)p(t n)
lim
0
ˆ (t ) f (t ) f
f ( ) (t ) d
2 .任意信号作用下的零状态响应
f (t)
根据h(t)的定义： 由时不变性：
LTI系统 零状态
yf(t) h ( t) h(t -τ) f (τ) h(t -τ)
对式(1)两端积分有
0
0
y' ' (t )dt 3 y' (t )dt 2 y(t )dt 2 (t )dt 6 (t )dt
0 0 0 0
0
0
0
0
由于积分在无穷小区间[0-，0+]进行的，且y(t)在t=0连续， 故 0 0 0 y(t )dt 0, 0 (t )dt 0
f (t ) f1 ( ) f 2 (t )d
为f1(t)与f2(t)的卷积积分，简称卷积；记为 f(t)= f1(t)*f2(t) 注意：积分是在虚设的变量τ下进行的，τ为积分变量， t为参变量。结果仍为t 的函数。
微分方程的多解与变号解
微分方程的多解与变号解
微分方程的多解与变号解是微积分的重要内容。
通过推导，可以知道一般地，一个非齐次线性微分方程的解可以是唯一的(特解)，也可以有多解，如果只有唯一的特解，则它的称为单解，如果有多解，则称为多解。
多解的存在是有条件的，它只出现在具有某种特殊条件的方程中。
一般地，一个非齐次线性微分方程有多解，当且仅当它具有部分恒定系数和部分给定函数。
即，当存在满足其中一些给定方程的连续解或连续解序列，其中某些解作为积分常数出现，多解存在。
即使多解存在，也可以将其转换为单解。
当微分方程有多解时，将其中的某些解作为积分常数出现，可以得到某个特解，即可以将多解变为单解，称为“变号解”。
变号解需要“类积分法”，其核心步骤是将同一非齐次线性微分方程系统多个解线组合分别当做积分因子，重新求解，最终变为单解。
因此，只要满足一些特殊条件，微分方程就可以具有多解，而当具有多解时，可以通过变号解将其变为单解，这是十分重要的技巧。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
微分方程的多解与变号解【摘要】：偏微分方程解的存在性与多重性是非线性分析的一个重要研究内容,有着广泛的背景,它来源于物理、生物工程、化学和医学等领域.近年来,许多学者对偏微分方程进行了研究,例如利用变分方法和临界点理论研究了各种Schrodinger方程解的存在性与多解性.这些研究都进一步促进了非线性分析的发展.本文利用变分方法、Morse 理论、临界点理论、拓扑度理论研究了几类偏微分方程的解与变号解.本文分为五章.第一章,我们介绍一些研究背景,国内外研究现状及本文的一些主要结果.第二章中,我们对经典椭圆方程Dirichlet边值问题进行了研究,其中Q(?)RN是具有光滑边界的有界区域,f∈C1(Ω×R1,R1)满足次临界增长条件|ft’(x,t)|≤C(1+|t|p-2),(x,t)∈Ω×R1,其中C0是一正常数,p∈(2,2’),如果N≥3,则2*=2N/(N-2)；如果N-1,2,则2*-∞.我们把拓扑度、临界群、不动点指数结合起来,得到它们之间的一些转化关系,给出一些假设条件,使得非线性项f在0点和∞共振和跨特征值.我们解决了仅仅使用不动点指数和拓扑度理论不能研究共振情形的问题,如文献[1](J.Math.Anal.Appl.314(2006)464-476).我们考虑在0点和∞都共振的情形,这是一般文章都没有考虑的情形.而且我们的共振条件去掉了[2](Math.Z.233(2000)655-677)中的有界性条件.我们还研究了单边共振的情形,这种共振条件又比我们常见的要弱,它去掉了极限存在的要求和增长性条件.我们得到的结论是：f在0点共振或跨特征值,在∞点共振或单边共振或跨特征值,那么上述边值问题都有七个非平凡解,其中两个正解、两个负解、三个变号解.在维数N=1时,我们可以进一步计算第七个解的临界群,从而得到第八个解,这个解是变号解,所以在维数为一的情形下,我们得到两个正解、两个负解、四个变号解.在维数N=1,方程是四阶时,就化为[1]中所研究的问题.我们不仅给出了零点和无穷远点共振的情形,而且当0点共振到奇特征值或跨奇特征值时,还可以通过计算临界群,得到七个非平凡解,其中三个变号解,这种情形利用[1]中的方法,不能考虑.而在0点和∞都共振到偶特征值时,得到六个非平凡解,这样的共振情形也是[1]中的方法所不能研究的.[1]中只能考虑在零点和无穷远点都是跨偶数个特征值的情形.从而我们推广了[1]的结论,而且得到比[1]中更好的结果.我们假设以下条件：(f1)f(x,t)t≥0,(x,t)∈Ω×R1；(f2)存在n0≥2满足λ。
0λn0+1,使得ft’(x,0)=λn0+1,xΩ.且存在δ0,使得f(x,t)t≤λn0+1t2,(x,t)∈Ω×[-δ,δ]；(f3)存在b0,使得|f(x,t)|b,(x,t)∈Ω×[-bc,bc],其中c=maxx∈Ωe(x),且e是边值问题：的解；(f4)存在满足λn1-1λn1的正整数n11及C10,α∈(0,1),使得(?)f(x,t)/t=λn1,对x∈Ω一致成立,|f(x,t)-λn1t|≤C1(1+|t|α),(x,t)∈Ω×R1,且(?)1/|t|2α(F(x,t)-1/2λn1t2)=(-1)n1+1∞对x∈Ω一致成立,其中F(x,t)=∫0tf(x,s)ds,(x,t)∈Ω×R1；(f5)存在n1,ε,R0,使得λ21-1+ε≤f(x,t)/t≤λ2n1,x∈Ω,|t|≥R,且(?)1/|t|(F(x,t)-1/2λn1t2)=-∞对x∈Ω一致成立；(f6)存在n1,ε,R0,使得λ2n1≤f(x,t)/t≤λ2n1+1-ε,x∈Ω,|t|≥R,且(?)1/|t|(F(x,t)-1/2λn1t2)=+∞对x∈Ω一致成立；(f7)存在n0≥2,使得λn0λn0+1且λn0ft’(x,0)λn0+1对x∈Ω一致成立；(f8)存在满足A2n1入2n1+1的n11,使得α(x)=(?)f(x,t)/t存在且λ2n1α(x)λ2n1+1对x∈Ω一致成立.定理2.1.1.设条件(f1)-(f3)成立,(f4),(f5),(f6),(f8)有一个成立,则边值问题(2.1.1)至少有七个非平凡解,其中两个正解,两个负解,三个变号解.定理2.1.2.设条件(f1),(f3)及(f7)成立,(f4),(f5),(f6),(f8)有一个成立,则边值问题(2.1.1)至少有七个非平凡解,其中两个正解,两个负解,三个变号解.在第三章中,我们考虑了RN上半线性椭圆方程-△u+u=f(x,x),u∈H1(RN)的解.在一些假设条件下,我们利用下降流线给出无穷多变号解的存在性.文献[3](Adv.Math.222(2009)2173-2195)中研究的是有界区域上的半线性椭圆方程的无穷多变号解,这里我们改进到无界区域.同时,我们也改进了文献[4](Comm.Math.Phys.55(1997)149-162)中只得到无穷多解,而没有确定出它们的符号的结论.本章中,我们给出以下条件：(A1)存在p∈(2,2*)及c0,使得|f(x,t)|≤c(|t|+|t|p-1),(x,t)∈RN×R1；(A2)存在α2,R0,使得αF(x,t)≤t,(x,t),(x,t)∈RN×R1,(?)F(x,t)0,其中F(x,t)=∫0t(x,s)ds,(x,t)∈RN×R1；(A3)limt→0f(x,t)/t=0在RN上一致成立；(A4)f(gx,t)=f(x,t),9∈O(N),(x,t)∈RN×R1；(A5)f(x,-t)=-f(x,t),tf(x,t)≥0,(x,t)∈RN×R1.那么我们有如下结果：定理3.1.1.假设(A1)-(A5)成立,则方程(3.1.1)有无穷多变号解.第四章中,我们研究Schrodinger-Poisson系统正径向解的存在性,其中f∈C(R1,R1).我们假设f满足limt→+∞f(t)/tp+∞,利用变分方法,当λ和p在不同范围时,给出一些解的存在性与不存在性结果.我们把[5](J.Funct.Anal.237(2006)655-674)中非线性项f(u)-up的结果推广到一般的非线性项.f,而且改进了文献[6](Ann.I.Poincare-AN,27(2010)779-791)中只对很小的参数得到一个正径向解的结论.事实上,我们假设以下条件：(H1)存在p∈(1,5),使得limsuPt→+∞f(t)/tp+∞；(H2)limt→0f(t)/t=0；(H3)f(t)t≥4F(t),t∈R1；(H4)存在q∈(2,5),使得liminft→+∞f(t)/tq0；(H5)limt→+∞f(t)/t=+∞.其中F(t)=∫0tf(s)ds.那么我们有下列结论.定理 4.1.1.如果条件(H1)中p∈(1,2),且(H2)及(H5)成立,则存在λ00,使得对所有的入∈(0,λ0),系统(4.1.1)至少有两个正径向解.定理4.1.2.如果条件(H1)中p∈(3,5),且(H2)-(H4)成立,则对所有的λ0,系统(4.1.1)至少有一个正径向解.定理4.1.3.如果条件(H1)中p∈[2,3],且(H2)及(H5)成立,则存在A00,使得对所有λ∈(0,入0),系统(4.1.1)至少有一个正径向解.定理4.1.4.如果条件(H1)中p∈(1,2],且(H2)成立,则存在λ00,使得对所有A∈(λ0,∞),系统(4.1.1)没有正解.我们的主要结果可以用下图表示.第五章中,我们研究了广义Kadomtsev-Petviashvili方程wt+wxxx+(f(w))x=Dx1-wyy的基态孤立波,其中我们在一些假设条件下,利用变分方法给出非平凡基态孤立波的存在性.我们去掉了[7](Appl.Math.Lett.15(2002)35-39)和[8](MinimaxTheorems,1996)中的条件：存在υ0∈Y:={gx:g∈C0∞(R2)},使得lims→+∞F(sυ0)/s2→+∞.这是研究广义Kadomtsev-Petviashvili方程经常会假设的条件.并且在本章中,相应的泛函不满足(PS)条件和(c)条件.我们假设(B1)f∈C(R1,R1),f(0)=0,且对某一p∈(3,6),(B2)limt→0,f(t)/t=0;(B3)存在θ≥1,使得θG(t)≥G(st),t∈R1,s∈[0,1],其中G(t)=f(t)t-2F(t).定理 5.1.1.假设(B1)-(B3)成立,那么问题(5.1.1)有一基态孤立波.另外,我们还研究了变系数的广义Kadomtsev-Petviashvili方程的基态孤立波：(-uxx+Dx-2uyy+a(x,y)u-f(u))x=0,(5.4.1)其中(x,y)∈R2.我们假设a∈C(R2,R1)且存在正数α,β,使得Oα≤a(x,y)≤β∞.且a(x,y)满足以下条件：a(x,y)a∞=(?)a(x,y)∞,(x,y)∈R2(a)这里我们同样也去掉了[9](J.Math.Anal.Appl.361(2010)48-58)中研究变系数p-Laplace方程所作的周期性假设.从而推广了[9]中的方法.则我们有下列结论：定理5.4.1.假设对p∈(3,4],(B1)-(B3)成立,且条件(a)成立,那么问题(5.4.1)有一基态孤立波.【关键词】：临界群拓扑度变号解Schr(o|¨)dinger-Poisson 系统广义Kadomtsev-Petviashvili方程基态孤立波【学位授予单位】：山西大学【学位级别】：博士【学位授予年份】：2011【分类号】：O175.2【目录】：中文摘要6-11ABSTRACT11-17第一章引言17-27第二章半线性椭圆方程的多解与变号解27-51§2.1临界群、拓扑度以及不动点指数的转化27-40§2.2变号解和解的存在性与多重性40-46§2.3维数为一的情形46-51第三章R~N上半线性椭圆方程的无穷多变号解51-63§3.1准备工作51-53§3.2一些必要的引理53-58§3.3无穷多变号解的存在性58-63第四章Schrodinger-Poisson系统的正解63-79§4.1主要结果63-64§4.2一些必要的引理64-75§4.3主要定理的证明75-79第五章广义Kadomtsev-Petviashvili方程的基态孤立波79-97§5.1广义Kadomtsev-Petviashvili方程孤立波的转化79-80§5.2一些必要的引理80-85§5.3常系数情形下主要定理的证明85-88§5.4变系数的情形88-97参考文献97-102攻读博士学位期间的主要研究成果102-103个人简介及联系方式103-104致谢104-106 本论文购买请联系页眉网站。