有用matlab做双谱图的高手？

xlyloveyou

作双谱对角切片，请各位多多帮忙　谢谢！！！！

[ 本帖最后由 zhangnan3509 于 2008-12-2 22:04 编辑 ]

simon21

有什么具体问题可以贴出来大家讨论

wangwlcm

n=1024;
%采样频率
fs=2000;
%输入信号
[DATAfile DATApath]=uigetfile('*.txt','输入信号');
FILENAME=[DATApath,DATAfile];
DATA=load(FILENAME)*0.000305+0.000153;
y=(DATA([10000:(n+9999)],2))
y1=mean(y);
y2=y-y1;
%双谱估计平面图
figure(1)
[bspec,waxis]=bispecd(y2,16)
title('双谱估计平面图')
%切片图，
figure(2)
t=50;
for k=-t:t
y_cum = cumest(y2,3,t,128,50,'biased',k,k)
end
%%取切片的正半部分
y3=y_cum([51:101],:);
plot(y3);
%时间单位
k2=(0:t)./fs;
figure(3)
plot(k2,y3)
%求切片图的fft
nfft=128;
yy=abs(psd(y3,nfft));
figure(4)
plot((0:nfft/2-1)/nfft*fs,yy(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('fft')

我做的上面的程序是不是那个地方做错了？为什么我从切片图中得不到信号的正确的频率？

[ 本帖最后由 VibInfo 于 2006-8-8 07:04 编辑 ]

aprillq

我也正好需要呢:@L

aprillq

cumest不是求累积量么？怎么成切片了？

tianyoume

我在这个网站上查到过双谱的程序, 你查查看有帮助没有:http://www.mathworks.com/matlabc ... nge/loadCategory.do

lovegirlmeijie

我也刚刚学习双谱，正在研究对角切片，不知道楼上是否已经解决这个问题？另外可不可以用matlab的高阶谱工具箱中求三阶累积量的函数cumest函数？？

真人

Warning: Integer operands are required for colon operator when used as index.
> In CUM3EST at 60
  In CUMEST at 71
  In readwaiquan1 at 27

cumest这个函数你用错了吧？

wyq-q

:loveliness: :loveliness: :loveliness: :loveliness: :loveliness:

xiaocai315

哥哥，有没有高阶谱的工具箱分享下，或则求三阶累积量，双谱的函数也行啊，在此谢过了！caiyihongts@126.com

200810100659

非常感谢大家的分享