MATLAB通信工具箱之Periodogram

Periodogram：Periodogrampower spectral density estimate

 

pxx = periodogram(x)：返回经过矩形窗口处理的输入信号x的PSD估计值。如果x是实信号，则pxx是单边PSD估计。如果x是复信号，则pxx是双边PSD估计。DFT中的点数nfft是最大值256或比信号长度大两倍的下一个功率

 

pxx = periodogram(x,window)：返回经过window处理的PSD值，且窗是与x相同长度的向量。

 

pxx = periodogram(x,window,nfft)：采用nfft点DFT变换。如果nfft大于信号长度，则x补零到长度nfft再做DFT。如果nfft小于信号长度，则对信号进行预包裹（wrapped modulo）成nfft长度，再对包裹（datawrap）求和。Forexample, the input signal [1 2 3 4 5 6 7 8] with nfft equal to 4results in the periodogram of sum([1 5; 2 6; 3 7; 4 8],2).

 

[pxx,w] = periodogram(___)：返回归一化频率向量，w。则w的范围为{nfft为偶则为[0,π]，奇则为[0,π）}；若pxx是双边的，则w的范围为 [0,2π）.

 

[pxx,f] = periodogram(___,fs)：以单位时间的周期数返回频率矢量f。采用频率fs为每单位时间内的采用数量。若单位时间的单位是秒，则f的单位是 cycles/sec(Hz)。对于实信号，f的范围为{nfft为偶则为[0,fs/2]，奇则为[0,fs/2）}；对于复信号，f的范围是[0,fs/2）。

 

[pxx,w] = periodogram(x,window,w)：返回双边周期图，频率归一化由矢量w指定。W，至少包含2个元素。

 

[pxx,f] = periodogram(x,window,f,fs)：返回双边周期图，频率归一化由矢量f指定。f，至少包含2个元素。f是每单位时间的周期数。采样频率fs是每单位时间的采样数。如果单位时间的单位是秒，则f的单位是 cycles/sec(Hz)。

 

[___] = periodogram(x,window,___,freqrange)：在freqrange指定的频率范围内返回周期图。freqrange的有效选项为：'oneided'，'twosided'或'centered'。

 

[___] = periodogram(x,window,___,spectrumtype)：如果spectrumtype值为“psd”则返回PSD（功率谱密度）；如果spectrumtype值为“power”，则返回功率谱。

 

[pxx,f,pxxc] = periodogram( ___,'ConfidenceLevel',probability)：返回对PSD做估计的概率置信区间，放置到pxxc中。

 

periodogram(___)：绘出没有输出参数的，当前数字窗口中每单位频率的PSD图（单位为dB）

 

Input arguments：（输入参数）

X：输入信号，为一个行或列向量。数据类型：单、双精度；支持复数。

 

Window：窗，跟输入信号一样长度的行或列向量。如果为空，则使用rectangular窗为默认窗。数据类型：单、双精度。

 

nfft：DFT运算点数，为正整数。对于实输入信号x的PSD评估，若nfft为偶，则pxx=(nfft/2+1)；若nfft为奇，则pxx=(nfft+1)2。对于复输入信号x的PSD评估，pxx的长度总是nfft。如果nfft为空，则使用默认值。数据类型：单、双精度。

 

fs：采样频率，为正变量。它是单位时间内的采用数。若单位时间的单位是秒，则它的单位是Hz。

 

w：Goertzel算法的归一化频率，至少有两个元素的行或列向量。归一化频率以弧度/采样为单位。如：w = [pi/4 pi/2]。数据类型：双精度。

 

f：Goertzel算法的周期性频率，至少有两个元素的行或列向量。它以周期/单位时间为单位。单位时间由采样频率fs指定。如果fs以采样/每秒为单位，那么f的单位是Hz。如：fs=1000；f=[100200]。数据类型：双精度。

 

Freqrange：PSD的频率范围，由'onesided','twosided', or'centered'之一指定。实信号默认'onesided'，复信号默认'twosided'。每个选项指定的范围如下：

 

'onesided'：返回实信号x的单边PSD。若nfft为偶，则pxx=(nfft/2+1)，计算范围为[0π]弧度/采样，f范围为[0,fs/2]周期/时间；若nfft为奇，则pxx=(nfft+1)2，计算范围为[0，π）弧度/采样，f范围为[0,fs/2)周期/时间。

 

'twosided'：返回实信号或复信号x的双边PSD。在这种情况下，pxx的长度为nfft，并在间隔[0,2π）弧度/样本上计算，f范围为[0,fs)周期/时间。

 

'centered'：返回实信号或复信号x的居中双边PSD。在这种情况下，pxx的长度为nfft；若nfft为偶，计算范围为(-ππ]弧度/采样，f范围为(-fs/2,fs/2]周期/时间；若nfft为奇，计算范围为(-π，π）弧度/采样，f范围为(-fs/2,fs/2)周期/时间。数据类型，字符串型。

 

Spectrumtype：功率谱形式，由“psd”、“power”。默认情况下或指定“psd”将返回功率谱密度。指定“power”经过窗的等效噪声带宽来对每个PSD值进行评估。使用“power”是可以获得每个频率功率估计值。

 

Probability：覆盖真实PSD的概率，指定为范围（0,1）内的标量。输出pxxc包括真实PSD的概率区间的下限和上限。

 

Output arguments：（输出参数）

pxx：PSD估计，为实数值、非负列向量。PSD的单位是每单位频率的时间序列的平方。

例如：如果输入是电压，则PSD是单位频率的单位电压的平方。

对于时间序列，如果假定电阻为1欧姆，并以Hz为单位指定采样频率，则PSD单位为瓦特/赫兹。数据类型：单、双精度

 

w：归一化频率，为实列向量。若pxx是单边的PSD，则若nfft为偶，w计算范围为[0，π]；则若nfft为奇，w计算范围为[0，π）；若pxx是双边的PSD，w计算范围为[0，π)。对于中心PSD，f的范围是（-pai,pai]弧度/采样（nfft偶）；f的范围是（-pai,pai）弧度/采样（nfft奇）。数据类型：双精度。

 

f：周期频率，实列向量。若pxx是单边的PSD，则若nfft为偶，f计算范围为[0，fs/2]；则若nfft为奇，w计算范围为[0，fs/2）；若pxx是双边的PSD，f计算范围为[0，fs)。对于中心PSD，f的范围是（-fs/2,fs/2]周期/时间（nfft偶）；f的范围是（-fs/2,fs/2）周期/时间（nfft奇）。数据类型：双精度。

 

pxxc:置信区间，Nx2矩阵实元素。矩阵的行维度等于PSD长度。第一列是置信下限，第二列是pxx行中相应PSD值的置信上限。置信区间的覆盖概率取决于输入概率值。数据类型：单、双精度。