基于Python的DFT计算余弦函数实部虚部相关问题咨询

1. 公式中12.5的参数含义 #

标准连续余弦信号的表达式为：

f(t) = cos(2πft + φ)
其中：

• f 是信号的频率，单位为Hz（次/秒）
• φ 是信号的初始相位，你给出的公式中默认φ=0
• 2πf 是角频率，单位为rad/s

所以你公式中的12.5就是信号的频率，含义为该余弦信号每秒完成12.5次周期性振动。结合你的采样参数：采样间隔0.01s，100个采样点的总时长为1s，刚好可以采集到12.5个完整的余弦周期。

2. 同类案例及实现示例 #

这类基础离散傅里叶变换计算是数字信号处理的入门级练习，大量开发者都做过同类实现。你可以直接用Python的numpy库完成计算，示例代码如下：

import numpy as np

# 定义参数
delta = 0.01  # 采样时间间隔，单位s
N = 100  # 采样点数
f0 = 12.5  # 余弦信号频率，单位Hz
Fs = 1 / delta  # 采样率，单位Hz

# 生成时间序列和余弦信号
t = np.arange(N) * delta
x = np.cos(2 * np.pi * f0 * t)

# 计算离散傅里叶变换
X = np.fft.fft(x)
# 提取DFT结果的实部和虚部
dft_real = X.real
dft_imag = X.imag

# 打印前20个频点的实部虚部示例
for i in range(20):
    print(f"频点索引{i}，实部：{dft_real[i]:.4f}，虚部：{dft_imag[i]:.4f}")

补充说明 #

• 你的余弦信号是实信号，本身没有虚部，你提到的实部虚部指的是该信号经过DFT变换后的频域结果的实部和虚部
• 本次采样的频率分辨率为Fs/N = 1Hz，12.5Hz刚好落在第12和13个频点中间，会出现频谱泄漏现象，因此DFT结果的虚部不会为0，如果换成整周期采样的12Hz或者13Hz信号，对应频点的虚部会趋近于0。

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