PCOV.

自回归功率谱密度估计协方差法

折叠所有页面

语法

pxx = pcov (x,顺序)

PXX = PCOV（x，订单，nfft）

[pxx，w] = pcov（＿＿＿）

[pxx，f] = pcov（＿＿＿，fs）

[PXX，W] = PCOV（x，订单，w）

[PXX，F] = PCOV（X，Order，F，FS）

［＿＿＿] = pcov (x,秩序,＿＿＿freqrange)

［＿＿＿pxxc] = pcov (＿＿＿、“ConfidenceLevel”概率)

pcov (＿＿＿）

描述

PXX.= PCOV（x，命令）返回功率谱密度(PSD)估计值，PXX.，离散时间信号的x，用协方差法得到。当x是一个矢量，它被视为单个通道。当x是一个矩阵，PSD独立计算每一列，并存储在PXX.．PXX.为单位频率功率分布。频率的单位为rad/sample。命令是用于产生PSD估计的自回归（AR）模型的顺序。

PXX.= PCOV（x，命令，nfft）使用nfft离散傅里叶变换(DFT)中的点。真实的x，PXX.长度(nfft/ 2 + 1）如果nfft为偶数，且(nfft+1）/ 2如果nfft是奇数。对于复杂的价值x，PXX.总长度nfft．如果你省略nfft，或将其指定为空，然后PCOV.使用默认的DFT长度256。

［PXX.，w] = PCOV（＿＿＿）返回归一化角频率的矢量，w，此时估计PSD。w单位为弧度/样本。为实值信号,w跨越间隔[0，π］什么时候nfft甚至和[0，π）什么时候nfft是奇数。为复值信号,w总是跨越这个区间[0,2π］．

例子

［PXX.，f] = PCOV（＿＿＿，FS.）返回频率向量，f，在每单位时间的周期中。采样频率,FS.，是每单位时间的样本数。如果时间单位是秒，那么f处于循环/秒（Hz）。为实值信号,f跨越区间[0，FS./ 2]什么时候nfft为偶数，[0，FS./ 2）何时nfft是奇数。为复值信号,f跨越区间[0，FS.）。

［PXX.，w] = PCOV（x，命令，w）返回在向量中指定的归一化频率下的双边AR PSD估计，w．矢量w必须包含至少两个元素，因为否则该功能将其解释为nfft．

［PXX.，f] = PCOV（x，命令，f，FS.）在向量中指定的频率下返回双面AR PSD估计，f．矢量f必须包含至少两个元素，因为否则该功能将其解释为nfft．频率f是单位时间内的周期。采样频率,FS.，是每单位时间的样本数。如果时间单位是秒，那么f处于循环/秒（Hz）。

［＿＿＿] = PCOV（x，命令，＿＿＿，freqrange）在指定的频率范围内返回AR PSD估计值freqrange．有效选项freqrange是:“单向的”，'twosiding',或“中心”．

［＿＿＿，pxxc] = PCOV（＿＿＿“ConfidenceLevel”,可能性）返回可能性×中PSD估计的100%置信区间pxxc．

例子

pcov (＿＿＿）在没有输出参数的情况下，在当前的图形窗口中以每单位频率的dB为单位绘制AR PSD估计图。

例子

全部收缩

协方差 - 方法PSD估计AR（4）过程

打开直播脚本

创建一个AR(4)广义平稳随机过程的实现。使用协方差方法估计PSD。将基于单一实现的PSD估计与随机过程的真实PSD进行比较。

创建AR（4）系统功能。获得频率响应并绘制系统的PSD。

a = [1 -2.7607 3.8106 -2.6535 0.9238];[h，f] = freqz（1，a，[]，1）;绘图（F，20 * log10（ABS（H）））XLabel（'频率（Hz）'）ylabel（“PSD (dB / Hz)”）

创建AR（4）随机过程的实现。将随机数生成器设置为可重复结果的默认设置。实现长度为1000个样本。假设1 Hz的采样频率。采用PCOV.来估计四阶过程的PSD。比较PSD估计与真实PSD。

RNG.默认的X = RANDN（1000,1）;y =过滤器（1，A，x）;[PXX，F] = PCOV（Y，4,1024,1）;抓住在情节(F, 10 * log10 (Pxx))传说(“真实功率谱密度”，“pcov PSD估计”）

多通道信号的协方差法PSD估计

打开直播脚本

创建由添加剂中的三个正弦曲线组成的多通道信号 $N （ 0 ， 1 ）$ 白色高斯噪音。正弦频率为100 Hz，200 Hz和300 Hz。采样频率为1 kHz，信号的持续时间为1秒。

FS = 1000;T = 0：1 / FS：1-1 / FS;f = [100; 200; 300];x = cos（2 * pi * f * t）'+ Randn（长度（t），3）;

估计信号的PSD使用协方差方法与12阶自回归模型。使用默认的DFT长度。情节的估计。

morder = 12;PCOV（x，morder，[]，fs）

输入参数

全部收缩

`x`- - - - - -输入信号
向量|矩阵

输入信号，指定为行或列向量或矩阵。如果x是矩阵，然后其列被视为独立信道。

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

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

数据类型：单|双倍的
复数支持：金宝app是的

`命令`- - - - - -自回归模型的顺序
正整数

自回归模型的顺序，指定为正整数。

数据类型：双倍的

`nfft`- - - - - -DFT点的个数
`256`(默认)|整数|`[]`

DFT点的数目，指定为正整数。对于实值输入信号，x，PSD估计，PXX.长度(nfft/ 2 + 1）如果nfft为偶数，且(nfft+1）/ 2如果nfft是奇数。对于复值输入信号，x，PSD估计总是长度nfft．如果nfft指定为空，则默认nfft用来。

数据类型：单|双倍的

`FS.`- - - - - -采样率
积极的标量

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

`w`- - - - - -归一化频率
向量

归一化频率，指定为具有至少两个元素的行或列向量。归一化频率以rad/sample为单位。

例子：w = [pi / 4 pi / 2]

数据类型：双倍的

`f`- - - - - -频率
向量

频率，指定为具有至少两个元素的行或列向量。频率以单位时间的周期为单位。单位时间由采样率确定，FS.．如果FS.有一个样本/秒，那么f有Hz的单位。

例子：fs = 1000;F = [100 200]

数据类型：双倍的

`freqrange`- - - - - -PSD估计的频率范围
`“单向的”`|`'twosiding'`|`“中心”`

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

“单向的”-返回实值输入信号的单边PSD估计，x．如果nfft甚至，PXX.长度nfft/ 2 + 1并计算在间隔内[0，π］Rad /样品。如果nfft是奇怪的，长度PXX.是(nfft+ 1)/2，区间是[0，π）Rad /样品。当FS.时，对应的间隔为[0，FS./2]周期/单位时间和[0，FS./2)偶数和奇数长度的周期/单位时间nfft分别。
'twosiding'- 返回实值或复数输入的双面PSD估计，x．在这种情况下，PXX.长度nfft并计算在间隔内[0,2π）Rad /样品。当FS.时，间隔为[0，FS.)周期/单位时间。
“中心”- 返回居住的双面PSD估计值，用于实值或复数输入的输入，x．在这种情况下，PXX.长度nfft并计算在间隔内（ -π，π］rad /样品均匀长度nfft和（ -π，π）Rad /样品奇数长度nfft．当FS.可选择指定，相应的间隔是（ -FS./ 2,FS./ 2]循环/单位时间和（ -FS./ 2,FS./2)偶数和奇数长度的周期/单位时间nfft分别。

`可能性`- - - - - -PSD估计的置信区间
0.95(默认)|量程(0,1)内的标量

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

输出参数

全部收缩

`PXX.`- PSD估计
向量|矩阵

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

数据类型：单|双倍的

`w`- 标准化频率
向量

归一化频率，作为实值列向量返回。如果PXX.是单面PSD估计，w跨越间隔[0，π］如果nfft甚至和[0，π）如果nfft是奇数。如果PXX.是一个双边PSD估计，w跨越间隔[0,2π）．对于以CC中心的PSD估计，w跨越间隔（ -π，π］甚至nfft和（ -π，π）为奇数nfft．

数据类型：双倍的

`f`——循环的频率
向量

循环频率，作为实值列向量返回。对于片面的PSD估计，f跨越区间[0，FS./ 2]什么时候nfft为偶数，[0，FS./ 2）何时nfft是奇数。对于双边PSD估计，f跨越区间[0，FS.）。对于以CC中心的PSD估计，f跨越间隔(-FS./ 2,FS./ 2]甚至长度的循环/单位时间nfft和(-FS./ 2,FS./2)奇数长度的周期/单位时间nfft．

数据类型：双倍的|单

`pxxc`- 信心界
矩阵

置信界限，返回为具有实值元素的矩阵。矩阵的行大小等于PSD估计的长度，PXX.．pxxc有两倍的列数PXX.．奇数列包含置信区间的下界，偶数列包含置信区间的上界。因此,pxxc (m, 2 * n - 1)是较低的置信度pxxc (m, 2 * n)是与估计相对应的上置信度pxx (m, n)．置信区间的覆盖概率由的值决定可能性输入。

数据类型：单|双倍的