丁宇堃，宋荊洲，尚志豪，韓亮亮

力反饋遙操作系統(tǒng)中的在線環(huán)境參數(shù)辨識(shí)研究

丁宇堃1，宋荊洲1，尚志豪1，韓亮亮2

（1.北京郵電大學(xué)自動(dòng)化學(xué)院，北京100876；2.上海宇航系統(tǒng)工程研究所，上海201108）

力反饋遙操作系統(tǒng)是載人航天工程中的重要技術(shù)組成部分，而環(huán)境的在線參數(shù)辨識(shí)是實(shí)現(xiàn)大時(shí)延下反饋力的本地實(shí)時(shí)預(yù)測(cè)的關(guān)鍵技術(shù)。使用交叉驗(yàn)證分析了不同建模方法的準(zhǔn)確性和適用條件；對(duì)常用的在線參數(shù)辨識(shí)算法進(jìn)行實(shí)驗(yàn)驗(yàn)證，從參數(shù)收斂性和對(duì)時(shí)變環(huán)境的跟蹤性等方面進(jìn)行對(duì)比分析。實(shí)驗(yàn)結(jié)果表明，Hunt-Crossley模型以高復(fù)雜獲得了更高的建模精度，適宜在數(shù)據(jù)充分的離線條件下采用；同時(shí)，在參數(shù)辨識(shí)算法方面，相比于同類算法，自擾動(dòng)遞歸最小二乘法在參數(shù)辨識(shí)的收斂性和對(duì)時(shí)變環(huán)境的跟蹤性上都有較好的表現(xiàn)，是實(shí)際應(yīng)用中的首選方法。該結(jié)果可為環(huán)境動(dòng)力學(xué)模型的在線參數(shù)辨識(shí)技術(shù)的工程應(yīng)用提供參考。

遙操作；環(huán)境模型；參數(shù)辨識(shí)

Abstract：Force feedback teleoperation system is an important part of the manned space flight project，and the on-line parameter identification of environment is the key technology to realize the local real-time prediction of the feedback force under large time delay.In this paper，the method of online environment parameter identification in force feedback teleoperation system was studied.The accuracy of different modeling methods and applicable conditions were analyzed through cross validation；then the commonly used online parameter identification algorithms were experimentally verified.The convergence of the environment parameters and the tracking of the time-varying environment were analyzed.The experimental results showed that the Hunt-Crossley model achieved higher modeling accuracy with high complexity and was suitable for use under sufficient offline conditions.At the same time， compared to similar parameter identification algorithms， the tracking of dynamic recursive least squares method had better performancein parameter identification of self convergence and time-varying environment and was the preferred method for practical application.This paper provides a reference for engineering application of technology for online parameter identification of dynamic model of environment.

Key words：teleoperation； environment model； parameter estimation

1 引言

隨著空間站建設(shè)和運(yùn)營(yíng)的深入，機(jī)器人天地遙操作裝配任務(wù)的需求越來(lái)越多。如果本地端操作設(shè)備能夠?qū)崟r(shí)提供準(zhǔn)確的反饋力信息，將有助于增強(qiáng)操作者對(duì)工作環(huán)境的感知，提高操作效率，改善操作者的操作體驗(yàn)。因而采用力預(yù)測(cè)的方法，構(gòu)造力預(yù)測(cè)模型，解算出與顯示的仿真環(huán)境同步的反饋力，通過(guò)力反饋手控器施加給操作者，使操作者感受到操作中的反饋力，能夠帶來(lái)更好的臨場(chǎng)感。

