声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: kevin19821

[小波] 我自己编的谐波小波程序(供大家参考)

[复制链接]
发表于 2009-9-18 23:03 | 显示全部楼层
程序好像有问题,fsai应该限定范围,另外,重构信号的幅值也不对。
期待楼主最新的成果。
回复 支持 反对
分享到:

使用道具 举报

发表于 2011-1-10 11:31 | 显示全部楼层
Newland D E. Harmonic wavelets in vibrations andacoustics 和 Newland D E. Wavelet analysis of vibration - part 1:theory 、Newland D E.Harmonic wavelet analysis这三篇文章哪位兄弟姐妹有啊,借我看看,本人不胜感激、本人邮箱:jinyan5407101@sina.com
发表于 2011-3-8 11:38 | 显示全部楼层
先看看!!!
发表于 2011-3-8 19:51 | 显示全部楼层
楼主牛人啊,谢谢分享,
发表于 2012-3-6 23:00 | 显示全部楼层
谢谢楼主分享啊!!
发表于 2012-3-6 23:10 | 显示全部楼层
学习了,谢谢!
发表于 2012-3-26 00:09 | 显示全部楼层
fs1=dt*fft(s)./(2*pi);
dt/2*pi为什么要有????????
发表于 2012-7-17 16:32 | 显示全部楼层
谢谢分享。不知楼主还在否?程序里面的Wd是什么意思啊?重构之后为啥分为实部和虚部?
发表于 2012-7-17 16:36 | 显示全部楼层
怎么看分了多少层啊?求指教,懂得请指教`~
发表于 2013-6-24 10:01 | 显示全部楼层
fs1=dt*fft(s)./(2*pi);是什么意思呀
发表于 2014-3-31 11:09 | 显示全部楼层
clc
clear all
fs=4096;
dt=1/fs;
N=2048;
f=fs*(0:N-1)./N;
t=0:dt:2047*dt;
Wd=hann(2048);
y1=3*sin(2*pi*50*t)+2*sin(2*pi*100*t);
y2=0.05*sin(2*pi*150*t);
s=y1+y2;
subplot(411)
plot(t,y2);
fs1=dt*fft(s)./(2*pi);
%j=7,s=127
figure(1)
subplot(412);
plot(f-fs/2,abs([fs1((N/2+1):N) fs1(1:N/2)]))
m=148;
n=152;
fsai=(exp(2*n*pi*t*i)-exp(2*pi*m*t*i))./(2*pi*(n-m)*t*i+1e-100).*Wd';
f1=dt*fft(fsai)./(2*pi);
subplot(413)
f=fs*(0:N-1)./N;
plot(f-fs/2,abs([f1((N/2+1):N) f1(1:N/2)]));
N=length(s);
W=f1.*fs1;
subplot(414)
plot(f-fs/2,abs([W((N/2+1):N) W(1:N/2)]));
w1=4*pi*2*pi*ifft(W)./dt;
figure(2)
subplot(211)
plot(t,real(w1))
subplot(212)
plot(t,imag(w1))

2.jpg
3.jpg
发表于 2014-5-28 16:13 | 显示全部楼层
为什么研究理论的比较多 而对于实际信号处理的反而比较少呢
发表于 2014-8-21 12:04 | 显示全部楼层
借鉴借鉴
发表于 2014-12-18 11:07 | 显示全部楼层
谢谢楼主啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 20:02 , Processed in 0.062910 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表