请编写一个MATLAB程序，功能为用一个序列的单个DFT计算两个偶数长度实序列的DFT，要求前者长度为后者的一半。

为了在MATLAB中编写一个函数，用于计算两个偶数长度实序列的DFT，其中第一个序列的长度是第二个序列长度的一半，我们可以利用MATLAB内置的`fft`函数。假设我们有两个序列`x`和`y`，`x`的长度是`y`长度的一半，我们可以这样做：

function [X, Y] = compute_half_length_dft(x, y)
    % 检查输入是否为实数向量，且长度为偶数
    if isnumeric(x) && isnumeric(y) && all(isreal(x)) && all(isreal(y)) && mod(length(x), 2) == 0 && mod(length(y), 2) == 0
        % 确保x的长度是y的一半
        assert(length(x) == floor(length(y)/2));

        % 计算两个序列的DFT
        X = fft(x);
        Y = fft(y);

        % 如果需要复数输出，可以取绝对值，对于实数序列，不需要这一步
        % X = abs(X);
        % Y = abs(Y);
    else
        error('Input vectors must be real and of even length.');
    end
end

你可以使用这个函数像下面那样调用：

% 示例数据
x = [1 2 3 4];   % 长度为4，因为是偶数，所以y长度将是8
y = [5 6 7 8 9 10 11 12];

[X, Y] = compute_half_length_dft(x, y);
disp("DFT of x:");
disp(X);
disp("DFT of y:");
disp(Y);