[发明专利]一种用于波形互相关双差定位算法的并行优化方法

权利要求书

1.一种用于波形互相关双差定位算法的并行优化方法，其特征在于，步骤如下：

步骤1，获取总任务数和总进程数，并将总任务平均分给每个进程，并得出每个进程的任务初始位置；

步骤2，计算获取初始匹配位置A_s1B_s2；

步骤3，计算获取结束匹配位置A_e1B_e2；

步骤4，根据获取的A_s1B_s2位置和A_e1B_e2位置读取对应数据进行计算；

步骤1具体包括：

步骤101，首先选定两个台站，台站A和台站B，每个台站都有N天的数据，A台站第i天的数据可表示为A_i，B台站第i天的数据可表示为B_i，台站A和台站B之间的数据需要进行匹配计算，匹配方式是对每个A_i，都需要和B_j(ji)进行匹配计算，匹配计算格式定义为A_iB_j(1≤ij≤N)，由此便可得到A和B匹配计算总的任务数为个，定义为PAIRS_LEN；

步骤102，将总进程个数定义为np，目标是求总任务平均分配给每个进程后，每个进程的任务初始匹配位置A_s1B_s2和结束匹配位置A_e1B_e2；

步骤103，将PAIRS_LEN整除np，得到每个进程要处理的任务数量，定义为task_num_k，k表示进程索引，再将PAIRS_LEN余除np，得到暂未被分配到进程的任务数量，定义为residul_day；将residul_day个任务数分配到总进程中前residul_day个进程，所以前residul_day个进程任务数task_num_k(k≤residul_day)自加1；至此完成进程任务数量分配；

步骤104，根据每个进程需要处理的任务数量task_num_k，累加计算得到每个进程的任务初始位置task_num_start_p，p表示进程索引，k＝1,2,3...p...；

步骤2具体包括：

步骤201，将A_i需要做匹配计算的个数定义为A_LEN_i，定义A_SUM_j为第j天之前需要做匹配计算数量；

步骤202，遍历整个A_SUM数组，找到第一个大于每个进程中task_num_start_p值的A_SUM_t，并保存t值；

步骤203，t-1即为开始匹配位置A_s1B_s2中的s1，再由t-1值得到A_SUM_t-1；

步骤204，将task_num_start_p-A_SUM_t-1即可得到每个进程在第t天的需要处理的任务数量

步骤205，再把步骤204中的task_num_start_p-A_SUM_t-1值加上t，即可得到每个进程的初始匹配位置A_s1B_s2中的s2，至此得到了每个进程的初始匹配位置A_s1B_s2；

步骤3具体包括：

步骤301，定义变量nday1用于记录每个进程处理的任务A_iB_j中i_max-i_min+1的值，nday1代表A_i中i的个数；

步骤302，定义变量nday2用于记录每个进程处理的任务A_iB_j中j_max-j_min+1的值，nday2代表B_j中j的个数；

步骤303，定义变量day1_min，day2_min用于记录每个进程处理的任务A_iB_j中最小的i和j；同时定义day1_start，day2_start表示步骤203、步骤205求出的A_s1B_s2中的s1和s2；

步骤304，将day1_min+nday1-1即可得到得到A_e1B_e2中的e1，将day2_start-day2_min+nday1即可得到A_e1B_e2中的e2；

步骤4具体包括：

步骤401，定义结构体数组st_k，k为进程索引，根据每个进程的初始匹配位置A_s1B_s2和结束匹配位置A_e1B_e2读入数据到结构体数组st_k中；

步骤402，对结构体数组st_k进行计算，公式为

t为A_i、B_j中数据的索引，至此得计算结果。