環(huán)境模型建模及參數(shù)辨識(shí)是實(shí)現(xiàn)接觸力預(yù)測(cè)的關(guān)鍵。通常將與機(jī)器人系統(tǒng)末端執(zhí)行器發(fā)生接觸的物體稱為“環(huán)境”，將環(huán)境的接觸動(dòng)力學(xué)特性建模為環(huán)境動(dòng)力學(xué)模型。在接觸階段，可以通過(guò)傳感器獲得從端執(zhí)行器和環(huán)境的相對(duì)位置、速度和接觸力等信息，從而對(duì)動(dòng)力學(xué)模型中在線參數(shù)辨識(shí)進(jìn)行修正。參數(shù)辨識(shí)結(jié)果可以輸入自適應(yīng)阻抗控制器或模型參考自適應(yīng)控制器來(lái)實(shí)現(xiàn)穩(wěn)定的接觸。此外，基于預(yù)測(cè)模型的力反饋遙操作系統(tǒng)也依賴在線環(huán)境參數(shù)辨識(shí)力修正主端預(yù)測(cè)模型的動(dòng)力學(xué)參數(shù)，為操作者提供準(zhǔn)確的力覺(jué)臨場(chǎng)感。

對(duì)于環(huán)境的幾何特征辨識(shí)，國(guó)內(nèi)外的研究者也開(kāi)展了相當(dāng)多的研究：Xu Xiao等基于飛行時(shí)間的立體相機(jī)，進(jìn)行距離測(cè)定從而進(jìn)行環(huán)境幾何參數(shù)辨識(shí)［1］。李會(huì)軍等使用力和位置信息的綜合處理來(lái)對(duì)環(huán)境的位置信息進(jìn)行辨識(shí)［2］，Velanas等使用遞歸最小二乘法對(duì)環(huán)境進(jìn)行曲率辨識(shí)［3］。Willaert等則研究了利用視覺(jué)信息辨識(shí)接觸點(diǎn)的位置和法線方向［4-5］。

相關(guān)研究應(yīng)用中采用的環(huán)境動(dòng)力學(xué)模型主要包括純剛度模型和Kelvin-Voigt模型［6］。近年來(lái)，學(xué)者開(kāi)始關(guān)注準(zhǔn)確性更高的Hunt-Crossley模型，然而它的非線性增加了對(duì)其進(jìn)行參數(shù)辨識(shí)的難度，導(dǎo)致它目前還難以實(shí)現(xiàn)實(shí)際應(yīng)用［7］。在不同模型的比較方面，Yamamoto等首先使用了交叉驗(yàn)證的方法比較了不同環(huán)境模型的建模精度［8］。

在環(huán)境的動(dòng)力學(xué)參數(shù)辨識(shí)方面，Haddadi等研究了基于梯度下降的最小二次均方算法（BLMS），并與基于窗函數(shù)的最小二乘法（BLS）做了比較［6-9］。Love等較早使用了遞歸最小二乘法來(lái)辨識(shí)環(huán)境剛度，并用于調(diào)整從端阻抗控制器的參數(shù)，改善接觸性能［10］。Colton等則使用了指數(shù)加權(quán)遞歸最小二乘法（EWRLS）進(jìn)行參數(shù)辨識(shí)［11］。Achhammer等使用了自擾動(dòng)遞歸最小二乘法（SPRLS）進(jìn)行了混合模型的參數(shù)辨識(shí)，并取得了較好的效果［12］。雖然研究者探索了多種參數(shù)辨識(shí)方法，遞歸最小二乘法類算法憑借其簡(jiǎn)單有效和魯棒性強(qiáng)的特點(diǎn)得到了最廣泛的應(yīng)用。

本文對(duì)機(jī)器人系統(tǒng)中的環(huán)境動(dòng)力學(xué)模型建模和參數(shù)辨識(shí)技術(shù)進(jìn)行比較驗(yàn)證研究，比較分析不同建模方法的準(zhǔn)確性和適用范圍，首次使用交叉驗(yàn)證的方法衡量在線參數(shù)辨識(shí)過(guò)程中模型準(zhǔn)確性的變化，以期為環(huán)境模型參數(shù)在線辨識(shí)技術(shù)在機(jī)器人空間裝配中的應(yīng)用提供指導(dǎo)。

