0

0

Matplotlib 绘图中因坐标轴范围过大导致曲线显示为扁平线的问题解析

霞舞

霞舞

发布时间:2026-01-10 13:55:02

|

519人浏览过

|

来源于php中文网

原创

Matplotlib 绘图中因坐标轴范围过大导致曲线显示为扁平线的问题解析

matplotlib 在绘制多条数据时,若某条数据存在极端异常值,会自动扩展 y 轴范围,致使其他正常数据的波动在视觉上被压缩成一条“扁平线”,实则并非数据错误,而是坐标尺度失衡所致。

该问题本质是 Matplotlib 自动设置 y 轴范围(ylim)时,被某列数据中的极大/极小异常值主导,从而掩盖了其余数据的真实形态。在您的代码中:

ax1.plot(FTBS_nonlin_out[0.0])  # 初始状态:幅度约 [-5, 5]
ax1.plot(FTBS_nonlin_out[4])    # t=4 时刻:含数值爆炸(如 -2e18)

FTBS_nonlin_out[4] 列因 FTBS 格式在非线性对流方程中缺乏耗散或稳定性控制(尤其当 c * dt / dx > 1 时),极易产生数值震荡与溢出,导致个别元素达到 ~ -1e18 量级。而 ax1.plot() 默认将所有传入数据的最大最小值作为绘图边界,因此 y 轴被拉伸至 [-2e18, 5] —— 此时幅度仅 ±5 的初始状态 FTBS_nonlin_out[0.0] 在像素级图像中几乎无法分辨,呈现为视觉上的“水平直线”。

验证方法:添加以下调试语句即可确认:

print("t=0.0 range:", FTBS_nonlin_out[0.0].min(), FTBS_nonlin_out[0.0].max())
print("t=4 range:  ", FTBS_nonlin_out[4].min(), FTBS_nonlin_out[4].max())

解决方案(推荐组合使用):

  • 显式限制 y 轴范围(快速诊断 & 可视化修复):

    ax1.set_ylim(-30, 30)  # 或 plt.ylim(-30, 30)
  • 过滤异常值再绘图(保障展示合理性):

    # 仅绘制有限数值(排除 inf / nan / 极端离群点)
    valid_0 = FTBS_nonlin_out[0.0].replace([np.inf, -np.inf], np.nan).dropna()
    valid_t4 = FTBS_nonlin_out[4].replace([np.inf, -np.inf], np.nan).dropna()
    ax1.plot(valid_0.index, valid_0, label='t=0.0')
    ax1.plot(valid_t4.index, valid_t4, label='t=4')
  • 根本性改进数值稳定性

    • 检查 CFL 条件:确保 c * dt / dx ≤ 1(当前 c=1, dt=0.5, dx=1.0 → CFL=0.5,合规,但非线性项可能引发隐式超限);
    • 改用更稳定格式(如 Lax-Wendroff、ENO/WENO)或添加人工粘性;
    • 对更新公式做截断保护(防御性编程):
      u_new = row[t] - row[t] * (dt/dx) * (row[t] - signal.loc[i-1,t])
      signal.loc[i, t1] = np.clip(u_new, -1e3, 1e3)  # 防止爆炸

? 总结:这不是 Matplotlib 的 bug,而是数值模拟中“结果失稳 → 数据污染绘图范围 → 可视化失效”的典型链式现象。务必养成先检查数据分布、再绘图的习惯;在科研可视化中,plt.ylim() 不是妥协,而是必要控制手段;而长期来看,算法稳定性永远比绘图技巧更重要。

相关标签:

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
页面置换算法
页面置换算法

页面置换算法是操作系统中用来决定在内存中哪些页面应该被换出以便为新的页面提供空间的算法。本专题为大家提供页面置换算法的相关文章,大家可以免费体验。

403

2023.08.14

Java JVM 原理与性能调优实战
Java JVM 原理与性能调优实战

本专题系统讲解 Java 虚拟机(JVM)的核心工作原理与性能调优方法,包括 JVM 内存结构、对象创建与回收流程、垃圾回收器(Serial、CMS、G1、ZGC)对比分析、常见内存泄漏与性能瓶颈排查,以及 JVM 参数调优与监控工具(jstat、jmap、jvisualvm)的实战使用。通过真实案例,帮助学习者掌握 Java 应用在生产环境中的性能分析与优化能力。

19

2026.01.20

PS使用蒙版相关教程
PS使用蒙版相关教程

本专题整合了ps使用蒙版相关教程,阅读专题下面的文章了解更多详细内容。

61

2026.01.19

java用途介绍
java用途介绍

本专题整合了java用途功能相关介绍,阅读专题下面的文章了解更多详细内容。

87

2026.01.19

java输出数组相关教程
java输出数组相关教程

本专题整合了java输出数组相关教程,阅读专题下面的文章了解更多详细内容。

39

2026.01.19

java接口相关教程
java接口相关教程

本专题整合了java接口相关内容,阅读专题下面的文章了解更多详细内容。

10

2026.01.19

xml格式相关教程
xml格式相关教程

本专题整合了xml格式相关教程汇总,阅读专题下面的文章了解更多详细内容。

13

2026.01.19

PHP WebSocket 实时通信开发
PHP WebSocket 实时通信开发

本专题系统讲解 PHP 在实时通信与长连接场景中的应用实践,涵盖 WebSocket 协议原理、服务端连接管理、消息推送机制、心跳检测、断线重连以及与前端的实时交互实现。通过聊天系统、实时通知等案例,帮助开发者掌握 使用 PHP 构建实时通信与推送服务的完整开发流程,适用于即时消息与高互动性应用场景。

19

2026.01.19

微信聊天记录删除恢复导出教程汇总
微信聊天记录删除恢复导出教程汇总

本专题整合了微信聊天记录相关教程大全,阅读专题下面的文章了解更多详细内容。

160

2026.01.18

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号