fft之后求模值和相位_Matlab FFT后求解相位差，结果与原论文不对（附程序）

原信号：sin(100*pi*t+初相角)是标准的50HZ的正弦信号，初相角从0度变化到14度

对该信号线性插值得到新序列：(t-t1)/(t0-t1)*sin(100*pi*t0+(q-1)/180*pi)+(t-t0)/(t1-t0)*sin(100*pi*t1+(q-1)/180*pi)

FFT后求解原序列与插值后新序列，对应50HZ信号的相位差，结果与原论文不对。请高人给予指点。

clc;

clear;

x1(:,1)=[0;0;0;0;0;  0;0;0;0;0;  0;0;0;0;0];

x2(:,1)=[0;0;0;0;0;  0;0;0;0;0;  0;0;0;0;0];

y1(1,:)=[0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0];

y2(1,:)=[0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0];

y1_max=[0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0]';

y2_max=[0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0;  0;0;0;0;0;0;0;0]';

N=24;%24个点，23个区间，n本来用N表示，为了好理解不用了，只是表示是24个点采样

fs=24/0.02;

N=24;

format long;

for q=1:15%(初相角从0度变换到14度)

for  m=1:15%(每两个采样点之间的区间15等份)

<