1.软件版本

matlab2013b

2.本算法理论知识

参考文献:

《广义S变换时频分析的应用研究》

《时频分布与地震信号谱分析研究》

《非平稳信号广义S不变换及其在SAR图像分析中的应用研究》

《S变换时变滤波在去噪处理中的应用研究》

《广义S变换域时频特征分析及微弱目标检测方法研究》

《广义S变换及其时频滤波》

《时频分析在地震信号处理中的研究及应用》

3.部分源码

clc;
clear;
close all;
warning off;
addpath 'func\S\'
addpath 'func\GS\'%产生线性调频信号
len         = 128;
t           = 0:len-1;
figure;
subplot(221);
y = cos(2*pi*(10+t/8).*t/len);
plot(y);
title('线性调频信号');%一维S变换
subplot(222);
[ST_matrix,ST_times,ST_frequencies] = func_ST(y);
contourf(ST_times,ST_frequencies,abs(ST_matrix));
title('一维S变换');%一维广义S变换
P=0.1;
R=1;
K=2;
subplot(223);
[GST_matrix,GST_times,GST_frequencies] = func_GST(y,P,R,K);
contourf(GST_times,GST_frequencies,abs(GST_matrix));
title('一维广义S变换');%一维广义S逆变换
subplot(224);
y2 = func_inv_GST(GST_matrix);
plot(y2);
title('线性调频信号');

4.仿真分析

注意,广义S和一般S的区别就是广义S具有调节参数。关于不同参数下广义S变换的窗函数的仿真的结果如下所示:

   基于广义S变换的时频滤波器。对于一个非平稳信号加上噪声,要把噪声滤除掉;而且要知道前后的信噪比。最后,再把设计得时频滤波器用到人工合成的地震信号和实际地震信号里。另外,要对这个信号用阈值滤波、小波软阈值滤波、时频滤波进行比较滤波前的信噪比和滤波后的信噪比。

这里主要是几个算法的滤波性能的对比。仿真结果如下所示:

         用二维广义S变换去除二维图像随机噪声,比较无噪声和加入噪声后去除噪声的图像比较及噪声比。程序要有二维广义S变换的算法实现。

这里,将使用二维图像的广义S变换处理:

注意,通过仿真可知,如果当图像像素大于50的时候,仿真速度及其缓慢,所以在实际进行二维图像滤波的时候哦,我们一般使用局部滤波的思路,即将图像分为多个局部图像,然后再将其二维的局部图转换为一维序列,然后完成滤波。

5.参考文献

[1]高静怀, 陈文超, 李幼铭,等. 广义S变换与薄互层地震响应分析[J]. 地球物理学报, 2003, 46(4):7.A28-09