2 動(dòng)力學(xué)模型建模

首先分析現(xiàn)有的環(huán)境動(dòng)力學(xué)模型建模方法，在大范圍實(shí)驗(yàn)條件下對(duì)其進(jìn)行實(shí)驗(yàn)驗(yàn)證，使用交叉驗(yàn)證的方法更完整地分析其建模準(zhǔn)確性。

2.1 接觸動(dòng)力學(xué)模型

純剛度模型（以下簡(jiǎn)稱K模型）只有一個(gè)參數(shù)，即環(huán)境剛度k，在接觸階段，接觸力f按環(huán)境剛度正比于侵入深度［7］，如公式（1）所示。

其中，x表示侵入物體簡(jiǎn)化的代理點(diǎn)位置，xe為環(huán)境位置。kK和fK的下標(biāo)表示對(duì)應(yīng)K模型的參數(shù)，以下類似。

Kelvin-Voigt模型（以下簡(jiǎn)稱KV模型）考慮了阻尼項(xiàng)bKV［7］，如公式（2）所示。

Hunt-Crossley模型（以下簡(jiǎn)稱HC模型）具有3個(gè)參數(shù)，通過(guò)阻尼項(xiàng)與位置的耦合，解決了KV模型的反向力問(wèn)題，更接近實(shí)際環(huán)境物理性質(zhì)。指數(shù)n的存在使其更好地描述接觸力隨侵入深度的非線性變化［8］，模型如公式（3）所示。

相同輸入信號(hào)下，典型的KV模型和HC模型的接觸力與侵入深度的磁滯回線如圖1所示。對(duì)于KV模型，當(dāng)（x-xe）→0而且x·＜0時(shí)，fKV為負(fù)值，而實(shí)際接觸時(shí)，接觸力始終為正值。為負(fù)值的接觸力會(huì)造成力反饋遙操作中的“黏附感”，影響操作者對(duì)環(huán)境的認(rèn)知和正確操作。而HC模型則不存在負(fù)值力的問(wèn)題，整個(gè)接觸過(guò)程中fHC始終為正值。

另一方面，KV模型的不連續(xù)問(wèn)題也給系統(tǒng)穩(wěn)定性帶來(lái)不利影響。由于直接相關(guān)，在接觸的開(kāi)始和結(jié)束階段，的突變會(huì)導(dǎo)致的突變，也導(dǎo)致了系統(tǒng)能量的突變。而實(shí)際的接觸中，即使考慮碰撞效應(yīng)，實(shí)際接觸力也在連續(xù)變化。

圖1 KV模型與HC模型的磁滯回線Fig.1 Hysteresis loops of KV model and HC model

2.2 模型準(zhǔn)確性分析

根據(jù)上述3種模型參數(shù)數(shù)量的依次遞增和實(shí)驗(yàn)結(jié)果驗(yàn)證，K、KV和HC模型的復(fù)雜度和準(zhǔn)確性依次遞增，通常在實(shí)際應(yīng)用中根據(jù)不同情況選取合適的建模方法。為了研究不同的模型在實(shí)際環(huán)境建模中的效果，首先比較不同模型在離線辨識(shí)中對(duì)不同剛度的環(huán)境的建模準(zhǔn)確性。實(shí)驗(yàn)構(gòu)造了剛度依次遞增的4種典型的接觸環(huán)境，分別為圖2中海綿、覆蓋于海綿上的橡膠墊、橡膠墊、金屬蓋板，并將其依次記為A、B、C、D。其剛度分別約為400 N／m、2000 N／m、8000 N／m和18 000 N／m。使用末端腕部裝有ATI mini40力傳感器的三自由度平動(dòng)機(jī)械臂與環(huán)境進(jìn)行接觸，機(jī)械臂末端位置和接觸力的采樣頻率為170 Hz。判斷發(fā)生接觸的力閾值設(shè)為0.05 N。

