谁能帮我提取一下这个信号的高频成分

16443

由于我能力有限，谁能帮我提取一下这个振动信号的高频成分。
该信号根据数值方法计算得到，不会存在噪声干扰，信号采样间隔Ts＝4E-4  
信号文件见附件

songzy41

可试试用小波分解的方法，只取高频部分。

16443

我对小波不太懂，弄了一晚上都没有结果 。
你能帮帮忙吗 ？

zhangnan3509

你的采样频率是多少，采样点是多少？这组数据是也就是几列，这些重要参数能不能写在帖子上，虽然这些都可以在处理中得到，但是毕竟不太方便

16443

采样点数N＝6250
采样频率Fs＝2500

zhangnan3509

应该出现在什么地方？在处理前你知道答案吗？

16443

高频出现在正半周期和负半周期过渡的短时间内。
这只是理论分析得到的

zhangnan3509

你做出来的结果是什么样啊！这个数据高频部分能量太小了，你想提取出来干什么？

[ 本帖最后由 zhangnan3509 于 2007-7-5 16:53 编辑 ]

wy558558558

打不开文件啊，晕！不过看你说ＦＳ才２５００，估计８００多是不是就算高频了啊？有２个方法首先，最好别用小波，用小波包更好一点，建议进行２到３层分解就可以．第二可以选择一个高频滤波器提取也可以．

songzy41

我这里用了6层分解，只取高频分量。分解层少的话，在低频分量中还包含有较多的高频成分，到分解6层时还稍满意。程序如下，其中用的是db3。笫1图中红色是原始信号，黑色是6层分解后的低频信号，把它们作一个比较；笫2图中便是高频信号。
load sig
N=length(sig);
T=4e-4;
nn=1:N;
[c,l]=wavedec(sig,6,'db3');  
y=wrcoef('a' ,c,l,'db3' ,6);
plot((nn-1)*T,sig,'r','linewidth',3); hold on;
plot((nn-1)*T,y,'k');
hold off; grid;
figure
z=wrcoef('d' ,c,l,'db3' ,6);
plot((nn-1)*T,z,'b'); grid;

zhangnan3509

N=6250;

f=2500;


s1=signal(zb,1)





%一维小波分解
[c,l]=wavedec(s1,6,'db2');
%重构第1-6层细节信号
d6=wrcoef('d',c,l,'db2',6);
d5=wrcoef('d',c,l,'db2',5);
d4=wrcoef('d',c,l,'db2',4);
d3=wrcoef('d',c,l,'db2',3);
d2=wrcoef('d',c,l,'db2',2);
d1=wrcoef('d',c,l,'db2',1);
figure(1)
subplot(6,1,1)
plot(d6,'linewidth',2);
ylabel('db6');
subplot(6,1,2)
plot(d5,'linewidth',2);
ylabel('db5');
subplot(6,1,3)
plot(d4,'linewidth',2);
ylabel('db4');
subplot(6,1,4)
plot(d3,'linewidth',2);
ylabel('db3');
subplot(6,1,5)
plot(d2,'linewidth',2);
ylabel('db2');
subplot(6,1,6)
plot(d1,'linewidth',2);
ylabel('db1');

[ 本帖最后由 zhangnan3509 于 2007-7-5 18:33 编辑 ]

zhangnan3509

上面就是一个简单的例子，不一定效果就好，你可以把它作为一个框架，很具自己具体需要修改就可以了

wy558558558

小波分析是无法得到准确高频信息的．你可以算一下你的高频的频率段是多少．小波包知道提取最后一层的所有信息就可以了．:@L

zhangnan3509

呵呵！能不能把小波包处理的效果图和程序传上来，比较一下啊。咱们版里好像小波包的实例不多呀，不吝赐教啊:handshake

16443

不好意思，下午有事情出去了，没能及时回复 。
我想看看高频随时间有怎样的变化 ，然后在根据变化情况做进一步的分析（机械结构结合部特性的研究）
我只用STFT分析出来结果，就是高频能量很低，看不出频率跟时间的关系。

[ 本帖最后由 zhlong 于 2007-7-5 22:04 编辑 ]