【广义S变换】一维广义S变换对非平稳信号处理的matlab仿真相关推荐

  1. 图像mnf正变换_基于MNF 变换的多元变化检测变化信息的集中

    展开全部 (一)基于MNF变换的变化信息集中的基本流程 MNF变换是一种多元62616964757a686964616fe58685e5aeb931333433616234线性统计变换方法,是针对一组 ...

  2. 广义线性模型GLM、GLMM、LMM、MLM、GMM、GEE、广义线性模型GLM和广义线性混合模型的GLMM区别

    R语言广义线性模型GLM.GLMM.LMM.MLM.GMM.GEE.广义线性模型GLM和广义线性混合模型的GLMM区别 目录

  3. 逆clarke变换_是clarke变换还是clark

    匿名用户 1级 2017-06-01 回答 是clarke变换还是clark 第三,关于abc坐标系,这是最重要的:四,:派克变换的初步推导:先推倒从abc系到静止的xyz系的过度矩阵M1,:这三个基 ...

  4. 理解离散傅立叶变换(一. 傅立叶变换的由来)

    理解离散傅立叶变换(一) ------傅立叶变换的由来 关于傅立叶变换,不管是书本还是在网上可以非常easy找到关于傅立叶变换的描写叙述,可是大都是些故弄玄虚的文章,太过抽象,尽是一些让人看了就望而生 ...

  5. K-L变换原理、图像压缩与人脸识别实战(附matlab代码)

    目录 摘要 K-L变换的基本原理 K-L变换在图像压缩中的应用 研究背景及意义 数字图像压缩与发展现状 K-L变换与图像压缩 K-L变换在人脸识别中的应用 人脸识别 K-L变换实现人脸识别 总结 参考 ...

  6. opencv-python:16_形态学处理【二】(开操作、闭操作、形态学梯度、顶帽变换、黑帽变换,去除皮肤镜中的毛发噪音、cv2.morphologyEx())

    形态学处理[二] 开操作.闭操作.形态学梯度.顶帽变换.黑帽变换 相关函数有:cv2.morphologyEx().cv2.getStructuringElement() 有趣的应用:去除皮肤镜中的毛 ...

  7. matlab电流检测,基于d-q变换的谐波电流检测方法与MATLAB仿真验证

    描述 引言 目前,谐波抑制的一个重要趋势是采用有源电力滤波器(Active Power Filter-APF).而该滤波器性能的好坏与它所采用的谐波电流检测方法有很大关系.因此,如何实时准确地检测出非 ...

  8. Matlab 伪彩色处理方法总结(密度分割法、灰度级变换法、频域变换法)

    伪彩色处理方法总结 伪彩色处理是将黑白图像转换为彩色图像,方法分为空域变换及频域变换.空域变换其基本原理是构建颜色映射函数,将灰度值转换为彩色值.因为人眼对彩色图像的分辨能力大于黑白图像,所以伪彩色处 ...

  9. MATLAB——拉氏变换及反其变换、Z变换及其反变换

    一.拉式正反变换 1.拉式变换 [时域转s域] laplace(F) 2.拉式反变换 [s域转时域] ilaplace(L) 二.Z正反变换 1.Z变换 [ztrans 时域转Z域] ztrans(f ...

  10. SLNR预编码:广义瑞丽熵与广义特征向量

    前言 这两天在调研MIMO预编码算法文献的时候,恰巧看到了一篇这篇paper: A Leakage-Based Precoding Scheme for Downlink Multi-User MIM ...

最新文章

  1. 汇编和python-编程语言与Python介绍
  2. ios开发 多人语音聊天_手游语音市场的现状、机遇与挑战
  3. Python办公自动化,对文件进行自由操作
  4. centos7通过阿里云配置docker加速镜像
  5. 戴尔G7 7588 Windows 10 未安装任何音频输出设备
  6. 带SN切换流程_一文详解贴片电阻生产工艺流程(15个步骤)
  7. 交换机不同vlan不同网段通过核心交换机配置VLANIF通信
  8. MarkDown语法学习--字体,分割线,删除线,下划线
  9. 感觉丧的时候,读一读曾国藩
  10. BMFont 制作字体时,无法导入图片
  11. c语言我喜欢你,【幻 仿】UC启动时的我喜欢你
  12. 迅雷协议分析–多链接资源获取
  13. python抢购茅台抢购描述
  14. 提供家用咖啡电子秤厨房秤方案芯片CS1238
  15. python 字符串和数字组合变量名
  16. 针对尚硅谷教学微服务硅谷课堂在线学习平台的部分功能的自我改进
  17. 俄罗斯起草法案建议设立离岸公司交易加密货币
  18. 计世网:IT人忧伤的五个阶段
  19. RHCE---请给openlab搭建web网站
  20. 百度、富士康、恒大跑去造车,蔚来、威马、理想、小鹏必须报团取暖?

热门文章

  1. 计算机二级考试vb内容,计算机二级VB考试考什么内容?(一)
  2. iphone ping_如何在iPhone上运行Ping(网络诊断)
  3. [备忘]方正字库中英文对照表
  4. linux ftp下载geo,GEO数据下载及处理详细过程
  5. python2.7中文手册下载_python 2.7.12中文文档下载|
  6. 计算机专业英语交换机,计算机专业英语词汇
  7. (二)网络传输协议:网络摄像机传输协议、流媒体传输协议
  8. SSE优化系列十:RGB和YUV格式图像互转SSE优化
  9. c语言sqlist结构体,数据结构的一个题目,有什么问题吗?问什么一直显示Sqlist结构体没定义...
  10. eova、easyui及普通DOM表单元素操作