圖2 實(shí)驗(yàn)接觸環(huán)境Fig.2 Contact environments in experiment

控制機(jī)械臂末端與不同環(huán)境發(fā)生接觸，并記錄接觸過(guò)程中機(jī)械臂末端位置和接觸力信息。每種接觸環(huán)境進(jìn)行5次接觸實(shí)驗(yàn)，共得到20組數(shù)據(jù)。每次接觸實(shí)驗(yàn)中，機(jī)械臂位置由操作者手動(dòng)控制，并在5次實(shí)驗(yàn)中采用不同的動(dòng)作模式來(lái)盡可能覆蓋可能的接觸情況，例如侵入深度偏大和偏小，速度較快和較慢。圖3為4種環(huán)境的典型力／位曲線圖。

使用最小二乘法根據(jù)得到的數(shù)據(jù)對(duì)各個(gè)模型選取最優(yōu)參數(shù)，比較其理論上的最高準(zhǔn)確度。對(duì)于線性模型［14］，預(yù)測(cè)力可根據(jù)公式（4）進(jìn)行計(jì)算。

其中，f（t）為t時(shí)刻的預(yù)測(cè)力，θ為未知參數(shù)，φ為回歸向量?？芍苯影凑帐剑?）計(jì)算θ。

其 中，X＝[φ1，φ2，…，φn]Τ，Y＝[f1，f2，…，fn]Τ。 對(duì)于非線性的HC模型，本文采用信賴域算法進(jìn)行迭代計(jì)算，當(dāng)目標(biāo)函數(shù)變化小于10-6時(shí)停止迭代。以環(huán)境B的第一組數(shù)據(jù)為例，根據(jù)所得參數(shù)計(jì)算得到的接觸力fp與實(shí)際接觸力fs的比較如圖4所示?？梢?jiàn)，各環(huán)境模型建模都基本準(zhǔn)確，K模型和KV模型具有明顯的力誤差，而HC模型的力誤差較小。

使用平均標(biāo)準(zhǔn)均方根誤差量化衡量模型準(zhǔn)確性。相比于平均力誤差，標(biāo)準(zhǔn)均方根誤差使用誤差的二次項(xiàng)，對(duì)較大的誤差更敏感，并考慮原始數(shù)據(jù)的變化范圍，能夠更好地衡量模型的準(zhǔn)確性。對(duì)誤差估計(jì)量y，標(biāo)準(zhǔn)均方根誤差的計(jì)算公式如式（6）：

實(shí)驗(yàn)結(jié)果如表1所示。由表可見(jiàn)，對(duì)于各種環(huán)境，HC模型都具有最低的力誤差，且對(duì)于中等剛度的對(duì)象優(yōu)勢(shì)較為明顯。KV模型相比于K模型，雖然具有更高的復(fù)雜度，但是準(zhǔn)確性只有微小的提高。

圖3 不同環(huán)境力／位曲線圖Fig.3 Force and position profiles of different environments

圖4 模型力誤差Fig.4 Force error of different models

表1 環(huán)境模型自驗(yàn)證誤差Table 1 Self-validation error of different models

2.3 模型準(zhǔn)確性交叉驗(yàn)證

常用的自驗(yàn)證方法是使用同一數(shù)據(jù)集進(jìn)行參數(shù)辨識(shí)和準(zhǔn)確性驗(yàn)證，而交叉驗(yàn)證則使用重復(fù)實(shí)驗(yàn)的多次數(shù)據(jù)分別進(jìn)行參數(shù)辨識(shí)和準(zhǔn)確性驗(yàn)證。相比于自驗(yàn)證，交叉驗(yàn)證的結(jié)果具有更強(qiáng)的通用性和魯棒性。為了更好地衡量辨識(shí)得到的模型對(duì)位置等狀態(tài)信息的魯棒性，綜合體現(xiàn)模型在實(shí)際應(yīng)用中的準(zhǔn)確性，使用交叉驗(yàn)證的方法對(duì)模型準(zhǔn)確性進(jìn)行分析。對(duì)每組數(shù)據(jù)所得參數(shù)，使用相同環(huán)境下的其它4組數(shù)據(jù)進(jìn)行準(zhǔn)確性檢驗(yàn)，并取平均標(biāo)準(zhǔn)均方根誤差。實(shí)驗(yàn)結(jié)果如表2所示。

