ddy41 发表于 2012-9-17 17:33
一定要截的!镜像法是最基本的抑制方法,效果一般
clear all;
close all;
load('IR014_1.mat','y')
load('zhengchang.mat','z')
subplot(211);plot(z);xlabel('点数');ylabel('幅值');title('正常信号');
subplot(212);plot(y);xlabel('点数');ylabel('幅值');title('故障信号');
fs=12000;
N=2400;
n=0:N-1;
t=n/fs;
f=1./(2400/12000)*n;
c=emd(y);
[r,s]=size(c);
figure;
for q=1:r
subplot(r,1,q);plot(t,c(q,:));xlabel('时间/t');ylabel('幅值');
end
for q=1:r-1
a(q)=kurtosis(c(q,101:2300));
end
峭度是这样写吗?怎么总感觉运行不对呀,请指点
还有我不会截取端点,怎么截呀?
|