说明：  计算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)