Spyder运行matplotlib简谐振子绘图代码弹窗空白无响应

这个故障和数值计算逻辑无关，是Matplotlib绘图调用逻辑、Tk后端事件循环和Spyder运行环境的冲突导致的：

• 代码没有为第一张位移曲线显式创建独立画布，所有循环内的plt.plot()默认绘制在隐式生成的第一个画布上，调用plt.figure()新建第二张收敛图时，第一张画布的渲染指令没有被完整注册到Tk事件队列，就会出现只有窗口框架、没有内容的空壳
• 代码中强制切换TkAgg后端的时机不对：Matplotlib后端必须在导入pyplot后第一时间切换，否则部分后端组件初始化不完整，加上Spyder默认不会自动接管外部Tk窗口的事件循环，渲染任务堆积就会导致第二个窗口无响应
• 循环内注释掉的plt.show(block=False)如果放开反而会加剧事件循环冲突，不要启用

1. 先配置Spyder环境：打开工具>偏好设置>IPython控制台>图形，将图形后端设置为TkAgg，保存后重启IPython内核，删掉代码里强制切换后端的plt.switch_backend('TkAgg')行，避免后端重复初始化
2. 绘图前显式声明两个独立的画布对象，分别对应位移曲线、收敛验证曲线，给每个画布添加必要的标签、标题等元素，明确画布边界
3. 所有绘图逻辑完成后统一调用一次plt.show()触发渲染，不要在循环中插入show调用

#!/usr/bin/env python3
# 蛙跳法求解简谐振子并验证收敛性
import matplotlib.pyplot as plt
import numpy

k      = 1.0 # 弹簧劲度系数
m      = 1.0 # 振子质量
cycles = 2.0 # 积分覆盖的振动周期数
x0     = 1.0 # 初始位移
v0     = 0.0 # 初始速度

def leapfrog( steps ):
    """用蛙跳法积分简谐运动方程，依赖全局定义的质量m、劲度系数k"""
    omega = (k/m)**0.5
    delta = 2.0*cycles*numpy.pi/omega/steps
    x     = numpy.empty( steps+1 )
    v     = numpy.empty( steps+1 )
    t     = numpy.empty( steps+1 )
    t[0]  = 0.0
    x[0]  = x0
    # 蛙跳法半步起步修正初始速度
    v[0]  = v0 + 0.5*delta*x[0]*omega**2.0
    for i in range(steps):
        t[i+1] = t[i] + delta
        v[i+1] = v[i] - omega**2*delta*x[i]
        x[i+1] = x[i] + delta*v[i+1]
    return t, x, v

def l2_error_norm( t , x ):
    """计算数值解相对于解析解的L2相对误差范数"""
    steps  = len( x ) - 1
    omega  = (k/m)**0.5
    l2_err = 0.0
    l2     = 0.0
    for i in range(steps):
        x_exact = x0*numpy.cos( omega*t[i] )
        l2_err += (x[i] - x_exact)**2.0
        l2     += x[i]**2.0
    return (l2_err/l2)**0.5

# 积分参数配置
n        = 14
steps    = 8
l2_error = numpy.empty( n )
delta    = numpy.empty( n )

# 绘制第一张图：不同步长下的位移随时间变化曲线
plt.figure(figsize=(8,4))
for i in range(0,n):
    t, x, v     = leapfrog( steps )
    delta[i]    = (k/m)**0.5*(t[1]-t[0])
    l2_error[i] = l2_error_norm( t , x )
    plt.plot( t , x, label=f'步长对应步数={steps}' )
    steps *= 2
plt.xlabel('时间 t (s)')
plt.ylabel('位移 x (m)')
plt.title('蛙跳法积分简谐振子：不同步长位移曲线')
plt.legend(fontsize=8, ncol=2)
plt.grid(alpha=0.3)

# 绘制第二张图：步长-误差双对数收敛验证曲线
plt.figure(figsize=(6,5))
plt.loglog( delta , l2_error , 'o', label='蛙跳法实际计算误差' )
plt.loglog( delta , l2_error[0]*(delta/delta[0])**1.0, '--', label='1阶收敛参考线' )
plt.loglog( delta , l2_error[0]*(delta/delta[0])**2.0, '--', label='2阶收敛参考线' )
plt.loglog( delta , l2_error[0]*(delta/delta[0])**3.0, '--', label='3阶收敛参考线' )
plt.xlabel('无量纲积分步长')
plt.ylabel('L2相对误差范数')
plt.title('蛙跳法收敛阶验证')
plt.legend()
plt.grid(alpha=0.3, which='both')

# 统一触发窗口渲染
plt.show()

如果修改后仍出现窗口无响应：

• 先打开任务管理器，结束所有残留的pythonw.exe进程，清理之前崩溃的Matplotlib窗口进程
• 确认Spyder使用的Python环境中安装的Matplotlib、Tkinter版本匹配，不要混用conda和pip安装的Matplotlib包
• 可以先把后端临时改成inline测试代码逻辑是否正常，如果inline模式下能正常出图，就说明问题完全出在Tk后端的事件循环配置上，和代码计算逻辑无关

内容的提问来源于stack exchange，提问作者Bríd Moloughney