用MATLAB实现，求信号x（n）=R(n)的离散时间傅立叶变换 DTFT，画出信号的幅度谱与相位谱。

在MATLAB中，你可以使用`fft`函数来计算离散时间傅立叶变换(DTFT)，对给定的信号`x(n)`。首先，你需要将信号存储在一个向量中，然后按照以下步骤操作：

1. **准备数据**:
假设`x`是一个长度为`N`的一维数组表示连续序列`x(n)`。

x = % 输入你的信号数据

2. **计算DTFT**:
使用`fft`函数计算X(k)，即DTFT的离散频谱，其中k是从0到N-1的频率指数。

X = fft(x);

3. **调整幅度谱和相位谱**:
- **幅度谱**: 取绝对值得到幅度响应，通常归一化。
- **相位谱**: 相位是弧度形式，需要转换成角度，并加上初始相位(通常是0)。

X_mag = abs(X);
X_phase = unwrap(angle(X)) * (180 / pi); % unwrapping for continuous phase

4. **绘制幅度谱和相位谱**:
对于幅度谱，可以使用`plot`函数；对于相位谱，你可能希望先创建一个新的频率轴（` linspace(0, Nyquist, N)`），然后在同一图上绘制两者。

% 创建频率轴
Fs = 1; % sample rate (假设是1)
f = Fs*(0:(N-1))/N; % frequency vector

% 绘制幅度谱
figure;
plot(f, X_mag, 'b', 'LineWidth', 2);
xlabel('Frequency [Hz]');
ylabel('Magnitude');
title('Amplitude Spectrum');

% 绘制相位谱
figure;
plot(f, X_phase, 'g', 'LineWidth', 2);
xlabel('Frequency [Hz]');
ylabel('Phase [°]');
title('Phase Spectrum');

记得替换`Fs`为你实际的采样率，并调整图表标题和标签内容以适应你的需求。运行上述代码后，你应该会看到两个图表，分别展示了信号`x(n)`的幅度谱和相位谱。