给出一个matlab代码，来计算给定数据文件的DFT，并绘制幅度和相位谱

在MATLAB中，计算离散傅立叶变换（Discrete Fourier Transform，DFT）并绘制幅度和相位谱通常会涉及到`fft`函数。假设你有一个名为`data.txt`的数据文件，其中包含一列实数数据。以下是一个简单的示例代码：

% 读取数据文件
data = textread('data.txt', '%f', 'delimiter', '\n'); % 假设每行是一个数据点

% 计算DFT
N = length(data);
dft = fft(data);

% DFT结果的幅度谱
abs_dft = abs(dft); % 获取复数的模
freqs = (0:N-1)*(1/(N)); % 创建频率向量

% 绘制幅度谱
figure;
plot(freqs, abs_dft);
title('Amplitude Spectrum');
xlabel('Frequency (Normalized)');
ylabel('Magnitude');

% 相位谱（仅保留第一半部分，因为第二半是对称的）
phase_dft = angle(dft); % 获取复数的相位
phase_dft = phase_dft(1:N/2 + 1)'; % 取前半部分并转置
figure;
plot(freqs(1:N/2), unwrap(phase_dft));
title('Phase Spectrum');
xlabel('Frequency (Normalized)');
ylabel('Phase (Radians)');

% 函数unwrap()用于消除相位跳跃

注意，这个例子假设数据是一维的，并且是偶数长度。如果你的数据是其他形式，例如二维数组或者奇数长度，代码需要相应调整。