ganhaohao
2009-03-16, 20:29
我有一组时间序列：
time=【12.000000000000 15.142982850836 18.285965701672 21.428948552508 24.571931403345 27.714914254181 30.857897105017 34.000879955853 37.143862806689 40.286845657525 43.429828508361 】
data=【16.0000000000000 21.2415061993540 34.6423808050926 28.5226576610713
36.0406016783245 29.8345124690840 28.7401134729748 27.5083562038488
26.8925311469514 24.6970909632913 22.9791261091705】
进行最大熵谱分析如下：
time=load('time.txt'); %时间序列
data=load('data.txt'); %原始信号
n=length(time);
dt=(time(end)-time(1))/(n-1);
Fs=1/dt; %采样频率
N=512;Nfft=256; %数据长度、分段长度
%采用最大熵MEM法，采用滤波器阶数15，估计功率谱
[Pxx1,f]=pmem(data,5,Nfft,Fs);
plot(f,10*log10(Pxx1)); %绘制功率谱
请问如何在生成的功率谱图中加一条显著性检验曲线，以直观的显示哪些频率是可信的？
求问，很着急，希望能得到帮助。
time=【12.000000000000 15.142982850836 18.285965701672 21.428948552508 24.571931403345 27.714914254181 30.857897105017 34.000879955853 37.143862806689 40.286845657525 43.429828508361 】
data=【16.0000000000000 21.2415061993540 34.6423808050926 28.5226576610713
36.0406016783245 29.8345124690840 28.7401134729748 27.5083562038488
26.8925311469514 24.6970909632913 22.9791261091705】
进行最大熵谱分析如下：
time=load('time.txt'); %时间序列
data=load('data.txt'); %原始信号
n=length(time);
dt=(time(end)-time(1))/(n-1);
Fs=1/dt; %采样频率
N=512;Nfft=256; %数据长度、分段长度
%采用最大熵MEM法，采用滤波器阶数15，估计功率谱
[Pxx1,f]=pmem(data,5,Nfft,Fs);
plot(f,10*log10(Pxx1)); %绘制功率谱
请问如何在生成的功率谱图中加一条显著性检验曲线，以直观的显示哪些频率是可信的？
求问，很着急，希望能得到帮助。