表2 環(huán)境模型交叉驗(yàn)證誤差Table 2 Cross-validation error of different models

由實(shí)驗(yàn)結(jié)果可見(jiàn)，HC模型在交叉驗(yàn)證中仍具有較小的預(yù)測(cè)誤差，體現(xiàn)出針對(duì)期望的實(shí)際物理特征具有更高的準(zhǔn)確性。但是相比于自驗(yàn)證誤差，其在交叉驗(yàn)證中相對(duì)于其他模型的優(yōu)勢(shì)有所降低。尤其是在與海綿的接觸實(shí)驗(yàn)中，誤差明顯超過(guò)了K和KV模型，其原因是包括海綿在內(nèi)的低剛度柔性環(huán)境的接觸動(dòng)力學(xué)特征具有較強(qiáng)的非線性，在辨識(shí)過(guò)程中，擬合能力較強(qiáng)的HC模型容易產(chǎn)生過(guò)擬合現(xiàn)象。過(guò)擬合的模型雖然在當(dāng)前數(shù)據(jù)下具有較小的誤差，但是整體模型準(zhǔn)確性卻明顯降低。

總體而言，相比于K模型，KV模型在準(zhǔn)確性上僅有微小提高，但是由于參數(shù)數(shù)量卻為K模型的2倍，所以在存儲(chǔ)空間、計(jì)算能力和通信帶寬等方面幾乎具有雙倍的占用。HC模型雖然具有更高的復(fù)雜度和資源占用，但是對(duì)中等剛度的環(huán)境的建模具有明顯提高的準(zhǔn)確性。同時(shí)，HC模型對(duì)低剛度柔性環(huán)境建模效果較差，對(duì)高剛度環(huán)境的建模沒(méi)有改善。

3 環(huán)境參數(shù)在線辨識(shí)

在環(huán)境模型建模的基礎(chǔ)上，模型參數(shù)決定了模型準(zhǔn)確性。由于往往不能事先獲得準(zhǔn)確的模型參數(shù)，需要在實(shí)際任務(wù)過(guò)程中，根據(jù)實(shí)時(shí)傳感信息，對(duì)環(huán)境參數(shù)進(jìn)行在線辨識(shí)。從參數(shù)收斂性和對(duì)模型變化的快速響應(yīng)特性等方面對(duì)最新的參數(shù)辨識(shí)算法進(jìn)行比較分析，首次使用交叉驗(yàn)證的方法衡量不同參數(shù)辨識(shí)算法在在線參數(shù)辨識(shí)過(guò)程中對(duì)模型準(zhǔn)確性的影響。得到參數(shù)辨識(shí)過(guò)程中模型準(zhǔn)確性的變化過(guò)程，更好地衡量參數(shù)辨識(shí)算法的實(shí)際效果。

3.1 環(huán)境參數(shù)的在線參數(shù)辨識(shí)方法

BLS算法取最近一段時(shí)間內(nèi)的數(shù)據(jù)求解具有最小二次誤差的參數(shù)。該方法的主要優(yōu)點(diǎn)在于簡(jiǎn)單有效，但缺點(diǎn)在于難以取得合適的窗口大小W，較大的W收斂性較好，但是對(duì)變化系統(tǒng)的響應(yīng)較慢，較小的W雖然能迅速反映系統(tǒng)變化，但是容易波動(dòng)劇烈，收斂性較差。

