7．4 交叉谱分析

7.4.1功能

计算两个不同的时间序列x1(t)，x2 (t)的交叉谱。

7.4.2方法说明

为不同的复谱，故交叉谱仍是复谱，它可以写成实部与虚部之和的形式，即：

7.4.3子程序语句

SUBROUTINE SPECTRUM2(N,M,X1,X2, R12,R21,P11,P22,P12,Q12,R212,THITA)

7.4.4哑元说明

N——整型变量，输入参数，序列长度。

M——整型变量，输入参数，最大落后长度。

X1——实型一维数组，输入参数，长度为N，存放第一要素观测序列值。

X2——实型一维数组，输入参数，长度为N，存放第二要素观测序列值。

R11——实型一维数组，输出参数，长度为M+1，存放X1的自相关系数。

R22——实型一维数组，输出参数，长度为M+1，存放X2的自相关系数。

R12——实型一维数组，输出参数，长度为M+1，存放落后相关系数（X2落后于X1）。 R21——实型一维数组，输出参数，长度为M+1，存放落后相关系数（X1落后于X2）。 P11——实型一维数组，输出参数，长度为M+1，存放X1的功率谱。

P22——实型一维数组，输出参数，长度为M+1，存放X2的功率谱。

P12——实型一维数组，输出参数，长度为M+1，存放X1和X2的协谱。

Q12——实型一维数组，输出参数，长度为M+1，存放X1和X2的正交谱。

R212——实型一维数组，输出参数，长度为M+1，存放X1和X2的凝聚谱。

THITA——实型一维数组，输出参数，长度为M+1，存放X1和X2的位相差谱。

7.4.5子程序（子程序名为：SPECTRUM2）

SUBROUTINE SPECTRUM2(N,M,X1,X2, R12,R21,P11,P22,P12,Q12,R212,THITA) REAL(4),DIMENSION(N)::X1,X2

REAL(4),DIMENSION(0:M)::R12,R21,R11,R22,P12,Q12,P22,P11,R212,THITA,B REAL(4):: X1BAR,X2BAR

INTEGER:: TAO,T

REAL(4),PARAMETER::PI=3.1415926

X1BAR=0

X2BAR=0

DO I=1,N

X1BAR=X1BAR+X1(I)

X2BAR=X2BAR+X2(I)

END DO

X1BAR=X1BAR/N

X2BAR=X2BAR/N

S1=0

S2=0

DO I=1,N

S1=S1+(X1(I)-X1BAR)**2

S2=S2+(X2(I)-X2BAR)**2

END DO

S1=SQRT(S1/N)