我自己编的谐波小波程序（供大家参考）

hhf_online · 发表于 2009-9-18 23:03

程序好像有问题，fsai应该限定范围，另外，重构信号的幅值也不对。
期待楼主最新的成果。

jinyan19860616 · 发表于 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

yonglets · 发表于 2011-3-8 11:38

先看看！！！

wildcranehust · 发表于 2011-3-8 19:51

楼主牛人啊，谢谢分享，

shoupi · 发表于 2012-3-6 23:00

谢谢楼主分享啊！！

shoupi · 发表于 2012-3-6 23:10

学习了，谢谢！

lragy · 发表于 2012-3-26 00:09

fs1=dt*fft(s)./(2*pi);
dt/2*pi为什么要有？？？？？？？？

kangyan139 · 发表于 2012-7-17 16:32

谢谢分享。不知楼主还在否？程序里面的Wd是什么意思啊？重构之后为啥分为实部和虚部？

kangyan139 · 发表于 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))

忍爱天下 · 发表于 2014-5-28 16:13

为什么研究理论的比较多而对于实际信号处理的反而比较少呢

cloudy · 发表于 2014-8-21 12:04

借鉴借鉴

danzhilan · 发表于 2014-12-18 11:07

谢谢楼主啊

账号		自动登录	找回密码
密码			我要加入

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