semilogx(f(1:N/2), 20*log10(mag(1:N/2))); % 绘制幅频特性曲线，Index exceeds the number of array elements. Index must not exceed 3.

这个错误提示表明程序试图访问一个超出数组范围的元素，数组的长度不足以满足程序的要求。

在这个问题中，错误出现在绘制幅频特性曲线的代码行中，可能是由于数组长度N过小导致的。具体来说，如果N的值小于等于3，那么在计算幅度谱和绘制幅频特性曲线时就会出现这个错误。

解决这个问题的方法是增加N的值，或者修改绘图代码以确保数组不会超出范围。例如，可以修改代码如下：

N = 1024; % 序列长度
Fs = 1000; % 采样频率
alpha = 0.8; % 衰减因子
f0 = 50; % 正弦波频率
phi = pi/4; % 正弦波相位

n = 0:N-1; % 时间序列
x = alpha.^n .* sin(2*pi*f0/Fs*n + phi); % 衰减正弦序列

X = fft(x); % 计算FFT
f = (0:N-1)/N*Fs; % 频率序列
mag = abs(X); % 幅度谱
phase = angle(X); % 相位谱

figure; % 新建一个图形窗口
semilogx(f(1:N/2+1), 20*log10(mag(1:N/2+1))); % 绘制幅频特性曲线
xlabel('Frequency (Hz)');
ylabel('Magnitude (dB)');

在这个修改后的代码中，将N设置为1024，可以保证序列长度足够大，不会导致数组超出范围。另外，为了防止出现数组超出范围的问题，将绘图代码修改为`semilogx(f(1:N/2+1), 20*log10(mag(1:N/2+1)))`，这样可以确保数组不会超出范围。