主要内容gydF4y2Ba

pmtmgydF4y2Ba

多锥功率谱密度估计gydF4y2Ba

折叠所有页面gydF4y2Ba

语法gydF4y2Ba

PXX = pmtm(x)gydF4y2Ba
pxx = pmtm(x，'Tapers'，tapertype)gydF4y2Ba
PXX = pmtm(x,nw)gydF4y2Ba
pxx = pmtm(x,m，'Tapers'，'sin ')gydF4y2Ba
PXX = pmtm(gydF4y2Ba___gydF4y2Banfft)gydF4y2Ba
[pxx,w] = pmtm(gydF4y2Ba___gydF4y2Ba）gydF4y2Ba
[pxx,f] = pmtm(gydF4y2Ba___gydF4y2Bafs)gydF4y2Ba
[pxx,w] = pmtm(x,nw,w)gydF4y2Ba
[pxx,w] = pmtm(x,m，'Tapers'，'sin '，w)gydF4y2Ba
[pxx,f] = pmtm(gydF4y2Ba___gydF4y2Ba, f, f)gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Bafreqrange)gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba，pxxc] = pmtm(gydF4y2Ba___gydF4y2Ba、“ConfidenceLevel”概率)gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba、“DropLastTaper”dropflag)gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba、方法)gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba= pmtm(x,e,v，gydF4y2Ba___gydF4y2Ba）gydF4y2Ba
［gydF4y2Ba___gydF4y2Ba] = pmtm(x,dpss_params，gydF4y2Ba___gydF4y2Ba）gydF4y2Ba
pmtm (gydF4y2Ba___gydF4y2Ba）gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

pxxgydF4y2Ba= pmtm (gydF4y2BaxgydF4y2Ba）gydF4y2Ba返回Thomson的多锥功率谱密度(PSD)估值，gydF4y2BapxxgydF4y2Ba为输入信号的gydF4y2BaxgydF4y2Ba使用gydF4y2Ba离散长球状(Slepian)序列gydF4y2Ba蜡烛。gydF4y2Ba

例子gydF4y2Ba

pxxgydF4y2Ba= pmtm (gydF4y2BaxgydF4y2Ba“蜡烛”,gydF4y2BatapertypegydF4y2Ba）gydF4y2Ba指定计算多锥度PSD估计时使用的锥度类型。您可以指定gydF4y2Ba“蜡烛”gydF4y2Ba，gydF4y2BatapertypegydF4y2Ba之后任意位置的名称-值对gydF4y2BaxgydF4y2Ba在函数调用中。gydF4y2Ba

例子gydF4y2Ba

pxxgydF4y2Ba= pmtm (gydF4y2BaxgydF4y2Ba，gydF4y2Ba西北gydF4y2Ba）gydF4y2Ba使用时半带宽乘积gydF4y2Ba西北gydF4y2Ba来控制频率分辨率时，计算PSD估计使用Slepian锥形。gydF4y2Ba

例子gydF4y2Ba

pxxgydF4y2Ba= pmtm (gydF4y2BaxgydF4y2Ba，gydF4y2Ba米gydF4y2Ba“蜡烛”,“正弦”)gydF4y2Ba指定计算PSD估计时使用的锥形数或平均权重gydF4y2Ba正弦蜡烛gydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

pxxgydF4y2Ba= pmtm (gydF4y2Ba___gydF4y2Ba，gydF4y2BanfftgydF4y2Ba）gydF4y2Ba使用gydF4y2BanfftgydF4y2Ba离散傅里叶变换(DFT)点结合任何以前的语法。如果gydF4y2BanfftgydF4y2Ba大于信号长度，gydF4y2BaxgydF4y2Ba是零填充长度gydF4y2BanfftgydF4y2Ba．如果gydF4y2BanfftgydF4y2Ba小于信号长度，信号被模包裹gydF4y2BanfftgydF4y2Ba．gydF4y2Ba

［gydF4y2BapxxgydF4y2Ba，gydF4y2BawgydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba）gydF4y2Ba返回一个具有归一化频率的向量gydF4y2BapxxgydF4y2Ba计算。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BapxxgydF4y2Ba，gydF4y2BafgydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Ba返回一个频率向量，gydF4y2BafgydF4y2Ba，单位为周期/单位时间。gydF4y2BafsgydF4y2Ba必须遵循gydF4y2BaxgydF4y2Ba，gydF4y2Ba西北gydF4y2Ba(或gydF4y2Ba米gydF4y2Ba对于正弦锥度)，和gydF4y2BanfftgydF4y2Ba在函数调用中。要输入一个采样率，并且仍然使用前面参数的默认值，请将这些参数指定为空，gydF4y2Ba［］gydF4y2Ba．gydF4y2Ba

［gydF4y2BapxxgydF4y2Ba，gydF4y2BawgydF4y2Ba] = pmtm(gydF4y2BaxgydF4y2Ba，gydF4y2Ba西北gydF4y2Ba，gydF4y2BawgydF4y2Ba）gydF4y2Ba中指定的归一化频率的Slepian序列计算的多锥PSD估计gydF4y2BawgydF4y2Ba．向量gydF4y2BawgydF4y2Ba必须包含至少两个元素。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BapxxgydF4y2Ba，gydF4y2BawgydF4y2Ba] = pmtm(gydF4y2BaxgydF4y2Ba，gydF4y2Ba米gydF4y2Ba“蜡烛”,“正弦”,gydF4y2BawgydF4y2Ba）gydF4y2Ba返回中指定的归一化频率上使用正弦锥度计算的多锥度PSD估计gydF4y2BawgydF4y2Ba．向量gydF4y2BawgydF4y2Ba必须包含至少两个元素。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BapxxgydF4y2Ba，gydF4y2BafgydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba，gydF4y2BafgydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Ba计算中指定频率的多锥PSD估计gydF4y2BafgydF4y2Ba．向量gydF4y2BafgydF4y2Ba必须包含至少两个与采样率相同单位的元素gydF4y2BafsgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba，gydF4y2BafreqrangegydF4y2Ba）gydF4y2Ba返回指定频率范围内的多锥PSD估计值gydF4y2BafreqrangegydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba，gydF4y2BapxxcgydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba“ConfidenceLevel”,gydF4y2Ba概率gydF4y2Ba）gydF4y2Ba返回gydF4y2Ba概率gydF4y2Ba中PSD估计的× 100%置信区间gydF4y2BapxxcgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba“DropLastTaper”,gydF4y2BadropflaggydF4y2Ba）gydF4y2Ba指定是否gydF4y2BapmtmgydF4y2Ba在计算多锥PSD估计时，下降最后一个Slepian锥度。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2Ba___gydF4y2Ba，gydF4y2Ba方法gydF4y2Ba）gydF4y2Ba中指定的方法组合各个锥形PSD估计值gydF4y2Ba方法gydF4y2Ba．此语法仅适用于Slepian锥形体。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2BaxgydF4y2Ba，gydF4y2BaegydF4y2Ba，gydF4y2BavgydF4y2Ba，gydF4y2Ba___gydF4y2Ba）gydF4y2Ba使用斯莱皮安锥形gydF4y2BaegydF4y2Ba特征值在gydF4y2BavgydF4y2Ba来计算PSD。使用gydF4y2Ba离散长gydF4y2Ba获得gydF4y2BaegydF4y2Ba而且gydF4y2BavgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

［gydF4y2Ba___gydF4y2Ba] = pmtm(gydF4y2BaxgydF4y2Ba，gydF4y2Badpss_paramsgydF4y2Ba，gydF4y2Ba___gydF4y2Ba）gydF4y2Ba使用单元格数组gydF4y2Badpss_paramsgydF4y2Ba将输入参数传递给gydF4y2Ba离散长gydF4y2Ba．此语法仅适用于Slepian锥形体。gydF4y2Ba

例子gydF4y2Ba

pmtm (gydF4y2Ba___gydF4y2Ba）gydF4y2Ba在没有输出参数的情况下，在当前图形窗口中绘制多锥PSD估计值。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开实时脚本gydF4y2Ba

获得由角频率为的离散时间正弦信号组成的输入信号的多锥PSD估计gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 加添加剂的Rad /samplegydF4y2BaNgydF4y2Ba(0,1)白噪声。gydF4y2Ba

创建一个角频率为的正弦波gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 加添加剂的Rad /samplegydF4y2BaNgydF4y2Ba(0,1)白噪声。信号长度为320个样本。使用缺省的时间-半带宽4和DFT长度乘积获得多锥PSD估计。DFT点的默认数量是512。由于信号是实值，PSD估计是片面的，在PSD估计中有512/2+1点。gydF4y2Ba

N = 0:319;X = cos(/4*n)+randn(size(n));PXX = pmtm(x);gydF4y2Ba

绘制多锥PSD估计值。gydF4y2Ba

pmtm (x)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

生成2048个嵌入加法的双通道信号样本gydF4y2BaNgydF4y2Ba(0,1)高斯白噪声。gydF4y2Ba

  • 第一个信道由两个归一化频率为的正弦波组成gydF4y2BaπgydF4y2Ba/ 3,gydF4y2BaπgydF4y2Ba/ 5 rad /样品。第一个正弦信号的振幅是第二个正弦信号的两倍。gydF4y2Ba

  • 第二个信道的归一化频率为gydF4y2BaπgydF4y2Ba/ 4 rad /样品。gydF4y2Ba

使用多锥方法估计信号的PSD在1024-样本区间从0.1开始gydF4y2BaπgydF4y2BaRad /sample到0.4gydF4y2BaπgydF4y2Barad /样品。使用13个等量的正弦锥度。gydF4y2Ba

N = (0:2047)';X = [sin(pi。/［3.5］．*n)*[2 1]' sin(pi/4*n)] + randn(length(n),2); w = linspace(0.1,0.4,1024); ntp = 13; pmtm(x,ntp,“蜡烛”gydF4y2Ba，gydF4y2Ba的正弦gydF4y2Baw *π)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的坐标轴对象包含2个类型为直线的对象。gydF4y2Ba

重复计算，但现在按线性降序对13个圆锥进行加权。您可以将gydF4y2Ba“蜡烛”、“正弦”gydF4y2Ba之后任意位置的名称-值对gydF4y2BaxgydF4y2Ba在函数调用中。gydF4y2Ba

pmtm (x,(国家结核控制规划:1:1)/(1:国家结核控制规划),和w *π,gydF4y2Ba“蜡烛”gydF4y2Ba，gydF4y2Ba的正弦gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的坐标轴对象包含2个类型为直线的对象。gydF4y2Ba

重复计算，但现在使用13个Slepian锥度，并指定7.5的时间-半带宽乘积。gydF4y2Ba

Nw = 7.5;pmtm (x, {nw,国家结核控制规划},w *π)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的坐标轴对象包含2个类型为直线的对象。gydF4y2Ba

重复计算，但现在指定采样率为2 kHz。gydF4y2Ba

Fs = 2e3;pmtm (x, {nw,国家结核控制规划},w * (fs / 2), fs)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的坐标轴对象包含2个类型为直线的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

使用指定的时间-半带宽乘积获得多锥PSD估计。gydF4y2Ba

创建一个角频率为的正弦波gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 加添加剂的Rad /samplegydF4y2BaNgydF4y2Ba(0,1)白噪声。信号长度为320个样本。在时间-半带宽乘积为2.5的情况下获得多锥PSD估计值。分辨率带宽为gydF4y2Ba ［gydF4y2Ba -gydF4y2Ba 2gydF4y2Ba ．gydF4y2Ba 5gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 3.gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba 2gydF4y2Ba ．gydF4y2Ba 5gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 3.gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba ］gydF4y2Ba rad /样品。DFT点的默认数量是512。由于信号是实值，PSD估计是片面的，在PSD估计中有512/2+1点。gydF4y2Ba

N = 0:319;X = cos(/4*n)+randn(size(n));pmtm (2.5 x))gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

获得由角频率为的离散时间正弦信号组成的输入信号的多锥PSD估计gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 加添加剂的Rad /samplegydF4y2BaNgydF4y2Ba(0,1)白噪声。使用与信号长度相等的DFT长度。gydF4y2Ba

创建一个角频率为的正弦波gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba 加添加剂的Rad /samplegydF4y2BaNgydF4y2Ba(0,1)白噪声。信号长度为320个样本。在时间-半带宽乘积为3和DFT长度等于信号长度的情况下，获得多锥PSD估计。由于信号是实值，默认情况下返回单边PSD估计，长度等于320/2+1。gydF4y2Ba

N = 0:319;X = cos(/4*n)+randn(size(n));pmtm (x 3长度(x))gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

获得在1khz采样的信号的多锥PSD估计。信号是一个100赫兹的正弦波加gydF4y2BaNgydF4y2Ba(0,1)白噪声。信号持续时间为2秒。使用时间-半带宽乘积3和DFT长度等于信号长度。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;X = cos(2* *100*t)+randn(size(t));[pxx,f] = pmtm(x,3,length(x)，fs);gydF4y2Ba

绘制多锥PSD估计值。gydF4y2Ba

pmtm (x 3长度(x), fs)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

获得一个多锥PSD估计，其中每个锥形的直接光谱估计在平均值中被赋予相等的权重。gydF4y2Ba

获得在1khz采样的信号的多锥PSD估计。信号是一个100赫兹的正弦波加gydF4y2BaNgydF4y2Ba(0,1)白噪声。信号持续时间为2秒。使用时间-半带宽乘积3和DFT长度等于信号长度。使用gydF4y2Ba“团结”gydF4y2Ba选择在平均中给予每个单独的锥形直接光谱估计相等的权重。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;X = cos(2* *100*t)+randn(size(t));[pxx,f] = pmtm(x,3,length(x)，fs，gydF4y2Ba“团结”gydF4y2Ba）;gydF4y2Ba

绘制多锥PSD估计值。gydF4y2Ba

pmtm (x 3长度(x), fs,gydF4y2Ba“团结”gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

这个例子检查了DPSS序列的频域浓度。该示例通过预计算Slepian序列并只选择那些超过99%的能量集中在分辨率带宽上的序列，产生输入信号的多锥PSD估计。gydF4y2Ba

信号是一个100赫兹的正弦波加gydF4y2BaNgydF4y2Ba(0,1)白噪声。信号持续时间为2秒。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;X = cos(2* *100*t)+randn(size(t));gydF4y2Ba

设置time-halfbandwidth积为3.5。对于2000个样本的信号长度和0.001秒的采样间隔，这将导致分辨率带宽为[-1.75,1.75]Hz。计算前10个Slepian序列，并检查它们在指定分辨率带宽内的频率集中。gydF4y2Ba

[e,v] = dpss(长度(x)，3.5,10);Lv =长度(v);阀杆(1:lv, v,gydF4y2Ba“填充”gydF4y2Baylim([0 1.2])gydF4y2Ba第k个Slepian序列[-w,w]的能量比例gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为Proportion of Energy in [-w,w] of k- Slepian Sequence的坐标轴对象包含一个stem类型的对象。gydF4y2Ba

确定能量浓度大于99%的Slepian序列的数量。利用选定的DPSS序列，得到多锥PSD估计。集gydF4y2Ba“DropLastTaper”gydF4y2Ba来gydF4y2Ba假gydF4y2Ba要使用所有选定的锥体。gydF4y2Ba

持有gydF4y2Ba在gydF4y2BaPlot ([1 lv]，0.99*[1 1]) hold住gydF4y2Ba从gydF4y2Ba

图中包含一个轴对象。k-th Slepian Sequence中标题为“Proportion of Energy in [-w,w]”的坐标轴对象包含2个类型为stem、line的对象。gydF4y2Ba

Idx = find(v>0.99,1，gydF4y2Ba“最后一次”gydF4y2Ba）gydF4y2Ba
Idx = 5gydF4y2Ba
[pxx f] = pmtm (x, e (:, 1: idx), v (1: idx),长度(x), fs,gydF4y2Ba“DropLastTaper”gydF4y2Ba、假);gydF4y2Ba

绘制多锥PSD估计值。gydF4y2Ba

pmtm (x, e (: 1: idx), v (1: idx),长度(x), fs,gydF4y2Ba“DropLastTaper”gydF4y2Ba假)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

获得100 Hz正弦波的多锥PSD估计gydF4y2BaNgydF4y2Ba(0,1)噪音。数据以1千赫采样。使用gydF4y2Ba“中心”gydF4y2Ba选项获取以dc为中心的PSD。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;X = cos(2* *100*t)+randn(size(t));[pxx,f] = pmtm(x,3.5,length(x)，fs，gydF4y2Ba“中心”gydF4y2Ba）;gydF4y2Ba

绘制以dc为中心的PSD估计。gydF4y2Ba

pmtm (x, 3.5,长度(x), fs,gydF4y2Ba“中心”gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为“功率谱密度”的axis对象包含一个类型为line的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

下面的例子说明了多锥PSD估计的置信界限的使用。虽然不是统计显著性的必要条件，但在多锥PSD估计中，当下置信界超过周围PSD估计的上置信界时，频率清楚地表明时间序列中存在显著的振荡。gydF4y2Ba

创建一个由100-Hz和150-Hz正弦波叠加组成的信号gydF4y2BaNgydF4y2Ba(0,1)噪音。两个正弦波的振幅是1。采样频率为1khz。信号持续时间为2秒。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;x = cos(2 *π* 100 * t) + cos(2 *π* 150 * t) + randn(大小(t));gydF4y2Ba

获得95%置信限的多锥PSD估计。绘制PSD估计值和置信区间，并放大100和150 Hz附近感兴趣的频率区域。gydF4y2Ba

[pxx,f,pxxc] = pmtm(x,3.5,length(x)，fs，gydF4y2Ba“ConfidenceLevel”gydF4y2Ba, 0.95);情节(f, 10 * log10 (pxx))gydF4y2Ba在gydF4y2Ba情节(f, 10 * log10 (pxxc),gydF4y2Bar -。gydF4y2Baxlim([85 175]) xlabel(gydF4y2Ba“赫兹”gydF4y2Ba) ylabel (gydF4y2Ba“数据库”gydF4y2Ba)标题(gydF4y2Ba95%置信限下的多锥PSD估计gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为Multitaper PSD Estimate with 95% confidence Bounds的axis对象包含3个类型为line的对象。gydF4y2Ba

在100和150 Hz附近的置信下限明显高于100和150 Hz附近以外的置信上限。gydF4y2Ba

打开实时脚本gydF4y2Ba

生成1024个由三个加性正弦信号组成的多通道信号样本gydF4y2Ba NgydF4y2Ba （gydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba 1gydF4y2Ba ）gydF4y2Ba 高斯白噪声。正弦信号的频率是gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ，gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 3.gydF4y2Ba ,gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba rad /样品。用汤姆逊多锥法估计信号的PSD并绘制图。gydF4y2Ba

N = 1024;n = 0: n -1;W = pi./[2;3;4];X = cos(w*n)' + randn(length(n)，3);pmtm (x)gydF4y2Ba

图中包含一个轴对象。标题为Thomson Multitaper功率谱密度估计的坐标轴对象包含3个类型为直线的对象。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

输入信号，指定为行或列向量或矩阵。如果gydF4y2BaxgydF4y2Ba是矩阵，则其列被视为独立的通道。gydF4y2Ba

例子:gydF4y2Ba因为(π/ 4 * (0:159))+ randn (1160)gydF4y2Ba是单通道行向量信号。gydF4y2Ba

例子:gydF4y2Ba因为(pi. / (4; 2) * (0:159)) ' + randn (160 2)gydF4y2Ba是一个双通道信号。gydF4y2Ba

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

锥度类型，指定为gydF4y2Ba“斯莱皮恩”gydF4y2Ba或gydF4y2Ba的正弦gydF4y2Ba．gydF4y2Ba

您可以指定gydF4y2Ba“蜡烛”gydF4y2Ba，gydF4y2BatapertypegydF4y2Ba之后任意位置的名称-值对gydF4y2BaxgydF4y2Ba在函数调用中。gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

时-半带宽积，指定为正标量。gydF4y2BapmtmgydF4y2Ba使用2 ×gydF4y2Ba西北gydF4y2BaPSD估计有1个Slepian锥度。典型的选择gydF4y2Ba西北gydF4y2Ba是gydF4y2Ba2gydF4y2Ba，gydF4y2Ba5/2gydF4y2Ba，gydF4y2Ba3.gydF4y2Ba，gydF4y2Ba7/2gydF4y2Ba,或gydF4y2Ba4gydF4y2Ba．gydF4y2Ba

在多锥谱估计中，用户可以指定多锥谱估计的分辨率带宽gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba在哪里gydF4y2BaWgydF4y2Ba＝gydF4y2BakgydF4y2Ba/gydF4y2BaNgydF4y2BaΔgydF4y2BatgydF4y2Ba对于一些小的gydF4y2BakgydF4y2Ba> 1gydF4y2Ba．同样,gydF4y2BaWgydF4y2Ba是DFT频率分辨率的一个小倍数。时间-半带宽乘积是分辨率-半带宽和输入信号中采样数的乘积，gydF4y2BaNgydF4y2Ba．其傅里叶变换集中的斯勒皮氏锥的数目gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba(特征值接近单位)为gydF4y2Ba2gydF4y2BaNgydF4y2BaWgydF4y2Ba- 1gydF4y2Ba．gydF4y2Ba

正弦锥度数或平均权重，指定为整数标量或矢量。gydF4y2Ba

  • 如果gydF4y2Ba米gydF4y2Ba为标量，它表示在计算PSD估计时用作数据窗口的正弦锥形数。正弦锥是均匀加权的。gydF4y2Ba

  • 如果gydF4y2Ba米gydF4y2Ba是一个向量，它表示在计算PSD估计时用于平均正弦锥度的权重。的长度gydF4y2Ba米gydF4y2Ba指示要使用的锥度数量。的要素gydF4y2Ba米gydF4y2Ba必须加1。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba

DFT点的数目，指定为正整数。对于一个实值输入信号，gydF4y2BaxgydF4y2Ba， PSD估计，gydF4y2BapxxgydF4y2Ba有长度(gydF4y2BanfftgydF4y2Ba/2 + 1) ifgydF4y2BanfftgydF4y2Ba是偶数，并且(gydF4y2BanfftgydF4y2Ba+ 1)/2如果gydF4y2BanfftgydF4y2Ba是奇数。对于复值输入信号，gydF4y2BaxgydF4y2Ba时，PSD估计总有长度gydF4y2BanfftgydF4y2Ba．如果gydF4y2BanfftgydF4y2Ba为空，默认gydF4y2BanfftgydF4y2Ba使用。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba

抽样率，指定为正标量。采样率是单位时间内采样的数量。如果时间的单位是秒，那么采样率的单位是Hz。gydF4y2Ba

规范化的频率，指定为至少有两个元素的行或列向量。归一化频率以弧度/样本为单位。gydF4y2Ba

例子:gydF4y2BaW = [pi/4 pi/2]gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

频率，指定为包含至少两个元素的行向量或列向量。频率是单位时间内的周期。单位时间由采样率指定，gydF4y2BafsgydF4y2Ba．如果gydF4y2BafsgydF4y2Ba有样本单位/秒吗gydF4y2BafgydF4y2Ba单位是Hz。gydF4y2Ba

例子:gydF4y2BaFs = 1000;F = [100 200]gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

标志，指示是删除还是保留作为逻辑指定的最后一个DPSS序列。默认为gydF4y2Ba真正的gydF4y2Ba而且gydF4y2BapmtmgydF4y2Ba落下最后一个锥度。在多重缩减估计中，前2gydF4y2Ba西北gydF4y2Ba- 1 DPSS序列的特征值接近统一。如果你用的少于2gydF4y2Ba西北gydF4y2Ba- 1个序列，很可能所有的圆锥都有接近1的特征值，你可以指定gydF4y2BadropflaggydF4y2Ba作为gydF4y2Ba假gydF4y2Ba保留最后一根。gydF4y2Ba

单个锥形PSD估计值的权重，指定为之一gydF4y2Ba“适应”gydF4y2Ba，gydF4y2Ba“特征”gydF4y2Ba,或gydF4y2Ba“团结”gydF4y2Ba．默认值是汤姆森自适应频率依赖权重，gydF4y2Ba“适应”gydF4y2Ba．这些权重的计算详见第368-370页gydF4y2Ba[２]gydF4y2Ba．的gydF4y2Ba“特征”gydF4y2Ba方法通过对应斯莱皮安锥度的特征值(频率浓度)对每个锥度PSD估计进行加权。的gydF4y2Ba“团结”gydF4y2Ba方法对每个锥形PSD估计值的权重相等。gydF4y2Ba

DPSS (Slepian)序列，指定为矩阵。如果gydF4y2BaxgydF4y2Ba长度gydF4y2BaNgydF4y2Ba,然后gydF4y2BaegydF4y2Ba有gydF4y2BaNgydF4y2Ba行。矩阵gydF4y2BaegydF4y2Ba的输出。gydF4y2Ba离散长gydF4y2Ba．gydF4y2Ba

DPSS (Slepian)序列的特征值，指定为列向量。DPSS序列的特征值表示序列能量集中在分辨率带宽中的比例，gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba．特征值的范围在区间内gydF4y2Ba(0,1)gydF4y2Ba一般来说，第一个gydF4y2Ba2gydF4y2BaNgydF4y2BaWgydF4y2Ba- 1gydF4y2Ba特征值趋近于1，然后递减到0。向量gydF4y2BavgydF4y2Ba的输出。gydF4y2Ba离散长gydF4y2Ba．gydF4y2Ba

的输入参数gydF4y2Ba离散长gydF4y2Ba，指定为单元格数组。的第一个输入参数gydF4y2Ba离散长gydF4y2Ba是DPSS序列的长度，从gydF4y2Badpss_paramsgydF4y2Ba因为它是由的长度得到的gydF4y2BaxgydF4y2Ba．gydF4y2Ba

例子:gydF4y2Bapmtm (randn (1000 1), {2.5, 3})gydF4y2Ba利用前3个时间-半带宽乘积为2.5的Slepian序列计算随机序列的PSD。gydF4y2Ba

PSD估计的频率范围，指定为其中之一gydF4y2Ba“单向的”gydF4y2Ba，gydF4y2Ba双侧的gydF4y2Ba,或gydF4y2Ba“中心”gydF4y2Ba．默认为gydF4y2Ba“单向的”gydF4y2Ba对于实值信号和gydF4y2Ba双侧的gydF4y2Ba对于复值信号。每个选项对应的频率范围为gydF4y2Ba

  • “单向的”gydF4y2Ba-返回一个实值输入信号的单边PSD估计，gydF4y2BaxgydF4y2Ba．如果gydF4y2BanfftgydF4y2Ba是偶数,gydF4y2BapxxgydF4y2Ba长度gydF4y2BanfftgydF4y2Ba/2 + 1，在区间内计算gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba］gydF4y2Barad /样品。如果gydF4y2BanfftgydF4y2Ba长度是奇数gydF4y2BapxxgydF4y2Ba是(gydF4y2BanfftgydF4y2Ba+ 1)/2，区间是gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba）gydF4y2Barad /样品。当gydF4y2BafsgydF4y2Ba，则对应的间隔为[0，gydF4y2BafsgydF4y2Ba/2]周期/单位时间和[0，gydF4y2BafsgydF4y2Ba/2)周期/偶数和奇数长度的单位时间gydF4y2BanfftgydF4y2Ba分别。gydF4y2Ba

    该函数将除0和奈奎斯特频率外的所有频率的功率乘以2以保持总功率。gydF4y2Ba

  • 双侧的gydF4y2Ba-返回实值或复值输入的双侧PSD估计，gydF4y2BaxgydF4y2Ba．在这种情况下，gydF4y2BapxxgydF4y2Ba长度gydF4y2BanfftgydF4y2Ba并且在区间内计算gydF4y2Ba[0, 2gydF4y2BaπgydF4y2Ba）gydF4y2Barad /样品。当gydF4y2BafsgydF4y2Ba，则间隔为[0，gydF4y2BafsgydF4y2Ba)周期/单位时间。gydF4y2Ba

  • “中心”gydF4y2Ba-返回实值或复值输入的居中双面PSD估计，gydF4y2BaxgydF4y2Ba．在这种情况下，gydF4y2BapxxgydF4y2Ba长度gydF4y2BanfftgydF4y2Ba并且在区间内计算gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Ba弧度/样品为均匀长度gydF4y2BanfftgydF4y2Ba而且gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba）gydF4y2Ba奇数长度的Rad /samplegydF4y2BanfftgydF4y2Ba．当gydF4y2BafsgydF4y2Ba，则对应的间隔为(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2]周期/单位时间和(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2)周期/偶数和奇数长度的单位时间gydF4y2BanfftgydF4y2Ba分别。gydF4y2Ba

真实PSD的覆盖概率，指定为范围(0,1)中的标量。输出,gydF4y2BapxxcgydF4y2Ba对象的下界和上界gydF4y2Ba概率gydF4y2Ba× 100%区间估计为真实PSD。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

PSD估计，返回为实值，非负列向量或矩阵。的每一列gydF4y2BapxxgydF4y2Ba对应列的PSD估计值是gydF4y2BaxgydF4y2Ba．PSD估计的单位是时间序列数据每单位频率的平方幅度单位。例如，如果输入数据以伏特为单位，则PSD估计的单位是每单位频率的平方伏特。对于以伏特为单位的时间序列，如果假设电阻为1 Ω，并指定以赫兹为单位的采样率，则PSD估计的单位为瓦特/赫兹。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba

归一化频率，作为实值列向量返回。如果gydF4y2BapxxgydF4y2Ba是片面的PSD估计，gydF4y2BawgydF4y2Ba跨越区间gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba］gydF4y2Ba如果gydF4y2BanfftgydF4y2Ba是偶数和gydF4y2Ba[0,gydF4y2BaπgydF4y2Ba）gydF4y2Ba如果gydF4y2BanfftgydF4y2Ba是奇数。如果gydF4y2BapxxgydF4y2Ba为双侧PSD估计，gydF4y2BawgydF4y2Ba跨越区间gydF4y2Ba[0, 2gydF4y2BaπgydF4y2Ba）gydF4y2Ba．对于以dc为中心的PSD估计，gydF4y2BawgydF4y2Ba跨越区间gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba］gydF4y2Ba甚至gydF4y2BanfftgydF4y2Ba而且gydF4y2Ba(-gydF4y2BaπgydF4y2Ba，gydF4y2BaπgydF4y2Ba）gydF4y2Ba为奇数gydF4y2BanfftgydF4y2Ba．gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

循环频率，作为实值列向量返回。对于片面的PSD估计，gydF4y2BafgydF4y2Ba跨度为区间[0，gydF4y2BafsgydF4y2Ba/ 2)当gydF4y2BanfftgydF4y2Ba是偶数，[0，gydF4y2BafsgydF4y2Ba/ 2)gydF4y2BanfftgydF4y2Ba是奇数。对于双面PSD估计，gydF4y2BafgydF4y2Ba跨度为区间[0，gydF4y2BafsgydF4y2Ba)．对于以dc为中心的PSD估计，gydF4y2BafgydF4y2Ba跨越区间(-)gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2]周期/单位时间为偶数长度gydF4y2BanfftgydF4y2Ba和(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2)周期/奇数长度的单位时间gydF4y2BanfftgydF4y2Ba．gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

置信边界，返回为具有实值元素的矩阵。矩阵的行长度等于PSD估计的长度，gydF4y2BapxxgydF4y2Ba．gydF4y2BapxxcgydF4y2Ba栏数是?的两倍gydF4y2BapxxgydF4y2Ba．奇数列包含置信区间的下界，偶数列包含上界。因此,gydF4y2Bapxxc (m, 2 * n - 1)gydF4y2Ba下置信界是和吗gydF4y2Bapxxc (m, 2 * n)gydF4y2Ba置信上限是否与估计相对应gydF4y2Bapxx (m, n)gydF4y2Ba．的值决定置信区间的覆盖概率gydF4y2Ba概率gydF4y2Ba输入。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

汤姆逊多锥谱估计gydF4y2Ba

周期图不是广义平稳过程真功率谱密度(PSD)的一致估计量。为了减少周期图中的变异性，从而产生PSD的一致估计，多锥方法对使用相互正交的窗口或或族获得的修改的周期图进行平均gydF4y2Ba蜡烛gydF4y2Ba．除了相互正交外，锥形还具有最佳的时频集中特性。锥面正交性和时频集中对多锥技术的成功与否至关重要。看到gydF4y2Ba离散长球状(Slepian)序列gydF4y2Ba以了解汤姆逊多锥方法中所使用的斯莱比安序列的简要描述。gydF4y2Ba

多锥法使用gydF4y2BaKgydF4y2Ba修改的周期图，每个周期图都使用不同的Slepian序列作为窗口。让gydF4y2Ba

年代gydF4y2Ba kgydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ΔgydF4y2Ba tgydF4y2Ba |gydF4y2Ba ∑gydF4y2Ba ngydF4y2Ba ＝gydF4y2Ba 0gydF4y2Ba NgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ggydF4y2Ba kgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba ngydF4y2Ba ΔgydF4y2Ba tgydF4y2Ba |gydF4y2Ba 2gydF4y2Ba

方法得到的修正周期图gydF4y2BakgydF4y2BaSlepian序列，gydF4y2BaggydF4y2BakgydF4y2Ba（gydF4y2BangydF4y2Ba)．在其最简单的形式，多锥方法只是平均gydF4y2BaKgydF4y2Ba修正周期图以产生多锥PSD估计:gydF4y2Ba

年代gydF4y2Ba （gydF4y2Ba 太gydF4y2Ba ）gydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba 1gydF4y2Ba KgydF4y2Ba ∑gydF4y2Ba kgydF4y2Ba ＝gydF4y2Ba 0gydF4y2Ba KgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba 年代gydF4y2Ba kgydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ．gydF4y2Ba

汤姆逊的多锥度方法，在gydF4y2Ba[4]gydF4y2Ba，类似于韦尔奇的重叠段平均方法，在这两个平均对PSD的近似不相关估计。然而，这两种方法在如何产生这些不相关的PSD估计方面有所不同。多锥方法在每个修正的周期图中使用整个信号。Slepian锥体的正交性揭示了不同的修正周期图。韦尔奇的方法使用每个修改的周期图中的信号片段，并且分割去关联不同的修改周期图。gydF4y2Ba

的方程gydF4y2Ba年代gydF4y2Ba(太)gydF4y2Ba（gydF4y2BafgydF4y2Ba）gydF4y2Ba对应于gydF4y2Ba“团结”gydF4y2Ba选项gydF4y2BapmtmgydF4y2Ba．然而，正如在gydF4y2Ba离散长球状(Slepian)序列gydF4y2Ba时，Slepian序列在感兴趣的频带内不具有相等的能量集中。序次越高，序能量在带内的集中程度越低gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba浓度由特征值给出。因此，使用特征值来加权是有益的gydF4y2BaKgydF4y2Ba平均前修改的周期图。这对应于gydF4y2Ba“特征”gydF4y2Ba选项gydF4y2BapmtmgydF4y2Ba．gydF4y2Ba

利用序列特征值产生修正周期图的加权平均，说明了Slepian序列的频率集中特性。但是，它没有考虑随机过程的功率谱密度与Slepian序列的频率浓度之间的相互作用。具体地说，在使用高阶Slepian序列的修正周期图中，随机过程功率较小的频率区域估计的可靠性较差。这证明了一个频率依赖的自适应过程，它不仅解释了Slepian序列的频率集中，而且还解释了时间序列中的功率分布。这种自适应加权对应于gydF4y2Ba“适应”gydF4y2Ba选项gydF4y2BapmtmgydF4y2Ba和是计算多锥估计的默认值。gydF4y2Ba

离散长球状(Slepian)序列gydF4y2Ba

Slepian序列的推导是从离散/连续频率集中问题出发的。对所有gydF4y2BaℓgydF4y2Ba2gydF4y2Ba序列索引受限到gydF4y2Ba0,1，…，gydF4y2BaNgydF4y2Ba- 1gydF4y2Ba时，该问题寻求在一个频带内具有最大能量集中的序列gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba与gydF4y2Ba|gydF4y2BaWgydF4y2Ba| < 1/2ΔgydF4y2BatgydF4y2Ba．gydF4y2Ba

这相当于找到an的特征值和对应的特征向量gydF4y2BaNgydF4y2Ba——- - - - - -gydF4y2BaNgydF4y2Ba自伴半正定算子。因此，特征值为实且非负，不同特征值对应的特征向量相互正交。在这个特殊的问题中，特征值以1为界，特征值是序列在频率区间内能量集中的度量gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba．gydF4y2Ba

特征值问题由gydF4y2Ba

∑gydF4y2Ba 米gydF4y2Ba ＝gydF4y2Ba 0gydF4y2Ba NgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba 罪gydF4y2Ba （gydF4y2Ba 2gydF4y2Ba πgydF4y2Ba WgydF4y2Ba （gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba ）gydF4y2Ba ）gydF4y2Ba πgydF4y2Ba （gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba ）gydF4y2Ba ggydF4y2Ba kgydF4y2Ba （gydF4y2Ba 米gydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba λgydF4y2Ba kgydF4y2Ba （gydF4y2Ba NgydF4y2Ba ，gydF4y2Ba WgydF4y2Ba ）gydF4y2Ba ggydF4y2Ba kgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba ，gydF4y2Ba ngydF4y2Ba ，gydF4y2Ba kgydF4y2Ba ＝gydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba 1gydF4y2Ba ，gydF4y2Ba 2gydF4y2Ba ，gydF4y2Ba .．.gydF4y2Ba ，gydF4y2Ba NgydF4y2Ba −gydF4y2Ba 1.gydF4y2Ba

0阶DPSS序列，gydF4y2BaggydF4y2Ba0gydF4y2Ba，是最大特征值对应的特征向量。一阶DPSS序列，gydF4y2BaggydF4y2Ba1gydF4y2Ba，为第二大特征值对应的特征向量，正交于零阶序列。二阶DPSS序列，gydF4y2BaggydF4y2Ba2gydF4y2Ba，为第三特征值对应的特征向量，与两个低阶DPSS序列正交。因为运算符是gydF4y2BaNgydF4y2Ba——- - - - - -gydF4y2BaNgydF4y2Ba，有gydF4y2BaNgydF4y2Ba特征向量。然而，对于给定的序列长度gydF4y2BaNgydF4y2Ba和指定的带宽gydF4y2Ba(- - - - - -gydF4y2BaWgydF4y2Ba，gydF4y2BaWgydF4y2Ba］gydF4y2Ba，大约有gydF4y2Ba2gydF4y2Ba西北gydF4y2Ba- 1gydF4y2Ba特征值非常接近单位的DPSS序列。使用gydF4y2Ba西北gydF4y2Ba指定gydF4y2Ba西北gydF4y2Ba．gydF4y2Ba

正弦蜡烛gydF4y2Ba

中提出的一种Slepian序列的替代方案——正弦锥gydF4y2Ba［3］gydF4y2Ba，定义为gydF4y2Ba

ggydF4y2Ba kgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba 2gydF4y2Ba NgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba 罪gydF4y2Ba πgydF4y2Ba kgydF4y2Ba ngydF4y2Ba NgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba ，gydF4y2Ba ngydF4y2Ba ，gydF4y2Ba kgydF4y2Ba ＝gydF4y2Ba 1gydF4y2Ba ，gydF4y2Ba 2gydF4y2Ba ，gydF4y2Ba .．.gydF4y2Ba ，gydF4y2Ba NgydF4y2Ba ．gydF4y2Ba

与Slepian序列不同，正弦锥可以直接计算，不需要建立和求解特征值方程。这使得正弦函数的计算速度更快。正弦锥具有接近于Slepian序列的光谱浓度，但不需要额外的参数来指定光谱带宽。使用正弦锥度计算的PSD估计带宽可以通过改变使用的锥度数量在本地进行调整gydF4y2Ba米gydF4y2Ba．gydF4y2Ba

比较Slepian和sin TapersgydF4y2Ba

生成前五个对应于时间-半带宽乘积3的斯莱皮安锥。指定锥度长度为1000。gydF4y2Ba

N = 1000;Nw = 3;Ns = 2*(nw)-1;tprs = dpss(N,nw,ns);磅=gydF4y2Ba“斯莱皮恩”gydF4y2Ba；gydF4y2Ba

生成前五个正弦锥度。gydF4y2Ba

n = 1: n;K = 1:ns;tprs(:，:，2) =√(2/(N+1))*sin(pi* N '*k/(N+1));磅(2)=gydF4y2Ba“正弦”gydF4y2Ba；gydF4y2Ba

画出两组锥形。gydF4y2Ba

为gydF4y2BaKj = 1:2 subplot(2,1, Kj) plot(tprs(:，:， Kj)) title(lbs(Kj)) legend(append(gydF4y2Ba'k = 'gydF4y2Ba字符串(k + kj-2)),gydF4y2Ba.．.gydF4y2Ba“定位”gydF4y2Ba，gydF4y2Ba“水平”gydF4y2Ba，gydF4y2Ba“位置”gydF4y2Ba，gydF4y2Ba“南”gydF4y2Ba)传说(gydF4y2Ba“boxoff”gydF4y2Ba) ylim([-0.09 0.07])gydF4y2Ba结束gydF4y2Ba

图中包含2个轴对象。标题为Slepian的Axes对象1包含5个类型为line的对象。这些对象分别代表k = 0, k = 1, k = 2, k = 3, k = 4。标题为sin的坐标轴对象2包含5个类型为line的对象。这些物体分别代表k = 1, k = 2, k = 3, k = 4, k = 5。gydF4y2Ba

参考文献gydF4y2Ba

[1] McCoy, Emma J.， Andrew T. Walden和Donald B. Percival。幂律过程的多锥谱估计gydF4y2BaIEEEgydF4y2Ba®gydF4y2Ba信号处理汇刊gydF4y2Ba46岁的没有。3(1998年3月):655-68。gydF4y2Bahttps://doi.org/10.1109/78.661333gydF4y2Ba．gydF4y2Ba

[2]珀西瓦尔，唐纳德·B和安德鲁·t·沃顿。gydF4y2Ba物理应用的光谱分析:多锥和常规单变量技术gydF4y2Ba．剑桥;美国纽约:剑桥大学出版社，1993。gydF4y2Ba

里德尔，库尔特·S，亚历山大·西多伦科。最小偏差多重锥度谱估计。gydF4y2BaIEEE信号处理汇刊gydF4y2Ba43岁的没有。1(1995年1月):188-95。gydF4y2Bahttps://doi.org/10.1109/78.365298gydF4y2Ba．gydF4y2Ba

[4]汤姆森，大卫J。“频谱估计和谐波分析。”gydF4y2BaIEEE论文集gydF4y2Ba70年,没有。9(1982): 1055-96。gydF4y2Bahttps://doi.org/10.1109/PROC.1982.12433gydF4y2Ba．gydF4y2Ba

版本历史gydF4y2Ba

R2006a之前介绍gydF4y2Ba

另请参阅gydF4y2Ba

|gydF4y2Ba|gydF4y2Ba

主题gydF4y2Ba