二阶系统的模拟实现以及时域分析基于Matlab&Pspice

本文主要包含以下内容：

一、典型二阶系统的参考模拟电路以及Pspice仿真
二、上述模拟线路对应结构框图以及传递函数推导
三、上述系统的时域分析

一、典型二阶系统的参考模拟电路以及Pspice仿真

下图为基于运放的积分线路和比例线路搭建的二阶系统：

调节C1R1、C2R3可以调节系统无阻尼的自然频率；调节Rfb/Rb可以调节系统的阻尼系数。
如下三张图片反应自然频率不变，阻尼系数增大的输出阶跃响应结果（Ri = 100K, Rfb = 50K:50K:150K)


如下三张图片反应阻尼系数不变，自然频率增大的输出阶跃响应结果（C:f分别取1uF，0.1uF，0.01uF)

二、上述模拟线路对应结构框图以及传递函数推导

由上述原理图整理得出如下结构框图：

用节点法列方程组求传递函数：
在Matlab运行如下脚本：

syms g1 g2 g3;
Q = [0 0 -1; g1 0 -g3; 0 g2 0];
P = [1; 0; 0];
W = inv(eye(3) - Q) * P;
W(3)

W(3）的值即为求得的化简后的传递函数：
(g1g2)/(g1g2 + g2*g3 + 1)
将G1 G2 G3的值带入上式,可以得到对应的传递函数:

备注传递函数有误，更正如下–2021/0929：
$$G(s)=\frac{\frac{1}{（RC{)}^2}}{s^2+\frac{Rf}{Ri}\ast \frac{1}{RC}s+\frac{1}{（RC{)}^2}}$$

对照二阶系统传递函数的标准形式可得:
自然频率 ωn = 1/RC
阻尼系数 ζ = Rfb/(2Ri)

三、上述系统的时域分析

有了以上化简后的传到函数,我们就可以在Matlab方便的进行时域分析了,
在Matlab中运行如下脚本,可得系统不改变自然频率 ωn,只改变阻尼系数 ζ 的阶跃响应曲线:

num = 100;
for sigma = 0: 0.25: 1.25
    den = [1, 2*sigma*10, 100];
    sys = tf(num, den);
    step(sys, 2)
    hold on
end

grid on
hold off

从图中可以方便看出,阻尼系数增加后,超调量会减小,但上升时间也会减小,稳定时间有一个相对最优的值，估值在0.75（理论值为0.707），我们回到Spice进行验证，设置Rfb = 150K ,得到如下响应波形：

同样，我们也可以设定改变自然频率 ωn,不改变阻尼系数 ζ 的响应曲线：

sigma = 0.5;
for i = 10: 10: 100
    den = [1, 2*sigma*i, i^2];
    sys = tf(i^2, den);
    step(sys, 2)
    hold on
end

grid on
hold off

可以发现，改变自然频率 ωn不影响超调量，只对上升时间和稳定时间由影响。