

















摘 要:針對低信噪比下神經網絡難以提取數字信號空間特征的問題,提出一種基于坐標注意力機制的數字信號識別方案。將8 種數字信號進行正交調制,根據其幅度、相位信息序列進行預編碼處理,在不同的訓練步長下,提取分析數字信號幅度和相位的關鍵特征,選取合適的神經網絡超參數,使網絡達到擬合面。坐標注意力機制將數字信號特征進行2 個一維特征編碼,分別沿縱向和橫向捕獲幅度和相位的遠程依賴關系;將生成的數字信號特征編碼為一對方向感知和位置敏感的權重系數,進行數字信號特征的重標定。仿真結果表明,8 種數字信號下,調制方式識別率高于95% 時,卷積神經網絡(Convolutional Neural Network,CNN) 中坐標注意力機制信噪比增益約為4 dB,殘差神經網絡中坐標注意力機制信噪比增益約為8 dB。坐標注意力機制取得了較高的識別率以及更好的信噪比增益,與通道注意力機制、空間注意力機制相比更適用于數字信號解調的應用。
關鍵詞:數字信號;調制識別;坐標注意力機制;權重系數
中圖分類號:TN911. 3 文獻標志碼:A 開放科學(資源服務)標識碼(OSID):
文章編號:1003-3106(2024)06-1398-09
0 引言
自動調制識別(Automatic Modulation Recognition,AMR)是在非坐標通信中從接收信號識別出調制類型的一種關鍵技術[1]。AMR 可以在電子戰、監視系統和威脅分析中截獲并識別對方的調制信號,并以該調制信號確定干擾信號的調制類型;同時在民用領域用于頻譜管理[2]。但在復雜的通信環境中,AMR的識別準確率因傳輸過程中的噪聲、多徑衰落等因素的影響,難以提升。傳統方法需要人工進行特征提取,存在著設計復雜度高、識別精度低等問題;在不同的通信系統中,分類的算法缺乏魯棒性[3]。
為了解決這方面問題,研究學者將神經網絡應用于AMR 中,并取得了更好的識別效果。文獻[4]研究了卷積神經網絡(Convolutional Neural Network,CNN)對無線電復值時域信號的適應性,證明在低信噪比情況下,深度CNN 對大型密集編碼時間序列進行盲學習的可行性。文獻[5]建立了新的數據集,拓展了通信信號的調制類型同時混雜了更多樣的信道損傷配置,以此設計了相應的CNN 和殘差網絡(Residual Network,ResNet)[6]結構,并可以有效地轉移到(Over The Air,OTA)數據集,識別精度損失大約為7% 。文獻[7]設計了一個8 層CNN 對信號的星座圖開展識別,獲取QAM 的調制階數,訓練好的網絡可以檢測到足夠的差異來區分16QAM 和64QAM。聯合參數調制識別方面,文獻[8]利用與樹結構相似的分層神經網絡在根節點分析接收信號是否連續,將不同性質的信號輸送到對應的網絡結構中,不同節點的網絡識別出調制方式和調制階數,最終將這些信息組合得出調制類型,但QAM 和8PSK 下的識別器并未完全收斂到最終的決策面。文獻[9]利用信號的異步延遲抽頭圖來訓練具有2 個隱藏層的深度神經網絡,實現了對調制速率和調制方式的聯合估計。文獻[10 ]通過平滑偽Wigner-Ville 分布和Born-jordan 分布變換將信號轉變為2 種時頻圖,并將其圖像特征和手工特征進行結合構成聯合特征。文獻[11]提出為每個神經網配置權重系數,通過壓縮感知增加權重系數的稀疏性,實現輕量化的AMR,但性能有所下降。因此,上述調制識別方案在復雜的網絡結構下能展現出良好的識別準確率,但計算復雜度也隨之提升。輕量化的神經網絡結構下,調制識別準確率還有一定的提升空間。
面對神經網絡無法有效提取不同調制類型信號的空間變換特征的問題,可將坐標注意力機制[12]與神經網絡相結合。本文將坐標注意力機制推廣到數字信號調制識別場景中,提出了基于坐標注意力機制的數字信號調制識別方案。首先,在不同的正交調制方式下,設置神經網絡的超參數、坐標注意力結構,計算不同坐標注意力機制下神經網絡基本單元的誤差項,利用適應性矩估計(Adaptive moment esti-mation,Adam)[13]優化算法,根據神經網絡基本單元的誤差項動態調整模型的權重和超參數,完成模型的自適應更新,從而選取較優的參數集使識別結果的精度最高并減小過擬合現象。其次,利用坐標注意力機制在空間維度上聚合幅度和相位特征,生成幅相通道描述符,根據識別誤差動態坐標注意力機制參數,進行幅度和相位特征的權重賦值并產生坐標注意力機制調制權值的集合,將其應用在特征映射上。最后,在不同網絡結構、步長以及不同調制數字信號下進行仿真分析。通過對數據集的訓練,取得了較高的識別率以及更好的信噪比增益,與通道注意力機制、空間注意力機制相比更適用于數字信號解調的應用。
1 數字信號調制識別系統模型
對寬帶接收信號完成信號檢測和初步的調制識別后,需要進一步對各信號的調制參數進行精確識別,以引導后續的解調譯碼和輻射源個體識別等處理過程。數字信號調制識別系統模型如圖1 所示。
在發射端,輸入信息i(t)經過信道編碼器,得到編碼后的信號c(t),隨后將其輸入調制器得到調制信號s(t)= S(c(t);θ),其中θ 為調制參數集合。s(t)經過多徑傳輸,路徑損耗和噪聲干擾后,到達接收端,接收信號為x(t):
式中:G 表示信道增益,fo 和θo 表示載波的頻率和相位偏移,s (τ)表示時間τ 處的發射信號樣本,p(·)表示脈沖整形,h(·)表示信道響應,εT 表示符號的定時誤差,dadd(t)表示額加性噪聲。
在接收端,調制識別的主要任務是針對接收信號x(t),設計識別算法以估計最佳的調制參數θ^。識別的調制參數θ^ 一方面可以用于對x(t)進行解調得到編碼信號c(t)的估計值c^(t),另一方面可以用于為輻射源個體識別模塊的信號預處理環節提供先驗信息。
2 注意力機制原理
2. 1 通道注意力機制原理
傳統的CNN 已經被證明可以出色地處理各種視覺任務[14-15]。卷積接收數據預處理后的二維數據的輸入,負責提取圖像的特征,池化層負責將卷積層提取的特征進行挑選,隨后全連接層將池化層的所有特征矩陣轉化成一維的特征向量,輸出模型學習到的結果。然而卷積只能構建局部關系,不能構建長期關系[16-17]。通道注意力機制是一種新的卷積處理方式[18],該方式通過引入自動學習,獲取特征圖每個通道的程度信息,然后將這個程度信息賦予每個特征圖一個權重值,從而師生神經網絡重點關注某些特征通道,增加對任務有用的特征通道。通道注意力機制原理如圖2所示。
圖2 中,W 和W′表示輸入特征寬度,H 和H′表示輸入特征高度,C 和C′表示輸入的通道數,Fsq 表示壓縮模塊,Fex 表示激勵模塊,Fscale 表示權重轉化。
通道注意力機制計算過程如下:
① Ftr:X→U,X∈CH′×W′×C′,U∈CH×W×C ,Ftr 視為卷積算子,V = [v1 ,v2 ,…,vc ]代表卷積核,Vc 代表卷積核參數,U = [u1 ,u2 ,…,uc]代表局部描述子集合,uc 表示為:
② Squeeze 負責全局信息嵌入,采用global average pooling 將一個通道上全局空間信息編碼為一個全局特征,其數學模型為:
式中:H 表示二維特征的高,W 表示二維特征的寬,Fsq 表示全局平均化。
③ Excitation 負責自適應調整來抓取通道之間的依賴關系,選用sigmoid 激活函數作為簡單的門控制單元:
s = Fex(z,W) = σ(g(z,W)) = σ(W2 δ(W1 z)), (4)
式中:δ 表示ReLU 激活函數[19],W1∈RC/r×C,W2 ∈RC×C/r。
最終的輸出X~ 由U 與激活權重轉換得到:
X~c = Fscale(uc ,sc ) = sc ·uc 。(5)
2. 2 基于坐標注意力機制的識別原理
通道注意力機制接收全局特征并編碼通道間的關系,為了利用這種表現特征的能力,根據數字調制信號幅度和相位特點,提出一種二次變換—坐標注意力機制,使用池化核的2 個空間范圍(h,1)和(1,w)沿水平方向和垂直方向對每個通道分別進行空間編碼,隨后將2 個方向的空間編碼鏈接起來,其數學模型表示為:
式中:F1 表示1 ×1 卷積采樣,[·,·]表示空間維度的拼接運算,δ 表示非線性激活函數,f∈CC / r×(H+W)表示對水平和垂直方向的空間信息編碼后的中繼特征,然后將f 沿著空間維度分為2 個獨立的張量f h ∈CC / r×H 、f w ∈CC / r×W ,隨后按式(7)和式(8)進行坐標通道注意力匹配:
gh = σ(Fh(f h )), (7)
gw = σ(Fw(f w ))。(8)
將上述得到的gh,gw 擴展并形成注意力權重,最終輸出yc:
yc(i,j) = xc(i,j)× ghc(i)× gwc(j)。(9)
當經過數字調制的IQ 信號輸入到網絡時,坐標注意力機制使用池化層的一種空間池化(h,1)沿著相位,另一種空間(1,w)沿著幅度對每個輸入信號進行精確編碼,沿著2 個方向聚合信號特征,產生一對方向感知的特征映射,隨后將此特征映射進行鏈接并激活,將信號的相位注意力和幅度注意力同時施加到信號張量上,更準確地提取到信號特征,提高識別精度。基于坐標注意力機制的殘差中坐標注意力機制原理如圖3 所示,其主要由輸入層、坐標注意力機制、殘差層和輸出層組成。在預測模型中,損失函數采用交叉熵損失(Cross Entropy Loss)函數,優化算法采用Adam 算法,并在PyTorch 框架下完成對模型的搭建測試,其模型的具體結構如下。
(1)數據預編碼。本文預編碼主要用于標識信號類型和增強信號。在譯碼過程中,模型根據所觀察到的狀態選擇類型。因此,應對設定參數所產生的數字調制信號進行編碼[20]。
① 標識信號編碼:標識正確是標識一位有效碼并區分不同調制信號的狀態,需要在一個向量中添加一個位來區分“信號調制類型”。本研究采用onehot 編碼,編碼中只有一位為1,其他位為0。對于長度為L+1 的二進制向量,前L 位與長度為X,寬度Y 的信號二維特征向量對應,L+1 位用于標識與之對應的調制信號。
② 增強信號編碼:A 表示AM 編碼器輸出的信號序列,采用批量反正切編碼;A~ 表示一個信號,由神經網絡的輸出來決定。
(2)輸入層。輸入層變量為經過onehot 編碼預處理后的8 種數字調制方式的IQ 數據。
(3)坐標注意力殘差層。基于坐標注意力的ResNet 預測模型中,每個坐標注意力殘差層包括多個殘差塊和多個坐標注意力塊。殘差塊的增加,可以解決梯度消失和梯度爆炸等問題,極大地提升有效訓練的網絡的深度;坐標注意力塊的增加,可以解決通道之間關鍵特征和空間關鍵特征提取等問題。同時,在坐標注意力殘差層中存在一個超參數學習率,用于決定識別最優調制信號類型的序列。本文通過數字信號數據集在預訓練下進行比較分析,篩選出信號識別效果較好的網絡預測模型。
(4)輸出層。輸出層變量為每一種數字調制信號的識別正確概率。
(5)損失計算。坐標注意力殘差模型一般采用交叉熵損失反映預測結果的準確度,其計算公式為:
Loss = - ylgy^ - (1 - y)lg(1 -y^ ), (10)
式中:y 表示真實值,y^ 表示預測值,Loss 表示損失值。
(6)優化算法。本文采用Adam 優化算法,根據學習率的大小與訓練輪次動態調整權重衰減系數。算法的參數設定如下。學習率因子a = 0. 01,a/λ為權重衰減系數,λ 為訓練輪次。
3 實驗部分及分析
為驗證所提出的坐標注意力方案在數字信號識別領域的性能,本文采用GNU Radio 軟件平臺仿真生成的IQ 信號數據集RML2016. 10A 進行相應的實驗,數據集在-20 ~ 18 dB 的信噪比環境下生成。并根據本文實驗選取其中8 種數字信號調制方式,每符號8 個I / Q 采樣點。每個樣本的數據格式為2×128,2 表示I、Q 兩路采樣信號,128 表示采樣點長度,經過數據預處理之后形成實驗所需要的數據集。本文實驗數字信號模擬參數數如表1 所示。部分數據集展示如圖4 所示,在不同的網絡結構和坐標注意力機制下進行了調制識別的系統仿真實驗。
針對不同的數字信號調制方式下,特征提取困難對信號調制識別系統造成不同程度的誤識別影響,需要增加網絡結構的泛化性,以應對各種復雜調制方式,從而建立識別性能更好的神經網絡模型。因此本文將數字信號預處理為(2,128)的IQ 數據,采集1. 6×105條數字信號序列樣本,將60% 的信號序列樣本劃分為訓練集,20% 劃分為驗證集,20% 劃分為測試集。
在該數據集上,本文重點測試了CLDNN、CLDNN_CA、ResNet、CA_ResNet 和ResNet_CA 五種神經網絡結構。實驗中,先進行數據預處理和數字信號特征選擇,學習率初始設置為0. 01(根據訓練誤差動態匹配),訓練的輪次(Epoch)設置為100,批量大小(batch size)設置為64。在不同的調制方式下,對于模型的性能進行討論與分析。
為了分析驗證坐標注意力機制對于網絡收斂速度的影響,得到如圖5 所示的4 種網絡的訓練損失。ResNet 在第5 輪的訓練后,訓練損失得益于殘差塊的作用,迅速下降,模型得到快速收斂。因為長短期記憶(Long Short-Term Memory,LSTM)的存在,其CLDNN 與CLDNN_CA 的訓練損失收斂得慢且整體訓練損失大于ResNet。在相同的網絡結構下,當加入坐標注意力機制時,隨著訓練輪次的增加,其訓練損失相比于傳統的網絡結構下降得更快,當訓練100 輪時,CLDNN_CA 的訓練損失相較于CLDNN 下降了2% ,ResNet_CA 的訓練損失則下降了18% ,且根據驗證集所得到的數據表明4 種網絡結構在100 輪的訓練下均能到達收斂。
為了驗證坐標注意力機制8 種調制信號的識別效果,本文將坐標注意力機制引入到ResNet 中。ResNet 和坐標注意力ResNet 的混淆矩陣分別如圖6 和圖7 所示,受限于殘差塊的作用,坐標注意力機制在GFSK 和QAM64 調制方式上提升了3% 和4%的識別準確率,但同時QAM16 和QAM64 調制方式的識別準確率較低。這是因為2 種調制信號具有較高的調制階數,64QAM 相當于16QAM 的對稱擴展,相似度較高。實驗表明,坐標注意力機制的引入仍然難以提升神經網絡識別不同階數的調制方式的精度,該調制中對幅度和相位分別提取特征,但對于幅度和相位聯合的特征提取較為困難。
為了充分說明坐標注意力機制能夠更好地提取數字信號幅度和相位的多層維度信息,將坐標注意力機制下的神經網絡模型分別對比文獻[5]設計的CLDNN 模型、文獻[20]采用的ResNet 模型,4 種網絡識別效果如圖8 所示。由圖8 可知,與CLDNN 相比,引入了坐標注意力機制的CLDNN_CA 在0 ~ 20 dB下的數字信號識別精度提升約為10% ,但在-10 ~-4 dB 下提升并不明顯。這是由于信號在低信噪比下,噪聲的相對強度更高,使得數字信號丟失了關鍵的信號特征。這一方面使得坐標注意力機制更多地提取到噪聲的特征;另一方面對于數字信號特征提取變得更加困難,造成數字調制信號識別準確率下降。由于沒有引入殘差塊,CLDNN 與CLDNN _CA的識別準確率都大幅低于ResNet,與文獻[21]中ResNet 相對比,在0 ~ 20 dB 下ResNet 和ResNet_CA識別準確率都在90% 以上。在高信噪比情況下,坐標注意力機制并沒有給ResNet 帶來額外的增益,但在-8 dB 下,ResNet _ CA 識別精度比ResNet 高17% ,表明在-8 dB 情況下,坐標注意力機制能夠更好地提高殘差結構的性能。
本文對坐標注意力機制是否影響數字信號的識別精度采取消融實驗,采取3 種網絡結構模型:ResNet 模型、CA_ResNet 模型、ResNet_CA 模型。部分網絡結構如圖9 和圖10 所示。
坐標注意力機制識別效果如圖11 所示。由圖11 可以看出,識別精度高于95% 時,卷積中坐標注意力機制識別性能相較于無坐標注意力機制增益為4 dB,殘差中坐標注意力機制信噪比增益大約為8 dB。在0 ~ 8 dB 情況下,ResNet 對于信號的識別準確率大于ResNet_CA,表明此時殘差塊中的坐標注意力機制對網絡的信號識別性能產生了抑制作用。信噪比為4 dB 時,坐標注意力機制對于ResNet識別正確率降低了1% 。當坐標注意力機制插入到ResNet 中的卷積部分時,坐標注意力機制能夠對ResNet 信號識別產生系統增益。信噪比為- 4 ~18 dB 時,CA_ResNet 對于信號的識別準確率均大于95% 。信噪比為10 dB 時,坐標注意力機制對于ResNet 的信號識別提高約2% 。以上結果表明,坐標注意力機制在卷積塊中相較于殘差塊能夠更好地提取信號的特征,且隨著信噪比的提高,信噪比增益效果更佳。
為分析與比較坐標注意力機制對于不同網絡結構算法的影響,列出不同已訓練完成的網絡結構中部分參數,如表2 所示。因CLDNN 結構簡單,CLDNN_CA 的信號識別時長明顯低于ResNet_CA,且參數大小同樣低于ResNet_CA。CLDNN 的平均絕對誤差(Mean Absolute Error,MAE)為1. 408 7,識別誤差效果最差。得益于坐標注意力機制的加入,CLDNN_CA 的誤差則降低到1. 181 6。同時,殘差中的添加坐標注意力機制同樣降低了平均絕對誤差。相較于未添加坐標注意力機制的神經網絡,坐標注意力機制的引入在降低平均絕對誤差的同時,參數量大小沒有明顯增加。但隨著網絡結構更復雜的同時,坐標注意力機制的引入會加長信號識別的時間。號數據集較小的情況下,ResNet_CA 和CA_ResNet 識別時間增加了大約1 s。
4 結束語
為了提高數字信號識別精度,本文設計了一種新的數字信號調制識別方案———基于坐標注意力機制的數字信號調制識別。通過仿真分析得到,在低信噪比情況下,坐標注意力機制對于數字信號調制識別能夠提升3% 左右的識別精度,識別性能也更加穩定。在此基礎上,本文進行了殘差塊和卷積塊中坐標注意力機制的對比實驗,結果表明卷積中的坐標注意力機制相較于殘差塊中平均增益為4 dB,卷積中的坐標注意力機制的識別準確率能夠達到95% 以上。當識別精度為95% 時,在CNN 中添加坐標注意力機制相較于無坐標注意力機制能夠得到4 dB 左右的信噪比增益,在殘差神經網絡中添加坐標注意力機制相較于無坐標注意力機制能夠得到8 dB 左右的信噪比增益。
實驗表明,相較于通道注意力機制、空間注意力機制,坐標注意力機制可根據數字信號的特點建立時頻域間的關系,同時精確地捕獲空間方向長期的依賴關系,提升數字信號調制識別精度。本文在QAM16 和QAM64 的調制方式識別上仍存在精度不足的問題,后續將繼續考慮提高對這2 類數字信號的識別精度。
參考文獻
[1] MENG F,CHEN P,WU L N,et al. Automatic ModulationClassification:A Deep Learning Enabled Approach [J].IEEE Transactions on Vehicular Technology,2018,67(11):10760-10772.
[2] KULIN M,KAZAZ T,MOERMAN I,et al. EndtoEndLearning from Spectrum Data:A Deep Learning Approachfor Wireless Signal Identification in Spectrum MonitoringApplications[J]. IEEE Access,2018,6:18484-18501.
[3] 李峰. 基于卷積神經網絡的調制信號識別算法研究[D]. 南京:南京信息工程大學,2021.
[4] O’SHEA T J,CORGAN J,CLANCY T C. ConvolutionalRadio Modulation Recognition Networks[C]∥EngineeringApplications of Neural Networks. Aberdeen:Springer,2016:213-226.
[5] O’SHEA T J,ROY T,CLANCY T C. Overtheair DeepLearning Based Radio Signal Classification [J ]. IEEEJournal of Selected Topics in Signal Processing,2018,12(1):168-179.
[6] LIU X Y,YANG D Y,GAMAL A E. Deep Neural NetworkArchitectures for Modulation Classification[C]∥2017 51stAsilomar Conference on Signals,Systems,and Computers.Pacific Grove:IEEE,2017:915-919.
[7] WANG Y,LIU M,YANG J,et al. Datadriven DeepLearning for Automatic Modulation Recognition inCognitive Radios [J]. IEEE Transactions on VehicularTechnology,2019,68(4):4074-4077.
[8] KARRA K,KUZDEBA S,PETERSEN J. Modulation Recognition Using Hierarchical Deep Neural Networks [C]∥2017 IEEE International Symposium on Dynamic SpectrumAccess Networks (Dyspan). Baltimore:IEEE,2017:1-3.
[9] KHAN F N,LU C,LAU A P T. Joint Modulation Format /Bitrate Classification and SignaltoNoise Ratio Estimationin Multipath Fading Channels Using Deep MachineLearning [J ]. Electronics Letters,2016,52 (14 ):1272-1274.
[10] ZHANG Z F,WANG C,GAN C Q,et al. Automatic Modulation Classification Using Convolutional Neural Networkwith Features Fusion of SPWVD and BJD [J ]. IEEETransactions on Signal and Information Processing overNetworks,2019,5(3):469-478.
[11] HOU Q B,ZHOU D Q,FENG J S. Coordinate Attentionfor Efficient Mobile Network Design [C]∥ 2021 IEEE /CVF Conference on Computer Vision and Pattern Recognition. Nashville:IEEE,2021:13708-13717.
[12] MANDT S,HOFFMAN M D,BLEI D M. Stochastic Gradient Descent as Approximate Bayesian Inference [EB /OL]. (2017 - 04 - 13)[2023 - 07 - 01]. https:∥ arxiv.org / pdf / 1704. 04289v1. pdf.
[13] SHELHAMER E,LONG J,DARRELL T. Fully Comvolutional Network for Semantic Segmentation [J ]. IEEETransactions on Pattern Analysis and MachineIntelligence,2017,39(4):640-651.
[14] KRIZHEVSKY A,SUTSKEVER I,HINTON G E.Imagenet Classification with Deep Convolutional NeuralNetworks[J]. Communications of the ACM,2017,60(6):84-90.
[15] ZHAO H S,SHI J P,QI X J,et al. Pyramid Scene ParsingNetwork[C]∥ Proceedings of the IEEE Conference onComputer Vision and Pattern Recognition. Honolulu:IEEE,2017:6230-6239.
[16] HOU Q B,ZHANG L,CHENG M M,et al. Strip Pooling:Rethinking Spatial Pooling for Scene Parsing[C]∥ Proceedings of the IEEE / CVF Conference on ComputerVision and Pattern Recognition. Seattle:IEEE,2020:4003-4012.
[17] HU J,SHEN L,SUN G. SqueezeandExcitation Networks[C]∥2018 IEEE / CVF Conference on Computer Visionand Pattern Recognition. Salt Lake City:IEEE,2018:7132-7141.
[18] NAIR V,HINTON G E. Rectified Linear Units ImproveRestricted Boltzmann Machines[C]∥ Proceedings of the27th International Conference on Machine Learning.Haifa:ACM,2010:807-814.
[19] 陶志勇,閆明豪,劉影,等. 基于AGCNN 的輕量級調制識別方法[J]. 電子測量與儀器學報,2022,36(4):241-249.
[20] 董聰,張傳武,高勇. 基于殘差神經網絡的通信混合信號識別[J]. 無線電工程,2020,50(9):727-731.
作者簡介
張 兢 女,(1965—),碩士,教授。主要研究方向:信號與信息處理、信號檢測與處理。
(*通信作者)蘭思源 男,(1996—),碩士研究生。主要研究方向:深度學習、信號調制識別。
曹 陽 男,(1977—),博士,教授。主要研究方向:現代數字信號處理、通信信道編碼理論。
彭小峰 男,(1980—),碩士,副教授。主要研究方向:無線激光通信技術、嵌入式與智能系統、物聯網與無線傳感器網絡。
基金項目:重慶市教委科學技術項目(KJQN201901125);重慶市基礎與前沿研究計劃項目(cstc2019jcymsxmX0233)