写一段matlab代码，使用傅里叶分析方法分析某个WAV文件的振幅、频率吧和初始相位

在MATLAB中，你可以使用傅里叶变换对音频信号进行频谱分析。以下是一段代码示例，用于读取一个WAV文件并计算其幅度、频率以及相位：

% 读取音频文件
[x1, fs_n] = audioread('test\Far_common.wav'); % 音频数据和采样率

% 傅立叶变换（包括幅度和频率）
x1_fd = fft(x1); % 傅里叶变换
n = 0:length(x1)-1;
f = n * fs_n / length(x1); % 频率轴

% 提取幅度
x1_abs = abs(x1_fd);

% 计算相位（仅取正半部分，因为负半部分是对称的）
ph = 2 * angle(x1_fd(1:length(x1)/2));
ph = ph * 180 / pi;

% 绘制幅值谱
plot(f(1:length(x1)/2), x1_abs(1:length(x1)/2), 'b', 'LineWidth', 1.5);
hold on;
% 绘制相位谱
plot(f(1:length(x1)/2), ph(1:length(x1)/2), 'r', 'LineWidth', 1.5);
hold off;

% 添加标签和标题
xlabel('频率 (Hz)');
ylabel('Amplitude');
title('Frequency Spectrum and Phase Spectrum');
legend('Amplitude', 'Phase');

% 初始相位通常指的是零频率处的相位，即第一个点的相位
initial_phase = ph(1); % 取第一点作为初始相位，注意这只是估计值，实际可能需根据具体算法确定

% 显示结果
disp(['Initial phase at DC (0 Hz): ', num2str(initial_phase)]);

这段代码首先读取音频文件，然后计算其傅立叶变换，提取幅度和频率信息。最后，它绘制出幅度谱和相位谱，并显示了直流分量（0 Hz）处的初始相位。