谁能帮我提取一下这个信号的高频成分
由于我能力有限，谁能帮我提取一下这个振动信号的高频成分。该信号根据数值方法计算得到，不会存在噪声干扰，信号采样间隔Ts＝4E-4
信号文件见附件 可试试用小波分解的方法，只取高频部分。
回复 #2 songzy41 的帖子
我对小波不太懂，弄了一晚上都没有结果 。你能帮帮忙吗 ？
回复 #3 16443 的帖子
你的采样频率是多少，采样点是多少？这组数据是也就是几列，这些重要参数能不能写在帖子上，虽然这些都可以在处理中得到，但是毕竟不太方便回复 #4 zhangnan3509 的帖子
采样点数N＝6250采样频率Fs＝2500
回复 #5 16443 的帖子
应该出现在什么地方？在处理前你知道答案吗？回复 #6 zhangnan3509 的帖子
高频出现在正半周期和负半周期过渡的短时间内。这只是理论分析得到的
回复 #7 16443 的帖子
你做出来的结果是什么样啊！这个数据高频部分能量太小了，你想提取出来干什么？[ 本帖最后由 zhangnan3509 于 2007-7-5 16:53 编辑 ] 打不开文件啊，晕！不过看你说ＦＳ才２５００，估计８００多是不是就算高频了啊？有２个方法首先，最好别用小波，用小波包更好一点，建议进行２到３层分解就可以．第二可以选择一个高频滤波器提取也可以． 我这里用了6层分解，只取高频分量。分解层少的话，在低频分量中还包含有较多的高频成分，到分解6层时还稍满意。程序如下，其中用的是db3。笫1图中红色是原始信号，黑色是6层分解后的低频信号，把它们作一个比较；笫2图中便是高频信号。
load sig
N=length(sig);
T=4e-4;
nn=1:N;
=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; N=6250;
f=2500;
s1=signal(zb,1)
%一维小波分解
=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 编辑 ] 上面就是一个简单的例子，不一定效果就好，你可以把它作为一个框架，很具自己具体需要修改就可以了
回复 #10 songzy41 的帖子
小波分析是无法得到准确高频信息的．你可以算一下你的高频的频率段是多少．小波包知道提取最后一层的所有信息就可以了．:@L回复 #13 wy558558558 的帖子
呵呵！能不能把小波包处理的效果图和程序传上来，比较一下啊。咱们版里好像小波包的实例不多呀，不吝赐教啊:handshake回复 #8 zhangnan3509 的帖子
不好意思，下午有事情出去了，没能及时回复 。我想看看高频随时间有怎样的变化 ，然后在根据变化情况做进一步的分析（机械结构结合部特性的研究）
我只用STFT分析出来结果，就是高频能量很低，看不出频率跟时间的关系。
[ 本帖最后由 zhlong 于 2007-7-5 22:04 编辑 ]