【摘 要】 目前,對(duì)于針灸手法的學(xué)習(xí)評(píng)價(jià),缺少針刺時(shí)施力大小、時(shí)長(zhǎng)和方向等可以進(jìn)行量化的物理參數(shù)。本文針對(duì)采集針刺過(guò)程中“捻轉(zhuǎn)”和“提插”動(dòng)作的觸覺(jué)參數(shù),設(shè)計(jì)了一種陣列式聚偏二氟乙烯膜(PVDF)觸覺(jué)傳感器;隨后,給出了面向觸覺(jué)信號(hào)的窗口分割方法,并提取窗口中的時(shí)域特征;最后,構(gòu)建了一種基于模糊 C均值聚類法(FCM)的針刺手法識(shí)別方法。通過(guò)實(shí)驗(yàn)證明該方法能夠有效識(shí)別捻轉(zhuǎn)補(bǔ)、捻轉(zhuǎn)瀉、提插補(bǔ)、提插瀉四種基本針刺動(dòng)作,有利于針刺手法的量化與傳播。
本文源自茍升異; 宿翀; 王磊; 趙亞楠; 陳捷, 針刺研究 發(fā)表時(shí)間:2021-06-25
【關(guān)鍵詞】 聚偏二氟乙烯膜傳感器;機(jī)器學(xué)習(xí);針刺手法分類
針刺作為一種傳統(tǒng)中醫(yī)療法,其教學(xué)大多是依靠經(jīng)驗(yàn)和課堂講解,這對(duì)于其傳承和推廣存在一定的制約。對(duì)于針刺的基本手法如捻轉(zhuǎn)補(bǔ)、捻轉(zhuǎn)瀉、提插補(bǔ)和提插瀉,由于個(gè)人理解的不同,同一針刺操作的時(shí)間、力量等方面均存在差異,從而在臨床療效上產(chǎn)生不同的結(jié)果,因此對(duì)針刺手法的量化研究就顯得極為重要。本文通過(guò)對(duì)針刺過(guò)程觸覺(jué)信息的采集,提取出針刺動(dòng)作的時(shí)域特征,完成對(duì)針刺操作的分類。量化的針刺手法將為針灸的診療和教學(xué)提供強(qiáng)有力的科學(xué)依據(jù),有助于針刺手法的推廣與傳承。
針刺手法量學(xué)的研究對(duì)針刺時(shí)施力大小、方向、時(shí)長(zhǎng)、頻率等要素進(jìn)行了界定。以捻轉(zhuǎn)和提插手法為基礎(chǔ),形成了捻轉(zhuǎn)補(bǔ)、捻轉(zhuǎn)瀉、提插補(bǔ)、提插瀉4種基本針刺手法[1]。為了提高針刺手法的規(guī)范性、可操作性,楊華元等[2]開(kāi)發(fā)了針刺手法參數(shù)測(cè)定儀,用其采集到的針刺手法仿真過(guò)程中的物理參數(shù)進(jìn)行分析[3]。王彩虹等[4]利用針刺手法測(cè)定儀采集并分析不同手法的參數(shù)和波形,進(jìn)而分析捻轉(zhuǎn)和提插兩種手法。丁光宏等[5]將基于光刻應(yīng)變硅片和壓電陶瓷的傳感器集成到針體上,實(shí)現(xiàn)了針刺過(guò)程中對(duì)針上作用力的測(cè)量。顧星等[6]開(kāi)發(fā)了針刺手法教學(xué)測(cè)量?jī)x,通過(guò)電磁感應(yīng)技術(shù)采集針刺過(guò)程信號(hào)。然而,上述針刺手法物理參數(shù)的提取系統(tǒng),大多架設(shè)于針柄上,鮮有從人體觸覺(jué)角度探討手法物理參數(shù)獲取與分析的方法。
聚 偏 二 氟 乙 烯 膜 (Polyvinylidenefluoride,PVDF)作為一種具有良好性能的柔性壓電材料,已經(jīng)被廣泛應(yīng)用于構(gòu)建柔性壓力傳感器[7]、心率和呼吸實(shí)時(shí)監(jiān)測(cè)傳感系統(tǒng)[8]、機(jī)械模具制造中的應(yīng)力分析等[9-10]各個(gè)方面的研究。近年來(lái),針對(duì) PVDF 應(yīng)用于人體動(dòng)作識(shí)別方面的 探 索,如 人 體 手 部 分割[11]、人體肌電信號(hào)提取與分類[12]等,逐步成為智能傳感系統(tǒng)研究的熱點(diǎn)。此外,諸如神經(jīng)網(wǎng)絡(luò)、聚類等機(jī)器學(xué)習(xí) 方 法[13]也 被 廣 泛 應(yīng) 用 于 人 體 連 續(xù) 動(dòng) 作的分類中,獲得豐碩的研究成果。
受上述研究啟發(fā),本研究設(shè)計(jì)了一款基于陣列式PVDF的觸覺(jué)傳感器,在獲取針灸醫(yī)師針刺過(guò)程中的觸覺(jué)壓電信號(hào)的同時(shí),構(gòu)建一種基于機(jī)器學(xué)習(xí)中模糊 C 均值聚類法(FuzzyC-Means,F(xiàn)CM)的 針刺手法識(shí)別方法。
1 針刺手法觸覺(jué)傳感器構(gòu)建及智能采集系統(tǒng)設(shè)計(jì)
PVDF壓 電 薄 膜 是 一 種 獨(dú) 特 的 高 分 子 傳 感 材料,外觀為半透明狀,具有分子鏈間排列緊密、不易燃、質(zhì)量輕、柔韌性好等特性;該材料電壓靈敏度高,能根據(jù)壓力的變化輸出電壓信號(hào),可加工成為高效可靠、低成本的柔性觸覺(jué)傳感器。
1.1 陣列式 PVDF觸覺(jué)傳感器的制備
區(qū)別于傳統(tǒng)單片 PVDF薄膜傳感器,陣列式傳感器采用多片 PVDF薄膜來(lái)構(gòu)造,其優(yōu)勢(shì)在于可以在同一時(shí)刻采集到不同位置的 PVDF壓電信號(hào),從而提取其時(shí)域特征。在進(jìn)行針刺時(shí),伴隨著指尖的相對(duì)移動(dòng),采用陣列式傳感器的設(shè)計(jì),可以增加有效的時(shí)域特征,為提升針刺手法分類的準(zhǔn)確率提供一種可行的途徑。本研究利用兩片尺寸相同、形狀完好的 PVDF 薄 膜,經(jīng) 過(guò) 形 狀 切 割、邊 緣 絕 緣 后 進(jìn) 行組裝同時(shí)引出導(dǎo)線,形成一個(gè)陣列式 PVDF壓電傳感器。(1)形 狀 切 割:分 別 在 兩 片 PVDF 成 品 上 畫(huà)出直徑為8 mm 的 圓,然 后 按 標(biāo) 記 線 條 進(jìn) 行 裁 剪。(2)邊緣絕緣:選用乙醇對(duì)薄膜邊緣進(jìn)行非金屬化的具體處理。(3)組裝引出導(dǎo)線:準(zhǔn)備透明塑料膜,在一片上涂抹少量導(dǎo)電銀膠,然后在導(dǎo)電銀膠中放入導(dǎo)線作為壓電傳感器的一極,再將 PVDF薄膜放在上面,按壓使其粘合,直至導(dǎo)電銀膠凝固。在 PVDF薄膜相對(duì)稱的另一邊同樣放置導(dǎo)電銀膠及導(dǎo)線作為壓電傳感器的另一極,再將透明薄膜放在導(dǎo)電銀膠之上并進(jìn)行按壓,直到最后凝固,形成一個(gè)“兩片分離式”的陣列式 PVDF壓電傳感器。見(jiàn)圖1。
1.2 針刺手法觸覺(jué)壓電信號(hào)采集系統(tǒng)
將圖1所示 PVDF壓電傳感器進(jìn)行“封裝”,形成可穿戴的“指套”觸覺(jué)傳感器。醫(yī)師的手指(一般為食指)固定住觸覺(jué)傳感器,在針刺過(guò)程中,手指對(duì)針體施加不同的壓力使指套產(chǎn)生形變,產(chǎn)生的壓電信號(hào)由放大器放大、采集模塊收集并發(fā)送至電腦保存,此時(shí)便完成了相應(yīng)數(shù)據(jù)的采集,如圖2所示,從左至右 依 次 為 PVDF 傳 感 器、信 號(hào) 放 大 器、信 號(hào) 采集模塊和電腦壓電曲線顯示。
2 針刺手法數(shù)據(jù)的采集和分類
2.1 針刺手法數(shù)據(jù)采集
如圖2所示,該系統(tǒng)采樣頻 率 設(shè) 置 為100 Hz,可以充分地捕獲針刺手法動(dòng)作所產(chǎn)生的數(shù)據(jù),數(shù)據(jù)的采集由醫(yī)師完成。將傳感器用醫(yī)用膠布固定在醫(yī)師的食指,實(shí)驗(yàn)中由醫(yī)師分別執(zhí)行捻轉(zhuǎn)補(bǔ)、捻轉(zhuǎn)瀉、提插補(bǔ)、提插瀉4個(gè)動(dòng)作。每個(gè)動(dòng)作重復(fù)2次,每次的持續(xù) 時(shí) 間 為80s。由 PVDF 壓 電 薄 膜 傳 感 器 采集到的壓電信號(hào)如圖3所示。
2.2 噪聲處理
移動(dòng)平均濾波器可以有效地減少隨機(jī)干擾的影響。本研究選 擇6 階 移 動(dòng) 平 均 濾 波 算 法 來(lái) 消 除 噪聲,去除信號(hào)毛刺。傳感器采集的原始加速度信號(hào)存在較多噪聲,在使用6階移動(dòng)平均濾波算法后,數(shù)據(jù)變得較為平滑。
2.3 數(shù)據(jù)窗口分割
在得到壓電數(shù)據(jù)后,無(wú)法對(duì)原始數(shù)據(jù)直接進(jìn)行分類。窗口分割的目的是將得到的壓電數(shù)據(jù)進(jìn)行劃分,將壓電數(shù)據(jù)流分割為若干的窗口,然后以一個(gè)窗口為基礎(chǔ),進(jìn)行特征的提取,只有有效提取出窗口的特征才能進(jìn)行分析,完成針刺動(dòng)作的分類。
在動(dòng)作的識(shí)別階段,需要對(duì)傳感器壓電信號(hào)數(shù)據(jù)流進(jìn)行數(shù)據(jù)窗口分割。常用的分割技術(shù)包括基于滑動(dòng)窗口的分割、基于事件定義的窗口分割以及基于動(dòng)作定義的窗口分割等3種方法[14-15]。滑動(dòng)窗口分割技術(shù)是指采用固定長(zhǎng)度的窗口對(duì)數(shù)據(jù)流進(jìn)行分割。滑動(dòng)窗口分割技術(shù)由于操作簡(jiǎn)單,因此被多數(shù)研究所采 用。本文中采用滑動(dòng)窗口分割的時(shí)間為1s。
2.4 特征提取
在濾除傳感器數(shù)據(jù)的噪聲后,用定長(zhǎng)度的滑動(dòng)窗口分割技術(shù)對(duì)傳感器數(shù)據(jù)流進(jìn)行劃分,窗口的長(zhǎng)度設(shè)置為1s。對(duì)于每個(gè)時(shí)間窗口,它所包含的數(shù)據(jù)可以用一個(gè) Nx1維的向量s=[s1,s2,...,sn]來(lái) 表示。由于本文 采 用“兩 片 分 離 式”陣 列 式 PVDF 壓電傳感器進(jìn)行針灸醫(yī)師手指觸覺(jué)信號(hào)的提取,因此從時(shí)域特征角度進(jìn)行分析,具體的特征類型如表1所示。
表1中所示特征的計(jì)算中,峰的位置信息反應(yīng)針體的位置,最大值反應(yīng)針灸動(dòng)作強(qiáng)度變化的幅度,均值反應(yīng)針灸動(dòng)作的平均強(qiáng)度;積分反應(yīng)針體旋轉(zhuǎn)與緊握程度,時(shí)間長(zhǎng)度反應(yīng)頻率大小。由于壓電信號(hào)無(wú)法直接獲取針體位置,通過(guò)對(duì)壓電信號(hào)的進(jìn)一步處理得到以下2個(gè)特征。特征ID9和10中所謂的“針體位置”,是指在一個(gè)窗口壓電數(shù)據(jù)的采集中,針體具體在哪個(gè) PVDF薄膜上滾動(dòng)。為獲得“針體位置”,我們將兩個(gè)薄膜上信號(hào)取最大值,最大值所在的薄膜序號(hào),即為針體所在位置。進(jìn)而求取一個(gè)壓電數(shù)據(jù) 窗 口 內(nèi) 針 體 位 置 的 平 均 值 和 方 差。公 式如下:
均值特征:!= 1N∑ Ni=1si (1)方差特征:σ2 = 1N∑ Ni=1(si - !)2 (2)采樣區(qū)間內(nèi)最大值:Smax = max (s) (3)采樣區(qū)間內(nèi)最小值:Smin = min (s) (4)采樣區(qū)間內(nèi)積分值:SS =∫ TnT1Stdt (5)
2.5 基于 FCM 的針刺手法分類器
本研究采用模糊C均值聚類算法對(duì)集合s進(jìn)行處理,即將大規(guī)模的訓(xùn)練樣本歸為數(shù)量相對(duì)較少的若干聚類,并基于聚類中心設(shè)計(jì)分類器。基于模糊C均值得到聚類中心:1)確定聚類數(shù)量 k和模糊參數(shù)b,設(shè)置迭代終止閾值ε;2)初始化每個(gè)聚類中心 m;3)計(jì)算各個(gè)聚類的隸屬度函數(shù)值;4)重新計(jì)算訓(xùn)練集的各類聚類中心;5)計(jì)算聚類損失函數(shù)值,如果與上次的函數(shù)值之差小于閾值,則退出迭代,否則返回步驟3重 復(fù) 運(yùn) 算,直 到 獲 得 穩(wěn) 定 聚 類 中 心 和 隸 屬度值。
每種針刺手法,例如捻轉(zhuǎn)補(bǔ),在得到針刺動(dòng)作壓電數(shù)據(jù)后對(duì)這些數(shù)據(jù)進(jìn)行窗口劃分,按照1s1個(gè)窗口,即100個(gè)采樣點(diǎn)1個(gè)窗口,但是這樣的滑動(dòng)窗口分割方法,無(wú)法保證這個(gè)窗口內(nèi)的壓電數(shù)據(jù)恰好反映了1個(gè)完整的針刺手法動(dòng)作。每次提取到的窗口數(shù)據(jù)所蘊(yùn)含的手法信息有如下6種可能:(1)1個(gè)完整的針刺動(dòng)作;(2)1個(gè)針刺動(dòng)作的前面部分;(3)1個(gè)針刺動(dòng)作的后面部分;(4)1個(gè)針刺動(dòng)作的后面部分和1個(gè)完整的針刺動(dòng)作;(5)1個(gè)完整的針刺動(dòng)作和1個(gè)針刺動(dòng)作的前面部分;(6)1個(gè)針刺動(dòng)作的后面部分和1個(gè)完整的針刺動(dòng)作和一個(gè)針刺動(dòng)作的前面部分。對(duì)于1個(gè)動(dòng)作,就會(huì)有這樣的6種情況的窗口,如果1個(gè)動(dòng)作只設(shè)立1個(gè)聚類中心,不符合窗口分割的實(shí)際情況,會(huì)影響分類的準(zhǔn)確率。這6種窗口,每種都有自己的聚類中心,所以1個(gè)動(dòng)作就有6個(gè)聚類中心,本文探討了4種針刺手法動(dòng)作,因此共有24 個(gè) 聚 類 中 心。 本 文 所 用 聚 類 模 型,如 圖 4所示。
在使用訓(xùn)練集訓(xùn)練得到24個(gè)聚類中心后,將待分類數(shù)據(jù),分別與24個(gè)聚類中心計(jì)算歐氏距離,然后同一類的6個(gè)聚類中心的歐氏距離相加,就得到4個(gè)歐氏距離之和,最小歐氏距離的對(duì)應(yīng)類別就 是分類結(jié)果。
3 實(shí)驗(yàn)設(shè)計(jì)與實(shí)驗(yàn)結(jié)果
3.1 實(shí)驗(yàn)設(shè)計(jì)與實(shí)驗(yàn)結(jié)果
本實(shí)驗(yàn)中,使用中國(guó)中醫(yī)科學(xué)院針灸研究所兩位醫(yī)師在佩戴本文中觸覺(jué)傳感器情況下進(jìn)行針刺操作中收集到的數(shù)據(jù)。要求醫(yī)師針對(duì)某個(gè)針刺手法動(dòng)作(例如:捻轉(zhuǎn)補(bǔ)、捻轉(zhuǎn)瀉、提插補(bǔ)和提插瀉)連續(xù)操作5min,由于本文提出的手法采集系統(tǒng)采樣頻率為100Hz,因此醫(yī)師的每種動(dòng)作可以采集到30000個(gè)采樣點(diǎn)。由于采樣初始和結(jié)束階段,存在粗大誤差,所以在初始和結(jié)束階段各剔除1000個(gè)采樣點(diǎn),最終保留28000個(gè)采樣數(shù)據(jù)。我們將此時(shí)采集成功的某個(gè)手法 動(dòng) 作 的24000個(gè) 采 樣 點(diǎn) 作 為 訓(xùn) 練 數(shù)據(jù),4000個(gè)采樣點(diǎn)作為測(cè)試數(shù)據(jù)。
循環(huán)執(zhí)行上述的采集動(dòng)作4次(針對(duì)4種基礎(chǔ)針刺手法動(dòng)作),其中:我們?cè)O(shè)置100個(gè)采樣點(diǎn)為一個(gè)窗口,對(duì)于訓(xùn)練數(shù)據(jù),每種手法動(dòng)作有240個(gè)窗口(共4種手 法 動(dòng) 作,總 計(jì) 960 個(gè) 窗 口);對(duì) 于 測(cè) 試 數(shù)據(jù),每種手法動(dòng)作40個(gè)窗口(共4種手法動(dòng)作,總計(jì)160個(gè)窗口)。
基于本文所提手法動(dòng)作的識(shí)別方法,得到針對(duì)測(cè)試集的識(shí)別準(zhǔn)確率如圖5所示。從圖中可以看出每一種針刺動(dòng)作的識(shí) 別情況,平 均 識(shí) 別 正 確 率 在88.22%,其中捻轉(zhuǎn)補(bǔ)與捻轉(zhuǎn)瀉的正確率較 高,均 在90%以上,而提插補(bǔ)、提插瀉則在80%左右。
以第一行第一列為例:醫(yī)師實(shí)際手法動(dòng)作為捻轉(zhuǎn)補(bǔ),且本文分類方法將該實(shí)際手法獲取的壓電數(shù)據(jù)分類為捻轉(zhuǎn)補(bǔ)的準(zhǔn)確率為92.5%;以第一行第二列為例,醫(yī)師實(shí)際手法動(dòng)作為捻轉(zhuǎn)補(bǔ),但是本文分類方法將該實(shí)際手法獲取的壓電數(shù)據(jù)分類為捻轉(zhuǎn)瀉的有7.5%。
圖6至 圖10是 訓(xùn) 練 樣 本 數(shù) 據(jù) 通 過(guò) FCM 聚 類后得到的數(shù)據(jù)分布情況。
捻轉(zhuǎn)補(bǔ)和捻 轉(zhuǎn) 瀉 的 壓 電 平 均 值 更 靠 近0,且 兩者 不易進(jìn)行區(qū)分,而提插補(bǔ)和提插瀉表現(xiàn)的更加分散,其中提插補(bǔ)前端的平均值比較集中,波動(dòng)不大,但是后端的表現(xiàn)就出現(xiàn)了較大的波動(dòng),從-1V 到1V都有分布。提 插 瀉 則 相 反,平 均 值 的 分 布 情 況在前端更加分散,而后端在0上下波動(dòng)。見(jiàn)圖6。
捻轉(zhuǎn)補(bǔ)和捻轉(zhuǎn)瀉的動(dòng)作壓電積分值更小,這是因?yàn)槟磙D(zhuǎn)動(dòng)作是一種滑動(dòng)動(dòng)作,在捻轉(zhuǎn)過(guò)程中,壓電的平均值較小,而且滑動(dòng)過(guò)程出現(xiàn)壓電信號(hào)的時(shí)間占比相對(duì)于提插過(guò)程較小,從而導(dǎo)致捻轉(zhuǎn)動(dòng)作的積分值更小。提插補(bǔ)動(dòng)作在后端的積分值較大,前端較小。而提插瀉則相反,在前端的積分值較大,后端的積分值較小。見(jiàn)圖7。
捻轉(zhuǎn)補(bǔ)和捻轉(zhuǎn)瀉這兩者在最大值上要小于提插補(bǔ)和提插瀉,因?yàn)槟磙D(zhuǎn)是滑動(dòng)動(dòng)作,需要使用的力較小,但是提插需要持握針體進(jìn)行上下動(dòng)作,此時(shí)手指施加的力較大,反映到壓電信號(hào)就出現(xiàn)了數(shù)值較大的最大值。捻轉(zhuǎn)補(bǔ)和捻轉(zhuǎn)瀉區(qū)分度不大,而提插補(bǔ)在后端有著更大的壓電最大值,提插瀉則在前端的壓電最大值更大,后端較小。見(jiàn)圖8。
捻轉(zhuǎn)補(bǔ)和捻轉(zhuǎn)瀉的最小值還是比較集中,不易進(jìn)行區(qū)分,而提插補(bǔ)的最小值在后端較為明顯,前端集中在0附近,提插瀉的最小值在前端較為明顯,后端集中在0附近。通過(guò)最大值和最小值的考慮,這兩者是相互對(duì)應(yīng)的,這是由于 PVDF壓電薄膜的壓電特性導(dǎo)致的,正向壓電越大,從而轉(zhuǎn)的負(fù)向壓電也就越大。所以對(duì)于壓電信號(hào)的最大值和最小值,這兩者是統(tǒng)一的關(guān)系。見(jiàn)圖9。
捻轉(zhuǎn)補(bǔ)的針體位置數(shù)值偏大,在后端薄膜附近的區(qū)域時(shí)間更長(zhǎng),而捻轉(zhuǎn)瀉則與之相反。提插補(bǔ)和提插瀉針體的位置則比較分散。見(jiàn)圖10。
3.2 所提方法的局限性
本研究中以觸覺(jué)傳感器測(cè)試醫(yī)師提插補(bǔ)和提插瀉兩種手法的結(jié)果顯示,提插動(dòng)作總是會(huì)在陣列式PVDF觸覺(jué)傳 感 器 其 中 的 一 片 區(qū) 域 產(chǎn) 生 更 大 的 壓力,且發(fā)力更加突然,瞬間的壓力更大,積分也就很大。與捻轉(zhuǎn)補(bǔ)與捻轉(zhuǎn)瀉兩種手法的識(shí)別率高達(dá)95%相比,本文所用方法對(duì)提插補(bǔ)和提插瀉兩種手法的識(shí)別準(zhǔn)確率還稍顯遜色。未來(lái),將嘗試引入頻域特征及能量熵特征,來(lái)衡量提插手法下不同操作的特征,以提升識(shí)別準(zhǔn)確率,也在理論上闡釋提插手法的作用機(jī)制。
4 結(jié)論
本研究構(gòu)建了一種陣列式 PVDF觸覺(jué)傳感器,可以有效地測(cè)量出針刺過(guò)程中指尖不同位置的壓電信號(hào),將力學(xué)信號(hào)轉(zhuǎn)化為電學(xué)信號(hào)進(jìn)行分析。隨后,本研究基于機(jī)器學(xué)習(xí)中的 FCM 方法對(duì)獲取的針刺手法壓電信號(hào)進(jìn)行聚類分析,構(gòu)建了一種針刺手法識(shí)別方法。經(jīng)過(guò)實(shí)驗(yàn)測(cè)試,此方法在捻轉(zhuǎn)補(bǔ)和捻轉(zhuǎn)瀉的動(dòng)作分類識(shí)別中有著 較高的準(zhǔn)確率,可 達(dá)95%;對(duì) 于 提 插 補(bǔ) 和 提 插 瀉 來(lái) 說(shuō),分 類 準(zhǔn) 確 率 可 達(dá)80%,還有較大的提升空間。在后續(xù)的工作中,計(jì)劃通過(guò)增加信號(hào)特征提取的種類并且構(gòu)建多特征融合下的識(shí)別方法,進(jìn)一步提高識(shí)別準(zhǔn)確率。該結(jié)果利用智能傳感 系 統(tǒng) 和 機(jī) 器 學(xué) 習(xí) 等 信 息 科 學(xué) 技 術(shù) 的 方法,實(shí)現(xiàn)了對(duì)針刺動(dòng)作特征的提取與分類,為針刺手法的研究提供科學(xué)依據(jù),對(duì)針刺手法的量化有著重要意義。
論文指導(dǎo) >
SCI期刊推薦 >
論文常見(jiàn)問(wèn)題 >
SCI常見(jiàn)問(wèn)題 >