xwvdgydF4y2Ba

交叉项能量分布和交叉平滑伪能量分布gydF4y2Ba

自从R2018bgydF4y2Ba

所有的页面崩溃gydF4y2Ba

语法gydF4y2Ba

d = xwvd (x, y)gydF4y2Ba

d = xwvd (x, y, fs)gydF4y2Ba

d = xwvd (x, y, ts)gydF4y2Ba

d = xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”)gydF4y2Ba

d = xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”,双胞胎,fwin)gydF4y2Ba

d = xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”NumFrequencyPoints = nf)gydF4y2Ba

d = xwvd (gydF4y2Ba___gydF4y2Ba,MinThreshold =打)gydF4y2Ba

[d, f t] = xwvd (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

xwvd (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba,gydF4y2BaygydF4y2Ba)gydF4y2Ba返回的交叉项能量分布gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba,gydF4y2BaygydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Ba当返回交叉项能量分布gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba采样率吗gydF4y2BafsgydF4y2Ba。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba,gydF4y2BaygydF4y2Ba,gydF4y2BatsgydF4y2Ba)gydF4y2Ba当返回交叉项能量分布gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba是采样时间间隔gydF4y2BatsgydF4y2Ba之间的样本。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”)gydF4y2Ba返回的交叉平滑伪能量分布gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba。输入信号的函数使用长度选择窗口的长度用于时间和频率平滑。这个语法可以包括从以前的语法输入参数的任意组合。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”,gydF4y2Ba双胞胎gydF4y2Ba,gydF4y2BafwingydF4y2Ba)gydF4y2Ba指定的时间窗口gydF4y2Ba双胞胎gydF4y2Ba和频率窗gydF4y2BafwingydF4y2Ba用于平滑。使用默认的窗口时间或频率平滑,指定相应的参数为空,gydF4y2Ba[]gydF4y2Ba。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba、“smoothedPseudo NumFrequencyPoints =gydF4y2BanfgydF4y2Ba)gydF4y2Ba计算交叉平滑伪能量使用gydF4y2BanfgydF4y2Ba频率点。您可以指定gydF4y2Ba双胞胎gydF4y2Ba和gydF4y2BafwingydF4y2Ba在这个语法,或者你可以忽略它们。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2BaMinThreshold =gydF4y2Ba打gydF4y2Ba)gydF4y2Ba这些元素设置为零gydF4y2BadgydF4y2Ba谁的振幅小于gydF4y2Ba打gydF4y2Ba。这个语法适用于交叉项能量和交叉平滑伪能量分布。gydF4y2Ba

