Matlab中信号经FFT变换后，如何求峰值点处的相位？

百度上的问题，看到后做了个实验如下：https://zhidao.baidu.com/question/511735134.html

用测试麦克风（UMIK-1)播放测试音频test.wav。先把喇叭的正负极接正确，然后在WIN10里用系统自带的录音机程序录制test1.m4a。再把喇叭的正负接反，录制test2.m4a。有了这两个录制文件后，可以用matlab观察喇叭接反后的变化。

修改后的代码如下：

%test1.m4a和test2.m4a是用同一音频文件播放时录制的两个文件，其中一个文件在录制时反接喇叭的正负极。

%分别得到最大幅值时的相位，然后判断相位差是否180度，即可判断是否反相。

clear;

clc;

[y_test1,fs]=audioread('d:\test1.m4a'); %将声音放于matlab中

N_test1=length(y_test1);

X_test1=fft(y_test1);

%DFT后的幅度，此为复数，包含幅度值和相位值

X_test1=X_test1/N_test1;

[y_test2,fs_test2]=audioread('d:\test2.m4a'); %将声音放于matlab中

N_test2=length(y_test2);

X_test2=fft(y_test2);

%DFT后的幅度，此为复数，包含幅度值和相位值

X_test2=X_test2/N_test2;

%归一化频率轴

f=(0:N_test1-1)/N_test1;

H=10*log10(abs(X_test1));%幅度，单位dB

P=angle(X_test1); %相位

figure();

subplot(4,1,1)

plot(f,H);

title('幅度/频率特性曲线');

subplot(4,1,2)

plot(f,P);

title('相位/频率特性曲线');

%画出test1峰值点出处的相位

X_max_test1=max(abs(X_test1));%找出最大峰值的数值

%找出最大值所在序列中的下标

flag_test1=find(abs(X_test1)==X_max_test1);

h_test1=angle(X_test1(flag_test1));%按着下标求出对应的相位

subplot(4,1,3)

plot(flag_test1,h_test1);

title('test1.m4a峰值处的相位')

%画出test2峰值点出处的相位

X_max_test2=max(abs(X_test2));%找出最大峰值的数值

%找出最大值所在序列中的下标

flag_test2=find(abs(X_test2)==X_max_test2);

h_test2=angle(X_test2(flag_test2));%按着下标求出对应的相位

subplot(4,1,4)

plot(flag_test2,h_test2);

title('test2.m4a峰值处的相位')

grid on

附件

结果如图：