基于梯度下降的算法利用目標(biāo)函數(shù)的梯度進(jìn)行迭代計(jì)算來(lái)搜索最優(yōu)參數(shù)。在BLMS算法中，目標(biāo)函數(shù)定義為預(yù)測(cè)誤差的二次函數(shù)如式（7）［13］：

未知參數(shù)的迭代公式如式（8）：

其中，λk為迭代步長(zhǎng)，為與梯度方向相反的單位向量，

EWRLS是環(huán)境參數(shù)辨識(shí)中一種常用的參數(shù)辨識(shí)算法，其迭代公式如式（9）所示［7］：

其中，P為協(xié)方差矩陣，η為遺忘因子。當(dāng)η＝1時(shí)，算法等效于標(biāo)準(zhǔn)的遞歸最小二乘法。通常取η為小于但接近1的數(shù)值來(lái)使算法對(duì)新數(shù)據(jù)賦予更高的權(quán)值。

上述參數(shù)辨識(shí)算法的缺點(diǎn)在于通常難以兼具收斂性和跟蹤性，即在實(shí)際環(huán)境不變時(shí)得到收斂的辨識(shí)結(jié)果，并在環(huán)境發(fā)生變化時(shí)及時(shí)反應(yīng)環(huán)境變化。

SPRLS的核心思想是在在線辨識(shí)的過(guò)程中檢測(cè)前向預(yù)測(cè)誤差，當(dāng)前向預(yù)測(cè)誤差超出預(yù)設(shè)閾值時(shí)，認(rèn)為待辨識(shí)系統(tǒng)發(fā)生了變化，并按前向預(yù)測(cè)誤差大小增大協(xié)方差矩陣，提高算法對(duì)新數(shù)據(jù)的敏感性，使算法具有快速響應(yīng)特性。同時(shí)，當(dāng)前向預(yù)測(cè)誤差預(yù)設(shè)閾值時(shí)，算法等效于標(biāo)準(zhǔn)的遞歸最小二乘法，具有良好的收斂性。因此，SPRLS兼具良好的收斂性和對(duì)變化系統(tǒng)的快速跟蹤性能。算法遞歸方程如式（10）所示［12］：

其中β和γ分別為擾動(dòng)增益和誤差增益，取整函數(shù)NINT（）的表達(dá)式如式（11）所示：

3.2 環(huán)境參數(shù)的在線參數(shù)辨識(shí)實(shí)驗(yàn)分析

對(duì)2.2節(jié)中設(shè)計(jì)的4種接觸環(huán)境各取2000個(gè)連續(xù)采樣數(shù)據(jù)，在K模型上應(yīng)用上述方法進(jìn)行參數(shù)辨識(shí)實(shí)驗(yàn)。使用試錯(cuò)法人工調(diào)節(jié)各算法的參數(shù)設(shè)置至最好的效果，實(shí)驗(yàn)結(jié)果如圖5所示。

圖5 辨識(shí)算法比較圖Fig.5 Comparison of different estimation algorithms

由圖5（a）可見(jiàn)，4種辨識(shí)算法總體上都能夠辨識(shí)出環(huán)境剛度并能及時(shí)反映環(huán)境的變化，但在收斂性和跟蹤速度上有差異。圖5中（b）為環(huán)境B切換到環(huán)境C的過(guò)渡階段的放大圖?？梢?jiàn)，SPRLS和BLMS算法對(duì)采樣編號(hào)為4000時(shí)發(fā)生的環(huán)境變化具有較快的反應(yīng)速度，而EWRLS和BLS則響應(yīng)較慢。圖5中（c）為環(huán)境C的平穩(wěn)階段的放大圖，可見(jiàn)，只有SPRLS表現(xiàn)出了良好的收斂性，其余3種算法都有一定程度的波動(dòng)?？梢?jiàn)，相比于其余算法，SPRLS兼顧了良好的參數(shù)收斂性和對(duì)環(huán)境變化的快速跟蹤特性。在此基礎(chǔ)上，使用交叉驗(yàn)證的方式衡量參數(shù)辨識(shí)過(guò)程中模型準(zhǔn)確性的變化。即使用4種接觸環(huán)境的五組采樣數(shù)據(jù)中的一組數(shù)據(jù)進(jìn)行參數(shù)辨識(shí)，然后取參數(shù)辨識(shí)過(guò)程中各個(gè)時(shí)刻的參數(shù)辨識(shí)值，并使用其余四組數(shù)據(jù)進(jìn)行交叉驗(yàn)證計(jì)算相應(yīng)的誤差。實(shí)驗(yàn)結(jié)果如圖6所示。

