求教：如何使用MATLAB计算两个周期函数的Fourier coefficients

云原生阿岚

2026-6-27

MATLAB计算周期函数傅里叶系数的实现方案

核心逻辑

傅里叶系数的计算核心是对周期函数在一个完整周期内做积分，针对分段周期函数，可通过符号积分或数值积分两种方式实现，以下是具体操作步骤：

步骤1：定义周期函数

以常见的分段周期函数（如方波）为例，先明确函数周期T及单周期内的表达式。比如周期为2π的方波，在[-π, π)区间内的定义：

syms t n
% 单周期内：|t|<π/2时f(t)=1，否则为0
f(t) = piecewise(abs(t) < pi/2, 1, abs(t) >= pi/2, 0);

步骤2：计算傅里叶系数

傅里叶系数的标准公式为：

直流分量a₀：a₀ = (1/T) * ∫(单周期) f(t) dt
余弦系数aₙ：aₙ = (2/T) * ∫(单周期) f(t)*cos(nωt) dt（ω=2π/T为角频率）
正弦系数bₙ：bₙ = (2/T) * ∫(单周期) f(t)*sin(nωt) dt

以周期T=2π（ω=1）为例，代码实现：

T = 2*pi;
omega = 2*pi/T;

% 计算直流分量a₀
a0 = (1/T)*int(f(t), t, -pi, pi);
a0 = simplify(a0);

% 计算余弦系数aₙ
an = (2/T)*int(f(t)*cos(n*omega*t), t, -pi, pi);
an = simplify(an);

% 计算正弦系数bₙ
bn = (2/T)*int(f(t)*sin(n*omega*t), t, -pi, pi);
bn = simplify(bn);

运行后可得到解析形式的系数，比如方波的aₙ为(2*sin(n*pi/2))/(n*pi)，bₙ=0，与示例结果一致。

步骤3：复杂函数的数值积分方案

若函数过于复杂无法得到解析解，可使用数值积分函数integral：

% 定义匿名函数形式的周期函数
f = @(t) (abs(t) < pi/2).*1 + (abs(t) >= pi/2).*0;
T = 2*pi;
omega = 2*pi/T;
n = 1; % 以n=1为例计算系数

a0_num = (1/T)*integral(f, -pi, pi);
an_num = (2/T)*integral(@(t) f(t).*cos(n*omega*t), -pi, pi);
bn_num = (2/T)*integral(@(t) f(t).*sin(n*omega*t), -pi, pi);