(gydF4y2BadgydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba)= xwvd (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba还返回一个向量的频率gydF4y2BafgydF4y2Ba和一个向量的时间gydF4y2BatgydF4y2Ba在这gydF4y2BadgydF4y2Ba计算。gydF4y2Ba

例子gydF4y2Ba

xwvd (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba没有输出参数块十字架的实数部分能量或交叉在当前图平滑伪能量。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

交叉信号的能量分布gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成两个信号采样1 kHz 1秒和嵌入在白噪声。一个信号的正弦信号频率150赫兹。另一个信号是频率啁啾的变化呈现正弦200赫兹到400赫兹。噪声的方差gydF4y2Ba $0gydF4y2Ba 。gydF4y2Ba {1gydF4y2Ba}^{2gydF4y2Ba}$ 。gydF4y2Ba

fs = 1000;t = (0:1 / fs: 1) ';x = cos(2 *π* t * 150) + 0.1 * randn(大小(t));y = vco (cos(3 *π* t) (200 400), fs) + 0.1 * randn(大小(t));gydF4y2Ba

计算的能量分布信号的总和。gydF4y2Ba

项(x + y, fs)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题项能量分布,包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

计算和情节的交叉项能量分布的信号。cross-distribution对应于能量分布的交叉项。gydF4y2Ba

fs xwvd (x, y)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题交叉项能量分布,包含时间(s), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

交叉啾啾的能量分布gydF4y2Ba

打开生活的脚本gydF4y2Ba

生成一个双通道信号,包括两个啾啾。信号采样3千赫一秒钟。第一个唧唧喳喳的初始频率400赫兹和达到800 Hz的抽样。第二啁啾起价500 Hz,最后达到1000 Hz。第二个唧唧喳喳的短促声波的两倍。gydF4y2Ba

fs = 3000;t = (0:1 / fs: 1 - 1 / fs) ';x1 =唧唧喳喳(800 t, 1400 t(结束),);x2 = 2 *唧唧喳喳(1000 t, 200 t(结束),);gydF4y2Ba

存储信号作为一个时间表。计算和情节的交叉项能量分布两个渠道。gydF4y2Ba

xt =时间表(秒(t), x1, x2);xwvd (xt (: 1), xt (:, 2))gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题交叉项能量分布,包含时间(ms), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

使用交叉项能量估计瞬时频率gydF4y2Ba

打开生活的脚本gydF4y2Ba

计算信号的瞬时频率通过使用一个已知的参考信号和交叉项能量分布。gydF4y2Ba

创建一个参考信号组成一个高斯原子1 kHz的采样1秒。一个高斯原子是一个正弦信号调制高斯。指定一个正弦信号的频率50赫兹。高斯集中在64毫秒的方差gydF4y2Ba $0gydF4y2Ba 。gydF4y2Ba {01gydF4y2Ba}^{2gydF4y2Ba}$ 。gydF4y2Ba

fs = 1 e3;t = (0:1 / fs: 1 - 1 / fs) ';μ= 0.064;σ= 0.01;傅氏国际= 50;(xr = exp () - tμ。^ 2 /(2 *σ^ 2))。*罪(2 *π*国际* t);gydF4y2Ba

创建“未知”的信号分析,组成的唧唧声。信号突然在0.4秒开始,突然半秒后结束。失误,唧唧喳喳的频率减少线性从400赫兹到100赫兹。gydF4y2Ba

f0 = 400;f1 = 100;xa = 0(大小(t));xa (t > 0.4 & t < = 0.9) =唧唧声((0:1 / fs: 0.5 - 1 / fs)”, f0, 0.5, f1);gydF4y2Ba

创建一个双组分组成的和未知的信号和参考信号。结果的平滑伪能量提供了一个“理想”的时频表示。gydF4y2Ba

计算并显示平滑伪能量分布。gydF4y2Ba

w =项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba);项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题平滑伪能量分布,包含时间(ms), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

计算的交叉项能量分布未知和参考信号。取绝对值的分布和设置为零的元素与振幅小于10。交叉项能量等于双组分信号的交叉项。gydF4y2Ba

情节十字架的实数部分能量分布。gydF4y2Ba

[c、fc、tc] = xwvd (xa xr, fs);c = abs (c);c (c < 10) = 0;xwvd (xa xr, fs)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象与标题交叉项能量分布,包含时间(ms), ylabel频率(赫兹)包含一个类型的对象的形象。gydF4y2Ba

提高能量交叉项通过添加交叉项能量分布的理想的时频表示。交叉项的能量发生中间参考信号和未知信号。gydF4y2Ba

d = w + c;d = abs(真正的(d));显示亮度图像(tc、fc、d)轴gydF4y2BaxygydF4y2BacolorbargydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象的形象。gydF4y2Ba

识别和相对应的高能脊交叉项的阴谋。孤立的山脊,找到cross-distribution有非零的时间值的能量。gydF4y2Ba

ff = tfridge (c, fc);电视=和(c) > 0;ff = ff(电视);tc = tc(电视);持有gydF4y2Ba在gydF4y2Ba情节(tc、ff,gydF4y2Ba“r——”gydF4y2Ba线宽= 2)gydF4y2Ba从gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含2图像类型的对象,线。gydF4y2Ba

重构未知信号的瞬时频率使用脊和参考功能。画出瞬时频率作为时间的函数。gydF4y2Ba

测试= 2 * tc -μ;节日= 2 * ff -国际;情节(测试、节日)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

`xgydF4y2Ba`,gydF4y2Ba`ygydF4y2Ba`- - - - - -gydF4y2Ba输入信号gydF4y2Ba
向量gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

输入信号,指定为向量或MATLABgydF4y2Ba^®gydF4y2Ba时间表每个包含一个向量变量。gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba都是向量或者两者都必须时间表和必须具有相同的长度。gydF4y2Ba

如果gydF4y2BaxgydF4y2Ba和gydF4y2BaygydF4y2Ba时间表,那么他们必须包含增加有限行乘以。gydF4y2Ba
如果时间表有缺失或重复的时间点,你可以修理它使用的技巧gydF4y2Ba清洁时间表与失踪、重复或不均匀gydF4y2Ba。gydF4y2Ba

如果输入信号有奇数长度,甚至函数附加一个零长度。gydF4y2Ba

例子:gydF4y2Ba因为(π/ 8 * (0:159))' + randn (160 1) / 10gydF4y2Ba指定一个正弦信号嵌入在白噪声。gydF4y2Ba

例子:gydF4y2Ba时间表(秒(0:5)”,兰德(6,1))gydF4y2Ba指定一个随机变量抽样1赫兹为4秒。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba
复数的支持:金宝appgydF4y2Ba是的gydF4y2Ba

`fsgydF4y2Ba`- - - - - -gydF4y2Ba采样率gydF4y2Ba
`2 *πgydF4y2Ba`(默认)|gydF4y2Ba积极的数字标量gydF4y2Ba

采样率,指定为一个积极的数字标量。gydF4y2Ba

`tsgydF4y2Ba`- - - - - -gydF4y2Ba样品时间gydF4y2Ba
`持续时间gydF4y2Ba`标量gydF4y2Ba

样品时间,指定为一个gydF4y2Ba持续时间gydF4y2Ba标量。gydF4y2Ba

`双胞胎gydF4y2Ba`,gydF4y2Ba`fwingydF4y2Ba`- - - - - -gydF4y2Ba时间和频率窗gydF4y2Ba
向量的奇数长度gydF4y2Ba

时间和频率窗口用于平滑,指定为奇数长度的向量。默认情况下,gydF4y2BaxwvdgydF4y2Ba使用Kaiser窗形状系数gydF4y2BaβgydF4y2Ba= 20gydF4y2Ba。gydF4y2Ba

默认的长度gydF4y2Ba双胞胎gydF4y2Ba是最小的奇数大于或等于gydF4y2Ba轮gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 10)gydF4y2Ba。gydF4y2Ba
默认的长度gydF4y2BafwingydF4y2Ba是最小的奇数大于或等于gydF4y2BanfgydF4y2Ba/ 4gydF4y2Ba。gydF4y2Ba

