【数字信号处理】基于matlab数字信号同步压缩变换【含Matlab源码 1534期】
一、获取代码方式
获取代码方式1:
完整代码已上传我的资源: 【数字信号处理】基于matlab数字信号同步压缩变换【含Matlab源码 1534期】
获取代码方式2:
通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码。
备注:订阅紫极神光博客付费专栏,可免费获得1份代码(有效期为订阅日起,三天内有效);
二、同步压缩变换简介
SST以WT为基础,利用同步压缩算子提高时频脊线在时频谱上的分辨率,实现了对瞬时频率的提取与重构。设ψ(b)为小波母函数,则信号x(t)的连续小波变换为:
式中:
x(t)——振动信号;
W(a,b)——x(t)的连续小波变换结果;
t——时间变量;
a——尺度因子;
b——平移因子;
通过分析可知,在小波域中(a,b)位置的瞬时频率信息为:
式中:
ωx(a,b)——瞬时频率;
j——虚数单位。
文献[10]发现,无论a取何值,W(a,b)在b上的振荡特性均指向初始频率Ω,因此:
根据定义的同步压缩变换,小波逆变换为:
式中:
x(b)——小波逆变换结果;
Cψ——相差系数;
ψ(aξ)——小波母函数。
对ωx(a,b)沿尺度a方向整合,归划到频域中ω=ωx(a,b)的位置上,则同步压缩变换定义为:
式中:
Sst (ω,b)——信号b的同步压缩函数;
ω——角频率。
通过式(6)的结果与相差系数Cψ,将信号的幅值归化到时频域中所在位置,最终获得高分辨率的时频谱。
三、部分源代码
% A numerical signal.
clear;
SampFreq = 100;t = 0 : 1/SampFreq : 14-1/SampFreq;Sig = [sin(2*pi*(25*t))];[m,n]=size(Sig);time=(1:n)/SampFreq;
fre=(SampFreq/2)/(n/2):(SampFreq/2)/(n/2):(SampFreq/2);Ts = SST(Sig',100);figure
imagesc(time,fre,abs(Ts));
axis xy
ylabel('Freq / Hz');
xlabel('Time / Sec');
title('SST');%signal reconstruction.
s=real(sum(Ts));%Reconstructed signal.
figure
plot(s);
title('Reconstructed signal');%Original signal
figure
plot(Sig,'r-');
title('Original signal');
function [Ts] = SST(x,hlength);
% Computes the SST (Ts) of the signal x.
% INPUT
% x : Signal needed to be column vector.
% hlength: The hlength of window function.
% OUTPUT
% Ts : The SST[xrow,xcol] = size(x);if (xcol~=1),error('X must be column vector');
end; if (nargin < 1),
error('At least 1 parameter is required');
end;if (nargin < 2),
hlength=round(xrow/5);
end;Siglength=xrow;
hlength=hlength+1-rem(hlength,2);
ht = linspace(-0.5,0.5,hlength);ht=ht';% Gaussian window
h = exp(-pi/0.32^2*ht.^2);
% derivative of window
dh = -2*pi/0.32^2*ht .* h; % g'[hrow,hcol]=size(h); Lh=(hrow-1)/2; N=xrow;
t=1:xrow;[trow,tcol] = size(t);tfr1= zeros (N,tcol) ;
tfr2= zeros (N,tcol) ; tfr= zeros (round(N/2),tcol) ;
Ts= zeros (round(N/2),tcol) ;
四、运行结果
五、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
[4]高彦彦,张晶,李莉,贾英茜.基于GUI的《数字信号处理》教学演示系统的设计[J].教育教学论坛. 2019,(48)
[5]李俊,张淑玲,帅晶.基于Matlab GUI界面的数字信号处理辅助教学系统[J].信息通信. 2020,(08)
[6]李亦佳,王静,王正方,隋青美.基于多重同步压缩变换的微震信号去噪方法研究[J].应用基础与工程科学学报
【数字信号处理】基于matlab数字信号同步压缩变换【含Matlab源码 1534期】相关推荐
- 【Matlab验证码识别】遗传算法和最大熵优化+大津法(OTSU)+自定义阈值数字验证码识别【含GUI源码 1694期】
一.代码运行视频(哔哩哔哩) [Matlab验证码识别]遗传算法和最大熵优化+大津法(OTSU)+自定义阈值数字验证码识别[含GUI源码 1694期] 二.matlab版本及参考文献 1 matlab ...
- 【Matlab人脸识别】KL变换人脸识别【含GUI源码 859期】
一.代码运行视频(哔哩哔哩) [Matlab人脸识别]KL变换人脸识别[含GUI源码 859期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 蔡利梅.MAT ...
- 【Matlab答题卡识别】hough变换答题卡判定与成绩统计【含GUI源码 752期】
一.代码运行视频(哔哩哔哩) [Matlab答题卡识别]hough变换答题卡判定与成绩统计[含GUI源码 752期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [ ...
- 【Matlab答题卡识别】hough变换答题卡判定与成绩统计(带面板)【含GUI源码 1017期】
一.代码运行视频(哔哩哔哩) [Matlab答题卡识别]hough变换答题卡判定与成绩统计(带面板)[含GUI源码 1017期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 ...
- 【Matlab肌电信号】肌电信号处理【含GUI源码 966期】
一.代码运行视频(哔哩哔哩) [Matlab肌电信号]肌电信号处理[含GUI源码 966期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 包子阳,余继周,杨 ...
- 【Matlab语音分析】语音信号分析【含GUI源码 1718期】
一.代码运行视频(哔哩哔哩) [Matlab语音分析]语音信号分析[含GUI源码 1718期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊,郑铁 ...
- 【Matlab语音加密】语音信号加密解密(带面板)【含GUI源码 181期】
一.代码运行视频(哔哩哔哩) [Matlab语音加密]语音信号加密解密(带面板)[含GUI源码 181期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆 ...
- 【Matlab语音处理】声音信号频谱分析仪【含GUI源码 325期】
一.代码运行视频(哔哩哔哩) [Matlab语音处理]声音信号频谱分析仪[含GUI源码 325期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊, ...
- 【Matlab心音信号】EMD心音信号特征提取【含GUI源码 1735期】
一.代码运行视频(哔哩哔哩) [Matlab心音信号]EMD心音信号特征提取[含GUI源码 1735期] 二.matlab版本及参考文献 1 matlab版本 2014a *2 参考文献 [1] 沈再 ...
- 【Matlab图像加密】正交拉丁方置乱算法图像加解密【含GUI源码 182期】
一.代码运行视频(哔哩哔哩) [Matlab图像加密]正交拉丁方置乱算法图像加解密[含GUI源码 182期] 二.matlab版本及参考文献 一.代码运行视频(哔哩哔哩) [Matlab图像处理]自动 ...
最新文章
- 第七篇:并发-恢复机制
- 高行健---江西赣州人
- 欢乐纪中A组赛【2019.8.17】
- php sql好处,关于php:使用Doctrine与SQL相比,使用DQL有什么好处?
- 谷歌浏览器32位安装包_谷歌团队新作!只需下载3M安装包,就能让你的手机浏览器跟踪眼球运动...
- c#使用webbrowser时,设定IE版本
- matlab 频域采样定理,信号时域和频域采样函数周期性与原信号的关系
- 主板检测卡c5_主板检测卡代码大全
- python禁用路径长度限制有啥影响吗_为什么Windows中存在260个字符的路径长度限制?...
- 搭建asp会议签到系统 第二章 生成会议签到二维码
- CompletableFuture学习
- ps 证件照(1,2寸)
- c语言中符号是什么作用是什么,c语言中的符号|=是什么意思?
- 精辟,16张图说透Modbus-RTU协议
- 不知道O2O产品的APP该如何运营推广?
- 前端开发_HTML5_布局-栅格布局
- 2037:【例5.4】约瑟夫问题
- Python+opencv调用摄像头实现拍照并保存
- 核心单词Word List 44
- python中引入sql的优点_引用sql-和引用sql相关的内容-阿里云开发者社区