使用SymPy绘制高频正弦曲线出现图形失真问题求助

问题原因 #

你遇到的曲线失真，本质是SymPy默认的绘图采样点数量不足。对于10kHz的高频正弦波，在(-0.001, 0.001)的区间内包含20个完整周期，但SymPy的plot函数默认的自适应采样或固定点数无法覆盖每个周期的足够细节，导致曲线呈现锯齿状失真。而缩小定义域/降低频率时，区间内的周期数减少，采样点相对足够，所以失真会减轻。

解决方案1：调整SymPy绘图参数 #

手动关闭自适应采样，并指定足够多的采样点，确保每个周期至少有10-20个采样点（高频信号需要更高密度）：

from sympy import *
%matplotlib inline

t = symbols('t')
f_carrier = 10000

carrier = cos(2*pi*f_carrier * t)
# 关闭自适应采样，设置2000个采样点（可根据需求调整）
plot(carrier, (t, -0.001, 0.001), adaptive=False, nb_of_points=2000)

解决方案2：用Matplotlib结合lambdify绘制（更可控） #

将SymPy的符号表达式转换为数值函数，用NumPy生成密集的时间点，再通过Matplotlib绘图，这种方式能完全控制采样密度：

from sympy import *
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

t = symbols('t')
f_carrier = 10000
carrier = cos(2*pi*f_carrier * t)

# 把符号函数转为支持NumPy计算的函数
carrier_func = lambdify(t, carrier, 'numpy')

# 生成步长远小于信号周期的时间点（周期=1/10000=0.0001秒，步长设为0.00001秒）
t_vals = np.arange(-0.001, 0.001, 0.00001)
y_vals = carrier_func(t_vals)

plt.plot(t_vals, y_vals)
plt.xlabel('t')
plt.ylabel('cos(2πf_c t)')
plt.title('高频载波正弦曲线')
plt.show()

内容的提问来源于stack exchange，提问作者hspil