每个窗口长度必须小于或等于gydF4y2Ba2 *gydF4y2Ba装天花板gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 2)gydF4y2Ba。gydF4y2Ba

例子:gydF4y2Ba凯撒gydF4y2Ba》(65,0.5)gydF4y2Ba指定一个65 -样本Kaiser窗形状系数为0.5。gydF4y2Ba

`nfgydF4y2Ba`- - - - - -gydF4y2Ba数量的频率点gydF4y2Ba
`2 *gydF4y2Ba装天花板gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 2)gydF4y2Ba`(默认)|gydF4y2Ba整数gydF4y2Ba

的频率点,指定为一个整数。该参数控制频率采样过密的程度。必须至少频率点的数量gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BafwingydF4y2Ba)+ 1)/ 2gydF4y2Ba,不能大于默认。gydF4y2Ba

`打gydF4y2Ba`- - - - - -gydF4y2Ba最小非零值gydF4y2Ba
`负gydF4y2Ba`(默认)|gydF4y2Ba真正的标量gydF4y2Ba

最小非零值,指定为一个真正的标量。零的那些元素的函数集gydF4y2BadgydF4y2Ba谁的振幅小于gydF4y2Ba打gydF4y2Ba。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

`dgydF4y2Ba`——交叉项能量gydF4y2Ba
矩阵gydF4y2Ba

交叉项能量,作为一个矩阵返回。时间增加的列gydF4y2BadgydF4y2Ba,和频率增加了行。矩阵的大小gydF4y2BaNgydF4y2Ba_fgydF4y2Ba×gydF4y2BaNgydF4y2Ba_tgydF4y2Ba,在那里gydF4y2BaNgydF4y2Ba_fgydF4y2Ba的长度是gydF4y2BafgydF4y2Ba和gydF4y2BaNgydF4y2Ba_tgydF4y2Ba的长度是gydF4y2BatgydF4y2Ba。gydF4y2Ba

`fgydF4y2Ba`——频率gydF4y2Ba
向量gydF4y2Ba

频率,作为一个向量返回。gydF4y2Ba

如果输入的时间信息gydF4y2BafgydF4y2Ba包含赫兹的频率。gydF4y2Ba
如果输入信息,没有时间gydF4y2BafgydF4y2Ba在rad /示例包含归一化频率表示。gydF4y2Ba

`tgydF4y2Ba`——时间瞬间gydF4y2Ba
向量gydF4y2Ba

瞬间的时间,作为一个向量返回。gydF4y2Ba

如果输入的时间信息gydF4y2BatgydF4y2Ba包含时间值表示在几秒钟内。gydF4y2Ba
如果输入信息,没有时间gydF4y2BatgydF4y2Ba包含样品数量。gydF4y2Ba

时间点的数量是固定的gydF4y2Ba4 *gydF4y2Ba装天花板gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 2)gydF4y2Ba。gydF4y2Ba

引用gydF4y2Ba

[1]科恩,莱昂。gydF4y2Ba时频分析:理论和应用程序gydF4y2Ba。恩格尔伍德悬崖,台北:普伦蒂斯·霍尔出版社,1995年。gydF4y2Ba

[2]Mallat,史蒂芬。gydF4y2Ba小波信号处理gydF4y2Ba。第二版。圣地亚哥,CA:学术出版社,1999年。gydF4y2Ba

一,[3]Malnar维克多Sucic, Boualem Boashash。“一个基于组件的方法的交叉项几何瞬时频率估计使用交叉项能量。”In11信息科学国际会议上,信号处理及其应用(ISSPA)gydF4y2Ba,1217 - 1222页。蒙特利尔:IEEEgydF4y2Ba^®gydF4y2Ba,2012年。gydF4y2Ba

