24–25 年英才班期末试题 (严军)
1. | 船行水上, 初速度为 v0, 水的阻力为 f=kv. 问船何时停? 总路程为多少? |
2. | 设连续函数 f(x) 在区间 (−∞,+∞) 上有界, 求证: 方程y′+y=f(x)在区间 (−∞,+∞) 内有且仅有一个有界解, 并求出这个有界解. 并进而证明: 当 f(x) 是以 ω 为周期的周期函数时, 这个有界解同时也是以 ω 为周期的周期函数. |
3. | 设 y 为关于 x 的连续可微函数, 已知 y′+a(x)y⩽0,∀x>0, 求证: y(x)⩽y(0)e−∫0xa(s)ds,∀x⩾0. |
4. | 设 x(t),y(t),z(t)∈C1(R), σ,ρ,β>0,ρ<1, 则如下方程: dtd⎝⎛xyz⎠⎞=⎝⎛−σρ0σ−1000−β⎠⎞⎝⎛xyz⎠⎞+⎝⎛0−xyxz⎠⎞的解渐近稳定. (提示: 考虑 Lyapunov 函数 V=x2+y2+z2.) |
5. | 求解偏微分方程(z−y)2∂x∂z−x∂y∂z=0过 x=0,y=z=s 的解. |
6. | 求波方程{x′′+λx=0x(0)=x(1), x′(0)=x′(1)的特征值和对应的特征函数. |
23–24 年英才班期末试题 (严军)
1. | (10 分) 设 y 为关于 x 的连续可微函数, 已知 y′+a(x)y⩽0,∀x>0, 求证: y(x)⩽y(0)e−∫0xa(s)ds,∀x⩾0. |
2. | (15 分) 设 φ1,⋯,φn 为一元 n 维实向量值连续函数, 且其 Wronsky 行列式 L 恒不等于 0. 求证: 存在一个 n×n 的矩阵值函数 A(t), 使得关于 n 维连续可微的向量值函数 y 的常微分方程dtdy=A(t)y的一组基本解为 φ1,⋯,φn. |
3. | (15 分) 设 x(t)∈C2(R), Q(t) 为 R 上的恒正函数, 则方程dt2d2x+Q(t)x=0的任何解都有零点. |
4. | (15 分) 设 x(t),y(t),z(t)∈C1(R), σ,ρ,β>0,ρ<1, 则如下方程: dtd⎝⎛xyz⎠⎞=⎝⎛−σρ0σ−1000−β⎠⎞⎝⎛xyz⎠⎞+⎝⎛0−xyxz⎠⎞的解全局渐近稳定. (提示: 考虑 Lyapunov 函数 V=x2+y2+z2.) |
5. | (15 分) 设 a(x),b(x) 在区间 I 上连续, y∈C1(R), 求证: 方程y′=a(x)y+b(x)的最大存在区间为 I. |
6. | (15 分) 求波方程{x′′+λx=0x(0)=x(1), x′(0)=x′(1)的特征值和对应的特征函数. |
7. | (15 分) 设 D 为 R2 上的单连通区域, h(x,y),P(x,y),Q(x,y)∈C1(D), 且 h 为恒正函数, ∂x∂(hP)+∂y∂(hQ)=0. 求证: 方程⎩⎨⎧dtdx=P(x,y)dtdy=Q(x,y)的解没有闭轨. |
22–23 年英才班期末试题 (严军)
1. | (10 分) 设连续函数 f(x) 在区间 (−∞,+∞) 上有界, 求证: 方程y′+y=f(x)在区间 (−∞,+∞) 内有且仅有一个有界解, 并求出这个有界解. 并进而证明: 当 f(x) 是以 ω 为周期的周期函数时, 这个有界解同时也是以 ω 为周期的周期函数. |
2. | (10 分) 求证: 齐次方程 P(x,y)dx+Q(x,y)dy 有如下积分因子: μ=xP+yQ1. |
3. | (20 分) 设 f 在 [0,+∞) 可微, f(0)=0, 且有常数 k>0 使得: ∣f′(x)∣⩽k∣f(x)∣,∀x∈[0,+∞).证明: f≡0. |
4. | (30 分) (1) 画出单摆方程dt2d2x+a2sinx=0,(a=0)的相图并判断讨论不动点的稳定性. (2) 判断耗散的单摆方程dt2d2x+a2sinx+bdtdx=0,(a=0,b>0)奇点的稳定性. |
5. | (20 分) 设 a(t),b(t) 在 0⩽t⩽1 上连续, 求证: x′′+a(t)x′+b(t)x=0满足初值条件 x(0)=0 的任意两个解都是线性相关的. |
20–21 年英才班期末试题 (袁小平)
1. | (10 分) 举例说明 Picard 存在唯一性定理中 Lipschitz 条件不可去. |
2. | (10 分) 解方程⎩⎨⎧dtdx=−3x−y,dtdy=4x+y. |
3. | (10 分) 求 dt2d2x=tx 的幂级数解. |
4. | (10 分) 证明 Hamilton 系统的极小值点为 Lyapunov 稳定的但不为渐近稳定的. |
5. | (10 分) 对自治系统 x˙=f(x) 的解 φt(ξ), 记 γ={φt(ξ)∣t∈R}. 令 ω(ξ) 为其 ω 极限集. 若 γ 为有界集, 证明: ω(ξ) 为非空闭集. |
6. | (15 分) 证明系统⎩⎨⎧dtdx=−y+x(x2+y2−1)+1001x3dtdy=x+y(x2+y2−1)有稳定的极限环. |
7. | (15 分) 证明边值问题⎩⎨⎧−dx2d2y=λy,0<x<1y(0)=0, y′(1)+y(1)=0的特征值均为正实数且有可列个, 记它们为 λ1<λ2<⋯<λk<⋯. 证明 λk∈((k−21)π,kπ). |
8. | (20 分) 证明 x˙=f(x) 的周期解 (非平衡点) 不可能为 Lyapunov 渐近稳定的. |
23–24 年期中测验试卷 (林伟)
一、 计算题 (本大题满分 40 分) 本大题共有 4 题, 每题满分 10 分.
1. | 求解微分方程 y′′+y′+y=e2xcosx. 解. 先求通解, 原方程的特征方程为 λ2+λ+1=0, 有特征根 λ1,2=2−1±3i, 因此原方程的实值通解为y=e−2x(C1cos(23x)+C2sin(23x)).再求特解, (D2+D+1)y=e2xcosx, 得: y=D2+D+11e2xcosx.转化得:y=D2+D+11e(i+2)x,由性质 4 得:y=D2+D+11e(i+2)x=(i+2)2+(i+2)+11e(i+2)x=5i+61e(i+2)x,从而有y=5i+61e(i+2)x=e2x(−615i+616)(cosx+isinx)=e2x(616sinx−615cosx)i+e2x(616cosx+615sinx).取实部得: y∗=e2x(616cosx+615sinx),也即原方程的解为y=e−2x(C1cos(23x)+C2sin(23x))+e2x(616cosx+615sinx).□ |
2. | 试找到所有可能的 λ (λ∈R), 使得微分方程 ⎩⎨⎧x′′(t)+λx(t)=0,x(0)=x(l)=0, l>0 存在非平凡解. 解. 原方程的特征方程可写为 r2+λ=0, 其解为 r=±−λ, 我们分三种情况讨论: (1) 当 λ<0 时, 方程的通解可写为: x(t)=C1e−λt+C2e−−λt,C1,C2∈R.要使其满足边界条件, 需要满足: ⎩⎨⎧C1+C2=0C1e−λl+C2e−−λl=0.由于 ∣∣1e−λl1e−−λl∣∣=0, 只有当 C1=C2=0, 所以此情况没有非平凡解; (2) 当 λ=0 时, 方程的通解可写为: x(t)=C1+C2t,C1,C2∈R.此种情况为满足边界条件也只能有 C1=C2=0, 无非平凡解; (3) 当 λ>0 时, 方程的通解可写为: x(t)=C1cosλt+C2sinλt,C1,C2∈R.由 x(0)=0 知 C1=0, 再由 x(l)=C2sinλl=0 可知, 为了使 C2=0, 我们必须有: λ=λk=l2k2π2 (k=1,2,…). 综上所述, λ 的取值为: l2π2,l24π2,l29π2,⋯,l2k2π2,⋯. |
3. | 求解初值问题: ⎩⎨⎧dtd⎝⎛xy⎠⎞⎝⎛x(0)y(0)⎠⎞=⎝⎛cos2t1+sintcost−1+sintcostsin2t⎠⎞⎝⎛xy⎠⎞+⎝⎛costsint⎠⎞,=⎝⎛01⎠⎞. 解. 事实上, 不难验证 ⎝⎛etcostetsint⎠⎞,⎝⎛−sintcost⎠⎞是上述齐次线性微分方程组 −∞<t<∞ 上的两个解; 而且它们的朗斯基行列式 W(t) 在 t=0 处的值为 W(0)=∣∣1001∣∣=1=0.故通解为 ⎝⎛xy⎠⎞=C1⎝⎛etcostetsint⎠⎞+C2⎝⎛−sintcost⎠⎞.这代表相应齐次线性微分方程组有一个基解矩阵 Φ(t)=⎝⎛etcostetsint−sintcost⎠⎞.容易求出 Φ−1(t)=⎝⎛e−tcost−sinte−tsintcost⎠⎞,Φ−1(0)=⎝⎛1001⎠⎞.利用公式, 就得到所求初值问题的解为 ⎝⎛xy⎠⎞=Φ(t)⎣⎡⎝⎛01⎠⎞+∫0t⎝⎛e−scoss−sinse−ssinscoss⎠⎞⎝⎛cosssins⎠⎞ds⎦⎤=⎝⎛(et−1)cost−sint(et−1)sint+cost⎠⎞. |
4. | 求解齐次线性方程组dtdx=⎝⎛3−441−1−800−2⎠⎞x. 解. 矩阵的 Jordan 标准型为: ⎝⎛−200010011⎠⎞, λ1=−2 对应的特征向量可取为 η1=⎝⎛001⎠⎞.对于二重特征根 λ2=1, 可以算出 (A−λ2E)2=⎝⎛2−441−2−800−3⎠⎞2=⎝⎛00280044009⎠⎞.因此, 方程 (A−λ2E)2r=0 有两个线性无关的解为 r10=⎝⎛11−70⎠⎞,r20=⎝⎛3−620⎠⎞.把它们分别代入, 并注意 ni=2, 就可得到 r11=⎝⎛2−441−2−800−3⎠⎞⎝⎛11−70⎠⎞=⎝⎛15−30100⎠⎞和 r21=⎝⎛2−441−2−800−3⎠⎞⎝⎛3−620⎠⎞=⎝⎛000⎠⎞.以上结果不唯一, 只要符合广义特征向量要求即可, 由此我们可得到一个基解矩阵 Φ(t)=⎝⎛00e−2t(11+15t)et(−7−30t)et100tet3et−6et20et⎠⎞x=Φ(t)c其中 c 为任意常数列向量. 这个通解也可以改写成下面更清晰的形式 x=C1⎝⎛001⎠⎞e−2t+C2⎝⎛11+15t−7−30t100t⎠⎞et+C3⎝⎛3−620⎠⎞et,其中 C1,C2 和 C3 是任意常数. |
二、 解答题 (本大题满分 60 分) 本大题共有 4 题, 解答下列各题并写出必要的步骤.
5. (本小题满分 10 分) 本题共有 2 个小题, 第 1 个小题满分 5 分, 第 2 个小题满分 4 分.
考虑自治系统 x˙=Ax,A∈Rn×n,
(1) | 若矩阵 A 是正定的, 也就是对于任意非零向量 x 都有 xTAx>0, 证明 ∥x∥A=xTAx 是一种范数 (注意, 正定矩阵不一定是对称矩阵) ; |
(2) | 若矩阵 A 是严格对角占优的, 也就是 ∣aii∣>j=1,…,n,j=i∑∣aij∣,i=1,…,n, 求出方程组的所有常值解. |
解.
(1) | (i) 正定性 (1 分) : 可由 A 的正定性直接推出; (ii) 齐次性 (1 分) : ∀x∈Rn,α∈R,∥αx∥A=(αx)TA(αx)=∣α∣⋅∥x∥A; (iii) 三角不等式: 我们需要说明, 对于任意 x,y∈Rn,∥x+y∥A≤∥x∥A+∥y∥A (1 分) . 当 A 为对称正定矩阵时 (1 分) , 可考虑采用谱分解或 Cholesky 分解等方法证明. 设 A=LTL,L 是下三角矩阵, 则 ∥x∥A=(Lx)T(Lx)=∥Lx∥2, 可推出 ∥Lx+Ly∥2≤∥Lx∥2+∥Ly∥2, 所以三角不等式得证; 当 A 不是对称阵时 (1 分) , 由于∥x∥A=xTAx=21xTAx+21xTATx=∥x∥2A+AT,由于 2A+AT 是对称正定矩阵, 所以此种情况也满足三角不等式. |
(2) | 证明严格对角占优矩阵非异即可 (例如戈式圆盘定理证明特征值严格位于虚轴右侧, 或直接证明 Ax=0 仅有 0 解, 都是常规的高等代数方法) . |
6. (本题满分 16 分) 本题共有 3 个小题, 第 1 个小题满分 4 分, 第 2 个小题满分 4 分, 第 3 个小题满分 6 分.
解的存在唯一定理是常微分方程中最基本的定理, 其最古典的证明中运用了逐次逼近法等数学思想. 现考虑矩形区域 R=[t0−a,t0+a]×[x0−b,x0+b] 上的柯西问题 ⎩⎨⎧dtdx=f(t,x)x(t0)=x0, 假设定理保证了微分方程在 I=[t0−h,t0+h] 上解的存在唯一,
(1) | 叙述定理中 f(t,x) 应在矩形区域 R 上满足的条件, 进一步验证当 f(t,x)={xln∣x∣,0,x=0x=0 时是否在包含 (t,0) 的矩阵区域上满足该条件, 讨论此时微分方程的解是否存在唯一; |
(2) | 若给定矩形区域 R1 和 R2, 定理保证了微分方程在 I1 和 I2 上的解存在唯一, 当 R1R2 时, 是否能推得 I1⊂I2? 请给出证明或举出反例; |
(3) | 运用逐次逼近法证明 x(t)=esint+∫0tarctan(t⋅τ)x(τ)dτ 的解在 t≥0 时存在唯一. |
解.
(1) | f(t,x) 应在矩形区域 R 上连续 (1 分) 且对 x 满足 Lipschitz 条件 (1 分) . 当 f(t,x)={xln∣x∣,0,x=0x=0 时, 在包含 (t,0) 的矩阵区域上不满足 Lipschitz 条件 (1 分) , 因为∂x∂f=1+ln∣x∣→−∞(x→0),可以验证不存在这样的 Lipschitz 常数 (1 分) . 但是, 原方程的解是存在唯一的 (1 分) , 可以由分离变量法求得原方程的解: ∫xln∣x∣dx=ln(ln∣x∣)(1 分) . (事实上是因为该方程满足 Osgood 条件, 因此存在唯一, 也可以通过分离变量法说明, 但因会漏解, 故严谨性略有不足. ) |
(2) | 不能推得 (2 分) , 可以举出如下反例 (2 分) (反例合理即可) : 考虑 ⎩⎨⎧dtdx=t2+x2x(0)=0, 考虑矩形区域 R1=[−1,1]×[−1,1] 和矩形区域 R2=[−2,2]×[−2,2], 我们易求得h1=min(1,21)=21,h2=min(2,82)=41,则I1=[−21,21],I2=[−41,41],此时 R1R2, 但 I2I1. |
(3) | 证明: (原试卷中误将 esint 写成 esinx, 只要作答了相应部分都给分. ) 存在性 (3 分) : 我们可以发现, esint 和 arctan(t⋅τ) 在定义域均为连续有界函数, 即esint≤e,arctan(t⋅τ)≤2π,构造序列⎩⎨⎧ϕ0(t)=esintϕn(t)=esint+∫0tarctan(t⋅τ)ϕn−1(τ)dτ,我们发现: ∣ϕ1(t)−ϕ0(t)∣=∣∣∫0tarctan(t⋅τ)ϕ0(τ)dτ∣∣≤∫0t∣arctan(t⋅τ)∣⋅∣esint∣dτ≤2π⋅e⋅t;∣ϕ2(t)−ϕ1(t)∣=∣∣∫0tarctan(t⋅τ)(ϕ1(τ)−ϕ0(τ))dτ∣∣≤∫0t∣arctan(t⋅τ)∣⋅∣ϕ1(τ)−ϕ0(τ)∣dτ≤2π∫0t∣ϕ1(τ)−ϕ0(τ)∣dτ≤(2π)2⋅e⋅2!t2;依此类推, ∣ϕn(t)−ϕn−1(t)∣=∣∣∫0tarctan(t⋅τ)(ϕn−1(τ)−ϕn−2(τ))dτ∣∣≤…≤(2π)n⋅e⋅n!tn.考虑级数 ϕ0(t)+n=1∑∞(ϕn(t)−ϕn−1(t)), 因为ϕ0(t)+n=1∑∞(ϕn(t)−ϕn−1(t))≤e+n=1∑∞(2π)n⋅e⋅n!tn=e⋅n=0∑∞(2πt)n⋅n!1=e2πt+1<+∞,所以 {ϕn(t)} 在 [0,t] 一致收敛. 所以对迭代序列两侧取极限, 可保证解的存在性. 唯一性 (3 分) : 设 ϕ1(t) 和 ϕ2(t) 是微分方程两个不同的解, 我们有: ⎩⎨⎧ϕ1(t)=esint+∫0tarctan(t⋅τ)ϕ1(τ)dτϕ2(t)=esint+∫0tarctan(t⋅τ)ϕ2(τ)dτ,由此我们有 ∣ϕ2(t)−ϕ1(t)∣=∣∣∫0tarctan(t⋅τ)(ϕ2(τ)−ϕ1(τ))dτ∣∣≤2πt∣ϕ2(t)−ϕ1(t)∣,不断在积分号内迭代, 可得 ∣ϕ2(t)−ϕ1(t)∣≤(2πt)nn!1∣ϕ2(t)−ϕ1(t)∣,由此可推得 ∣ϕ2(t)−ϕ1(t)∣=0, 从而唯一性得证. |
7. (本题满分 16 分) 本题共有 3 个小题, 第 1 个小题满分 4 分, 第 2 个小题满分 6 分, 第 3 个小题满分 6 分.
已知光滑映射 Q:R→O(3) 满足 Q(0)=I3 为三阶单位矩阵, 这里 O(3) 为 R3×3 上的正交矩阵组成的集合. 若一个在 R3 中运动质点轨迹满足 x(t)=Q(t)b, 其中 b 为 R3 中的常值向量.
(1) | 证明质点的运动轨迹满足线性微分方程 x˙=A(t)x, 且系数矩阵 A(t) 是反对称阵. |
(2) | 若 (1) 中的系数矩阵 A(t)≡A, 证明该质点在 R3 中绕一根过原点的定轴旋转, 即存在非零向量 u∈R3, 使得 (a) 质点运动速度 x˙ 与 u 垂直, (b) ∥x−u∥2 关于 t 不变. |
(3) | 若 (1) 中 A(t)≡⎣⎡0−3−4300400⎦⎤, 试确定光滑映射 Q(t). |
解.
(1) | 由于 Q(t) 是正交矩阵, 对于任意 t, 我们由正交矩阵的性质有: Q(t)QT(t)=I3, 两边对 t 求导有: Q˙(t)QT(t)+Q(t)Q˙T(t)=O,于是Q˙(t)QT(t)=−Q(t)Q˙T(t).由 x(t)=Q(t)b, 求导得 x˙(t)=Q˙(t)b, 由此得x˙(t)=Q˙(t)QT(t)b=−Q(t)Q˙T(t)b,所以A(t)=Q˙(t)QT(t)=−Q(t)Q˙T(t),所以 A(t) 为反对称矩阵. |
(2) | (a) 设A=⎝⎛0c−b−c0ab−a0⎠⎞,x=⎝⎛x1x2x3⎠⎞,u=⎝⎛u1u2u3⎠⎞,若 uT⋅x˙=0, uT⋅Ax=0, uT⋅Ax=(u1 u2 u3)⋅⎝⎛−cx2+bx3cx1−ax3−bx1+ax2⎠⎞=u1(−cx2+bx3)+u2(cx1−ax3)+u3(−bx1+ax2)=x1(cu2−bu3)+x2(−cu1+au3)+x3(bu1−au2)=0.由 x1,x2,x3 的任意性可知⎩⎨⎧cu2−bu3=0−cu1+au3=0bu1−au2=0,解得 u=⎝⎛abc⎠⎞, 此时 x˙ 关于 u 垂直. (b) 下面验证 (a) 中的 u 关于 t 不变, 等价于说明 ∥x−u∥22 对 t 求导数为 0. ∂t∂∥x−u∥22=∂t∂(xTx−2uTx+uTu)=2xTx˙−2uTx˙=2(xTAx+uTAx),因为 xTAx=uTAx=0, 所以 ∂t∂∥x−u∥22=0. |
(3) | 由矩阵指数求解可得 x(t)=eAtc, 满足初值条件 x(0)=Q(0)b=b, 则 c=b. 所以Q(t)=eAt=⎝⎛0−4/315i/43/41−5i/43/41⎠⎞⋅⎝⎛1000e−5it000e5it⎠⎞⋅⎝⎛0−4/315i/43/41−5i/43/41⎠⎞−1=⎝⎛5cos5t−3sin5t/5−4sin5t/53sin5t/59cos25t/25+16/2512cos5t/25−12/254sin5t/512cos5t/25−12/2516cos5t/25+9/25⎠⎞. |
8. (本题满分 18 分) 本题共有 4 个小题, 第 1 小题满分 5 分, 第 2 小题满分 4 分, 第 3 小题满分 4 分, 第 4 小题满分 5 分.
考虑 n 维线性周期常微分方程: x˙=A(t)x, 其中 A(t) 是 R 上以 T>0 为周期的连续矩阵值函数.
(1) | 设上述方程的一个基本解方阵为 Φ(t), 证明存在一个可逆的常数矩阵 B, 使得 Φ(t+T)≡Φ(t)B 对 t∈R 恒成立. |
(2) | 设 B 的特征值为 ρ1,ρ2,…,ρn, 则称这 n 个值为原方程的 Floquet 特征值. 请证明上述定义是良定的, 也即一个线性周期常微分方程的 Floquet 特征值与基本解方阵的选取无关. |
(3) | 证明: 对于任意一个 Floquet 特征值 ρ, 存在原方程的一个解 x=ϕ(t), 使得 ϕ(t+T)=ρϕ(t), 并以此分析原方程的 Floquet 特征值两两不同的情况下所有解都有界的充要条件. |
(4) | 对于一个非线性的自治系统 x˙=f(x), f=(f1,…,fn),x∈Rn, 假设其有 T-周期解 x=ψ(t). 考虑周期解 ψ(t) 的摄动 x=ψ+v (∥v∥≪∥ψ∥, 且 x 满足 x˙=f(x)). 那么 v 满足系统的变分方程 v˙=(∂xj∂fi)1≤i,j≤n∣∣ψ(t)v, 其中 (∂xj∂fi)1≤i,j≤n∣∣ψ(t) 为原系统的 Jacobi 矩阵, 请证明此线性方程有一个 Floquet 特征值为 1. |
证明.
(1) | 先考虑 Φ(t+T), 记为 Y(t), 则有Y′(t)=Φ′(t+T)=A(t+T)Φ(t+T)=A(t)Φ(t+T)=A(t)Y(t),因此 Φ(t+T) 依然为原方程的一个基解矩阵, 由基解矩阵性质知, 存在一个可逆的常数矩阵 B, 使得 Φ(t+T)≡Φ(t)B 对 t∈R 恒成立. |
(2) | 假设原方程的另一基解矩阵为 Ψ(t), 则存在可逆常矩阵 C 使得 Ψ(t)=Φ(t)C, 从而由 (1) 知存在可逆常数矩阵 Bˉ, 使得 Ψ(t+T)≡Ψ(t)Bˉ, 那么有: Bˉ=Ψ−1(t)Ψ(t+T)=C−1Φ−1(t)Φ(t+T)C=C−1BC,也即 Bˉ 与 B 相似, 自然拥有相同的特征值. 这代表 Floquet 特征值的定义是良定的. |
(3) | 设 b 为 B 对应于特征值 ρ 的特征向量, 考虑原方程的解 ϕ(t)=Φ(t)b, 我们有ϕ(t+T)=Φ(t+T)b=Φ(t)Bb=Φ(t)ρb=ρϕ(t),ϕ(t−T)=ρ1ϕ(t).由此即知, 若 Floquet 特征值两两不同, 那么原方程所有解有界的充要条件是 ∣ρ∣=1,∀ρ (∣ρ∣≤1 也得全分) . |
(4) | 对于摄动 x=ϕ+v, 代回原方程有v˙=f(ϕ+v)−f(ϕ)=(∂xj∂fi)i,jv+o(v).忽略高阶项的意义下, v 满足线性周期常微分方程 v˙=(∂xj∂fi)i,jv. 注意到当 v=ϕ′(t) 时 (线性 ODE 的求解仅与 “方向” 相关, 在这个解上除以一个足够大的数即可使得 ∥v∥≪∥ϕ∥ 成立) , 有下式成立: v˙=ϕ′′(t)=dtdf(ϕ(t))=(∂xj∂fi)i,jϕ′(t)=(∂xj∂fi)i,jv.这代表周期解 ϕ′(t) 为该 (变分/线性) 方程的解, 从而这个方程有一个 Floquet 特征值为 1, 证毕. |
23–24 年平时测验试卷 (林伟)
一、 解下列各微分方程 (组) (本大题满分 60 分) 本大题共有 6 题, 每题满分 10 分.
1. | dt4d4x=x+a3t3+a2t2+a1t+a0, 其中 a0,a1,a2,a3∈R. |
2. | (axy)dx+(x2+y2−1)dy=0, 其中 a∈R,a=0. |
3. | ⎩⎨⎧dtdx=yt(x2+y2),dtdy=xt(x2+y2). |
4. | dt3d3x+t1dt2d2x−t24dtdx=t3. |
5. | x=t(dtdx+1+(dtdx)2). |
6. | x+tdtdx=t−1tx−3+tan(t−1tx−3t). [答案未必表示为 x=x(t) 的显式形式.] |
二、 解答题 (本大题满分 40 分) 本大题共有 3 题, 解答下列各题并写出必要的步骤.
7. | (本题满分 12 分) 本题共有 2 个小题, 第 1 小题满分 6 分, 第 2 小题满分 6 分. 设 g(t) 是以 T∈R+ 为周期的连续函数, 考虑以下非齐次线性方程: dtdx=−x+g(t).(1) 证明: 若这个方程有解 φ(t), 且满足 φ(T)=φ(0), 则 φ(t) 是以 T 为周期的周期函数; (2) 试给出一个使 (1) 中的周期解 φ(t) 唯一存在的充分条件, 并说明理由. |
8. | (本题满分 14 分) 本题共有 2 个小题, 第 1 小题满分 6 分, 第 2 小题满分 8 分. 考虑微分方程 dtdx+ax=f(t), 其中 a∈R, f(t) 为 [0,+∞) 上的连续函数, 试证明: (1) 当 a=1, 且 f(t) 有界时, 微分方程一切解 x=x(t) 均在 [0,+∞) 上有界; (2) 当 a>0, 且 t→+∞limf(t)=b∈R 时, 微分方程一切解 x=x(t) 均有 t→+∞limx(t)=ab. |
9. | (本题满分 14 分) 本题共有 3 个小题, 第 1 小题满分 4 分, 第 2 小题满分 4 分, 第 3 小题满分 6 分. 微分方程的具体形式通常反映了系统的结构特征. 例如, 单变量一阶线性常微分方程 dtdx=ax,a∈R 与系统结构有如下对应关系, (i) | a>0⇔x◯←x◯;
| (ii) | a<0⇔x◯ ⊢x◯;
| (iii) | a=0⇔x◯ x◯. |
这样的关系可以自然地推广至由多变量线性常微分方程组描述的系统. 现考虑方程组⎩⎨⎧dtdx=ay,dtdy=bx+cy,(a,b,c∈R, ab=0).(1) 证明: 若方程组的一切解 (x(t),y(t)) 均满足 t→+∞lim(x(t),y(t))=(0,0), 则系统必有 “y◯⊢y◯” 结构; (2) 证明: 若方程组除零解外的一切解在 t∈(0,+∞) 均无界, 则系统必有 “y◯←y◯” 结构; (3) 画出所有能使方程组存在非常值周期解的系统结构, 并说明理由. |
22–23 年期末考试试卷 (林伟)
考生注意:
1. | 这份试卷共 9 道试题, 满分 100 分, 共 8 页. 考试时间: 2022 年 12 月 28 日 8:30~10:30. |
2. | 请在自备的纸张上, 依次写清题号 (阿拉伯数字) 并答题. 考试结束后, 请扫描答卷并上传 “姓名+学号+学院.pdf” 文件. |
3. | 线上考试, 请考生依旧秉持诚实守信宗旨, 谨守考试纪律, 摒弃考试作弊. 学生如有违反学校考试纪律的行为, 将按《复旦大学学生纪律处分条例》规定予以严肃处理. |
4. | 本课程学期总评分 (满分 100 分) 由平时成绩 (满分 50 分, 包含期中, 作业以及小论文成绩) 与本试卷成绩的 50% 组合计算 (方案 A) ; 若考生确实仅选择计算本试卷成绩的 90% 与小论文成绩 (满分 10 分) 作为学期总评分 (方案 B) , 请在答卷最开始注明 “选择方案 B”. |
一、 填空题 (本大题满分 27 分) 本大题共有 3 题, 将答案直接写在题后的横线上或者绘制在对应的坐标系中, 每题满分 9 分.
1. | 考虑矩阵值函数: Ψ(t)=⎝⎛e−tsinte−tcoste−tcost−e−tsint⎠⎞ 和 Φ(t)=⎝⎛e−tcoste−tsinte−tsinte−tcost⎠⎞. 当上述矩阵值函数是二维自治系统 [∗]x˙=Ax 的基本解方阵时, 分别计算 ∥A∥, det[eA⊤+5A] 的值, 并求二次型 V(x)=⟨x,Hx⟩ 的具体形式, 使其沿着上述自治系统的全导数 dtdV∣∣[∗](x)=−∥x∥2 (以上的 ∥⋅∥ 表示 R2 中的向量 2-范数或由该向量范数诱导的矩阵范数) : . |
2. | 设常数 aij (i,j=1,2) 中至少有两个不为 0, 且 f 在一定的区域中为足够光滑的函数. 若方程组 x˙=a11x+a12y+x2,y˙=−a21x−a22y 存在形如 f(x)=C 的首次积分, 则该微分方程奇点 O(0,0) 的类型是: . |
3. | 考虑极坐标 (r,θ) 下的平面系统: r˙=100r(−6+r2),θ˙=1+100c(−6+r2), 其中 c 为参数. 试在对应的欧氏坐标系 x–O–y 中, 分别绘制参数 c<0,c=0,c>0 三种情形下系统轨线的大致图形 (x=rcosθ,y=rsinθ): |
二、 辨析题 (本大题满分 17 分) 本大题共有 2 题, 每题都给出了若干个结论, 分别判断这些结论是否正确, 正确的在结论前的方括号内打 “”, 否则打 “×”, 同时根据每题的要求完成证明或说理.
4. | (本题满分 8 分) 设 φ(t) 是向量场充分光滑的 n 维自治系统 x˙=f(x) 一特解. 【】 | 若 ∥x∥→+∞limsup∥x∥∥f(x)∥=2023, 那么上述系统的任何一个解的最大存在区间是 R; | 【】 | 若存在满足 k→∞limtk=+∞ 的数列 {tk}k=0∞, 成立 k→∞limφ(tk)=x∗, 那么 x∗ 为上述系统的奇点或周期轨道上的点; | 【】 | 若 n≥3 且 φ(t) 是非平凡周期解, 则该解一定不是 Lyapunov 意义下渐近吸引的; | 【】 | 若 Jacobian 矩阵 Df(φ(t)) 所有的特征值在 [0,+∞) 上都小于 −2023, 则上述系统的解 φ(t) 在 Lyapunov 意义下是渐近稳定的. |
在以上结论中, 任意选择一个你认为正确的结论, 证明之. [证明] |
5. | (本题满分 9 分) 考虑二阶微分方程 x¨+a(t)x+b(t)=f(t,x˙), 其中 a(t) 与 b(t) 是 R 上的连续函数, f(t,y) 是 R×R 上充分光滑的二元函数. 【】 | 存在上述函数 a,b 及 f, 使得上述方程任意两个不同的解在平面 x–O–x˙ 上对应的轨线于有限时刻相交, 相交也不重合; | 【】 | 存在上述函数 a,b 及 f, 使得 η(t)=et−1,ξ(t)=tcos2t 都为上述方程的解; | 【】 | 对于任意给定的上述函数 a,b 及 f, 以上方程一切解的最大存在区间为 R; | 【】 | 若 f(t,y)≡2t, 存在上述函数 a 和 b, 使得上述方程的任意一解关于初值在 [0,+∞) 上是连续依赖的; | 【】 | 若 a(t)≡1,b(t)≡0 且 f(t,y)≡f(y), 存在这样的 f 使得上述微分方程的零解是渐近稳定的. |
在以上结论中, 任意选择一个你认为错误的结论, 说明错误的理由. [理由] |
三、 解答题 (本大题满分 56 分) 本大题共有 4 题, 解答下列各题须写明推理过程.
6. | (本题满分 8 分) 叙述并证明: 自治系统的零解在 Lyapunov 意义下是渐近稳定的判定定理. (6) [叙述定理] (6) [证明] |
7. | (本题满分 21 分) 本题共有 3 个小题, 第 1 小题满分 10 分, 第 2 小题满分 6 分, 第 3 小题满分 5 分. 设 Φ(x,y) 为 R2 上充分光滑的函数, 考虑含有不全为零的参数 a,b 的平面系统: x˙=∂y∂Φ−a∂x∂Φ,y˙=−∂x∂Φ−b∂y∂Φ. (1) | 试用压缩映象原理证明: 对于给定的 h>0, 存在 T<h, 使得以上系统满足条件 x(h)=0,y(h)=1 的解在 [T,h] 上存在; | (2) | 若 ab=0 且 Φ(x,y)=x6+y4+x2y2, 证明: 上述系统在相平面中没有闭轨线; | (3) | 若 a>0,b>0 且原点 O(0,0) 是 Φ(x,y) 孤立的极小值点, 证明: 原点 O 对应的零解在 Lyapunov 意义下是稳定的. |
(7-1) [证明] (7-2) [证明] (7-3) [证明] |
8. | (本题满分 20 分) 本题共有 3 个小题, 第 1 小题满分 7 分, 第 2 小题满分 5 分, 第 3 小题满分 8 分. 考虑如下二维平面自治系统: {x˙=αx−βxy,y˙=−γy+ζxy,其中, 参数 α,β,γ,ζ 均大于零. (1) | 求上述系统在相平面 x–O–y 中坐标轴上及第一象限内的奇点, 并确定奇点类型; | (2) | 证明: 系统在第一象限内出发的轨线永远停留在第一象限中; | (3) | 若 x>0 表示蝗虫数量, y>0 表示益鸟数量. 如果对于蝗虫使用杀虫剂, 同时这个杀虫剂对于益鸟也有同样的抑制作用. 若视这样的抑制作用为对于参数 α 和 γ 的扰动, 试讨论蝗虫、益鸟在一个周期内平均数量如何分别与这样的扰动连续依赖. |
(8-1) [解] (8-2) [证明] |
9. | (本题满分 7 分) 设常数 τ>0, 且实值函数 x(t) 满足: x∈C([−τ,+∞))∩C1(∪n=0∞Tn),x(0)=0,其中区间 Tn≜(nτ−τ,nτ) (n=0,1,2,⋯). 若 x(t) 在 [0,+∞) 中的可导处成立: dtdx(t)=a⋅x(t)+b⋅h(t)⋅x(t−τ), 其中 h(t)={0,1,t∈T2m+1,t∈T2m+2,m=0,1,2,⋯.试确定参数 a 与 b 的范围使得 t→+∞limx(t)=0. (9) [解] |
13–14 年期末考试试卷 (林伟)
考生注意:
1. | 这份试卷共 8 道试题, 满分 100 分, 共 8 页. |
2. | 考试时间: 2014 年 1 月 16 日 13:00~15:00. |
3. | 请将每张试卷装订线内的学号、姓名及专业填写清楚. |
4. | 学期总评分 (满分 100 分) 由平时成绩 (满分 50 分) 与本试卷成绩的 50% 组合计算. |
一、 填空题 (本大题满分 14 分) 本大题共有 2 题, 将答案直接写在题后的横线上, 每题满分 7 分.
1. | 设参数 a,b,c∈R, 而平面自治系统 x˙=cy+x2,y˙=ax+by+xy3 具有一个连续的首次积分 f(x,y)=C, 且其一次近似系统的系数矩阵行列式不等于零. 那么该系统的奇点 (0,0) 的类型是 . |
2. | 设参数 λ∈R, 3 阶矩阵 A 满足: Ap1+λp1=0,Ap2+λp2=0,Ap3+λp3=p2, 其中 p1,2,3 为线性无关的向量. 则微分方程组 x˙=Ax+[0,cost,e2t]⊤ 的任意一解在 Lyapunov 意义下的稳定性是 . |
二、 辨析题 (本大题满分 24 分) 本大题共有 2 题, 每题都给出了五个结论, 分别判断这些结论是否正确, 正确的在结论前的方括号内打 “✓”, 否则打 “×”, 同时根据要求完成说理或证明, 每题满分 12 分.
3. | 考虑 2 阶微分方程 x¨+p(t)x=q(t), 其中实值函数 p(t),q(t) 在 R 上连续. 【】 | 该方程的任意一解的存在区间均为 R; | 【】 | 该方程至多有 3 个线性无关的解; | 【】 | 若 p(t)≡−1, q(t) 是周期为 π 的函数, 则该方程有且仅有一个周期为 π 的周期解. | 【】 | 若 q(t)≡0, 且该方程有惟一的满足条件 x(0)=x(1)=0 的解, 则该方程对应的齐次方程满足条件 x(0)=x(1)=0 的只有零解; | 【】 | 若 q(t)≡0, 则该方程的零解在 Lyapunov 意义下可以是渐近稳定的; |
在以上结论中, 任意选择一个你认为错误的结论, 说明错误的理由. [理由] |
4. | 考虑右端在区域 D⊂Rn 上具有连续可微的向量场的 n 维自治系统 x˙=f(x). 【】 | 系统 x˙=f(x)/[∥f(x)∥+1] 任意一解的最大存在区间为 R; | 【】 | 若维数 n=1, 则该系统可以有非平凡的周期解; | 【】 | 记上述系统从 (0,x0) 出发的解为 x(t,x0), 则对于充分小的任意 t1,t2 成立 x(t1+t2,x0)=x(t2,x(t1,x0)); | 【】 | 若 x=0 是上述系统的零解, 且对于初始时刻 t0, 该解是 Lyapunov 意义下稳定的, 则对于任意其它的初始时刻 t^=t0, 该解亦是 Lyapunov 意义下稳定的; | 【】 | 若 x0 是上述系统的奇点, 并且一次近似线性系统 x˙=Df(x0)x 的零解在 Lyapunov 意义下是不稳定的, 则奇点 x0 在 Lyapunov 意义下是不稳定的. |
在以上结论中, 任意选择一个你认为正确的结论, 证明之. [证明] |
三、 解答题 (本大题满分 62 分) 本大题共有 4 题, 解答下列各题须写明推理过程.
5. (本题满分 10 分) 从以下两个小题中选择一个完成.
(1) | 叙述并证明: 自治系统的零解在 Lyapunov 意义下渐近稳定的判定定理. |
(2) | 叙述并证明: 平面自治系统同宿轨不存在的 Bendixson 判定定理 (其中, 同宿轨指相平面中从奇点出发又回到该奇点的奇异闭轨线). |
(5) [叙述定理]
(5) [证明]
6. (本题满分 25 分) 本题共有 3 个小题, 第 1 小题满分 12 分, 第 2 小题满分 7 分, 第 3 小题满分 6 分.
考虑以 μ,ξ 为参数的平面系统 {x˙=μx−2y+ξx(x2+y2),y˙=2x+μy+ξy(x2+y2).
(1) | 利用压缩映象原理证明: 存在 h>0, 上述系统满足初值条件 x(0)=x0,y(0)=y0 的解在区间 [0,h] 上存在且惟一; |
(2) | 若 ξ=1, 则上述系统的解是否关于参数 μ 连续依赖? 论证你的结论; |
(3) | 若 ξ=−1, 分别绘制当 μ<0,μ=0,μ>0 时, 上述系统在相平面上轨线的大致图形. |
(6-1) [证明]
(6-2) [结论与论证]
(6-3) [作图]
7. (本题满分 20 分) 本题共有 3 个小题, 第 1 小题满分 7 分, 第 2 小题满分 7 分, 第 3 小题满分 6 分.
设 Φ(x,y) 是定义在 R×R 上无穷次可微的二元函数, 考虑如下平面系统: dtdx=∂y∂Φ(x,y),dtdy=−∂x∂Φ(x,y).(⋯⋯[⋇])
(1) | 若设 Φ(x,y)=ax2+bxy+cy2+y4, 试确定非零参数 a,b,c 满足的条件, 使得系统 [⋇] 的零解在 Lyapunov 意义下是不稳定的; |
(2) | 若 (0,0) 是函数 Φ(x,y) 孤立的极小值点, 证明: 系统 [⋇] 的零解是 Lyapunov 意义下稳定的; |
(3) | 证明: 系统 [⋇] 任意的平衡点对应的解都不是 Lyapunov 意义下渐近稳定的. |
(7-1) [解]
(7-2) [证明]
8. (本题满分 7 分)
若连续函数 p(t)>2, 考虑平面系统: x˙=2y,y˙=−2x−p(t)y3. 试证明: 存在 x∗∈R 使得 x(t)→x∗,y(t)→0 (t→+∞).
(8) [证明]
12–13 年期末考试试卷 (林伟)
考生注意:
1. | 这份试卷共 8 道试题, 满分 100 分, 共 8 页. |
2. | 考试时间: 2013 年 1 月 18 日 8:30~10:30. |
3. | 请将每张试卷装订线内的学号、姓名及专业填写清楚. |
4. | 学期总评分 (满分 100 分) 由平时成绩 (满分 50 分) 与本试卷成绩的 50% 组合计算. |
一、 填空题 (本大题满分 14 分) 本大题共有 2 题, 将答案直接写在题后的横线上, 每题满分 7 分.
1. | 设参数 λ∈R, 3 阶矩阵 A 满足: Ap1=λp1,Ap2=λp2,Ap3=λp3+p2, 其中 p1,2,3 为线性无关的向量. 则微分方程组 x˙=Ax+[sint,0,e−t]⊤ 的任意一解在 Lyapunov 意义下的稳定性是 . |
2. | 设参数 a,b,c∈R, 平面自治系统 x˙=ax+by+x3y,y˙=cx+y2 具有连续的首次积分 f(x,y)=C, 且其一次近似系统的系数矩阵行列式不等于零. 那么该系统的奇点 (0,0) 的类型是 . |
二、 辨析题 (本大题满分 24 分) 本大题共有 2 题, 每题都给出了五个结论, 分别判断这些结论是否正确, 正确的在结论前的方括号内打 “”, 否则打 “×”, 同时根据要求完成说理或证明, 每题满分 12 分.
3. | 考虑 2 阶微分方程 x¨+p(t)x=q(t), 其中实值函数 p(t),q(t) 在 R 上连续. 【】 | 该方程的任意一解的存在区间均为 R; | 【】 | 该方程至多有 3 个线性无关的解; | 【】 | 若 q(t)≡0, 且该方程有惟一的满足条件 x(0)=x(1)=0 的解, 则该方程对应的齐次方程满足条件 x(0)=x(1)=0 的只有零解; | 【】 | 若 q(t)≡0, 则该方程的零解在 Lyapunov 意义可以是渐近稳定的; | 【】 | 若 p(t)≡−1, q(t) 是周期为 π 的函数, 则该方程有且仅有一个周期为 π 的周期解. |
在以上结论中, 任意选择一个你认为错误的结论, 说明错误的理由. [理由] |
4. | 考虑右端在 D⊂Rn 上具有连续可微的向量场的 n 维自治系统 x˙=f(x). 【】 | 系统 x˙=f(x)/[∥f(x)∥+1] 任意一解的最大存在区间为 R; | 【】 | 若维数 n=1, 则该系统可以有非平凡的周期解; | 【】 | 若 t→+∞limx(t;0,x1)=x0 且 x0=x1, 则 x0 是系统的奇点; | 【】 | 若 x0 是上述系统的奇点, 并且线性系统 x˙=Df(x0)x 的零解在 Lyapunov 意义下是不稳定的, 则奇点 x0 在 Lyapunov 意义下是不稳定的; | 【】 | 若上述系统具有稳定的极限圈, 则该极限圈对应的解在 Lyapunov 意义下是渐近稳定的. |
在以上结论中, 任意选择一个你认为正确的结论, 证明之. [证明] |
三、 解答题 (本大题满分 62 分) 本大题共有 4 题, 解答下列各题须写明推理过程.
5. | (本题满分 11 分) 从以下两个小题中选择一个完成. (1) | 叙述并证明: 自治系统的零解在 Lyapunov 意义下渐近稳定的判定定理. | (2) | 叙述并证明: 平面自治系统闭轨线不存在的 Bendixson 判别法. |
|
6. | (本题满分 24 分) 本题共有 3 个小题, 第 1 小题满分 12 分, 第 2 小题满分 6 分, 第 3 小题满分 6 分. 考虑以 μ,ξ 为参数的平面系统 {x˙=μx−y+ξx(x2+y2),y˙=x+μy+ξy(x2+y2). (1) | 利用压缩映象原理证明: 存在 h>0, 上述系统满足初值条件 x(0)=x0,y(0)=y0 的解在区间 [0,h] 上存在且惟一; | (2) | 若 ξ=1, 则上述系统的解是否关于参数 μ 连续依赖? 论证你的结论; | (3) | 若 ξ=−1, 分别绘制当 μ<0,μ=0,μ>0 时, 上述系统在相平面上轨线的大致图形. |
(6-1) [证明] (6-2) [结论与论证] (6-3) [作图] |
7. | (本题满分 19 分) 本题共有 3 个小题, 第 1 小题满分 5 分, 第 2 小题满分 5 分, 第 3 小题满分 9 分. 二元函数 Φ(p,q)=λp2+q4 (λ=0), 考虑以 λ,ξ 为参数的平面系统: dtdx=2y3+λξxy4,dtdy=−λx−3x2y,(⋯⋯[ξ]) (1) | 如果 Φ(x,y)=C (C 为任意正常数) 是方程组 [ξ0] 的一个首次积分, 试确定 ξ0 的值; | (2) | 若 ξ=ξ0, 证明: 在相平面中不存在系统 [ξ] 的闭轨线; | (3) | 若 λ>0,ξ⩽ξ0, 试利用 李雅普诺夫直接方法 确定系统 [ξ] 的零解的稳定性. |
(7-1) [解] (7-2) [证明] (7-3) [解] |
8. | (本题满分 8 分) 若连续函数 p(t)>2, 考虑平面系统: x˙=2y,y˙=−2x−p(t)y3. 试证明: 存在 x∗∈R 使得 x(t)→x∗,y(t)→0 (t→+∞). (8) [解] |