MATLAB中if-else语句设置数组n=0时的值问题求助

你的代码存在三个核心问题：

1. 数组逻辑判断错误：if n == 0是对整个数组做布尔判断，只有当所有元素全为0时才会执行if分支，完全不符合你单独处理n=0元素的需求。
2. 运算顺序错误：原公式为 $\frac{\sin(0.1\pi n)}{\pi n}$，但你的代码写成(1./n*pi).*(sin(0.1*pi.*n))，等价于 $\frac{\pi \cdot \sin(0.1\pi n)}{n}$，和原公式结构完全相反，会导致计算结果错误。
3. 初始计算产生无效值：直接对包含0的n数组执行除法，会得到NaN（0/0）或Inf，后续再属于多余步骤。

修正后的代码方案 #

方案一：索引定位处理（推荐，高效清晰） #

先初始化数组，分别处理非0和0元素的位置：

% 替换为你实际使用的奇数N值
N = 5;
n = -(N-1)/2:(N-1)/2;
hlp = zeros(size(n)); % 初始化与n同尺寸的0数组

% 计算所有n≠0的元素
non_zero_pos = n ~= 0;
hlp(non_zero_pos) = sin(0.1*pi.*n(non_zero_pos))./(pi.*n(non_zero_pos));

% 将n=0的位置设为0.1
zero_pos = n == 0;
hlp(zero_pos) = 0.1;

方案二：利用极限特性简化 #

根据极限 $\lim_{x \to 0} \frac{\sin(x)}{x} = 1$，当n→0时，$\frac{\sin(0.1\pi n)}{\pi n} = 0.1 \cdot \frac{\sin(0.1\pi n)}{0.1\pi n} \to 0.1$，因此可以先直接计算整个数组，再替换无效值：

N = 5;
n = -(N-1)/2:(N-1)/2;
hlp = sin(0.1*pi.*n)./(pi.*n);
% 替换n=0时产生的NaN值
hlp(isnan(hlp)) = 0.1;

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