常方敏 李素敏

摘? 要： 針對(duì)最小二乘算法在處理坐標(biāo)轉(zhuǎn)換存在準(zhǔn)則單一這一問(wèn)題，本文提出一種思想，將和聲搜索算法加入最小二乘坐標(biāo)轉(zhuǎn)換模型中，從而使傳統(tǒng)坐標(biāo)轉(zhuǎn)換這一P類(lèi)問(wèn)題，轉(zhuǎn)換為NP類(lèi)問(wèn)題。通過(guò)不斷迭代使其結(jié)果趨近于所設(shè)定的函數(shù)解，從而實(shí)現(xiàn)解的多樣性。本文以方向精度最優(yōu)為例，分別采用模擬數(shù)據(jù)和實(shí)測(cè)點(diǎn)云數(shù)據(jù)來(lái)驗(yàn)證本文思想。結(jié)果表明，在求某一方向精度時(shí)將和聲搜索算法加入模型中會(huì)提高方向精度，加權(quán)后也可以保證其垂直方向一定的精度。

關(guān)鍵詞： 和聲搜索算法;坐標(biāo)轉(zhuǎn)換;方向精度;坐標(biāo)轉(zhuǎn)換準(zhǔn)則

中圖分類(lèi)號(hào)： P228? ? 文獻(xiàn)標(biāo)識(shí)碼： A? ? DOI：10.3969/j.issn.1003-6970.2020.10.027

本文著錄格式：常方敏，李素敏. 和聲搜索算法在坐標(biāo)轉(zhuǎn)換中的應(yīng)用[J]. 軟件，2020，41（10）：109113

【Abstract】： Aiming at the problem that the least squares algorithm has a single criterion when dealing with coordinate transformation， this paper proposes an idea to add an intelligent optimization algorithm to the least square coordinate transformation model， so that the traditional polynomial problem of coordinate transformation is converted to the Non-polynomial problem. Through continuous iteration， the results approach the set function solution， so as to realize the diversity of the solution. In this paper， taking the best direction accuracy as an example， the simulation data and the measured point cloud data are used to verify the idea of this paper. The results show that HS algorithm is added to the model， the direction accuracy can be improved， and the vertical direction accuracy can be guaranteed after weighting.

【Key words】： Harmony search algorithm; Coordinate transformation; Direction accuracy; Coordinate transformation criterion

0? 引言

在2001年，Zong W G等[1]提出了一種新穎的啟發(fā)式智能優(yōu)化算法——和聲搜索（Harmony Search，HS）算法，由于該算法操作簡(jiǎn)單，收斂性與初始值無(wú)關(guān)，具備一定魯棒性并且可以避免陷入局部最優(yōu)等優(yōu)點(diǎn)，因此近年來(lái)HS應(yīng)用在各個(gè)領(lǐng)域，例如：穩(wěn)健大地測(cè)量[2]、地下水位的長(zhǎng)期預(yù)測(cè)[3]、斜坡堤護(hù)面塊體[4]等等。

傳統(tǒng)的最小二乘（Least Square，LS）算法在用于坐標(biāo)轉(zhuǎn)換模型存在一定的缺陷[5]。文獻(xiàn)[6]在坐標(biāo)轉(zhuǎn)換模型中引入基準(zhǔn)旋轉(zhuǎn)中心，解決轉(zhuǎn)換系數(shù)矩陣病態(tài)的問(wèn)題。文獻(xiàn)[7]提出基于單位四元數(shù)的任意旋轉(zhuǎn)角度的三維坐標(biāo)轉(zhuǎn)換，解決了大角度轉(zhuǎn)換的問(wèn)題。文獻(xiàn)[8]采用選權(quán)迭代的總體最小二乘算法，通過(guò)降低粗差的權(quán)重，從而得到更高精度轉(zhuǎn)換參數(shù)。文獻(xiàn)[9]在加權(quán)總體最小二乘的坐標(biāo)轉(zhuǎn)換模型中加入驗(yàn)后估計(jì)，解決了加權(quán)總體最小二乘單位權(quán)方差不同的問(wèn)題。以上文獻(xiàn)分別從LS坐標(biāo)轉(zhuǎn)換模型不同角度的缺陷出發(fā)，解決了一些基于LS坐標(biāo)轉(zhuǎn)換的問(wèn)題。

本文從另一角度出發(fā)，引入啟發(fā)式智能優(yōu)化算法，將坐標(biāo)轉(zhuǎn)換這一確定性問(wèn)題即P類(lèi)問(wèn)題轉(zhuǎn)換為NP類(lèi)問(wèn)題，從而實(shí)現(xiàn)不同的優(yōu)化準(zhǔn)則，即解決坐標(biāo)轉(zhuǎn)換模型轉(zhuǎn)換準(zhǔn)則單一這一問(wèn)題。本文將和聲搜索算法，這一智能算法加入坐標(biāo)轉(zhuǎn)換模型中，以某一方向精度最優(yōu)為例，通過(guò)模擬數(shù)據(jù)和點(diǎn)云數(shù)據(jù)進(jìn)行對(duì)比實(shí)驗(yàn)從而驗(yàn)證本文思想。

1? 基本理論與流程

1.1? 和聲最小二乘算法基本理論

HS算法是類(lèi)比樂(lè)師反復(fù)調(diào)整不同樂(lè)器達(dá)到最優(yōu)美和聲這一過(guò)程。和聲創(chuàng)作過(guò)程的音樂(lè)元素與對(duì)應(yīng)數(shù)學(xué)關(guān)系如表1所示，即將最優(yōu)美和聲類(lèi)比于最優(yōu)解。每一個(gè)樂(lè)器的聲音Ri（i=1，2，…，n）類(lèi)比于每個(gè)解向量Xi （i=1，2，…，n），樂(lè)器的和聲類(lèi)比于函數(shù)解的值，其最優(yōu)美和聲對(duì)應(yīng)的是需要求得的最優(yōu)解。其中和聲搜索算法包含以下幾個(gè)初始變量[1]：和聲記憶庫(kù)（HM）、和聲記憶庫(kù)的大?。℉MS）、和聲記憶庫(kù)的取值概率（HMCR）、音調(diào)微調(diào)概率（PAR）、音調(diào)微調(diào)步長(zhǎng)（bw）。關(guān)于標(biāo)準(zhǔn)和聲搜索算法的詳細(xì)操作流程，請(qǐng)參見(jiàn)文獻(xiàn)[1]。

1.2? 本文算法的基本理論

基于LS算法布爾莎七參數(shù)模型方程為：

其中為目標(biāo)坐標(biāo)系下的坐標(biāo);為原坐標(biāo)系下的坐標(biāo);為平移參數(shù);、、分別為坐標(biāo)系繞X、Y、Z軸旋轉(zhuǎn)的矩陣;m為原坐標(biāo)系到目標(biāo)坐標(biāo)系的尺度參數(shù)，通過(guò)模型方程，基于最小二乘計(jì)算準(zhǔn)則便可得到七參數(shù)，進(jìn)而通過(guò)七參數(shù)實(shí)現(xiàn)坐標(biāo)轉(zhuǎn)換。

HS并不能直接加入坐標(biāo)轉(zhuǎn)換模型中，因此本文對(duì)傳統(tǒng)的HS算法進(jìn)形了適當(dāng)?shù)恼{(diào)整，從而使其可以結(jié)合LS算法應(yīng)用于坐標(biāo)轉(zhuǎn)換模型中。對(duì)HS的調(diào)整如下：