第十部分 数组高级运算函数(高等数学与统计)【参数数组为正序,返回数组也为正序，适合做中间变量计算用】

头文件声明变量：

1.         //第十部分 数组高级运算函数(高等数学与统计)
   
2.         vector<double>SeriesSine(vector<double> pc); //数组元素进行正弦转化
   
3.         vector<double>SeriesPearsoncorrelationcoefficient(vector<double> pc1, vector<double>pc2, int num); //皮尔逊相关系数计算公式

复制代码

公式模块增加：

1. //第十部分 数组高级运算函数(高等数学与统计)【参数数组为正序,返回数组也为正序】
   
2. //皮尔逊相关系数计算公式【参数数组为正序,返回数组也为正序】
   
3. vector<double> test::SeriesPearsoncorrelationcoefficient(vector<double> pc1, vector<double> pc2, int num)
   
4. {
   
5.         //皮尔逊相关系数（Pearson correlation coefficient）设计
   
6.         //声明两个vector变量存入X,Y两组数据
   
7.         vector<double>X, Y;//两个容器变量存入相对应数据元素
   
8.         vector<double>avgX, avgY;//X,Y的平均值
   
9.         vector<double>covXY;//X,Y的协方差(covariance)
   
10.         vector<double>stdX, stdY;//X,Y标准差
    
11.         vector<double>Rxy;//相关系数
    
12. 
    
13.         //采集X数据元素
    
14.         X = pc1;
    
15. 
    
16.         //采集Y数据元素(Y与X数据不同品种合约,但是相同周期)
    
17.         Y = pc2;
    
18. 
    
19.         //求length周期X,Y的平均值
    
20.         avgX = SeriesAverage(X, num);//调用数组计算公式计算
    
21.         avgY = SeriesAverage(Y, num);//计算Y在length周期的平均数
    
22. 
    
23.         vector<double>x_avgx, y_avgy, covxy;
    
24.         x_avgx = SeriesMinus(X, avgX);//计算X[i] - avgX[i]
    
25.         y_avgy = SeriesMinus(Y, avgY);//计算Y[i] - avgY[i]
    
26.         covxy = SeriesTimes(x_avgx, y_avgy);//计算（x-avgx）*(y-avgy)
    
27.         covXY = SeriesSummation(covxy, num);//sumcovxy
    
28. 
    
29.         //计算X,Y各自的方差
    
30.         vector<double>x_avgxpow, y_avgypow;
    
31.         vector<double>sum_x_avgx, sum_y_avgy;
    
32.         vector<double>std_x_avgx, std_y_avgy, std_xy;
    
33.         x_avgxpow = SeriesPowerN(x_avgx, 2);//(X[i] - avgX)^2
    
34.         y_avgypow = SeriesPowerN(y_avgy, 2);//(Y[i] - avgY)^2
    
35.         sum_x_avgx = SeriesSummation(x_avgxpow, num);
    
36.         sum_y_avgy = SeriesSummation(y_avgypow, num);
    
37.         std_x_avgx = Seriessquareroot(sum_x_avgx);
    
38.         std_y_avgy = Seriessquareroot(sum_y_avgy);
    
39.         std_xy = SeriesTimes(std_x_avgx, std_y_avgy);
    
40.         //计算相关性
    
41.         Rxy = SeriesDivide(covXY, std_xy);//计算相关性
    
42.         return Rxy;
    
43. }
    
44. //Vector正弦值数组计算【参数数组为正序,返回数组也为正序】
    
45. vector<double> test::SeriesSine(vector<double> pc)
    
46. {
    
47.         vector<double>newvector;//声明新变量容器
    
48.         for (size_t i = 0; i < pc.size(); i++)
    
49.         {
    
50.                 newvector.push_back(sin(pc[i]));//sin(x)
    
51.         }
    
52.         return newvector;
    
53. }
    
54. //第十部分 数组高级运算函数(高等数学与统计)【参数数组为正序,返回数组也为正序】

复制代码