|
|
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
本帖最后由 wdhd 于 2016-9-1 13:55 编辑
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')
我做的上面的程序是不是那个地方做错了?为什么我从切片图中得不到信号的正确的频率?
[ 本帖最后由 zhlong 于 2007-6-4 21:51 编辑 ]
|
|
发表于 2006-6-9 17:00
提升卡
变色卡
显身卡
楼主
发表于 2006-6-9 18:59
