关于双谱切片的程序，错在哪？请指教！
本帖最后由 wdhd 于 2016-9-1 13:55 编辑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')
我做的上面的程序是不是那个地方做错了?为什么我从切片图中得不到信号的正确的频率?
[ 本帖最后由 zhlong 于 2007-6-4 21:51 编辑 ]
双谱
请大家帮帮忙啊！QQ：47629765，wangwlcm@sohu.com 请注意：不要在多个版内重复发贴！谢谢！
页:
[1]