自动控制原理实验四：基于MATLAB/Simulink的系统频率特性分析与可视化

1. 实验背景与核心概念频率特性分析是自动控制领域最实用的工具之一它就像给系统做心电图——通过不同频率的输入信号观察系统的心跳反应。我在工业现场调试时经常用这种方法快速判断系统稳定性。这次我们要用MATLAB/Simulink这个工程师的瑞士军刀对二阶系统进行完整的频率特性体检。幅频特性和相频特性是频率特性的双胞胎指标。举个生活中的例子当你用不同频率摇晃弹簧玩具时玩具的摆动幅度幅频特性和滞后时间相频特性都会变化。实验中传递函数G(s)500/(s²10s500)描述的就是这样一个弹簧系统其中分母的s²项代表惯性s项代表阻尼常数项代表刚度。2. 实验环境搭建2.1 Simulink建模技巧打开MATLAB后在命令行输入simulink召唤建模界面。我习惯先建立空白模型CtrlN然后从库浏览器拖拽这些关键模块Sine Wave藏在Simulink/Sources里这是我们的信号发生器Transfer Fcn在Continuous分类下系统的心脏ScopeSinks分类里的示波器Bus CreatorSignal Routing分类下的信号合成器连接模块时有个小技巧按住Ctrl键点击模块端口可以自动连线。配置传递函数参数时直接在分子分母输入框填写[500]和[1 10 500]注意方括号表示多项式系数从高次到低次。2.2 参数调试经验设置正弦波频率时建议先扫频测试从1rad/s开始逐步提高到40rad/s。我在调试工业伺服系统时发现临界频率点往往出现在系统自然频率附近本例约22.36rad/s。具体操作双击Sine Wave模块频率参数设置为w*1后续通过变量w控制采样时间设为0.001秒保证波形光滑3. 三大频率特性图谱实战3.1 波特图绘制技巧波特图就像系统的体检报告单由幅频曲线和相频曲线组成。新版MATLAB2020以后的操作更智能用Inport和Outport替换原有信号接口点击APP → Linear Analysis Tool在Analysis I/Os选择Root level inports and outports点击Bode按钮瞬间生成读图要点幅频曲线纵轴单位是dB计算公式20log10(A)相频曲线纵轴是角度-180°是稳定性警戒线转折频率处会出现明显拐点3.2 奈奎斯特图判据奈氏图是判断稳定性的照妖镜绘制步骤与波特图类似最后选择Nyquist图表。重点观察曲线是否包围(-1,j0)点相位裕度是否大于30°幅值裕度是否大于6dB实测案例中当R3电阻从500kΩ降到200kΩ时奈氏曲线会明显向右扩张这时就需要调整控制器参数了。3.3 李沙育图观测方法李沙育图(XY图)是传统实验室的土办法但在现场调试中依然实用添加XY Graph模块在Sinks分类输入接系统输入输出信号调整x-y轴比例至1:1当输入输出同频时图形呈椭圆。椭圆倾斜方向反映相位差右倾滞后左倾超前长短轴比值就是幅值比。我在电机测试中常用这种方法快速估算相位差。4. 深度分析与工程应用4.1 频率特性物理意义改变输入频率w时可以观察到三个典型区域低频段wwn输出几乎无衰减相位滞后可忽略谐振区w≈wn幅值出现峰值相位快速变化高频段wwn幅值急剧衰减相位滞后接近180°这解释了为什么音响系统播放低音时喇叭振幅大而高音时几乎看不到振动。工业中的振动筛也是利用谐振原理工作。4.2 参数影响分析通过修改传递函数参数可以直观看到阻尼系数增大谐振峰值降低系统变得更迟钝自然频率提高整个曲线向右平移增益变化幅频曲线整体上下移动在PID调参时我常用这种方法观察参数调整效果。比如增大积分时间会使相位滞后增加这在波特图上会表现为相频曲线下移。5. 常见问题排查5.1 波形异常处理如果Scope显示波形失真检查正弦波模块的采样时间是否为0确认仿真步长设为auto或足够小尝试调整求解器为ode45遇到奈氏图不闭合的情况通常是频率范围不够宽建议将频率上限提高到10倍自然频率。5.2 精度提升技巧要提高测量精度在谐振区附近加密频率测试点使用Zoom FFT工具局部放大分析对数据取多次测量平均值有个小技巧在命令行输入ltiview可以调出更专业的线性系统分析工具里面包含更多高级分析功能。