简介面向自动控制原理教学与工程分析场景Matlab平台提供了一套伯德图坐标纸定制生成工具帮助自动化专业学生与工程师规范化绘制系统频率响应曲线用于幅频特性与相频特性的直观呈现。压缩包共4个文件包含3个m脚本对应1x1、2x1、2x2等不同比例的坐标纸绘制与1份Readme.txt使用说明文档整体仅5KB结构精简、即取即用。脚本以bodeplot函数为核心支持自定义频率范围、采样点数量与轴刻度样式可借助subplot灵活组织幅值图和相位图的版面布局同时覆盖传递函数建模tf/ss/zpk、图线样式调整线型、颜色、标签、外部数据读取及图像保存导出等完整作图流程。已有1034人学习下载适合正在学习频域分析法或需要制作伯德图坐标纸以完成实验报告、课程设计的人员也可作为控制理论学习中的实用绘图辅助。1. 打印伯德图手绘纸为什么还要自己写脚本有次做带阻滤波器的频域分析想按实验报告要求手绘一张伯德图结果找了半天现成坐标纸有的频率范围只到 1kHz带宽的转折区根本画不下有的自带半页水印打印机一出来灰蒙蒙一片还有的纵轴步长是 3dB 一格算好的 -20dB/dec 渐近线对不上格子。最后只好自己写脚本生成一张 A4 横向的半对数坐标纸频率范围、每十倍频程分几格、纵轴每格多少分贝全部自己定导出 PDF 打印出来线条干净、尺寸可控。Bode_PrintPaper.zip 这类打包工具的实质就是帮你在半小时内得到一张“合身”的伯德图手绘专用纸适合自动控制原理实验、频域分析课设和滤波器草图绘制的人。2. 半对数坐标纸的设计要点decade 宽度与网格分层一张能用的伯德图坐标纸不是简单画几条网格线就行。频率轴必须按对数展开纵轴按 dB 或角度线性展开网格线的疏密还要匹配手绘时读斜率的习惯。先把这个几何逻辑想清楚生成的纸才不会画到一半发现格子不对。2.1 频率轴必须对数展开Bode 坐标纸的几何原理伯德图的核心是频率响应随频率变化的趋势。一个系统可能要从 1Hz 看到 100kHz跨 5 个数量级用线性坐标画根本装不下低频段。对数坐标把频率相乘的关系变成纸上等距的位置每一段“十倍频程”decade在纸上占据固定宽度。这就是为什么坐标纸的横轴起始频率不能是 0而必须是一个正数比如 0.1、1、10。幅频特性纸的纵轴是 dB 值线性分布相频特性纸的纵轴是角度值也是线性分布。这样设计带来一个直接好处一阶惯性环节的幅频渐近线是 -20dB/dec在纸上是一条固定斜率的直线二阶振荡环节的渐近线是 -40dB/dec斜率恰好是 -20dB/dec 的两倍。如果你手头的坐标纸纵横比例不对或者频率轴被等距刻度冒充对数刻度这些斜率关系就全乱了。2.2 三个关键参数决定一张纸是否好用第一个参数是 decade 个数。拿 A4 横向纸举例去掉打印边距后有效画图宽度大约 270 毫米如果题目范围是 1Hz 到 100kHz正好 5 个 decade每段约 54 毫米这个宽度画渐近线很舒服如果范围只有 0.1Hz 到 1kHz4 个 decade每段可以放到 67 毫米。第二个参数是每个 decade 内的分格方式。教科书手绘常用 1-2-5 分格也就是十倍频程里只画 2、5 两条次要网格线显得干净如果需要在图上精确读出频率值比如做校正网络设计就用 2-9 全分格代价是线多了会显得密。第三个参数是纵轴步长。幅频纸常用每格 10dB 或 20dB相频纸常用每格 45° 或 30°。步长定得太大画 -40dB/dec 的线手抖一点就偏定得太小网格过密打印出来全是线。2.3 现成坐标纸 vs 脚本生成为什么选后者对比三种常见来源在线网页生成的坐标纸快是快但分辨率往往只有屏幕级打印出来线迹发虚还常带站点标识别人共享的 PDF 坐标纸省事但频率范围、纸张方向、纵轴范围都是定死的遇到非典型题目就尴尬自己用脚本生成全部参数可控输出 PDF 是矢量格式打印不会糊还能把生成脚本固化到自己的工具箱里下次改个范围直接再跑一次。时间成本上脚本第一次写半小时之后每次用时几秒钟长期看最划算。来源频率范围可控网格密度可控打印清晰度时间成本在线网页部分基本不可控一般低共享 PDF固定固定取决于源文件低脚本生成完全可控完全可控高矢量输出首次较高3. 用 Python 从零生成 Bode 坐标纸最小脚本与参数调优理解了设计逻辑落地就快了。我一般用 Python 的 matplotlib 来画因为它能精确控制对数坐标轴、自定义刻度位置和网格线型。下面这个脚本是我常用版本的骨架幅频纸和相频纸统一放在一个文件里跑一次命令行就能输出 PDF。3.1 生成脚本幅频纸与相频纸在一个文件里搞定import argparse import numpy as np import matplotlib.pyplot as plt from matplotlib.ticker import LogLocator, MultipleLocator # 中文字体修复避免轴标签变成方块 plt.rcParams[font.sans-serif] [SimHei, Microsoft YaHei] plt.rcParams[axes.unicode_minus] False def draw_paper(xlim, ylim, db_step, freq_divs, title, save_path): fig, ax plt.subplots(figsize(11.69, 8.27)) # A4 横向单位英寸 ax.set_xscale(log) ax.set_xlim(*xlim) ax.set_ylim(*ylim) # 主刻度每个 decade 的起点如 1, 10, 100, 1000 ax.xaxis.set_major_locator(LogLocator(base10, subs[1.0])) # 次刻度decade 内的细分格freq_divs 可传 [2,5] 或 2~9 ax.xaxis.set_minor_locator(LogLocator(base10, subsfreq_divs)) # 纵轴主次刻度 ax.yaxis.set_major_locator(MultipleLocator(db_step)) ax.yaxis.set_minor_locator(MultipleLocator(db_step / 2)) # 网格分层主网格粗实线次网格细虚线 ax.grid(True, whichmajor, linewidth0.8, colorblack) ax.grid(True, whichminor, linewidth0.4, colorblack, linestyle(0, (3, 2))) ax.set_title(title) plt.savefig(save_path, formatpdf, bbox_inchestight, pad_inches0.1) plt.close(fig) if __name__ __main__: parser argparse.ArgumentParser(descriptionGenerate Bode plot paper) parser.add_argument(--start, typefloat, default1.0) parser.add_argument(--stop, typefloat, default10000.0) parser.add_argument(--type, choices[magnitude, phase], defaultmagnitude) parser.add_argument(--db-step, typefloat, default10.0) parser.add_argument(--subs, typestr, default2,3,4,5,6,7,8,9, helpdecade 内分格逗号分隔如 2,5) parser.add_argument(--out, typestr, defaultbode_paper.pdf) args parser.parse_args() freq_divs [float(x) for x in args.subs.split(,)] ylim (-20.0, 60.0) if args.type magnitude else (-180.0, 0.0) title 幅频特性坐标纸 if args.type magnitude else 相频特性坐标纸 draw_paper((args.start, args.stop), ylim, args.db_step, freq_divs, title, args.out)这个脚本的核心有两点一是把横轴设置为对数刻度再分别用 LogLocator 指定主刻度和次刻度位置二是把网格线按主次分层主网格对应 1、10、100 这类 decade 边界次网格对应内部的 2、3、4 等分格。纵轴用 MultipleLocator 按固定步长布点主步长由 db_step 控制次步长自动取一半。这样画出来的纸主次分明不会一打印就糊成一团。参数方面--start和--stop直接决定十年频程的个数比如从 1 到 10000 正好 4 个 decade--subs控制 decade 内的网格密度传2,5就是 1-2-5 分格传2,3,4,5,6,7,8,9就是全分格--db-step控制纵轴主网格线的间隔。跑通的最小命令是python bode_paper.py --start 1 --stop 100000 --type magnitude --subs 2,5 --out mag_paper.pdf3.2 输出尺寸、分格数和线宽的推荐取值A4 横向的 figsize 用 (11.69, 8.27) 英寸对应 297×210 毫米这是最常用的报告用纸。纵向坐标纸在频率范围超过 6 个 decade 时才考虑否则横向排版更适合手写标注。线宽取值上主网格线 0.8pt、次网格线 0.4pt 是我试下来打印最稳的组合如果打印机分辨率偏低把次网格线提到 0.5pt虚线间隔用 (3, 2) 而不是默认的短点线保证碳粉能均匀附着。分格数选择有个经验写课程报告建议全分格方便老师检查你读出的转折频率自己工程草图用 1-2-5 分格图面清爽画渐近线更快。--db-step则要看纵轴范围如果增益范围是 80dB用 20dB 一格会太疏画出来的 -20dB/dec 斜率几乎是平线用 10dB 一格刚刚好。3.3 从 PNG 到 PDF打印场景的正确导出方式同样的图导出成什么格式差别很大。PNG 是像素图放大打印时线条边缘会发虚PDF 是矢量格式文字和线条由数学描述打印设备直接渲染任意缩放都清晰。在 matplotlib 里只要 savefig 的路径后缀是.pdf就会自动走矢量后端。如果一定要在报告里插入位图比如学校模板只认图片那就把 dpi 提到 600 再导出 PNG网格线才不至于断线。4. Bode 坐标纸打印避坑5 条让图纸白做的典型问题脚本能跑通只是第一步真正把纸打印出来、拿尺子量一量才会碰到几个让人抓狂的问题。下面这几条都是实操里反复遇到的按“现象、原因、解决”列出来。4.1 打印后主次网格糊成一片现象在屏幕上明明主网格粗、次网格细打印出来却全是同样粗细的黑线甚至虚线变成了实线。原因次网格线宽设得太细比如低于 0.3pt打印机碳粉或喷墨扩散后和主网格线分辨不开另外虚线样式如果间隔太密比如 (1, 0)打印分辨率不够时会被渲染成连续线段。解决把次网格线宽调到 0.4 到 0.5pt虚线样式用 (3, 2)即 3pt 实线加 2pt 空白的节奏并提前用打印预览放大检查主次是否可区分。4.2 首尾 decade 宽度不一致现象生成的纸上 1-10 和 10-100 两段宽度不一样频率轴像被拉伸过。原因matplotlib 在 log 刻度下如果坐标轴留了 margin 或者 xlim 设置后又被自动扩展首尾会多出多余的空白区导致相邻 decade 不等宽。解决设置 xlim 之后再调用一次ax.margins(0)去掉自动留白频率上限和下限也尽量取 10 的整数幂比如 1 和 10000不要取 2.7 这种奇数边界。4.3 图上的中文标签全部显示为方块现象标题写的是“幅频特性坐标纸”生成 PDF 后变成一排豆腐块。原因matplotlib 默认字体不含中文字形系统也没有自动回退到中文字体。解决在脚本开头加两行指定 SimHei 或 Microsoft YaHei同时设置axes.unicode_minus为 False否则纵轴的负号也会变成方块。注意这行配置必须在创建 figure 之前执行。4.4 频率刻度被“手动修正”成等距现象有人拿普通坐标纸脚本硬改把 1 到 10 的刻度等距排开看起来频率轴很匀称实际读频率完全错误。原因对数坐标的本质就是低频疏、高频密同一 decade 内 1 到 2 的距离大于 9 到 10 的距离。这个特性不是 bug也不应该被修正。解决认准 LogLocator 生成的刻度不要手动传等距数组给 set_xticks。检查方法很简单看任意一个 decade 内部2 的位置应该离 1 较远离 10 较近这是判断坐标纸真伪的快速手段。4.5 打印缩放让坐标纸的比例尺完全失真现象屏幕上看坐标纸长宽正好打印出来却变小了用尺子量刻度间距和设计值对不上画的斜率线自然也不对。原因打印对话框默认选了“适合页面”把 A4 内容缩放到实际纸张边缘坐标纸被整体缩放。解决打印时在属性里选“实际大小”或“缩放 100%”不要选“适应边框”另外可以在坐标纸角落加一个 10mm 长的比例尺线段打印后先拿尺子核对长度不对就知道打印设置出了问题。5. 把生成脚本变成管用的小工具参数化、自校验和打印验证脚本跑通后下一步是让它更顺手。我把生成命令按使用场景拆成两条幅频纸用--type magnitude相频纸用--type phase频率范围从 0.1Hz 到 100kHz 这种跨 6 个 decade 的场景也只需换参数。另外加了一个自检函数生成之后自动验证相邻 decade 在纸上的像素宽度是否一致避免 margin 干扰导致图纸失真。def check_decade_width(fig, ax, xlim): fig.canvas.draw() p1, p2, p3 ax.transData.transform([(xlim[0], 0), (xlim[0] * 10, 0), (xlim[0] * 100, 0)]) d1 p2[0] - p1[0] d2 p3[0] - p2[0] return d1, d2 # 调用示例在生成 PDF 前检查 fig, ax plt.subplots(figsize(11.69, 8.27)) # ... 各种坐标轴设置 ... w1, w2 check_decade_width(fig, ax, (1, 10000)) assert abs(w1 - w2) 1.0, decade width mismatch!这个自检函数利用 transData 把数据坐标映射到显示坐标算出相邻两个 decade 的像素间距。正常对数坐标下两者应完全相等如果差超过 1 个像素说明页面上有额外留白或缩放需要回头检查 xlim 和 margins。我习惯在保存 PDF 前跑一次这个断言保证输出文件不会带着隐性比例错误。实际使用时还有一个小技巧生成坐标纸后不要急着打印先在 PDF 阅读器里用测量工具量一段 1-2 之间的宽度记下数值打印后用直尺再量一次两次一致才算校准通过。我现在每次画手绘伯德图都走这套流程坐标纸固定用脚本生成的一版配合 45° 三角板画渐近线转折频率和斜率都能对得上题目要求。这套习惯帮我少返工了很多次希望帮到你。本文还有配套的精品资源点击获取