MATLAB傅里叶级数绘图报错：矩阵维度不匹配求助

常见错误原因 #

你遇到的维度不匹配问题，核心是数组运算维度不一致：

• 原始的*、/是矩阵运算，要求操作数维度严格符合矩阵乘法规则；改成.*、./元素级运算后，若参与运算的数组长度不同，依然会报错。
• 多数情况是初始化的近似值数组与x轴数组维度不匹配，或循环中生成的级数项维度错误。

可运行示例代码（以周期2π的方波为例，n=4） #

% 生成x轴数据，覆盖一个周期（[-π, π]），确保足够密集
x = linspace(-pi, pi, 1000);
% 初始化近似值数组，与x维度完全一致
f_approx = zeros(size(x));

% 计算n=4的傅里叶级数项（方波为奇函数，仅含正弦项）
% n=4表示取前4项中的有效奇数项：k=1,3
for k = 1:2:4
    % 每一项都是与x同维度的数组
    term = (4/(pi*k)) * sin(k*x);
    f_approx = f_approx + term;
end

% 绘制近似曲线与原函数
plot(x, f_approx, 'LineWidth', 1.5);
hold on;
% 绘制原方波函数
f_original = sign(x);
f_original(x==0) = 1; % 处理x=0的间断点
plot(x, f_original, '--', 'LineWidth', 1.5);
legend('n=4傅里叶近似','原方波');
xlabel('x');
ylabel('f(x)');
title('傅里叶级数近似(n=4)');
hold off;

关键修改要点 #

1. 维度对齐：初始化f_approx时必须用zeros(size(x))，确保其与x长度一致，后续累加的每一项term也会自动与x同维度。
2. 级数项范围：明确n=4对应的级数项数量——若函数是奇函数（如方波），仅含奇次谐波，因此k取1、3（到4为止的奇数）；若为偶函数则取余弦项，需调整循环逻辑。
3. 元素级运算：确保所有逐元素计算的操作使用.*、./，标量与数组相乘时MATLAB会自动广播，无需额外处理。

排查维度问题的通用步骤 #

• 运行代码前，用length(x)、length(f_approx)检查核心数组的长度是否一致。
• 循环中单独输出size(term)，确认每一项的维度与x匹配。
• 核对傅里叶级数的数学表达式，避免因公式错误（如遗漏系数、错误的项数范围）导致运算维度异常。

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