交叉小波分析matlab,[转载]Matlab 小波分析及物理含义
一、为什么要进行小波变换
最初的原因很简单,傅里叶变换没有时间信息,也就是说,我们不知道傅里叶变换结果频率出现在什么时间,所以出现了时频分析,就是在一张结果图上同时表明信号的频率及其出现的时间,
这样做更有利于瞬时信号的研究。小波分析就是时频分析的一种。
二、Matlab连续小波变换实例
% cwt_test
clear all
close all
clc
%%
Ts = 0.001;
Fs = 1/Ts;
f1 = 20;
f2 = 50;
f3 = 100;
dt = 0.2;
t1 = (0:Ts:dt-Ts) + 0;
t2 = (0:Ts:dt-Ts) + dt;
t3 = (0:Ts:dt-Ts) + 2*dt;
y1 = sin(2*pi*f1*t1);
y2 = sin(2*pi*f2*t2);
y3 = sin(2*pi*f3*t3);
t = [t1 t2 t3];
y = [y1 y2 y3];
figure
plot(t,y)
xlim([t(1) t(end)])
ylim([min(y) max(y)])
xlabel('时间t')
ylabel('信号y(t)')
title('原始信号')
%%
scale = 1:50;
cw2 = cwt(y,scale,'morl');
figure
subplot(1,3,[2,3]) % 频率轴化为频率
[X,Y] = meshgrid(t,5/(2*pi)./scale*Fs);
mesh(X,Y,abs(cw2))
view(0,90)
title('时频图')
xlabel('时间')
ylabel('频率')
xlim([t(1) t(end)])
set(gca,'ylim',[0,max(max(Y))])
set(gca,'YScale','log')
set(gca,'YTick',[1:9,10:10:90,100:100:900,1000,2000])
三、物理含义
上图是用mesh函数绘制的三维小波变换图。
1. 横轴x是时间T
2.
纵轴y是小波的缩放scale,通过转换成为了频率,这个应该是跟选取的小波有关,这里的转换公式如代码所示,小波工具箱提供了求解的函数,具体见http://blog.sina.com.cn/s/blog_84024a4a01019ms1.html
3. Z轴是小波变换系数,代表了与小波与原信号之间的相关程度,越大表示相关度越高
4.
从图中可以看出,0~0.2s,主要频率是20Hz左右,0.2~0.4s,主要频率是50Hz左右,0.4~0.6s,主要频率在100Hz左右。跟原信号比较可知,结果还是很准确的。
上面的详细原理可以参考第五点。
四、小波变换部分更为规范的写法
wavename='morl';
totalscal= length(t);
Fc=centfrq(wavename); % 小波的中心频率
c=2*Fc*totalscal;
scals=c./(1:totalscal);
f=scal2frq(scals,wavename,1/Fs); % 将尺度转换为频率
cw2=cwt(y,scals,wavename); % 求连续小波系数
五、小波变换原理简单介绍
摘自《小波分析完美教程经典》
信号分析一般是为了获得时间和频率域之间的相互关系。傅里叶变换提供了有关频率域的信息,但时间方面的局部化信息却基本丢失。与傅里叶变换不同,小波变换通过平移母小波(mother
wavelet)可获得信号的时间信息,而通过缩放小波的宽度(或者叫做尺度)可以获得信号的频率特性。对母小波的缩放和平移操作是为了计算小波的系数,这些系数代表小波和局部信号之间的相互关系。本节将介绍小波分析中常用的三个基本概念:连续小波变换、离散小波变换和小波重构。
1. 连续小波变换
傅里叶分析是把一个信号分解成各种不同频率的正弦波,因此正弦波是其基函数。同样,小波分析是把一个信号分解成将原始小波经过移位和缩放之后的一系列小波,因此小波同样可以用作表示一些函数的基函数。可以说,凡是能够用傅里叶分析的函数都可以用小波分析,因此小波变换也可以理解为用经过缩放和平移的一系列函数代替傅里叶变换的正弦波。
数学上,傅里叶表示如下
同样,连续小波变换(continuous wavelet transform, CWT)用小式表示
这个式子的含义就是,小波变换是信号f(t)与被缩放和平移的小波函数(Psi)之积在信号存在的整个周期间求和。CWT变换的结果是许多小波系数C,这些系数是缩放因子(scale)和位置因子(position)的函数。
CWT的变换过程可分为如下5个步骤:
步骤1:把小波(Psi)和原始信号f(t)的开始部分进行比较
步骤2:计算系数C。该系数表示该部分信号与小波的近似程度。系数C的值越高表示信号与小波越相似,因此系数C可以反映这种波形的相关程度。
步骤3:把小波向右移,距离为k,得到的小波函数为(Psi(t -
k)),然后重复步骤1和2。再把小波向右移,得到小波(Psi(t -
2k)),重复步骤1和2,按照上述步骤一直进行下去,知道信号f(t)结束。
步骤4:扩展小波(Psi(t)),例如开展一倍,得到的小波函数为(Psi(t/2))。
步骤5:重复步骤1~4。
CWT的整个变换过程如图所示
小波变换完成之后得到的系数是在不同的缩放因子下由信号的不同部分产生的。这些小波系数、缩放因子和时间之间的关系和它们的含义可以用下图表示,该图是用Matlab软件绘制的。图(a)是用二维图像表示的小波分析图,x轴表示沿信号的时间方向上的位置,y轴表示缩放因子,每个x-y点的颜色表示小波系数C的幅度大小。图(b)是用三维图像表示的小波分析图,z轴表示小波变换之后的系数。
小波的缩放因子与信号频率之间的关系可以这样来理解。缩放因子小,表示小波比较窄,度量的是信号细节,表示频率(omega)较高;相反,缩放因子大,表示小波比较宽,度量的是信号的粗糙程度,表示频率(omega)比较低。
交叉小波分析matlab,[转载]Matlab 小波分析及物理含义相关推荐
- parrallel for matlab,[转载]MATLAB并行计算工具箱 -- Paralleln
龙泉居士基于文档原创,转载请注明 Parallel Computing Toolbox是一个matlab2011开始提供的组件,用于提供交互式的并行计算功能 一.运用的场合 很多应用程序中包含多个重复 ...
- 天空之城 matlab,[转载]matlab演奏《天空之城》代码
%matlab演奏<天空之城>代码 l_dao=262; %将"l_dao"宏定义为低音"1"的频率262Hz l_re =286; %将" ...
- mul matlab,[转载]Matlab boxplot for Mul
由于 命令 格式如下 :产生矩阵X的每一列的盒图和"须"图,"须"是从盒的尾部延伸出来,并表示盒外数据长度的线,如果"须"的外面没有数据,则 ...
- 求噪声功率 MATLAB,[转载]matlab中噪声功率、噪声方差关系
以matlab中awgn函数为例说明: 在matlab中无论是wgn还是awgn函数,实质都是由randn函数产生的噪声.即:wgn函数中调用了randn函数,而awgn函数中调用了wgn函数. 根据 ...
- 混沌pso matlab,[转载]MATLAB-PSO工具箱
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%% PSO工具箱函数 %%%%%%%%%%%%%%% pso PSO主程序 pso ...
- 小波变换的MATLAB原理,[转载]Matlab 小波分析及物理含义
一.为什么要进行小波变换 最初的原因很简单,傅里叶变换没有时间信息,也就是说,我们不知道傅里叶变换结果频率出现在什么时间,所以出现了时频分析,就是在一张结果图上同时表明信号的频率及其出现的时间, 这样 ...
- matlab作动态函数曲线图,[转载]Matlab中使用Plot函数动态画图方法总结
本帖最后由 sonictl 于 2012-12-31 12:18 编辑 请删除我 清楚超靠靠靠 没办法,一会儿限制这不能发表,那不能发表的.... [转载]Matlab中使用Plot函数动态画图方法总 ...
- matlab制作动画并保存,[转载]Matlab中动画的实现、制作和保存(转载Matlab技术论坛)...
Matlab中动画的实现.制作和保存(转载Matlab技术论坛) 本帖由Matlab技术论坛原创,原帖参见 http://www.matlabsky.com/thread-592-1-1.html 光 ...
- 结构matlab,MATLAB做晶体结构图(固体物理)
MATLAB做晶体结构图(固体物理).md 写在前面 最近在复习考研复试<固体物理>这一门课,去年学的内容已经忘干净了,所以就翻开前几页.突然看到了面心立方和体心立方结构图,想到了去年室友 ...
- 阿英讲频率f,角频率Ω和数字频率w的物理含义--附MATLAB仿真
http://anony3721.blog.163.com/blog/static/51197420111129503233/ 即以此功德,庄严佛净土.上报四重恩,下救三道苦.惟愿见闻者,悉发菩提心. ...
最新文章
- glance系列一:glance基础
- 安装Jenkins及配置自由风格自动打包部署
- 小波阈值去噪法基础理论
- 深入理解 KVC\KVO 实现机制 — KVC
- 如何改变对话或窗体视窗的背景颜色
- beeline安装_hive使用beeline配置远程连接
- 面试python说会为什么_Python面试必须要看的15个问题
- Qt5Widgetsd.lib(Qt5Widgetsd.dll) : fatalerror LNK1112: 模块计算机类型“x64”与目标计算机类型“X86”冲突
- centos5安装PHP5时遇到问题
- JAVA 实验报告总结
- 【JY】YJK前处理参数详解及常见问题分析(六):地震信息
- 74LS 系列与 74HC,74HCT,CD系列的区别
- 计算机键盘上删除,电脑哪个是删除键
- poi批量导入导出Excel(三、通过模板导出Excel并能实现导入多个sheet表)
- 【论坛】交通需求管理政策与实践——中国城市交通发展论坛第十一次研讨会讨论精选...
- c++程序设计报告总结
- 输入一个字符串转换成十进制整数
- 哪些“惯性”在2021可能被打破?(20210106).PDF
- 怎么学好数据结构?北大计算机博士亲授。
- 有趣的110道LeetCode算法题小抄,势如破竹的面试
热门文章
- x线计算机体层成像设备教案,医学影像设备学第4章-数字X线设备ppt课件
- CentOS mysql重置密码
- SpringBoot+Quartz+数据库存储
- 微信小程序使用iconfont在线icon
- ELK-日志收集工具nxlog
- 【Java】PAT乙级真题全记录(二)21到40题
- 电力拖动计算机控制系统讲什么,电力拖动自动控制系统
- linux 系统调试工具,Linux 系统调试...
- 可变模糊集理论matlab程序,可变模糊集理论.pdf
- Appium:ios自动化环境搭建