帮忙看看我这个相位-频率图有什么问题吗？

lzh19840809

clc
clear
close all
SampleData=load('3.txt')
x=SampleData(1:51200)
x=x-mean(x)
subplot(311)
plot(x)
title('51200个样本点，fs=100')
grid on
%---------------------------------------
%FFT变换
%---------------------------------------
fs=100
N=51200
n=0:N-1;
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=n*fs/length(y);%进行对应的频率转换
subplot(312);
plot(f(1:N/2),mag(1:N/2));%做频谱图
title('幅值-频率曲线图')
grid on
Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(y(i)); %计算相位
Pyy(i)=Pyy(i)*180/pi; %换算为角度
end;
subplot(313);
plot(f(1:N/2),Pyy(1:N/2));    %显示相位图
title('相位-频率曲线图');
grid on

tisang

第一个图横坐标不是频率，第二个图一般表示成对数，第三个图更看不出什什么

lzh19840809

第一个横坐标是51200个样本点。
第三个图我怀疑有问题，重点是第三个的问题

gotlanbo

同问 ，第3个图有什么问题 ，应该怎样改 ？？

hcharlie

FFT只是一个数学方法，不要忘了你想解决什么问题。
从图1看，你采集的大概是随机信号。随机信号只有“谱估计”没有“谱计算”，所以图2只是一个谱估计值，而从图2来看，看出来幅值的一些趋势，但这些值的精度和置信度很低，也就是不可信的。
LZ不知道什么是随机信号的功率谱统计平均？
正确的做法是可以将这51200个数据分成比如50组1024点数据，每一组要加（海宁）窗，FFT，再做这50组数据的（功率）谱平均，再画出谱的幅值-频率图，置信度就高多了。
随机信号没有相位问题。只有周期信号才有相位问题。所以图3是没有物理意义的。

[ 本帖最后由 hcharlie 于 2009-11-4 07:53 编辑 ]

gotlanbo

楼上回答的太经典了。可是，如果那个不是完全的随机信号，而是周期信号和噪声的组合呢，也不能用图3的相位图表示出来么？
还想问一下，如果是周期信号，如何用MATLAB的程序实现相位求法呢？

[ 本帖最后由 gotlanbo 于 2009-11-4 09:04 编辑 ]

hcharlie

如果是周期信号和噪声的组合，先用随机法处理，在周期信号的频率附近有频谱的尖峰。
为什么对相位如此关心。相位与采集起始点和另一个参考信号是相关联而不是绝对独立的。