写出用MATLAB验证吉布斯现象的代码

下面是一个简单的MATLAB代码，用于验证吉布斯现象：

% 生成一个矩形波信号
x = zeros(1, 100);
x(20:40) = 1;
x(60:80) = -1;

% 进行傅里叶级数展开
N = length(x);
X = fft(x);
X = X(1:N/2+1);
X(2:end-1) = 2*X(2:end-1);

% 用ifft计算重建信号
xn = ifft(X);
xn = real(xn);

% 绘制原始信号和重建信号
subplot(2,1,1);
plot(x);
title('原始信号');
subplot(2,1,2);
plot(xn);
title('重建信号');

上述代码生成了一个矩形波信号，并进行了傅里叶级数展开和重建。通过绘制原始信号和重建信号的图像，可以观察到在信号的边缘处出现了明显的震荡，这就是吉布斯现象的表现。