小虾米

在谈傅立叶变换时，大家关注的往往是其幅度谱，那么位相谱的作用呢？

查阅了一些资料，说位相谱表示为信号的个频率成分的初始位相。那么我在matlab做了一个计算，代码如下：
t = 0:0.001:0.6;
x = sin(2*pi*50*t)+sin(2*pi*120*t);
%y = x + 2*randn(size(t));
y=x;
figure;
plot(1000*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')

Y = fft(y,512);
Pyy = Y.* conj(Y) / 512;
phase=angle(Y)*180./pi;
f = 1000*(0:256)/512;
figure;
plot(f,Pyy(1:257));
figure;
plot(f,phase(1:257));
title('Frequency content of y');
xlabel('frequency (Hz)');

这是matlab帮助文件中的一个例子，我稍加改动，增加了计算计算傅立叶变换后计算位相谱的语句及相应的作图。

从理论上讲，信号由两个初位相为零的正弦函数组成，那么画出的位相谱和频率的曲线至少在频率为50和120处为零，但发现位相值不为零，而且曲线很不规则。这是为什么？


请各位大侠指教，小女子先谢过了。