[发明专利]一种基于MATLAB的典型相关分析方法及装置

权利要求书

1.一种基于MATLAB的典型相关分析方法，其特征在于，在MATLAB中设定典型相关分析算法，还包括：

获取待分析的原始数据；

将获取的所述原始数据组织为原始矩阵；

根据所述原始矩阵，将所述原始矩阵转换为MATLAB中的原始变量；

根据MATLAB中的预先设定的典型相关分析算法，对所述原始变量进行典型相关分析。

2.根据权利要求1所述的方法，其特征在于，所述根据所述原始矩阵，将所述原始矩阵转换为MATLAB中的原始变量，包括：

将所述原始矩阵导入并存储至MATLAB，存储格式为data.mat；

为所述原始矩阵中的每个原始数据赋予变量名，形成所述原始变量。

3.根据权利要求1所述的方法，其特征在于，对所述原始变量进行典型相关分析，包括：

计算所述原始矩阵的协方差矩阵DXY：

DXY＝cov{[zscore(X)zscore(Y)]}；

其中，第一原始矩阵X＝[x1x2x3x4]，第二原始矩阵Y＝[y1y2y3y4]；

根据所述协方差矩阵，设定下述矩阵A、矩阵B：

A＝inv[DXY(1:4,1:4)]*DXY(1:4,5:8)*inv[DXY(5:8,5:8)]*DXY(5:8,1:4)；

B＝inv[DXY(5:8,5:8)]*DXY(5:8,1:4)*inv[DXY(1:4,1:4)]*DXY(1:4,5:8)；

分别计算矩阵A的特征值D_A和特征向量U_A和矩阵B的特征值D_B和特征向量U_B；

在矩阵A的特征值D_A和矩阵B的特征值D_B相等时，计算典型变量u(:,i)、v(:,i)：

u(:,i)＝U_A(:,nA2(4,4-i+1))./sqrtm((U_A(:,nA2(4,4-i+1)))'*

DXY(1:4,1:4)*U_A(:,nA2(4,4-i+1)))；

v(:,i)＝inv(DXY(5:8,5:8))*DXY(5:8,1:4)*u(:,i)./sqrtm(mA2(4,4-i+1))；

[mA1,nA1]＝sort(D_A,1)；

[mA2,nA2]＝sort(mA1,2)；

其中，i＝1，2，3，4；mA1是针对D_A排序时的中间矩阵，nA1是排序的第一索引矩阵；mA2是针对D_A排序完成后的矩阵，nA2是排序的第二索引矩阵。

4.根据权利要求3所述的方法，其特征在于，进一步包括：

通过调用MATLAB中的eig()函数，分别计算矩阵A的特征值D_A和特征向量U_A和矩阵B的特征值D_B和特征向量U_B：

[U_A,D_A]＝eig(A)

[U_B,D_B]＝eig(B)

通过调用MATLAB中的sqrtm()函数，计算典型变量：u(:,i)、v(:,i)。

5.根据权利要求3或4所述的方法，其特征在于，在所述对所述原始变量进行典型相关分析之后，进一步包括：

根据计算的所述典型变量，确定所述原始变量在所述典型变量上的负荷矩阵，并输出所述负荷矩阵作为所述典型相关分析的结果；其中，所述原始变量在所述典型变量上的负荷矩阵ans通过下式计算：

ans1＝inv(u(:,i))；

ans2＝inv(v(:,i))。

6.一种基于MATLAB的典型相关分析装置，其特征在于，包括：

设定单元，用于在MATLAB中设定典型相关分析算法；

获取单元，用于获取待分析的原始数据；

组织单元，用于将获取的所述原始数据组织为原始矩阵；

转换单元，用于根据所述原始矩阵，将所述原始矩阵转换为MATLAB中的原始变量；

分析单元，用于根据MATLAB中的预先设定的典型相关分析算法，对所述原始变量进行典型相关分析。

7.根据权利要求6所述的装置，其特征在于，所述转换单元，用于将所述原始矩阵导入并存储至MATLAB，存储格式为data.mat；为所述原始矩阵中的每个原始数据赋予变量名，形成所述原始变量。