📄 matlab_wavelet.txt
字号:
第 18 楼发表于 2006-4-26 16:18 资料 文集 短消息
请问楼主有对称小波的程序吗?谢谢了
[em04]
q2384
应门书童
精华: 0
积分: 7
帖子: 12
威望: 6 点
振动币: 2698 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2005-11-23
第 19 楼发表于 2006-4-26 19:17 资料 文集 短消息
谢谢
zhawenhua
应门书童
精华: 0
积分: 2
帖子: 6
威望: 1 点
振动币: 360 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-4-26
第 20 楼发表于 2006-4-26 21:07 资料 文集 短消息
辛苦了
谢谢!
ruola
应门书童
精华: 0
积分: 1
帖子: 2
威望: 1 点
振动币: 163 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-4-30
第 21 楼发表于 2006-4-30 10:50 资料 文集 短消息
请问有没有小波系数模极大值的程序啊!
arlay
布衣书生
精华: 0
积分: 3
帖子: 24
威望: 1 点
振动币: 711 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-3-13
第 22 楼发表于 2006-4-30 14:34 资料 文集 短消息
谢谢分享!^_^
simon21
新科状元
精华: 6
积分: 400
帖子: 670
威望: 321 点
振动币: 23552 个
管理积分: 0 点
人气指数: 0
阅读权限: 80
注册: 2005-7-24
第 23 楼发表于 2006-5-23 14:44 资料 个人空间 短消息
回复:(simon21)[分享]个人收集的一些关于小波分析的...
消失矩作用的程序
clear;clc;
f=50;
T=0.001;
n=1:50;
y=sin(2*pi*f*n*T);
noise=[zeros(1,17),0.2*randn(1,15),zeros(1,18)];
y_noise=y+noise;
figure(1);
subplot(2,1,1);
plot(y);
title('signal');
subplot(2,1,2);
plot(y_noise);
title('noise & signal');
[yl2,yh2]=dwt(y,'db2');
[yl10,yh10]=dwt(y,'db10');
figure(2);
subplot(2,1,1);
plot(yl2);
title('db2 low frequency signal');
subplot(2,1,2);
plot(yh2);
title('db2 high frequency signal');
figure(3);
subplot(2,1,1);
plot(yl10);
title('db10 low frequency signal');
subplot(2,1,2);
plot(yh10);
title('db10 high frequency signal');
[yl2,yh2]=dwt(y_noise,'db2');
[yl10,yh10]=dwt(y_noise,'db10');
figure(4);
subplot(2,1,1);
plot(yl2);
title('db2 low frequency noise & signal');
subplot(2,1,2);
plot(yh2);
title('db2 high frequency noise & signal');
figure(5);
subplot(2,1,1);
plot(yl10);
title('db10 low frequency noise & signal');
subplot(2,1,2);
plot(yh10);
title('db10 high frequency noise & signal');
simon21
新科状元
精华: 6
积分: 400
帖子: 670
威望: 321 点
振动币: 23552 个
管理积分: 0 点
人气指数: 0
阅读权限: 80
注册: 2005-7-24
第 24 楼发表于 2006-5-23 14:45 资料 个人空间 短消息
回复:(simon21)[分享]个人收集的一些关于小波分析的...
平移变换平移法(cycle_spinning)消除gibbs效应
clear;
clc; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 1.原始信号 f=50; % 信号频率
fs=800; % 采样频率
N=128; % 采样点 % 信号赋值
n=1:N;
y=sin(2*pi*f*n/fs); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 2.噪声
noise=0.4*rand(1,128); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 3.染噪信号
y_noise=y+noise; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 4.硬消噪采用cycle_spinning技术 % 累加量
z5=zeros(1,N); % 平移变换频移法
for i=1:N;
z=circshift(y_noise.',i-1).'; % 源信号右平移
[z1,z2]=lwt(z,'db3'); % 小波正变换
z2=zeros(1,N/2); % 高频分量全部为零(主要噪声,硬消噪)
z3=ilwt(z1,z2,'db3'); % 小波反变换
z4=circshift(z3.',-(i-1)).'; % 变换后信号左平移
z5=z5+z4/N; % 平均
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 5.显示 error=norm(y-z5)/norm(y); % 相对误差 figure(1);
subplot(2,1,1)
plot(y,'r');
legend('源信号');
subplot(2,1,2);
plot(y_noise);
legend('染噪信号');
figure(2);
subplot(2,1,1)
plot(y,'r');
legend('源信号');
title(error);
subplot(2,1,2);
plot(z5);
legend('消噪后信号');
yushouzhuo
应门书童
精华: 0
积分: 2
帖子: 10
威望: 1 点
振动币: 654 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2005-11-15
第 25 楼发表于 2006-5-23 16:39 资料 文集 短消息
哇!真是好东西啊
cinda
应门书童
精华: 0
积分: 1
帖子: 1
威望: 1 点
振动币: 120 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-5-23
第 26 楼发表于 2006-5-24 16:08 资料 文集 短消息
请教一个问题,关于三楼那个程序(2代小波示意程序 )
2代小波示意程序 (三楼)中
……
f1(:,T+1)=f1(:,1); % 补列
f2(T+1,:)=f2(1,:); % 补行
……
补行是把分裂的偶数的第一行补到偶数的最后,但是补列
是怎么进行的?比如图像是256×256,分裂后奇数为128
×256,偶数也为128×256,
进行f1(:,T+1)=f1(:,1); 补列之后,相当于把原来的第129列
替换成第1列
这样补行和补列后,奇数和偶数的分裂结果为128×256和
129×256
还请大侠不吝赐教,小生在此谢过了
[em24]
[此贴子已经被作者于2006-5-24 16:09:37编辑过]
crosstree
应门书童
精华: 0
积分: 1
帖子: 2
威望: 1 点
振动币: 136 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-6-10
第 27 楼发表于 2006-6-10 13:17 资料 文集 短消息
多谢楼主,都是好东西~
灰色轨迹
应门书童
精华: 0
积分: 2
帖子: 11
威望: 1 点
振动币: 152 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-6-10
第 28 楼发表于 2006-6-11 08:34 资料 个人空间 短消息
经典,谢谢
wj118
应门书童
精华: 0
积分: 2
帖子: 7
威望: 1 点
振动币: 359 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-6-8
第 29 楼发表于 2006-6-13 19:13 资料 文集 短消息
thanks
LZ是个大好人,谢谢了.收藏.
cike5
布衣书生
精华: 0
积分: 4
帖子: 26
威望: 1 点
振动币: 528 个
管理积分: 0 点
人气指数: 0
阅读权限: 10
注册: 2006-8-13
第 30 楼发表于 2006-8-14 19:55 资料 文集 短消息
真的不错
其实我一直在努力,但是水平有限,还是先向各位学习,谢谢大家!
53 2/4 ??1234?? 投票
交易
悬赏
活动
振动论坛
家园建设
> 论坛公告
> 论坛建设
> 版面建设、版主申请
> 专家讲坛
> 专家讲坛资料
> 新手上路
> 新人送积分活动
振动基础
> 基础力学理论
> 一般力学和振动理论
> 数学与非线性科学
> 算法原理及编程
> 结构可靠性及优化设计
> 控制理论及应用
> NVH及振动利用工程
CFD专区
> CFD基础理论与实践
> Fluent讨论区
> 其他CFD软件讨论区
信号处理
> 振动实验与测试技术
> 信号处理方法
> 故障诊断及健康监测
噪声分析
> 声学基础及噪声治理
> 声学分析软件
> 声学测量
数学软件
> Matlab讨论区
> Simulink仿真专区
> 虚拟仪器技术
> 其他数学软件
工程软件
> Ansys讨论区
> MSC系列软件专区
> ABAQUS讨论区
> CAD软件讨论区
> 其他软件及相关讨论
综合讨论
> 电脑安全维护
> 科学研究方法
> 信息广场
休闲娱乐
> 心情故事
> 图行天下
管理专区
当前时区 GMT+8, 现在时间是 2007-3-31 16:04
黑ICP备06005726号
Powered by Discuz! 5.0.0 ? 2001-2006 Comsenz Inc.
Processed in 0.100503 second(s), 9 queries , Gzip enabled TOP
清除 Cookies - 联系我们 - 中国振动联盟 - Archiver - WAP
社区银行
论坛勋章
虚拟股市
宠物中心
宣传中心
人气排行
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -