计算ARMA(p，q)模型的功率谱密度。

形参说明：

b——双精度实型一维数组，长度为(q+1),存放ARMA(p，q)模型的滑动平均系数。

a——双精度实型一维数组，长度为(p+1),存放ARMA(p,q)模型的自回归系数。

q——整型变量,ARMA(p,q)模型的滑动平均阶数。

p——整型变量,ARMA(p,q)模型的自回归阶数。

sigma2——双精度实型变量,ARMA(p,q)模型白噪声激励的方差。

fs——双精度实型变量,采样频率(Hz)。

x——双精度实型一维数组,长度为len。当sign＝0时，存放功率谱密度；当sign＝

1时，存放用分贝表示的功率谱密度。

freq——双精度实型一维数组，长度为len。存放功率谱密度所对应的频率。

len——整型变量，功率谱密度的数据点数。

sign——整型变量,当sign＝0时，计算功率谱密度；当sign＝1时，计算用分贝表

示的功率谱密度。

-calculated ARMA (p, q) model of the power spectral density.- Participation : b-- double-precision real one-dimensional arrays, length (q 1), storage ARMA (p, q) model of the average coefficient of sliding. A-- double-precision real one-dimensional arrays, length (p 1), storage ARMA (p, q) model of regression coefficients. Q-- integer variables, ARMA (p, q) moving average model of order. P-- integer variables, ARMA (p, q) model of order since the reunification. Sigma2-- double precision real-type variables, ARMA (p, q) model white noise variance. Fs-- double precision real-type variables, sampling frequency (Hz). X-- double-precision real one-dimensional array, the length of len. When the sign = 0, storage power spectral density; When the sign = 1:00, storage decibel levels, said th