xlyloveyou 发表于 2006-5-21 10:49

有用matlab做双谱图的高手？

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

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

simon21 发表于 2006-6-7 15:37

回复：（xlyloveyou）有用matlab做双谱图的高手？

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

wangwlcm 发表于 2006-6-9 16:53

双谱切片分析

n=1024;
%采样频率
fs=2000;
%输入信号
=uigetfile('*.txt','输入信号');
FILENAME=;
DATA=load(FILENAME)*0.000305+0.000153;
y=(DATA(,2))
y1=mean(y);
y2=y-y1;
%双谱估计平面图
figure(1)
=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(,:);
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 发表于 2007-11-5 17:18

求助

我也正好需要呢:@L

aprillq 发表于 2008-3-13 17:57

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

tianyoume 发表于 2008-3-14 11:43

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

lovegirlmeijie 发表于 2008-12-2 10:54

回复 板凳 wangwlcm 的帖子

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

真人 发表于 2008-12-21 17:17

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 发表于 2009-2-13 15:40

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

xiaocai315 发表于 2010-7-7 15:38

回复 板凳 wangwlcm 的帖子

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

200810100659 发表于 2010-12-6 22:58

{:{27}:}非常感谢大家的分享

402144999 发表于 2014-7-28 17:49

页: [1]
查看完整版本: 有用matlab做双谱图的高手？