圖6中，各條曲線在采樣編號(hào)為2000、4000和6000處的尖峰表明環(huán)境的突變導(dǎo)致環(huán)境模型的誤差突然增大，但是各個(gè)參數(shù)辨識(shí)算法都能夠通過(guò)在線辨識(shí)修正模型參數(shù)，有效降低模型誤差，即獲得了較高的模型準(zhǔn)確性。其中，SPRLS算法在環(huán)境突變后既能夠較快地降低誤差，又能在誤差降低后保持很好的穩(wěn)定性。各個(gè)算法的模型誤差曲線與圖5中的模型參數(shù)變化規(guī)律相一致。SPRLS算法在采樣編號(hào)為2000的環(huán)境突變過(guò)程中保持了一小段時(shí)間的較高誤差之后才迅速降低誤差，其原因是該時(shí)刻主端位置較小，環(huán)境突變并沒(méi)有使預(yù)測(cè)力產(chǎn)生明顯誤差，沒(méi)有觸發(fā)SPRLS中的自擾動(dòng)項(xiàng)。它表明SPRLS算法的快速跟蹤特性依賴于明顯的前向預(yù)測(cè)誤差，在缺乏特定輸入條件的情況下可能會(huì)導(dǎo)致參數(shù)收斂速度過(guò)慢。

4 結(jié)論

本文通過(guò)理論分析和實(shí)驗(yàn)驗(yàn)證對(duì)力反饋遙操作系統(tǒng)中環(huán)境動(dòng)力學(xué)模型的建模和在線參數(shù)辨識(shí)方法進(jìn)行了深入的比較分析，得到如下結(jié)論：

1）交叉驗(yàn)證的結(jié)果說(shuō)明：純剛度模型能以較低的復(fù)雜度得到一般的建模精度，Kelvin-Voigt模型相比于純剛度模型并沒(méi)有明顯優(yōu)勢(shì)；Hunt-Crossley模型以高復(fù)雜度獲得了更高的建模精度，適宜在數(shù)據(jù)充分的離線條件下采用。

2）在參數(shù)辨識(shí)算法方面，相比于同類算法，自擾動(dòng)遞歸最小二乘法在參數(shù)辨識(shí)的收斂性和對(duì)時(shí)變環(huán)境的跟蹤性上都有較好的表現(xiàn)，是實(shí)際應(yīng)用中的首選方法，但是需要注意其快速跟蹤特性對(duì)較大前向誤差的依賴問(wèn)題。

圖6 模型準(zhǔn)確性變化圖Fig.6 Model accuracy profile during online estimation

（References）

［1］ Xu X， Cizmeci B， Al-Nuaimi A， et al.Point cloud-based model-mediated teleoperation with dynamic andperceptionbased model updating［J］.Instrumentation and Measurement，IEEE Transactions on， 2014， 63（11）： 2558-2569.

［2］ Liu H J， Aiguo S.Virtual-environment modeling and correction for force-reflecting teleoperation with time delay［J］.Industrial Electronics， IEEE Transactions on， 2007， 54（2）：1227-1233.

