写一段matlab程序计算给定函数lambda(x)在一个特定区间[a, b]内的离散傅立叶级数（Discrete Fourier Series, DFT），并展示其幅度谱、实部、虚部和相位谱。

在MATLAB中，计算离散傅里叶级数（DFT）通常使用`fft`函数。以下是一个基本的示例程序，它假设你有一个名为`lambda`的函数定义在区间`[a, b]`上，并且你想对其进行DFT分析：

% 定义函数
function y = lambda(x)
    % 这里替换成你的实际函数定义
    % 例如，y = sin(2*pi*x); 或 y = x.^2;

end

% 定义区间和采样点
a = 0;           % 区间开始
b = 1;           % 区间结束
n_points = 100;   % 采样点数
dx = (b - a) / n_points; % 间隔

% 创建x向量
x = linspace(a, b, n_points + 1);
y = lambda(x(1:end-1));  % 减去最后一个点避免重复

% 计算离散傅立叶变换
Y = fft(y);

% 获取幅度谱、实部、虚部和相位谱
abs_Y = abs(Y);          % 幅度谱
real_Y = real(Y);        % 实部
imag_Y = imag(Y);        % 虚部
phase_Y = angle(Y);      % 相位谱 (以弧度表示)

% 显示结果
figure;
subplot(2, 2, 1), plot(x, y, 'LineWidth', 1.5), title('Original Function');
subplot(2, 2, 2), plot(abs_Y, 'LineWidth', 1.5), title('Magnitude Spectrum');
subplot(2, 2, 3), plot(real_Y, 'LineWidth', 1.5), title('Real Part');
subplot(2, 2, 4), plot(phase_Y, 'LineWidth', 1.5), title('Phase Spectrum');

% 如果需要，也可以将相位转换为角度范围(-π, π]
phase_degrees = unwrap(angle(Y)) * 180 / pi;