Matlab_Marr小波尺度变换图
多尺度的连续小波变换分解后时间尺度图
墨西哥Marr小波函数:
%delta 小波变换尺度
%N 小波函数的长度
%s 原始信号
%g 原始信号某个尺度下的小波变化系数
function g=Singularity_Detection(delta,N,s);
%保持信号长度
n=length(s);
% 构造墨西哥小帽小波函数
for index_x=1:N;
x=index_x-(N+1)/2;
phi_x(index_x)=((pi^(-1/4))*(2/sqrt(3)))*(1-x.*x/(delta^2))*exp(-(x.*x)/(2*delta^2));
end;
phi_x=phi_x/norm(phi_x); %能量归一化
g=conv(s,phi_x); %卷积
g=wkeep(g,n); %保持信号长度
绘图:
clc;clear
load vonkoch
vonkoch=vonkoch(1:510);
S_Min=1;
S_Max=32;
index=0;
for scale=S_Max:-1:S_Min;
index=index+1;
cwt_coef(index,:)=Singularity_Detection(scale,32*(scale),vonkoch);
end
cwtcoef_abs=abs(cwt_coef);
for index=S_Min:S_Max
max_coef=max(cwtcoef_abs(index,:));
min_coef=min(cwtcoef_abs(index,:));
ext=max_coef-min_coef;
cwtcoef_abs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext;
end
figure(1)
subplot(2,1,1);
plot(vonkoch);
xlabel('Time');
ylabel('Amplitude');
title('Sigal');
axis([1 510 0 0.02])
subplot(2,1,2)
colormap(pink(64));
image(cwtcoef_abs)
set(gca,'YTick',2:3:32)
set(gca,'YTickLabel',32:-3:2)
title('连续小波变换时间尺度图')
xlabel('时间')
ylabel('尺度')
Matlab_Marr小波尺度变换图相关推荐
- Matlab绘制小波尺度函数和小波函数
%小波尺度函数和小波函数 clear,clc,close all; figure,[phi1,psi1,xval1] = wavefun('haar',8,'plot'); figure,[psi2, ...
- matlab小波 尺度,matlab常用小波基的尺度函数和小波函数
[实例简介] matlab自带小波基的尺度函数和小波函数 [实例截图] [核心代码] matlab小波基的尺度函数和小波函数 ├── bior1.1小波基的尺度函数和小波函数.png ├── bior ...
- 小波基函数的选择和小波构造
滤波器组完美重构与小波快速算法. 前面的分析可以知道Vj相当于在j分辨率的逼近,Vj-1相当于j-1分辨率的逼近,这样Wj-1相当于两个分辨率逼近的差.在高分辨率下,我们可以用f在(2^j*t)的采样 ...
- 认识小波——MATLAB中国的讲座(笔记)
第一部分:什么是小波 https://www.bilibili.com/video/BV1xu411R7yK 现实世界中的数据或信号经常表现出缓慢变化的趋势或者有瞬变的振荡 Real world da ...
- matlab中daubechie小波,[转载]小波滤波器
滤波器组完美重构与小波快速算法. 前面的分析可以知道Vj相当于在j分辨率的逼近,Vj-1相当于j-1分辨率的逼近,这样Wj-1相当于两个分辨率逼近的差.在高分辨率下,我们可以用f在(2^j*t)的采样 ...
- 【数字信号去噪】基于matlab小波阙值数字信号去噪和求信噪比【含Matlab源码 2191期】
⛄一.小波阈值法去噪概述 电能质量扰动信号的噪声大多以高斯白噪声的形式存在,利用小波变换对信号进行多分辨率分解,由于小波变换具有去除数据相关性的特点,故可以将有用信号与噪声的能量分离开来.信号中有效的 ...
- 【心电信号】基于matlab小波阙值心电信号去噪【含Matlab源码 2188期】
⛄一.小波阈值法去噪概述 电能质量扰动信号的噪声大多以高斯白噪声的形式存在,利用小波变换对信号进行多分辨率分解,由于小波变换具有去除数据相关性的特点,故可以将有用信号与噪声的能量分离开来.信号中有效的 ...
- 基于小波Elman神经网络的短期风电功率预测
风力发电在全球范围内快速发展,装机容量逐年增加,截止2013 年底,中国风电新增装机容量约 16.1GW,较 2012 年的 12.96GW 大幅提高了 24%,中国风电累计装机已超过 90GW. 风 ...
- 常用数据特征提取,时域特征、频域特征、小波特征提取汇总;特征提取;有效matlab代码
clc;clear %% 导入数据 load('ct.mat') Fs=12800; a=[];c=[];w=[]; %% 1时域特征提取 for i=2:8 y=ct(:,i); a(1,i) = ...
最新文章
- 过滤输入内容中是否含有特殊字符与表情
- linux下基本命令
- attribute 'downsample' of type 'NoneType' is not usable in a script method
- WebStorm 提示快捷键
- Android持久化存储(1)文件存储
- 弹性式分布数据集RDD——Pyspark基础 (二)
- javafx 调用接口_JavaFX技巧3:使用回调接口
- update和delete操作忘加where条件导致全表更新的处理方法
- POJ 1155 TELE 树形DP
- PyCharm安装配置教程
- Unity 3D 实战核心技术详解
- 多变量微积分(4)——多重积分之三重积分
- 学习pathon的几大步骤
- 音程示范歌曲,lilypond example
- 为什么在线性模型中相互作用的变量要相乘
- word文档中的毕业论文的页眉与页脚设置
- es拼音分词 大帅哥_elasticsearch实现中文分词和拼音分词混合查询+CompletionSuggestion...
- 数据同步工具Sqoop
- 【收藏】最全计算机网络基础思维导图
- SEO诊断,如何查询网站是否被百度降权