python复数转化为相位角,Scipy FFT-如何获得相位角

I'm having trouble getting the phase of a simple sine curve using the scipy fft module in python. I followed this tutorial closely and converted the matlab code to python. However, no matter what phase I use for the input, the graph always shows 3. What am I missing?

import numpy as np

import matplotlib.pyplot as plt

import scipy.fftpack

import cmath

A=10

fc = 10

phase=60

fs=32#Sampling frequency with oversampling factor of 32

t = np.arange(0,2,1/fs)

#Convert the phase shift to radians from degrees.

phi = phase*np.pi/180

x=A*np.cos(2*np.pi*fc*t+phi)

N=256

X = scipy.fftpack.fftshift(scipy.fftpack.fft(x,N))/N

df=fs/N #Frequency resolution.

sampleindex = np.arange(-N/2,N/2,1) #Ordered index for FFT plot.

f &