Y=fft(y,512);
F =10*f*[0:256]/512;
fp=2*sqrt(Y.*conj(Y));％幅度譜
xp=angle(Y); ％相位譜
gl=abs(Y).^2; %功率譜
magif=ifft2(abs(f2));%幅度重構
pha=angle(f2);%取相位
phaif=ifft2(exp(j*pha));%相位重構
數據是x(i)，***N個點，采樣頻率是fsample
dt = 1/fsample
df = 1/N/dt
f = 0:df:(N-1)*df
X=fft(fftshift(x))
plot(f,real(X))
plot(f,imag(X))
擴展資料：
相位譜是調整聲音相位的，最容易理解的就是左右聲道的位置調整，實際上相位還決定著其他很多聲音的屬性。
對於壹個系統，能夠通過其相位譜來判斷該系統是否為線性相位系統。線性相位系統故名思義，看相位是否隨頻率線性變化。但相位譜的作用不僅限於此，奧本海姆在壹篇經典文獻中認為信號的相位包含的信息大於幅度，實際上從最初的最小相位系統，倒譜分析，到系統辨識，高階譜估計等理論都是以相位譜為突破口。
百度百科-相位譜