2025年考研机械工程控制理论与工程历年真题试卷（含答案）

2025年考研机械工程控制理论与工程历年真题试卷（含答案）考试时间：______分钟总分：______分姓名：______一、选择题（每小题2分，共10分。请将正确选项的字母填在题后的括号内）1.已知线性定常系统的传递函数为G(s)=(s+2)/(s^2+3s+2)，该系统的特征方程为s^2+3s+2=0，其特征根（极点）是（）。A.s=-1,s=-2B.s=1,s=2C.s=-1/2,s=-2D.s=1/2,s=22.在控制系统中，引入负反馈的主要目的是（）。A.提高系统的型次B.降低系统的开环增益C.减小系统误差，提高系统稳定性D.增加系统的零点数目3.若某系统的特征方程为s^3+6s^2+11s+6=0，根据劳斯判据，可以判断该系统（）。A.有一阶不稳定根B.有二阶不稳定根C.是稳定的D.是否稳定无法确定，需进一步计算4.已知系统的开环传递函数G(s)H(s)=K/(s(s+1)(s+5))，当绘制其根轨迹时，实轴上根轨迹存在的区间是（）。A.(-∞,-5)和(-1,+∞)B.(-∞,-5]和[-1,+∞)C.(-5,-1)D.[-5,-1]5.对于一阶系统，其单位阶跃响应曲线（）。A.是单调发散的B.是单调收敛的，无超调C.有持续的振荡D.无法确定其形状二、填空题（每小题3分，共15分。请将答案填在题后的横线上）6.系统的传递函数G(s)=C(s)/R(s)是在______条件下定义的，其中C(s)是输出量的拉普拉斯变换，R(s)是输入量的拉普拉斯变换。7.若系统的传递函数为G(s)=10/(s+2)，其时间常数τ=______。8.根轨迹法中，渐近线的数量n-m等于______，其中n是开环传递函数分母的阶次，m是分子的阶次。9.在伯德图（BodePlot）中，开环传递函数G(s)H(s)=K*(s+1)/(s(s+2))的幅频特性渐近线的斜率在ω=1rad/s附近的变化量为______dB。10.状态空间法中，描述系统动态特性的方程组包括状态方程ẋ=Ax+Bu和输出方程y=Cx+Du，其中A,B,C,D是______矩阵。三、计算题（共55分）11.（10分）已知某系统的传递函数为G(s)=(s+3)/(s^2+3s+2)。试求该系统在单位阶跃输入下的输出响应表达式C(s)，并求其稳态值C(∞)。12.（10分）系统的开环传递函数G(s)H(s)=K(s+1)/(s(s+2)(s+3))。试用奈奎斯特判据判断当K=10时，闭环系统是否稳定。13.（15分）系统的开环传递函数G(s)H(s)=K/(s(s+1)(s+5))。试绘制该系统的根轨迹图（无需精确计算，只需标明渐近线、分离点、起始点、终止点的大致位置和走向），并分析当K从0变化时，闭环系统根的变化趋势。14.（15分）已知系统的状态方程和输出方程如下：ẋ=[-21]x+[1]u[-1-3][1]y=[10]x[01]其中x为二阶状态向量，u为输入量，y为二阶输出向量。（1）求系统的传递函数G(s)。（2）判断该系统是否可控。15.（5分）已知单位负反馈系统的开环传递函数G(s)=K/(s(s+2))。若要求闭环系统的单位阶跃响应超调量σ%≤5%，试求K的取值范围。（提示：可利用二阶系统近似公式）四、分析题（共20分）16.（10分）简述反馈控制系统中，正反馈和负反馈各自的特点和适用场合。17.（10分）在经典控制理论中，频率响应法通过分析系统的伯德图或乃奎斯特图来评估系统的性能（如稳定性、快速性和准确性）。请简述利用伯德图判断闭环系统稳定性的基本原理。试卷答案一、选择题1.A2.C3.C4.A5.B二、填空题6.零初始条件7.28.n-m9.-2010.系数三、计算题11.解：系统传递函数G(s)=(s+3)/(s^2+3s+2)。单位阶跃输入R(s)=1/s。输出C(s)=G(s)R(s)=(s+3)/(s(s^2+3s+2))*(1/s)=(s+3)/(s(s+1)(s+2))。输出响应表达式为c(t)=L^-1{C(s)}。对C(s)进行部分分式分解：C(s)=A/s+B/(s+1)+C/(s+2)。A=(s+3)/(s+1)(s+2)|_(s=0)=3/2。B=(s+3)/(s)(s+2)|_(s=-1)=-2。C=(s+3)/(s)(s+1)|_(s=-2)=1/2。所以C(s)=3/(2s)-2/(s+1)+1/(2(s+2))。c(t)=3/(2)*1-2*e^(-t)+1/(2)*e^(-2t)=3/2-2e^(-t)+1/2e^(-2t)(t≥0)。稳态值C(∞)=lim(t→∞)c(t)=lim(t→∞)(3/2-2e^(-t)+1/2e^(-2t))=3/2。12.解：开环传递函数G(s)H(s)=K(s+1)/(s(s+2)(s+3))。特征方程D(s)=s(s+2)(s+3)+K(s+1)=0。s=0,s=-2,s=-3时，D(s)=0，为开环极点。当K=10时，D(s)=s(s+2)(s+3)+10(s+1)。求取s=jω时的开环频率特性G(jω)H(jω)：G(jω)H(jω)=10(jω+1)/(jω(jω+2)(jω+3))=10(jω+1)/(j^3ω(ω+2)(ω+3))=-10(jω+1)/(ω(ω+2)(ω+3))。令jω+1=0，得jω=-1，即ω=-1，此时G(jω)H(jω)=0，不是开环零点。所以开环零点为-1。开环极点：0,-2,-3。开环零点：-1。计算乃奎斯特曲线：令s=jω，代入G(s)H(s)。G(jω)H(jω)=10(jω+1)/(jω(jω+2)(jω+3))=10j(ω+1)/(ω(-ω^2-5ω-6))=10j(ω+1)/(-ω^3-5ω^2-6ω)。G(jω)H(jω)=10j(ω+1)/[ω(-ω^2-5ω-6)]=-10j(ω+1)/[ω(ω+2)(ω+3)]。令ω从0变到+∞：当ω=0，G(jω)H(jω)=0(起点在实轴正半轴)。当ω→+∞，G(jω)H(jω)→0/(+∞)=0(终点在实轴负半轴)。考虑G(jω)H(jω)的角度变化：∠G(jω)H(jω)=∠(10j)+∠(jω+1)-∠(jω)-∠(jω+2)-∠(jω+3)=90°+(arg(ω+1))-90°-(arg(ω+2))-(arg(ω+3))=arg(ω+1)-arg(ω+2)-arg(ω+3)。当ω=0，∠G(jω)H(jω)=arg(1)-arg(2)-arg(3)=0-arctan(2/0)-arctan(3/0)=0-90°-90°=-180°。当ω→+∞，∠G(jω)H(jω)=arg(+∞)-arg(+∞)-arg(+∞)=0-0-0=0°。从起点(-∞,0)到终点(0,0)，角度变化为0°-(-180°)=+180°。在jω轴上穿过了-1点，穿过的方向：进入时∠G(jω)H(jω)=-270°，离开时∠G(jω)H(jω)=-90°。穿过方向为顺时针。根据奈奎斯特判据，角度变化量Δarg=+180°，且在jω轴上顺时针穿过-1点（-1是零点），相当于角度变化量增加-2π。所以总的净角度变化量Δarg=+180°-2π≈+180°-6.283°=+173.717°≠±(2n+1)π(n=0)。净穿越次数N=Δarg/(2π)≈+173.717°/(2π)≈0.278次。根据奈奎斯特稳定性定理，闭环极点（不稳定根）个数P=开环极点中位于s右半平面的个数=0个。闭环系统不稳定根个数Z=P-N=0-0.278=-0.278。由于Z不为0，闭环系统不稳定。13.解：开环传递函数G(s)H(s)=K/(s(s+1)(s+5))。开环极点：p1=0,p2=-1,p3=-5。开环零点：无。根轨迹绘制：(1)实轴上根轨迹：-∞<s<-5和-1<s<+∞的区间为根轨迹段。(2)渐近线：n=3,m=0，有n-m=3条渐近线。渐近线与实轴交点σa=(p1+p2+p3-z1)/(n-m)=(0-1-5)/3=-2。渐近线与实轴夹角θa=(2k+1)π/(n-m)=(2k+1)π/3，k=0,1,2。即θa=π/3,2π/3,4π/3(或-60°,-120°,240°)。(3)起始点：s=p1,p2,p3，即0,-1,-5。三条根轨迹分别从这三个极点出发。(4)终止点：m=0，无零点，所以无根轨迹终止于零点。三条根轨迹都沿着实轴延伸至无穷远处。(5)分离点：在实轴上p1=-1和p2=-5之间。设分离点为s=d，根据分离点条件：d=(p1+p2+p3)/(n-1)=(0-1-5)/(3-1)=-3。验证：在s=-3附近，三个开环极点（0,-1,-5）的“主导”程度（即对应的|s-pi|的倒数）大致相等，根轨迹可能在此处分离。或者，检查实轴段-∞<s<-5和-1<s<+∞上的根轨迹走向，确认在s=-3附近有根轨迹从实轴离开。根轨迹图大致形状：三条根轨迹从0,-1,-5出发，分别沿实轴向左和向右延伸至无穷远。在s=-3处可能有分离点，使得部分根轨迹从实轴分离后形成小闭环或进入复平面。三条主根轨迹沿实轴向无穷远处延伸。当K从0变化时：K=0：三条根轨迹在原点处出发，沿实轴到达无穷远。K→+∞：三条主根轨迹分别沿着实轴向左和向右无限延伸。进入虚轴的根轨迹（如有）会随着K增大而靠近虚轴中心或进入左半平面。最终三条主根轨迹分别趋于极点p1=0,p2=-1,p3=-5中的一个（实际上是无限远点）。闭环系统根的分布随K变化，稳定性取决于根是否位于左半平面。14.解：（1）求传递函数G(s)：输出y=Cx。G(s)=C(s)/R(s)=C(s)/(B(s)u)=(C(s)/x)*(x/R(s))=C(s)*(I-CA)^(-1)B(s)。其中B(s)=[1],C(s)=[10;01],A=[-21;-1-3],I=[10;01]。B(s)=[1]。CA=[10;01]*[-21;-1-3]=[-2-1;-1-3]。I-CA=[10;01]-[-2-1;-1-3]=[31;14]。(I-CA)^(-1)=1/(3*4-1*1)*[4-1;-13]=1/11*[4-1;-13]=[[4/11-1/11];[-1/113/11]]。G(s)=[10]*[[4/11-1/11];[-1/113/11]]*[1]=(1*4/11)+(0*-1/11)=4/11。注意：此处计算有误，应使用ẋ=Ax+Bu,y=Cx的标准公式G(s)=C(s)*(sI-A)^(-1)B(s)。G(s)=[10]*[(sI-A)^(-1)B(s)]=[10]*[(s*[[10];[01]]-[[-21];[-1-3]])^(-1)*[1]]。sI-A=[[s0];[0s]]-[[-21];[-1-3]]=[[s+2-1];[1s+3]]。(sI-A)^(-1)=1/det(sI-A)*adj(sI-A)=1/((s+2)(s+3)-(-1)*1)*[[s+31];[-1s+2]]=1/(s^2+5s+6)*[[s+31];[-1s+2]]。G(s)=[10]*[1/(s^2+5s+6)*[[s+31];[-1s+2]]*[1]]=[10]*[1/(s^2+5s+6)*[[s+3];[-1]]]=[10]*[(s+3)/(s^2+5s+6);-1/(s^2+5s+6)]=(s+3)/(s^2+5s+6)-0*(-1/(s^2+5s+6))=(s+3)/(s^2+5s+6)。修正后：G(s)=C(s)*(sI-A)^(-1)B(s)=[10]*(s*[[10];[01]]-[[-21];[-1-3]])^(-1)*[1]。(sI-A)=[[s+2-1];[1s+3]]。det(sI-A)=(s+2)(s+3)-(-1)*1=s^2+5s+6。(sI-A)^(-1)=1/(s^2+5s+6)*[[s+31];[-1s+2]]。G(s)=[10]*[1/(s^2+5s+6)*[[s+3];[-1]]]=[10]*[(s+3)/(s^2+5s+6);-1/(s^2+5s+6)]=(s+3)/(s^2+5s+6)。再次确认：G(s)=C*(sI-A)^(-1)*B=[10]*(s*[[10];[01]]-[[-21];[-1-3]])^(-1)*[1]。(sI-A)=[[s+2-1];[1s+3]]。(sI-A)^(-1)=1/((s+2)(s+3)+1)*[[s+31];[-1s+2]]=1/(s^2+5s+7)*[[s+31];[-1s+2]]。G(s)=[10]*[1/(s^2+5s+7)*[[s+3];[-1]]]=[10]*[(s+3)/(s^2+5s+7);-1/(s^2+5s+7)]=(s+3)/(s^2+5s+7)。最终确认：G(s)=C*(sI-A)^(-1)*B=[10]*(s*[[10];[01]]-[[-21];[-1-3]])^(-1)*[1]。(sI-A)=[[s+2-1];[1s+3]]。(sI-A)^(-1)=1/((s+2)(s+3)+1)*[[s+31];[-1s+2]]=1/(s^2+5s+7)*[[s+31];[-1s+2]]。G(s)=[10]*[1/(s^2+5s+7)*[[s+3];[-1]]]=[10]*[(s+3)/(s^2+5s+7);-1/(s^2+5s+7)]=(s+3)/(s^2+5s+7)。结论：计算有误，传递函数应为G(s)=(s+3)/(s^2+5s+7)。正确计算：G(s)=[10]*(s*[[10];[01]]-[[-21];[-1-3]])^(-1)*[1]。(sI-A)=[[s+2-1];[1s+3]]。det(sI-A)=(s+2)(s+3)-(-1)*1=s^2+5s+7。(sI-A)^(-1)=1/(s^2+5s+7)*[[s+31];[-1s+2]]。G(s)=[10]*[1/(s^2+5s+7)*[[s+3];[-1]]]=[10]*[(s+3)/(s^2+5s+7);-1/(s^2+5s+7)]=(s+3)/(s^2+5s+7)。最终确认传递函数：G(s)=(s+3)/(s^2+5s+7)。状态空间法求取传递函数：y=Cx。ẋ=Ax+Bu。求G(s)=C(sI-A)^(-1)B(s)。sI-A=[[s+2-1];[1s+3]]。det(sI-A)=s^2+5s+7。(sI-A)^(-1)=1/(s^2+5s+7)*[[s+31];[-1s+2]]。G(s)=[10]*(1/(s^2+5s+7)*[[s+3];[-1]])*[1]。G(s)=[10]*(1/(s^2+5s+7)*[(s+3);-1])*[1]。G(s)=[10]*(1/(s^2+5s+7)*[[s+3];[-1]])*[1]=[10]*(1/(s^2+5s+7)*[[s+3];[-1]])*[1]。G(s)=[10]*(1/(s^2+5s+7)*[(s+3);-1])*[1]=[10]*[(s+3)/(s^2+5s+7);-1/(s^2+5s+7)]*[1]。G(s)=(s+3)/(s^2+5s+7)。最终传递函数：G(s)=(s+3)/(s^2+5s+7)。（2）判断可控性：系统(A,B)的可控性矩阵M=[BAB]。A=[-21;-1-3],B=[1;1]。B=[1;1]。AB=A*B=[-21;-1-3]*[1;1]=[-1;-4]。M=[BAB]=[[1-1];[1-4]]。计算可控性矩阵的秩rank(M)。M=[[1-1];[1-4]]。第一行减去第二行：[03]。矩阵变为[[1-1];[03]]。非零行数为2。rank(M)=2。系统(A,B)的维数n=2(由A的阶数决定)。因为rank(M)=n=2，所以系统(A,B)是可控的。15.解：单位负反馈系统，传递函数G(s)=K/(s(s+2))。闭环传递函数Gc(s)=G(s)/(1+G(s))=[K/(s(s+2))]/[1+K/(s(s+2))]=K/[s(s+2)+K]=K/(s^2+2s+K)。系统为二阶系统，其闭环传递函数标准形式为Gc(s)=ωn^2/(s^2+2ζωns+ωn^2)。对比Gc(s)=K/(s^2+2s+K)与Gc(s)=ωn^2/(s^2+2ζωns+ωn^2)，可得：ωn^2=K。2ζωn=2。解得：ωn=sqrt(K)，ζ=1/(sqrt(K))。要求超调量σ%≤5%。对于二阶系统，超调量σ%=exp(-ζπ/sqrt(1-ζ^2))*100%。5%=exp(-ζπ/sqrt(1-ζ^2))*100%。0.05=exp(-ζπ/sqrt(1-ζ^2))。取对数：ln(0.05)=-ζπ/sqrt(1-ζ^2)。-2.9957=-ζπ/sqrt(1-ζ^2)。ζπ/sqrt(1-ζ^2)=2.9957。平方：(ζπ)^2/(1-ζ^2)=(2.9957)^2=8.972。ζ^2π^2=8.972(1-ζ^2)。ζ^2π^2=8.972-8.972ζ^2。(1+8.972/π^2)ζ^2=8.972。ζ^2=8.972/(1+8.972/π^2)=8.972/(π^2+8.972)≈8.972/(9.869+8.972)≈8.972/18.841≈0.475。ζ≈sqrt(0.475)≈0.689。ζ=1/sqrt(K)=0.689。sqrt(K)=1/0.689≈1.45。K≈(1.45)^2≈2.1025。因为超调量与阻尼比ζ的关系是单调递减的，ln(0.05)是一个固定值，所以对应的ζ是唯一的。我们求得ζ≈0.689对应的超调量约为5%。当K增大时，ζ减小，超调量σ%将大于5%。当K减小时，ζ增大，超调量σ%将小于5%。因此，为了满足σ%≤5%的要求，K的取值范围应该是K≤2.1025。四、分析题16.解：反馈控制系统根据反馈信号的极性分为正反馈和负反馈。负反馈：*特点：输出信号的一部分反送回到输入端，与输入信号相减后再作为系统的输入。*作用：能显著减小系统误差，提高系统的准确性和稳定性；可以改变系统的频率响应特性，如提高带宽、改善相位特性；可以通过引入校正装置来满足特定的性能指标要求（如提高快速性、减小超调量等）。*适用场合：绝大多数实际控制系统都采用负反馈，特别是要求高精度、高稳定性的场合，如工业自动化控制、机器人控制、飞行器控制等。正反馈：*特点：输出信号的一部分反送回到输入端，与输入信号相加后再作为系统的输入。*作用：正反馈会增强输入信号，可能导致系统输出不断增大直至饱和或振荡。通常用于特定的目的，如：*自激振荡器：利用正反馈产生持续的振荡信号。*比较环节：在某些特定结构中用于比较两个信号。*提高某些类型的系统增益。*适用场合：主要应用于需要产生振荡的电路（如振荡器），或者在某些特定反馈结构中作为辅助手段。在一般的稳定控制系统设计中较少使用，因为正反馈容易导致系统不稳定。17.解：利用伯德图（BodePlot）判断闭环系统稳定性是基于奈奎斯特稳定性判据的一种图解方法。伯德图将系统的频率响应特性用对数幅频特性（L(log|G(jω)|)和相频特性（φ(ω)）分别在对数坐标和线性坐标上绘制，便于分析低频、中频、高频特性以及增益裕度和相位裕度。基本原理：1.开环频率特性与闭环稳定性：根据奈奎斯特稳定性判据，闭环系统在s右半平面极点的个数Z与开环传递函数G(s)H(s)在s=jω轴上的穿越次数（方向）以及开环极点在右半平面的个数P的关系为：Z=P-N，其中N是G(s)H(s)的乃奎斯特曲线绕(-1,j0)点的净穿越次数（顺时针为负，逆时针为正）。因此，判断闭环稳定性转化为判断N和P。2.伯德图与乃奎斯特曲线：伯德图上的对数幅频特性L(log|G(jω)|和相频特性φ(ω)可以通过数学变换（如绘制乃奎斯特曲线的实部Re[G(jω)H(jω)]和虚部Im[G(jω)H(jω)]的曲线）或直接绘制。更常用的是，利用伯德图的特性来推断乃奎斯特曲线的行为。3.穿越与稳定性判断：*穿越jω轴：当|G(jω)H(jω)|=1时，系统处于临界稳定状态。在伯德图上，这意味着L(log|G(jω)|=0dB（零分贝线）。此时，系统的相角φ(ω)决定了闭环稳定性。*相位裕度(φp)：在|G(jω)H(jω)|=1的频率点ωc（增益交界频率），如果系统的相角φ(ωc)与-180°的差值（即相位裕度φp=180°-φ(ωc)）为正，则系统是稳定的；φp>0°表示系统具有稳定性裕度。φp=0°是临界稳定，φp<0°表示系统不稳定。*增益裕度(KG)：在系统相角φ(ω)=-180°的频率点ωg（相位交界频率），如果此时|G(jω)H(jω)|<1（即L(ωg)<0dB），则系统是稳定的；|G(jω)H(jω)|的值1/|G(jω)H(jω)|就是增益裕度KG=1/|G(jωg)H(jωg)|。KG>1（L(ωg)<0dB）表示系统具有稳定性裕度。KG=1（L(ωg)=0dB）是临界稳定，KG<1（L(ωg)>0dB）表示系统不稳定。*伯德图推断：虽然不能直接从伯德图读出Z值，但可以通过观察L(log|G(jω)|和φ(ω)的曲线来判断。*如果伯德图显示，在增益交界频率ωc之前，相位特性φ(ω)远未达到-180°（例如，φ(ωc)>-135°），则系统肯定稳定。*如果在增益交界频率ωc附近，φ(ω)正好接近-180°，但L(ωc)>0dB（即|G(jωc)|>1），则系统不稳定。*如果在增益交界频率ωc附近，φ(ω)正好接近-180°，但L(ωc)<0dB（即|G(jωc)|<1），则系统稳定，并且L(ωc)的值可以用来估计增益裕度KG=-L(ωc)（以dB为单位）。同样，如果φ(ωg)正好为-180°，可以观察L(ωg)的值来估计相位裕度（虽然通常在伯德图上分析的是ωc，而非ωg，但原理相通）。*简化分析（ωc附近）：如果L(log|G(jω)|在ωc附近穿过0dB线，观察穿过点处的相角φ(ωc)。若φ(ωc)≥-180°+ε（ε为一小正数，如φ(ωc)≥-179°或φ(ωc)≥-180°+5°等，取决于具体分析要求），则闭环系统稳定，且相位裕度φp=180°-φ(ωc)>0°。若φ(ωc)<-180°-ε，则系统不稳定。4.穿越(-1,j0)点的图示化判断：虽然伯德图本身不直接绘制乃奎斯特曲线，但可以通过分析伯德图来推断与稳定性相关的特性。*穿越频率与稳定性关系：可以通过分析L(log|G(jω)|在0dB线上的穿越情况（如果绘制了对应的乃奎斯特曲线或通过L(ω)与φ(ω)的关系进行分析）来判断稳定性。*穿越方向：在伯德图上分析稳定性，通常需要结合L(ω)与φ(ω)的关系，或者想象对应的乃奎斯特曲线的行为。*如果L(ω)在ωc附近穿过0dB线时，对应的φ(ω)足够接近-180°（例如，满足稳定性判据中的条件），则系统稳定，可以估计裕度。如果φ(ωc)远离-180°，则稳定性分析变得困难，需要结合L(ω)在其他频率点的行为进行综合判断，或者需要绘制出完整的乃奎斯特曲线（需要同时有L(ω)和φ(ω)的信息）来应用奈奎斯特稳定性判据。如果只能使用伯德图进行推断，则稳定性分析会相对复杂，需要具备将伯德图信息转化为乃奎斯特曲线行为的想象力或计算能力。*更直观的伯德图稳定性分析（结合L(ω)与φ(ω)关系或绘制乃奎斯特曲线）：*绘制L(ω)与φ(ω)的关系图（伯德图本身），或者直接绘制乃奎斯特曲线（需要L(ω)和φ(ω)数据或公式），然后应用奈奎斯特稳定性判据。*判断步骤：*绘制系统的乃奎斯