设x(n)=RN(n)，用matlab求x(n)的幅频特性和相频特性

在MATLAB中，要计算离散时间序列x(n) = Rn(n)的幅频特性（Magnitude Frequency Response, MFR）和相频特性（Phase Frequency Response, PFR），通常需要使用`fft`函数来获取频谱分析结果，然后分别计算幅值和相位。假设R是一个给定的实数序列，n是序列索引。

首先，我们需要对x(n)进行快速傅立叶变换（FFT）。这里有一个基本步骤：

% 假设R是输入的实数序列
R = ...; % 你需要提供具体的数值或生成一个示例

% 创建等长的全零序列，用于FFT
N = length(R); % 获取序列长度
n = (0:N-1)';
x = R .* exp(-1i * 2*pi * n/N); % x(n) = Rn(n)

% 对x进行FFT
X = fft(x);

% 幅频特性
X_mag = abs(X);
f = (0:N-1) / N; % 频率向量，单位一般为Hz
MFR = X_mag ./ max(X_mag); % 归一化到最大值，便于比较

% 相频特性
X_phase = angle(X);
PFR = unwrap(angle(X)) / (2*pi); % 将相位从弧度转换到角度，并去unwrap周期性

% 可视化
figure;
subplot(2,1,1)
plot(f, MFR, 'LineWidth', 1.5)
title('幅频特性')
xlabel('频率 (Hz)')
ylabel('幅度')

subplot(2,1,2)
plot(f, PFR, 'LineWidth', 1.5)
title('相频特性')
xlabel('频率 (Hz)')
ylabel('相位 (度)')