信号的时频分析MATLAB,[转载]时频特性分析(Matlab)
最近瞅了一些关于时频分析工具箱的matlab函数使用方法,总结一下吧.
我使用的是2011a的matlab,貌似没有自带的时频分析工具箱,可以到网上下载,google一搜就能搜到,安装后就可以使用了(所谓安装就是把工具箱的目录包含到matlab工作目录中即可).
下面说一些时频工具箱函数的用法(由于我下的工具箱没有html版的帮助,所以只能使用help
funname的方法查看帮助信息).
AMGAUSS
Generate gaussian amplitude modulation.
Y=AMGAUSS(N,T0,T) generates a gaussian amplitude modulation
centered on a time T0, and with a spread proportional to T.
This modulation is scaled such that Y(T0)=1
and Y(T0+T/2) and Y(T0-T/2) are approximately equal to 0.5 .
产生高斯幅值调制信号,其以信号点数为N,中心为T0,传播0.5的比例到T,即Y(T0)=1,Y(T0+T/2) =
Y(T0-T/2) = 0.5。 T0默认为N/2,T默认为2*sqrt(N)。
如z=amgauss(160,90,40); plot(z);
FMCONST
Signal with constant frequency modulation.
[Y,IFLAW] = FMCONST(N,FNORM,T0) generates a frequency
modulation
with a constant frequency fnorm.
The phase of this modulation is such that y(t0)=1.
产生一个固定频率的频率调制信号。N为产生信号的点数,FNORM为标准化频率(默认为0.25),T0表示此时刻为正幅值(默认为round(N/2)),相当于规定了相位。IFLAW为设置的频率变化情况,这里就是一个常数。
如z=fmconst(128,0.05,50); plot(real(z));
FMLIN
Signal with linear frequency modulation.
[Y,IFLAW]=FMLIN(N,FNORMI,FNORMF,T0) generates a linear
frequency
modulation.
The phase of this modulation is such that Y(T0)=1.
产生一个线性调频信号。N为信号的点数,FNORMI为开始标准频率(默认为0.0),FNORMF为结束标准频率(默认为0.5),相位由T0规定,即Y(T0)=1。IFLAW为设置的频率变化情况,这里就是一个线性变化的直线。
如[z, f]=fmlin(128,0.05,0.3,50); plot(real(z));figure,
plot(f)
以上是比较常用的产生信号的函数。
下面说一些时频分析常用的分析方法,主要有短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD等,其他还有很多分析方法,都有相关函数,具体可以查看帮助。
短时傅里叶变换STFT
功能:计算时间序列的短时离散傅里叶变换,得到瞬时频率。
格式:
[tfr, t, f] = tfrstft(x) %
计算时间序列x的短时傅里叶变换,参数tfr为短时傅里叶变换系数,t为系数tfr对应的时刻,f为归一化频率向量
[tfr, t, f] = tfrstft(x, t) % 计算对应时刻t的短时傅里叶变换
[tfr, t, f] = tfrstft(x, t, n) % 计算n点对应时刻t的短时傅里叶变换
[tfr, t, f] = tfrstft(x, t, n, h) % 参数h为归一化频率平滑窗
[tfr, t, f] = tfrstft(x t, n, h, trace) % trace显示算法进程
说明:
x--信号
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
h--频率滑窗,h归一化为单位能量(缺省值为hamming(n/4))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
如
sig=[fmconst(128,0.2);
fmconst(128,0.4)]; % 产生由两个常值调频信号(即正弦信号)的组合信号
tfr=tfrstft(sig);
subplot(211); imagesc(abs(tfr)); % 瞬时频率
subplot(212); imagesc(angle(tfr)); % 瞬时相位
Wigner-Ville时频分布图
功能:计算时间序列的Wigner-Ville时频分布图,得到瞬时频率
格式:
[tfr, t, f] = tfrwv(x)
[tfr, t, f] = tfrwv(x, t)
[tfr, t, f] = tfrwv(x, t, n)
[tfr, t, f] = tfrwv(x, t, n, trace)
说明:
x--信号
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
如
sig=fmlin(128,0.1,0.4);
tfrwv(sig);
伪Wigner-Ville时频分布图
功能:计算时间序列的伪Wigner-Ville时频分布图,得到瞬时频率
格式:
[tfr, t, f] = tfrpwv(x)
[tfr, t, f] = tfrpwv(x, t)
[tfr, t, f] = tfrpwv(x, t, n)
[tfr, t, f] = tfrpwv(x, t, n, trace)
说明:
x--信号
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
f--归一化频率
如
sig=fmlin(128,0.1,0.4);
tfrpwv(sig);
还可以设置成三维显示,如下
其他时频分析
fmt
计算时间序列的快速梅林变换
ifmt
计算快速梅林逆变换。
tfrbj
生成时间序列的Bord-Jondan时频分布图,得到瞬时频率。
tfrbud
生成时间序列的Butterworth时频分布图,得到瞬时频率。
tfrcw
生成时间序列的Choi-Williams时频分布图,得到瞬时频率。
tfrgrd
生成时间序列的广义矩形时频分布图,得到瞬时频率。
tfrmh
生成时间序列的Margenau-Hill时频分布图,得到瞬时频率。
tfrmhs
生成时间序列的Margenau-HillP谱图时频分布图,得到瞬时频率。
tfrsp
计算时间序列的频谱图分布,得到瞬时频率。
tfrscalo
计算时间序列的尺度图时频表示,即连续小波变换的幅值的平方。
tfrrsp
计算时间序列的时频分布和它的重排形式,得到瞬时频率。
friedman
计算时间序列的瞬时频率密度。
ridges
从时间序列的重排时频分布中提取瘠。
tfrideal
计算给定时间序列的瞬时频率规律的理想时频表示。
以上给了很多函数,但常用的基本就是短时傅里叶变换STFT、Wigner-Ville分布WVD、伪Wigner-Ville分布PWVD和提取瘠ridges、得到理想时频表示tfrideal。
信号的时频分析MATLAB,[转载]时频特性分析(Matlab)相关推荐
- union matlab,[转载]intersect,unique,union在matlab中的用法(I lo
Matlab中集合的交运算函数intersect intersect Set intersection of two vectors Syntax c = intersect(A,B) c = int ...
- 时频特性分析(Matlab)
最近瞅了一些关于时频分析工具箱的matlab函数使用方法,总结一下吧. 我使用的是2011a的matlab,貌似没有自带的时频分析工具箱,可以到网上下载,google一搜就能搜到,安装后就可以使用了( ...
- 地震 位错模型 matlab,8解(EEMD)的汶川Ms8[1].0强震动记录时频特性分析
第34卷第3期地震学报V01.34,No.32012年5月(350-362)ACTASEISMOLOGICASINICAMay,201Z李大虎,赖敏,何强,马新欣,顾勤平.2012.基于聚类经验模态分 ...
- 基于MATLAB的语音信号的时域特性分析(一)——分帧、窗函数
语音信号是一个非稳态的.时变的信号.但在"短时间"范围内可以认为语音信号是稳态的.时不变的.这个短时间一般值10~30ms.由于有这个特性,故常把语音信号称为"准稳态 ...
- matlab仿真时三相电流设置,三相输入电流波形与三相整流MATLAB仿真
三相输入电流波形与三相整流MATLAB仿真 三相全桥整流输入电流与MATLAB仿真 本文简要介绍三相全桥整流,并建立MATLAB仿真模型.针对输入电流波形做一定的分析,简要介绍Multimeter 和 ...
- 同频干扰信号 仿真 matlab,直接序列扩频系统抗同频干扰的MATLAB仿真8.29(可靠性分析12月)...
直接序列扩频系统抗同频干扰的MATLAB仿真Simulation of DSSS System with Co-Channel Interference by MATLAB 陈吉文,郭伟(西安电子科技 ...
- 添加布林带择时策略有多便捷!股票量化分析工具QTYX-V2.4.7
前言 布林带通道(Bollinger Bands)是非常经典的技术指标,常用于研判市场中长期运动趋势. 比如我们以[350, 2,2] 这组长线参数来绘制恒瑞医药.贵州茅台10年行情走势的布林带通道, ...
- 【通信仿真】基于matlab空时自适应处理降维仿真【含Matlab源码 1957期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[通信仿真]基于matlab空时自适应处理降维仿真[含Matlab源码 1957期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式 ...
- matlab信号分析实验报告,信号分析及处理实验报告(基于MATLAB).doc
WORD 格式 整理 PAGE 学习 参考 资料 分享 武 汉 工 程 大 学 电气信息学院 专业班级 14通信3班 实验时间 12-16周周二上午 学生学号 1404200529 实验地点 4B31 ...
最新文章
- 这38篇原创文章,带我入门深度学习!
- php umount强制,linux mount挂载与umount 卸载及“Device is busy”问题 | 璞玉(POOY)
- POJ 1661 Help Jimmy(递推DP)
- DataFrame/Series获取列名以及更改列名(转)
- 微服务介绍及Asp.net Core实战项目系列之微服务介绍
- java svn 版本号_eclipse中的Java文件自动根据svn版本号生成注释
- 将字符串或数字倒序输出
- 用户控件中图片路径问题(用户控件、图片路径) ,ResolveUrl
- STM32F407 硬件IIC驱动MCP4017 数字电位器
- java毕业设计大数据在线考试系统在线阅卷系统及大数据统计分析源码+系统+数据库+lw文档+调试运行
- 简繁体计算机术语对照表
- I2C电路故障排除---边沿时间与杜邦线
- java实现微信小程序获取手机号(htts接口实现)
- 孙鑫老师VC笔记(转)
- Python-docx 读写 Word 文档:读取正文、表格文本信息、段落格式、字体格式等
- 计算机控制运行内存,运行内存
- NLP 基础应用研究方向简介
- 消除oracle 中 job 的时间误差?
- 百度地图API支持HTTPS
- 软件工程复习13:软件发布
热门文章
- S4HANA事务码MM01创建物料的调试
- where and when is Angular scope initialized
- Why Opportunity list is empty
- 步骤4 - 微服务提供者接收请求,提供服务并传回给Orchestra
- KPI tile在launchpad里正常显示number,但是进不了detail view
- ABAP并发计算的一个实例
- UDO report generate DDIC table
- 后台nodejs程序如何主动推送数据给浏览器
- 为什么SAP UI5框架在应用整个生命周期只调用onBeforeRendering一次
- JDBC + SAP云平台 = 运行在云端的数据库应用