傅里叶级数的3D表示 包括源码

傅里叶级数的3D表示 包括源码

flyfish
傅里叶级数的基本形式 : y ( t ) = ∑ n = 1 , 3 , 5 , ... N 4 A n π sin ⁡ ( n π T t ) y(t) = \sum_{n=1,3,5,\ldots}^{N} \frac{4A}{n\pi} \sin\left(\frac{n\pi}{T} t\right) y(t)=n=1,3,5,...∑Nnπ4Asin(Tnπt)

其中, A A A 是振幅, T T T 是周期, n n n 是傅里叶级数的项数。

傅里叶级数中每一项的形式
y n ( t ) = 4 A n π sin ⁡ ( n π T t ) y_n(t) = \frac{4A}{n\pi} \sin\left(\frac{n\pi}{T} t\right) yn(t)=nπ4Asin(Tnπt)

这个公式表示傅里叶级数中每一项的贡献。

解释

t (时间) : 横坐标表示时间。
y (幅度) : 纵坐标表示傅里叶级数在该时间点的值。
f (频率) : 深度坐标表示傅里叶级数的项数对应的频率。

py 复制代码
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
# 定义傅里叶级数的参数
T = 2 * np.pi  # 周期
A = 1  # 振幅
N = 11  # 级数的最大项数

# 定义傅里叶级数的函数
def fourier_series(t, N, A, T):
    y = np.zeros_like(t)
    for n in range(1, N + 1, 2):  # 只考虑奇数项
        y += (4 * A / (np.pi * n)) * np.sin((n * np.pi / T) * t)
    return y

# 创建时间域和频率域的网格数据
t = np.linspace(0, 2 * T, 1000)
f = np.array([1 / T, 3 / T, 5 / T, 7 / T, 9 / T, 11 / T])
T_grid, F_grid = np.meshgrid(t, f)

# 创建3D图
fig = plt.figure(figsize=(12, 8))
ax = fig.add_subplot(111, projection='3d')

# 绘制傅里叶级数的时域图像
for n in range(1, N + 1, 2):
    y = (4 * A / (np.pi * n)) * np.sin((n * np.pi / T) * t)
    ax.plot(t, n / T, zs=y, label=f'n={n}')

# 设置坐标轴标签和刻度
ax.set_xlabel('t (时间)')
ax.set_ylabel('f (频率)')
ax.set_zlabel('y (幅度)')

# 设置坐标轴的刻度
ax.set_xticks([0, np.pi, 2 * np.pi])
ax.set_xticklabels(['0', r'$\pi$', r'$2\pi$'])
ax.set_yticks(f)
ax.set_yticklabels([r'$\frac{1}{T}$', r'$\frac{3}{T}$', r'$\frac{5}{T}$', r'$\frac{7}{T}$', r'$\frac{9}{T}$', r'$\frac{11}{T}$'])
ax.set_zticks([-A, 0, A])
ax.set_zticklabels([f'-{A}', '0', f'{A}'])

plt.title('傅里叶级数的3D表示')
plt.legend()
plt.show()
相关推荐
huoyingcg8 小时前
3D Mapping秀制作:沉浸式光影盛宴 3D mapping show
科技·3d·动画·虚拟现实
luoganttcc1 天前
FastPillars:一种易于部署的基于支柱的 3D 探测器
3d
工业3D_大熊1 天前
3D Web轻量化引擎HOOPS Communicator在装配件管理上的具体优势
3d·3d web轻量化·3d渲染·3d模型可视化·工业3d·web端3d可视化·3d复杂模型轻量化
在下胡三汉1 天前
3dmax批量转glb/gltf/fbx/osgb/stl/3ds/dae/obj/skp格式导出转换插件,无需一个个打开max,材质贴图在
3d·材质·贴图
xhload3d1 天前
智能网联汽车云控平台 | 图扑数字孪生
3d·gis·智慧城市·html5·webgl·数字孪生·可视化·工业互联网·车联网·智慧交通·智能网联·汽车云控
木木黄木木1 天前
使用HTML5和CSS3实现炫酷的3D立方体动画
3d·css3·html5
lb29171 天前
CSS 3D变换,transform:translateZ()
前端·css·3d
淮北4942 天前
ros调试工具foxglove使用指南三:在3d空间写写画画(Panel->3D ->Scene entity)
python·学习·3d·机器人
大有数据可视化2 天前
为什么可视化大屏越来越多应用3D元素呢?
3d
千野竹之卫2 天前
2025最新云渲染网渲100渲染农场使用方法,渲染100邀请码1a12
开发语言·前端·javascript·数码相机·3d·3dsmax