二阶微分方程的教学课件

二阶微分方程从基础概念到数值与解析求解的完整教学Contents课程目录二阶微分方程的系统学习路径，从基础概念到物理应用。01基础概念与物理背景02转化为一阶方程组03欧拉法数值求解04特征值与解析解法05相图分析与解读06物理应用案例07综合练习与课程总结Chapter01基础概念与物理背景理解二阶微分方程的定义、分类与物理意义SECONDORDERODE什么是二阶微分方程二阶微分方程是最高阶导数为二阶的微分方程，广泛描述加速度相关的物理过程。理解其一般形式与线性/非线性、齐次/非齐次的分类，是选择正确求解策略的前提。01一般形式为F(x,y,y',y'')=0，其中y''是未知函数y(x)对自变量x的二阶导数，方程阶数由最高阶导数决定02线性二阶ODE标准形式：y''+p(x)y'+q(x)y=g(x)，当p、q为常数时称为常系数线性方程，是最常见且可解析求解的类型03按g(x)分类：g(x)=0时为齐次方程，解具有叠加性；g(x)≠0时为非齐次方程，通解等于齐次通解加特解04非线性二阶方程如y''+sin(y)=0（单摆方程）通常无法获得闭式解，需依赖数值方法或线性化近似处理教学场景·黑板上的微分方程公式SECOND-ORDERODE·物理建模物理世界的二阶微分方程牛顿第二定律F=ma本质上是二阶微分方程，因为加速度是位移的二阶导数。弹簧振子、单摆、RLC电路等经典物理模型都可以用二阶线性常系数ODE精确描述，这是学习二阶方程的核心物理动机。01弹簧-质量-阻尼系统：mx″+cx′+kx=0，其中m为质量、c为阻尼系数、k为弹簧刚度，完整描述受迫振动的动力学行为02单摆运动方程：θ″+(g/L)sin(θ)=0，小角度近似简化为简谐振动，周期T=2π√(L/g)03RLC串联电路：Lq″+Rq′+(1/C)q=V(t)，与弹簧系统数学结构完全相同，体现不同物理系统的数学统一性04共同特征：系统状态由位置和速度两个变量完全确定，需要y(0)和y′(0)两个初始条件确定唯一解弹簧-质量-阻尼系统物理实验装置COMPARISON·ODE一阶与二阶微分方程的关键差异二阶微分方程的通解包含两个独立的任意常数，需要两个初始条件才能确定唯一解。与一阶方程相比，二阶方程的解空间维度更高，能描述振荡等更复杂的动力学行为，同时也要求更高级的求解工具。2C解的结构差异一阶ODE通解含1个任意常数C，二阶ODE通解含2个独立常数C₁、C₂，解空间从一维扩展到二维y₀+y₀′初始条件需求一阶方程仅需y(x₀)=y₀即可定解，二阶方程需要同时给出y(x₀)=y₀和y′(x₀)=y₀′两个条件OSC动力学行为丰富度一阶自治系统只能单调趋近或远离平衡点，二阶系统可产生振荡、螺旋收敛、极限环等复杂行为λ²求解方法升级一阶的分离变量法和积分因子法不直接适用于二阶，需引入特征方程法、待定系数法、参数变易法等CharacteristicEquationMethod特征方程法：常系数齐次方程的解析解特征方程法是求解常系数齐次二阶线性ODE的核心方法。通过代入试探解将微分方程转化为代数方程，根据判别式分为三种情况。Δ>0两个不等实根01通解y=C₁er₁x+C₂er₂x，解表现为纯指数增长或衰减02物理对应"过阻尼"：物体缓慢回到平衡位置而不振荡，如液压缓冲器过阻尼Δ=0重根01通解y=(C₁+C₂x)erx，第二个线性无关解通过乘以x获得02物理对应"临界阻尼"：系统以最快且不振荡的方式回到平衡态临界阻尼Δ<0共轭复根01通解y=eαx(C₁cosβx+C₂sinβx)，实部控制衰减，虚部控制振荡频率02物理对应"欠阻尼"：弹簧系统在小阻尼下做衰减振荡，能量逐步耗散欠阻尼Chapter02转化为一阶方程组将二阶ODE降阶为一阶系统，架起微积分与线性代数的桥梁Second-OrderODE·Reduction降阶转化：从二阶到一阶系统通过引入新变量u=y和v=y'，任何二阶微分方程都可以等价转化为一阶方程组。这种降阶技巧不仅统一了高阶ODE的处理框架，还为数值方法和矩阵分析提供了标准输入格式。01核心变量替换对二阶方程y''=f(x,y,y')，令u=y、v=y'，则原方程等价于一阶系统u'=v,v'=f(x,u,v)，维度从1个二阶变为2个一阶1→202具体示例y''+2y'+5y=0转化为u'=v,v'=-5u-2v，初始条件y(0)=1,y'(0)=0对应u(0)=1,v(0)=0ODE03推广到n阶方程令x₁=y,x₂=y',…,xₙ=y⁽ⁿ⁻¹⁾，则n阶ODE转化为n维一阶系统x'=F(x)，为统一处理任意阶方程奠定框架n维04物理直觉u和v分别对应系统的"位置"和"速度"，二者共同构成完整状态描述，与经典力学中的相空间概念一致PhaseSpace线性代数方法矩阵形式：x'=Ax常系数线性二阶ODE转化为一阶系统后，可以写成紧凑的矩阵形式x'=Ax。系数矩阵A的特征值和特征向量直接决定了系统解的结构和动力学行为，从而将微分方程求解问题完全转化为线性代数问题。线性代数课堂教学场景01矩阵形式构建：对系统u'=v,v'=-5u-2v，状态向量x=[u,v]ᵀ，系数矩阵A=[[0,1],[-5,-2]]，方程简化为x'=Ax02矩阵指数的解：若x'=Ax，则通解为x(t)=e^(At)x₀，矩阵指数e^(At)的计算依赖于A的特征分解A=PDP⁻¹03特征值的物理意义：A的特征值λ的实部决定解的增减性（稳定/不稳定），虚部决定振荡频率，与特征方程法中根的分析完全一致04特征向量的几何意义：特征向量p指明了相空间中的"特殊方向"，系统沿这些方向的运动是纯指数形式，所有其他运动都是这些基本模式的线性组合WorkedExample完整示例：从方程到矩阵求解以y''−3y'+2y=0为例，完整演示二阶ODE到通解的五步标准化流程，适用于所有常系数线性二阶方程。01降阶与矩阵化：令u=y,v=y'得系统u'=v,v'=−2u+3v，矩阵A=[[0,1],[−2,3]]，初始向量x₀=[1,0]ᵀu'=v,v'=−2u+3v02求特征值：det(A−λI)=λ²−3λ+2=(λ−1)(λ−2)=0，得λ₁=1（不稳定增长）和λ₂=2（更快增长模式）λ₁=1,λ₂=203求特征向量：λ₁=1→p₁=[1,1]ᵀ，λ₂=2→p₂=[1,2]ᵀ，两个特征向量线性无关，构成相空间的基p₁=[1,1]ᵀp₂=[1,2]ᵀ04组合通解并定初值：x(t)=C₁eᵗ[1,1]ᵀ+C₂e²ᵗ[1,2]ᵀ，代入x(0)=[1,0]ᵀ解得C₁=2,C₂=−1y(t)=2eᵗ−e²ᵗNUMERICALMETHODS·EULER欧拉法计算示例：y''+y=0以简谐振动方程y''+y=0为例（精确解为cost），用h=0.2的欧拉法逐步计算y(1)的近似值。通过对比数值解与精确解，可以直观看到欧拉法的精度特征和误差积累规律。欧拉法逐步迭代计算表（h=0.2）tₙuₙ≈y(tₙ)vₙ≈y'(tₙ)精确解cos(tₙ)绝对误差0.01.00000.00001.00000.00000.21.0000-0.20000.98010.01990.40.9600-0.40000.92110.03890.60.8800-0.59200.82530.05470.80.7616-0.76800.69670.06491.00.6080-0.92030.54030.0677随迭代步数增加，欧拉法的累积误差逐步增大，在t=1.0时绝对误差达到约0.068问题设定y''+y=0，y(0)=1，y'(0)=0，转化为一阶方程组u'=v，v'=-uh=0.2迭代过程逐步推进：(1,0)→(1,-0.2)→(0.96,-0.4)→…5步计算结果y(1)≈0.54vs精确解cos(1)≈0.5403，误差约0.7%≈0.7%误差分析折线逼近走过圆弧外侧，振幅缓慢增长，能量不守恒数值不稳定NumericalMethods欧拉法的局限性与改进方法欧拉法作为一阶方法，精度有限且在振荡系统中存在数值不稳定问题。改进欧拉法和四阶龙格-库塔法通过多点采样显著提升了精度和稳定性，是工程实践中的主流选择。Limitations核心局限性01一阶精度限制：全局误差O(h)，步长减半仅使误差减半，要达到10⁻⁶精度可能需要h<10⁻⁶，计算量不可接受02振荡系统能量漂移：对保守系统（如简谐振动），数值解振幅逐步增长，违反能量守恒，长时间模拟严重失真03刚性方程失效：特征值量级差异很大时，欧拉法要求步长极小才能稳定，计算效率极低Improvements改进方法01改进欧拉法（Heun）：先用欧拉法预测，再用终点斜率修正，精度提升到O(h²)，计算量仅翻倍但精度提高一个数量级02四阶龙格-库塔法（RK4）：每步取4个采样点加权平均，精度达O(h⁴)，是工程计算和科学模拟中最常用的单步方法03自适应步长策略：根据局部误差估计动态调整h，在解变化剧烈处加密步长，平缓处放松步长，兼顾精度与效率CHAPTER04特征值与解析解法利用线性代数的特征值分解工具求解二阶微分方程的精确解MATRIXEXPONENTIAL矩阵指数e^(At)：解析解的核心工具矩阵方程x'=Ax的精确解为x(t)=e^(At)x₀。矩阵指数通过特征分解A=PDP⁻¹简化为Pe^(Dt)P⁻¹的计算，其中e^(Dt)只需对角元素取指数。这一方法将微分方程求解完全转化为线性代数运算。级数定义e^(At)=I+At+(At)²/2!+(At)³/3!+…，与标量指数级数形式完全一致，保证了微分方程x'=Ax解的解析性无穷级数对角化计算若A可对角化为PDP⁻¹（P为特征向量矩阵，D=diag(λ₁,λ₂)），则e^(At)=P·diag(e^(λ₁t),e^(λ₂t))·P⁻¹，计算量大幅降低特征分解Jordan标准形当A有重特征值且几何重数小于代数重数时，需用A=PJP⁻¹，此时e^(Jt)含te^(λt)项，对应微分方程中的xe^(rx)解广义特征物理意义e^(At)是系统的状态转移矩阵（基本解矩阵），描述初始状态x₀在时间t后的演化结果，是系统动力学的完整描述状态转移SolutionMethod特征值方法求解五步法特征值方法将二阶ODE求解标准化为五个清晰步骤，特征值的类型直接决定了通解的函数形式和系统的动力学行为。01矩阵化将二阶ODE转化为一阶系统x'＝Ax，确定系数矩阵A和初始状态向量x₀x'＝Ax02特征分解求解det(A−λI)＝0得特征值，再解(A−λᵢI)pᵢ＝0得特征向量λ·Eigenvector03构造通解实根叠加指数解；复根提取实虚部构造三角函数解；重根引入teλt项实根·复根·重根04定初始条件将x(0)＝x₀代入通解，解线性方程组确定C₁、C₂得唯一特解C₁·C₂05验证与解释将解代回原方程验证正确性，分析特征值实部符号判断系统稳定性稳定性·验证CompleteWorkedExample完整计算示例：阻尼振荡的解析解以y''+4y'+13y=0为例，特征值为共轭复根-2±3i，对应物理上的阻尼振荡。通过特征向量分解和实虚部提取，得到精确解y(t)=e−2t(cos3t+⅔sin3t)，清晰展示了衰减因子和振荡频率的来源。Step01矩阵化与特征值A=[[0,1],[-13,-4]]，det(A-λI)=λ²+4λ+13=0，解得λ=-2±3i（共轭复根，实部-2控制衰减，虚部3控制振荡频率）λ=-2±3iStep02特征向量计算λ₁=-2+3i时解(A-λ₁I)p=0，得p₁=[1,-2+3i]ᵀ，实部Re(p₁)=[1,-2]ᵀ，虚部Im(p₁)=[0,3]ᵀRe(p₁)=[1,-2]ᵀIm(p₁)=[0,3]ᵀStep03构造实值通解x(t)=e-2t{C₁(cos3t·[1,-2]ᵀ−sin3t·[0,3]ᵀ)+C₂(sin3t·[1,-2]ᵀ+cos3t·[0,3]ᵀ)}C₁·cos+C₂·sinStep04定初始条件x(0)=[C₁,-2C₁]ᵀ=[1,0]ᵀ→C₁=1，再由-2C₁+3C₂=0→C₂=2/3，最终得精确解y(t)=e−2t(cos3t+⅔sin3t)MethodologyComparison两种方法的等价性：特征方程vs特征值分解传统特征方程法与矩阵特征值方法在数学上完全等价——两者的特征多项式相同，根的分类和通解结构一致。矩阵方法的额外价值在于提供相空间的几何视角，并能自然推广到高维系统和数值计算。特征方程法视角01方法对ay″+by′+cy=0代入y=erx，得特征方程ar²+br+c=0，根据根的三种情况直接写出y的通解表达式02优势计算过程简短直接，不需要线性代数知识，适合手算和考试场景03局限难以推广到高阶方程组和非自治系统，缺乏几何直觉矩阵特征值视角01等价性A=[[0,1],[−c/a,−b/a]]的特征方程det(A−λI)=λ²+(b/a)λ+c/a=0与原特征方程仅差常数因子a，根完全相同02优势自然延伸到n维系统x′=Ax，特征向量给出相空间中的特殊方向，便于可视化相图和稳定性分析03关联P矩阵的列就是特征向量，P⁻¹x₀给出展开系数C₁,C₂，与特征方程法中代入初始条件求常数本质相同SolvingStrategy非齐次方程的求解策略非齐次方程的通解=齐次通解+特解。待定系数法适用于g(x)为多项式、指数或三角函数的情况，通过猜测特解形式并代入确定系数；参数变易法更通用，将齐次解中的常数变为函数后求解。两种方法互补构成完整的非齐次求解工具箱。待定系数法Presumptiong(x)为多项式、eαx、sin/cos(βx)及其乘积时，根据g(x)形式预设yp的结构并代入方程确定未知系数Resonance预设形式与齐次解重合时，需乘以x提升阶数（yp=Axe2x），物理上对应受迫振动的共振现象参数变易法GeneralMethod令C₁→u₁(x),C₂→u₂(x)，由辅助条件u₁'y₁+u₂'y₂=0和u₁'y₁'+u₂'y₂'=g(x)解出u₁',u₂'，再积分得特解MatrixForm对x'=Ax+f(t)，特解为xp(t)=Φ(t)∫Φ−1(s)f(s)ds，其中Φ(t)是基本解矩阵，与线性代数变参数公式一致CHAPTER05相图分析与解读用相平面上的轨迹图形直观理解二阶系统的全局动力学行为PhasePortrait相图的基本概念相图是在位置-速度平面（u-v平面）上描绘系统所有可能运动轨迹的图形工具。每条相轨迹对应一组初始条件，轨迹的几何形状直接反映系统的动力学特征。01相平面定义以u=y（位置）为横轴、v=y'（速度）为纵轴的二维平面，每个点(u,v)代表系统的一个完整状态，系统演化即该点在相平面上的运动。02向量场在相平面每个点处，系统有确定的速度向量(u',v')=(v,f(u,v))，所有向量构成向量场，相轨迹就是沿向量场方向的积分曲线。03相轨迹的物理读法闭合轨道→周期运动（能量守恒）；向内螺旋→衰减振荡（能量耗散）；向外螺旋→发散振荡（能量输入）；趋向固定点→系统趋于平衡。04关键优势无需解出y(t)的显式表达式，仅通过相图的几何结构就能判断系统的稳定性、周期性和长期行为，特别适合非线性系统的定性分析。PhasePortraitClassification特征值类型与相图形态的对应关系矩阵特征值的类型完整决定了相图的几何形态：实根对应结点（过阻尼）、鞍点（不稳定平衡）；复根对应焦点（阻尼振荡）、中心（保守振荡）。特征值实部的正负决定稳定性（收敛/发散），虚部决定是否有振荡。实根情形01两负实根（稳定结点）：所有轨迹单调趋向原点，无振荡，对应过阻尼系统如浸入蜂蜜中的弹簧振子λ₁,λ₂<002两正实根（不稳定结点）：所有轨迹从原点单调发散，系统在任何微扰下都会远离平衡态λ₁,λ₂>003异号实根（鞍点）：存在两条稳定流形和两条不稳定流形，大多数轨迹先靠近后远离原点，如倒立摆的不稳定平衡λ₁·λ₂<0复根情形01负实部复根（稳定焦点）：轨迹螺旋收敛至原点，振荡幅度指数衰减，对应欠阻尼系统如汽车减震器正常工作Re(λ)<0,Im(λ)≠002正实部复根（不稳定焦点）：轨迹螺旋发散，振荡幅度指数增长，如负阻尼电路中的自激振荡初始阶段Re(λ)>0,Im(λ)≠003纯虚根（中心）：轨迹为闭合椭圆，系统做等幅周期振荡，能量守恒，对应理想无阻尼弹簧或LC电路Re(λ)=0,Im(λ)≠0PhasePortrait手绘相图的实用步骤手绘相图遵循"定点→定性→定向→描迹"四步流程：先找平衡点并确定类型，再画特征向量方向，然后补充向量场箭头，最后描绘光滑轨迹。01找平衡点—解u'=0,v'=0的方程组；线性系统中det(A)≠0则原点为唯一平衡点02确定类型—计算特征值，根据实根/复根/重根及正负号确定结点、焦点、中心或鞍点03画特征方向—特征向量给出直线轨迹方向；"慢"特征值方向是轨迹趋向的主导方向04描绘轨迹—沿向量场方向画光滑曲线；焦点螺旋收敛，鞍点沿双曲线型流形运动黑板手绘相图与向量场教学示意PhasePortraitAnalysis相图绘制示例与参数变化通过三个典型系统展示相图从"中心→焦点→结点"的连续演变：无阻尼产生闭合椭圆，添加阻尼后螺旋收敛，继续增大则直接趋向原点。稳定结点A=[[-1,0],[0,-2]]特征值λ₁=-1,λ₂=-2，均为负实根，所有轨迹单调趋向原点，无振荡行为慢方向为x轴，轨迹最终沿x轴趋于原点；快方向y轴快速衰减λ₁=-1·λ₂=-2中心A=[[0,1],[-1,0]]特征值λ=±i（纯虚根），相轨迹为闭合椭圆，系统做等幅周期振荡，能量守恒旋转方向为逆时针：第一象限向右运动，第四象限向左运动λ=±i稳定焦点A=[[0,1],[-1,-0.5]]λ=-0.25±0.97i，负实部使振幅指数衰减，虚部产生旋转，轨迹为向内螺旋线阻尼参数增大时，相图经历中心→焦点→结点的分岔过程中心→焦点→结点CHAPTER06物理应用案例将二阶微分方程工具应用于振动、电路和共振等经典物理问题SECOND-ORDERODE·DAMPEDOSCILLATION弹簧-质量-阻尼系统的三种运动模式弹簧系统的运动方程x″+2ζω₀x′+ω₀²x=0中，阻尼比ζ完全决定了系统的运动模式。ζ>1为过阻尼（缓慢回位）、ζ=1为临界阻尼（最快回位）、0<ζ<1为欠阻尼（衰减振荡），这三种模式对应特征方程的两种实根、重根和共轭复根。01物理参数标准化：令ω₀=√(k/m)（固有频率）、ζ=c/(2√(mk))（阻尼比），方程化为x″+2ζω₀x′+ω₀²x=0，物理意义更清晰02过阻尼ζ>1：特征根为两个负实数，x(t)=C₁e^(r₁t)+C₂e^(r₂t)，系统缓慢回到平衡位置不振荡，如液压缓冲器、门吸03临界阻尼ζ=1：特征根为重根−ω₀，x(t)=(C₁+C₂t)e^(−ω₀t)，系统以最快且不超调的方式回到平衡态，精密仪器常设计为此模式04欠阻尼0<ζ<1：特征根为−ζω₀±iω₀√(1−ζ²)，x(t)=Ae^(−ζω₀t)cos(ωdt−φ)，衰减率由ζω₀控制，ωd=ω₀√(1−ζ²)为阻尼频率弹簧-质量-阻尼器实物结构SECOND-ORDERODE·CIRCUITDYNAMICSRLC串联电路中的二阶动力学RLC串联电路的电荷方程Lq''+Rq'+q/C=V(t)与弹簧系统在数学上完全同构：电感L↔质量m、电阻R↔阻尼c、1/C↔刚度k。这种跨领域的数学统一性使得二阶ODE的分析工具可以同时服务于力学和电学问题。基尔霍夫电压定律Lq''+Rq'+(1/C)q=V(t)，令ω₀=1/√(LC)、ζ=R/(2√(L/C))，化为与弹簧系统完全相同的标准形式ω₀=1/√(LC)物理量对应电感L（磁能储存）↔质量m（动能储存）；电容C（电能储存）↔弹簧1/k（势能储存）；电阻R（焦耳热耗散）↔阻尼c（摩擦耗散）L↔m·C↔1/k·R↔c收音机调谐通过调节C改变固有频率ω₀=1/√(LC)，使电路与目标电台频率共振，从混合信号中选出特定频率频率共振选台品质因数Q=(1/R)√(L/C)：Q值越高，振荡衰减越慢，频率选择性越好；高Q电路在通信和滤波器设计中至关重要Q=(1/R)√(L/C)ForcedVibration&Resonance受迫振动与共振现象受迫振动系统mx''+cx'+kx=F₀cos(ωt)的稳态响应振幅在驱动频率接近固有频率时急剧增大，产生共振。共振峰值由阻尼比ζ控制：ζ越小，共振峰越尖锐，振幅放大倍数Q≈1/(2ζ)越大。TacomaNarrows大桥因风致共振坍塌，1940SteadyState特解xp(t)=A·cos(ωt-δ)，振幅A(ω)=F₀/√((k-mω²)²+(cω)²)，相位差tan(δ)=cω/(k-mω²)，均取决于驱动频率ωResonancedA/dω=0时振幅最大，共振频率ωr=ω₀√(1-2ζ²)，峰值Amax=F₀/(2mζω₀²√(1-ζ²))QFactorQ=1/(2ζ)衡量共振峰尖锐度，共振时振幅为静态位移的Q倍；半功率带宽Δω≈2ζω₀，Q越高带宽越窄EngineeringTacomaNarrows大桥为经典教训；MRI、微波加热、无线电接收等技术则利用共振实现信号放大与能量传递Cross-Validation数值解与解析解的交叉验证在实际问题中，解析方法和数值方法应当互补使用：解析解提供精确基准和物理洞察，数值解处理复杂模型并验证理论预测。验证数值精度对已知解析解的线性系统，用欧拉法和RK4分别计算，对比不同步长下的误差收敛阶，验证算法实现的正确性。EulervsRK4长时间行为对比欧拉法在保守系统中破坏能量守恒导致振幅漂移，而辛方法（如Verlet积分）能精确保持系统能量。Verlet积分非线性系统大角度单摆θ″+sin(θ)=0无闭式解，必须依赖数值方法；对比小角度近似解和数值解，可量化线性化误差。θ″+sin(θ)=0变系数方程Mathieu方程y″+(a+b·cost)y=0出现在参数共振问题中，解析解仅在某些参数区间存在，数值方法可全面扫描参数空间。Mathieu方程CHAPTER07综合练习与课程总结通过练习巩固核心方法，回顾课程知识体系Practice·课堂练习课堂练习题三道递进难度的练习题覆盖课程核心知识点：基础题检验特征方程法，进阶题考察矩阵方法与相图分析，挑战题测试数值计算与误差分析。基础题：特征方程法1求y″−5y′+6y=0的通解，并在y(0)=1,y′(0)=0下求特解（特征根r=2,3）2验证解满足原方程和初始条件，判断过阻尼、临界阻尼或欠阻尼并解释物理含义提示：先因式分解特征方程，再代入初始条件确定常数r²−5r+6=0进阶题：矩阵方法与相图1将y″+6y′+25y=0化为x′=Ax，求特征值和特征向量，构造实值通解2在u-v相平面上画相图草图，标注旋转方向、衰减趋势和特征向量方向提示：复特征值对应螺旋相图，实部决定收敛或发散x′=Ax·相图挑战题：欧拉法数值计算1用h=0.1欧拉法计算y″+2y′+y=0在t=0.5的近似值，完成5步迭代表2求精确解（重根r=−1），对比数值解误差，分析步长h对精度的影响提示：注意欧拉法的一阶精度特性，误差随步长线性减小h=0.1·5步EXERCISEREVIEW练习题解答要点特征方程实根、矩阵复根与欧拉数值法构成三类解题模板，掌握后可举一反三。基础题解答特征根r₁=2,r₂=3（两不等正实根），通解y=C₁e^(2t)+C₂e^(3t)；代入初始条件得C₁=3,C₂=−2，特解y=3e^(2t)−2e^(3t)。两个正特征根意味着系统不稳定（如倒立摆），物理上对应"负阻尼"或能量输入系统，解指数增长。r₁=2,r₂=3进阶题解答A=[[0,1],[−25,−6]]，特征值