如何使用MATLAB的bode函数进行伯德图分析？
一、bode函数基础作用
bode是Control System Toolbox核心频域分析函数，用于绘制线性系统幅频特性、相频特性两张子图，横轴为对数频率，纵轴分别是分贝幅值、相位角度。
通过伯德图可直观判断系统稳定裕度、谐振峰值、带宽、低频稳态误差、高频衰减特性，是校正装置设计、稳定性判据（奈奎斯特、对数稳定判据）的必备工具。
前提：必须用tf、ss、zpk生成LTI线性系统模型，缺少控制系统工具箱会提示函数未定义。
二、系统建模两种常用写法
1. 分子分母系数向量建模（课程实验常用）
传递函数G(s)=10s2+2s+10G(s)=\frac{10}{s^2+2s+10}G(s)=s2+2s+1010
num=[10];den=[1210];G=tf(num,den);2. s符号建模，复杂分式更简洁
s=tf('s');G=10/(s^2+2*s+10);多因式相乘用conv拼接分母
典型一阶积分加惯性环节G(s)=3s(s+5)G(s)=\frac{3}{s(s+5)}G(s)=s(s+5)3
num=[3];den=conv([10],[15]);G=tf(num,den);三、bode基础绘图语法
1. 最简调用，自动生成频率区间
clear;clc;s=tf('s');G=16/(s^2+4*s+16);figure;bode(G);grid on;title('二阶振荡系统伯德图');图像分为上下两幅：
上图：幅值（dB），下图：相位（°）；横轴自动对数刻度。
右键画布选择「特性」，可读取指定频率下幅值、相位数值。
2. 自定义仿真频率范围
默认频率由MATLAB自适应，需要限定起止频率时，构造对数频率向量w：
% 频率区间0.1~100 rad/s，对数等分w=logspace(-1,2,1000);bode(G,w);grid on;logspace(a,b,n)：生成10^a ~ 10^b 之间n个对数均匀分布频率点，点数越多曲线越平滑。
四、高频进阶实操场景
场景1：多系统伯德图同图对比
对比校正前后、不同阻尼、不同增益系统频域特性，自定义线型颜色区分：
clear;clc;s=tf('s');G1=10/(s^2+2*s+10);G2=25/(s^2+2*s+25);w=logspace(-1,2,800);figure;bode(G1,w,'r-',G2,w,'b--');grid on;legend('G1 阻尼大','G2 阻尼小','Location','best');title('两组二阶系统伯德图对比');场景2：仅提取幅值、相位数据，不绘图
需要手动计算增益裕度、相位裕度、带宽时，返回频率wout、幅值mag、相位phase：
s=tf('s');G=16/(s^2+4*s+16);[mag,phase,wout]=bode(G);% mag为三维数组，提取一维向量mag_dB=20*log10(squeeze(mag));phase_deg=squeeze(phase);% 手动绘制伯德图subplot(2,1,1);semilogx(wout,mag_dB);grid on;ylabel('幅值 / dB');title('手动绘制幅频曲线');subplot(2,1,2);semilogx(wout,phase_deg);grid on;xlabel('频率 ω / (rad/s)');ylabel('相位 / °');title('手动绘制相频曲线');场景3：闭环系统伯德图分析
工程中常分析闭环系统带宽、谐振峰值，先用feedback构造闭环再绘图：
clear;clc;s=tf('s');G_open=5/(s*(s+3));% 开环传递函数G_close=feedback(G_open,1);% 单位负反馈闭环w=logspace(-1,3,1000);figure;bode(G_open,w,'r',G_close,w,'b');grid on;legend('开环伯德图','闭环伯德图');title('开环与闭环系统伯德图对比');场景4：离散系统伯德图
离散传递函数指定采样时间Ts，tf末尾传入采样周期：
num=[0.2];den=[1-0.70.1];Ts=0.1;Gz=tf(num,den,Ts);figure;bode(Gz);grid on;title('离散控制系统伯德图');五、搭配裕度函数快速判定稳定性
bode仅绘图，margin可直接在伯德图标注增益裕度GM、相位裕度PM，工程调试最常用：
s=tf('s');G=3/(s*(s+1)*(s+5));figure;margin(G);grid on;title('带稳定裕度标注的伯德图');运行后图像自动标出0dB穿越频率、-180°穿越频率，同时命令行输出裕度数值。
六、常见报错与实操注意事项
- 未定义函数’bode’
未安装Control System Toolbox，命令行输入ver control校验，附加功能管理器安装对应工具箱。 - 曲线出现尖峰、折线粗糙
频率点数量太少，logspace第三个参数调大至800~1500，曲线会更顺滑。 - 多系统绘图线条重叠
手动指定线型：实线-、虚线--、点线:，搭配r/g/b/k区分颜色。 - 分不清开环/闭环伯德图用途
开环伯德图用于稳定裕度、校正装置设计；闭环伯德图用于分析系统响应速度、谐振、带宽。 - 幅值数值不是dB
bode输出原始幅值mag，转为分贝必须用公式20*log10(mag)。
七、完整可直接运行综合示例脚本
% bode伯德图完整仿真代码clear;clc;close all;% 1. 建立三阶开环传递函数s=tf('s');G_open=4/(s*(s+2)*(s+6));% 2. 自定义频率区间绘制伯德图w=logspace(-2,3,1200);figure(1);bode(G_open,w);grid on;title('三阶系统开环伯德图');% 3. 绘制带稳定裕度的伯德图figure(2);margin(G_open);grid on;title('标注增益裕度、相位裕度伯德图');% 4. 构造闭环系统对比频域特性G_close=feedback(G_open,1);figure(3);bode(G_open,w,'r-',G_close,w,'b--');grid on;legend('开环','闭环');title('开环、闭环伯德图对比');% 5. 提取频域数据[mag,phase,w_data]=bode(G_open);mag_dB=20*log10(squeeze(mag));disp('前10组频率-幅值数据：');disp([w_data(1:10),mag_dB(1:10)]);