［3］ Velanas S V，Tzafestas C S.Model-mediated telehaptic perception of delayed curvature［C］／／RO-MAN， 2012 IEEE.IEEE，2012：941-947.

［4］ Willaert B， Bohg J， Van Brussel H， et al.Towards multi-DOF model mediated teleoperation：using vision to augment feedback［C］ ／／Haptic Audio Visual Environments and Games（HAVE）， 2012 IEEE International Workshop on.IEEE，2012：25-31.

［5］ 徐效農(nóng)，宋愛(ài)國(guó)，朱澄澄，等.基于環(huán)境建模與修正的視覺(jué)／力覺(jué)輔助遙操作系統(tǒng)［J］.載人航天，2016，22（1）：55-61.XU Xioanong， SONG Aiguo， ZHU Chengcheng， et al.Visualhaptic aid teleoperation system based on environment modeling and updating［J］.Manned Spaceflight， 2016， 22（1）：55-61.（in Chinese）

［6］ Haddadi A，Hashtrudi-Zaad K.Online contact impedance identification for robotic systems［C］／／Haddadi A.2008 IEEE／RSJ International Conference on Intelligent Robots and Systems.Piscataway， NJ： IEEE， 2008：974-980.

［7］ Haddadi A， Hashtrudi-Zaad K.Real-time identification of Hunt-Crossley dynamic models of contactenvironments［J］.Robotics， IEEE Transactions on， 2012， 28（3）： 555-566.

［8］ Yamamoto T，Vagvolgyi B，Balaji K，et al.Tissue property estimation and graphical display for teleoperatedrobot-assisted surgery［C］ ／／Yamamoto T.IEEE International Conference on Robotics and Automation， 2009.Piscataway， NJ： IEEE，2009：4239-4245.

［9］ Haddadi A，Hashtrudi-Zaad K.A new fast online identification method for linear time-varying systems［C］ ／／Haddadi A.2008 American Control Conference.Piscataway， NJ： IEEE，2008：1322-1328.

［10］ Love L J， Book W J.Environment estimation for enhanced impedance control［C］／／Love L J.1995 IEEEInternational Conference on Robotics and Automation.Piscataway， NJ：IEEE，1995，2：1854-1859.

［11］ Colton M B， Hollerbach J M.Identification of nonlinear passive devices for haptic simulations［C］／／Colton M B.2005 First Joint Eurohaptics Conference and Symposium on Haptic Interfaces for Virtual Environment andTeleoperator Systems.］ ／／Piscataway， NJ： IEEE， 2005： 363-368.

［12］ Achhammer A， Weber C， Peer A， et al.Improvement of model-mediated teleoperation using a new hybridenvironment estimation technique［C］ ／／Achhammer A.2010 IEEE International Conference on Robotics andAutomation.Piscataway，NJ： IEEE，2010： 5358-5363.

［13］ Haddadi A.Stability， performance， and implementation issues in bilateral teleoperation control and haptic simulation systems［J］.2012.

［14］ Yamamoto T， Bernhardt M， Peer A， et al.Techniques for environment parameterestimation duringtelemanipulation［C］／／IEEE Ras＆ Embs International Conference on Biomedical Robotics and Biomechatronics.IEEE， 2008：217-223.

（責(zé)任編輯：康金蘭）

Online Environment Parameter Estimation for Force Feedback Teleoperation System

DING Yukun1， SONG Jingzhou1，SHANG Zhihao1，HAN Liangliang2
（1.Automation School in Beijing University of Posts and Telecommunications， Beijing 100876， China；2.Aerospace System Engineering Shanghai， Shanghai 201108， China）

V11

A

1674-5825（2017）05-0704-07

2017-03-06；

2017-08-05

載人航天預(yù)先研究項(xiàng)目（030401）

丁宇堃，男，碩士研究生，研究方向?yàn)闄C(jī)器人遙操作。E-mail：dingyukun1994＠126.com