扩展功能gydF4y2Ba

C / c++代码生成gydF4y2Ba
生成C和c++代码使用MATLAB®编码器™。gydF4y2Ba

版本历史gydF4y2Ba

介绍了R2018bgydF4y2Ba

全部展开gydF4y2Ba

R2023a:gydF4y2Ba时间表对代码生成的支金宝app持gydF4y2Ba

的gydF4y2BaxwvdgydF4y2Ba函数支持代码生成的金宝app时间表的输入。gydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

运行cpsdgydF4y2Ba|gydF4y2BafsstgydF4y2Ba|gydF4y2BapspectrumgydF4y2Ba|gydF4y2BaxspectrogramgydF4y2Ba|gydF4y2Ba项gydF4y2Ba

主题gydF4y2Ba

时频画廊gydF4y2Ba

xwvdgydF4y2Ba

语法gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

交叉信号的能量分布gydF4y2Ba

交叉啾啾的能量分布gydF4y2Ba

使用交叉项能量估计瞬时频率gydF4y2Ba

输入参数gydF4y2Ba

xgydF4y2Ba,gydF4y2BaygydF4y2Ba- - - - - -gydF4y2Ba输入信号gydF4y2Ba向量gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

fsgydF4y2Ba- - - - - -gydF4y2Ba采样率gydF4y2Ba2 *πgydF4y2Ba(默认)|gydF4y2Ba积极的数字标量gydF4y2Ba

tsgydF4y2Ba- - - - - -gydF4y2Ba样品时间gydF4y2Ba持续时间gydF4y2Ba标量gydF4y2Ba

双胞胎gydF4y2Ba,gydF4y2BafwingydF4y2Ba- - - - - -gydF4y2Ba时间和频率窗gydF4y2Ba向量的奇数长度gydF4y2Ba

nfgydF4y2Ba- - - - - -gydF4y2Ba数量的频率点gydF4y2Ba2 *gydF4y2Ba装天花板gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 2)gydF4y2Ba(默认)|gydF4y2Ba整数gydF4y2Ba

打gydF4y2Ba- - - - - -gydF4y2Ba最小非零值gydF4y2Ba负gydF4y2Ba(默认)|gydF4y2Ba真正的标量gydF4y2Ba

输出参数gydF4y2Ba

dgydF4y2Ba——交叉项能量gydF4y2Ba矩阵gydF4y2Ba

fgydF4y2Ba——频率gydF4y2Ba向量gydF4y2Ba

tgydF4y2Ba——时间瞬间gydF4y2Ba向量gydF4y2Ba

更多关于gydF4y2Ba

交叉项能量gydF4y2Ba

引用gydF4y2Ba

扩展功能gydF4y2Ba

C / c++代码生成gydF4y2Ba生成C和c++代码使用MATLAB®编码器™。gydF4y2Ba

版本历史gydF4y2Ba

R2023a:gydF4y2Ba时间表对代码生成的支金宝app持gydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

`xgydF4y2Ba`,gydF4y2Ba`ygydF4y2Ba`- - - - - -gydF4y2Ba输入信号gydF4y2Ba
向量gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

`fsgydF4y2Ba`- - - - - -gydF4y2Ba采样率gydF4y2Ba
`2 *πgydF4y2Ba`(默认)|gydF4y2Ba积极的数字标量gydF4y2Ba

`tsgydF4y2Ba`- - - - - -gydF4y2Ba样品时间gydF4y2Ba
`持续时间gydF4y2Ba`标量gydF4y2Ba

`双胞胎gydF4y2Ba`,gydF4y2Ba`fwingydF4y2Ba`- - - - - -gydF4y2Ba时间和频率窗gydF4y2Ba
向量的奇数长度gydF4y2Ba

`nfgydF4y2Ba`- - - - - -gydF4y2Ba数量的频率点gydF4y2Ba
`2 *gydF4y2Ba装天花板gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba)/ 2)gydF4y2Ba`(默认)|gydF4y2Ba整数gydF4y2Ba

`打gydF4y2Ba`- - - - - -gydF4y2Ba最小非零值gydF4y2Ba
`负gydF4y2Ba`(默认)|gydF4y2Ba真正的标量gydF4y2Ba

`dgydF4y2Ba`——交叉项能量gydF4y2Ba
矩阵gydF4y2Ba

`fgydF4y2Ba`——频率gydF4y2Ba
向量gydF4y2Ba

`tgydF4y2Ba`——时间瞬间gydF4y2Ba
向量gydF4y2Ba

C / c++代码生成gydF4y2Ba
生成C和c++代码使用MATLAB®编码器™。gydF4y2Ba