matlab已知幅度怎么求相位,用FFT求已知波形的幅值和相角

运行以下的程序得到的结果如下：

频率               幅值         初始相角(弧度)  初始相角(度)

59.996548      6.999343      -1.569867      -89.946751

300.004871   0.999988       -1.640208      -93.976971

780.018408   0.300000         2.078221     119.073274

xx=load('test_signal-21.txt');

time=xx(:,1);

dt=time(2)-time(1);

fs=1/dt;

x=xx(:,2);

subplot 211; plot(time,x,'r','linewidth',3);

N=length(time);

rad=pi/180;

X=fft(x);

n=1:(N+1)/2;

freq=(n-1)*fs/N;

subplot 212; plot(freq,20*log10(abs(X(n)))); grid;

Z1=Specorr(x,fs,N,40,70,2);

Z2=Specorr(x,fs,N,280,320,2);

Z3=Specorr(x,fs,N,760,800,2);

fprintf('%5.6f   %5.6f   %5.6f   %5.6f\n',Z1,Z1(3)/rad);

fprintf('%5.6f   %5.6f   %5.6f   %5.6f\n',Z2,Z2(3)/rad