|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clc
clear all
load 72HrStress_102_1
s=chanvals(:,1)';
s=(s-mean(s))/std(s,1);
fs=102400;
[c,l]=wavedec(s,6,'db10');
for i=1:6
d=wrcoef('d',c,l,'db10',i);
y=abs(hilbert(d));
Sper=abs(fft(y));
Sper=Sper*2/length(Sper);
fn=(0:length(Sper)-1)*fs/length(Sper);
figure(1),
subplot(6,1,i),
plot(fn(1:length(fn)/2),Sper(1:length(Sper)/2))
axis([0,1200,0,0.1]);
end |
|