(19)國(guó)家知識(shí)產(chǎn)權(quán)局(10)申請(qǐng)公布號(hào)CN120222972A(71)申請(qǐng)人陜西東方華通微波科技有限公司(74)專(zhuān)利代理機(jī)構(gòu)深圳齊茂專(zhuān)利代理事務(wù)所(普通合伙)44396(54)發(fā)明名稱(chēng)一種射頻微波信號(hào)調(diào)制方法及系統(tǒng)本發(fā)明涉及頻率變換技術(shù)領(lǐng)域，具體為一種射頻微波信號(hào)調(diào)制方法及系統(tǒng)，包括以下步驟：通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與所述時(shí)域峰值間距的差過(guò)峰值間距與主邊帶能量比值分析提取周期特征生成能量耦合因子，主頻點(diǎn)偏移量與能量?jī)A斜向量融合生成諧波傾向參數(shù)，結(jié)合頻率軌跡向量定位干擾源路徑，聚類(lèi)能量與諧波參數(shù)劃分?jǐn)_動(dòng)等級(jí)，灰色預(yù)測(cè)主邊帶頻率差值趨勢(shì)預(yù)判沖突方時(shí)域與頻域特征提取時(shí)域與頻域特征提取主頻點(diǎn)偏移與能量?jī)A斜處理參數(shù)融合與K均值聚類(lèi)2S1:通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與所述時(shí)域峰值間距的差分化運(yùn)算結(jié)果，通過(guò)差異化窗口偏移量下的皮爾遜相關(guān)系數(shù)提取同步性S2:通過(guò)計(jì)算當(dāng)前調(diào)制周期內(nèi)的主頻點(diǎn)偏移量，基于離散傅里葉變換的基頻分量相位差計(jì)算諧波漂移量，結(jié)合能量?jī)A斜向量的歸一化處理結(jié)果，調(diào)用協(xié)方差矩陣對(duì)諧波漂移量與能量?jī)A斜度進(jìn)行加權(quán)融合，輸出諧波影響傾向參數(shù)，同時(shí)基于連續(xù)周期內(nèi)主邊帶中心頻率的浮動(dòng)位置生成主邊帶頻率軌跡向量序列；S3:基于所述能量耦合因子與所述諧波影響傾向參數(shù)，將兩者進(jìn)行逐點(diǎn)乘積運(yùn)算并按樣本編號(hào)升序排列為64維向量，采用K均值聚類(lèi)算法對(duì)兩者的乘性融合結(jié)果進(jìn)行分類(lèi)，劃分諧波擾動(dòng)狀態(tài)等級(jí)，生成矩陣控制項(xiàng)并觸發(fā)對(duì)應(yīng)的通帶壓制行為決策指令。2.根據(jù)權(quán)利要求1所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述能量耦合因子具體為梯度差值、峰值間隔差值、同步因子，所述諧波影響傾向參數(shù)包括加權(quán)偏移量、能量?jī)A斜系數(shù)、協(xié)方差權(quán)重，所述主邊帶頻率軌跡向量序列具體為中心頻率坐標(biāo)、遷移速率、周期索引，所述諧波擾動(dòng)狀態(tài)等級(jí)具體為高干擾交疊狀態(tài)、低干擾孤立狀態(tài)、中干擾過(guò)渡狀3.根據(jù)權(quán)利要求2所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述差異化窗口偏移量為滑動(dòng)信號(hào)相對(duì)于基準(zhǔn)信號(hào)的時(shí)延步進(jìn)量，步進(jìn)量范圍為當(dāng)前調(diào)制周期的1/10至1/5;所述諧波漂移量為基頻分量與基準(zhǔn)頻率的相位差對(duì)應(yīng)的頻率偏移量，通過(guò)離散傅里葉變換提取基頻幅值波動(dòng)量，并以線性差值計(jì)算得到。4.根據(jù)權(quán)利要求3所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述能量耦合因子的獲取步驟具體為：S101:檢測(cè)射頻信號(hào)時(shí)域波形的局部振幅極值點(diǎn)，采用滑動(dòng)窗口極值檢測(cè)法遍歷波形數(shù)據(jù)，標(biāo)記振幅高于前后相鄰采樣點(diǎn)的位置，記錄極值點(diǎn)對(duì)應(yīng)時(shí)間戳，計(jì)算相鄰極值點(diǎn)的時(shí)間間隔，構(gòu)建時(shí)域峰值間距序列并歸一化至相同時(shí)間軸，對(duì)頻域主邊帶范圍進(jìn)行頻點(diǎn)能量累加，并計(jì)算與次邊帶頻點(diǎn)能量總和的比值，生成主邊帶能量比值；S102:基于所述時(shí)域峰值間距序列與所述主邊帶能量比值，將序列數(shù)據(jù)歸一化至相同時(shí)間軸，以主邊帶能量比值為基準(zhǔn)信號(hào)，時(shí)域間距序列為滑動(dòng)信號(hào)，計(jì)算兩信號(hào)在差異化窗口偏移量下的皮爾遜相關(guān)系數(shù)，提取相關(guān)系數(shù)曲線的全局最大值，生成同步性相關(guān)系數(shù)；S103:對(duì)頻域能量密度分布曲線進(jìn)行逐點(diǎn)差分運(yùn)算，計(jì)算相鄰頻點(diǎn)能量差值，生成能量密度梯度值，調(diào)用所述同步性相關(guān)系數(shù)，并將絕對(duì)值與所述能量密度梯度值輸入減法器，計(jì)算兩者差值，疊加時(shí)域峰值間距序列的均方根值，生成能量耦合因子。5.根據(jù)權(quán)利要求4所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述主邊帶頻率軌跡向量序列的獲取步驟具體為：S201:基于當(dāng)前調(diào)制周期的主頻點(diǎn)數(shù)據(jù)，采用離散傅里葉變換計(jì)算基頻分量與基準(zhǔn)頻率的相位差，提取基頻幅值波動(dòng)量，通過(guò)線性差值計(jì)算得到主頻偏移量，對(duì)能量?jī)A斜向量的幅值分量進(jìn)行極差歸一化處理，將多分量壓縮至0-1區(qū)間，生成歸一化能量?jī)A斜向量；S202:將所述主頻偏移量映射至諧波漂移量維度，基于協(xié)方差矩陣的對(duì)角線元素權(quán)重，3將諧波漂移量與所述歸一化能量?jī)A斜向量的幅值乘積作為融合因子，采用標(biāo)準(zhǔn)差動(dòng)態(tài)修正融合因子權(quán)重系數(shù)，構(gòu)建諧波漂移量與能量?jī)A斜度的線性組合，輸出諧波影響傾向參數(shù)；S203:提取連續(xù)10個(gè)調(diào)制周期的主邊帶中心頻率均值，將所述諧波影響傾向參數(shù)作為修正因子，對(duì)相鄰周期均值差值進(jìn)行增益補(bǔ)償，基于三次樣條插值函數(shù)對(duì)補(bǔ)償后的離散點(diǎn)序列進(jìn)行軌跡擬合，并設(shè)定三次樣條插值函數(shù)在相鄰節(jié)點(diǎn)處一階導(dǎo)數(shù)連續(xù)，生成主邊帶頻率軌跡向量序列。6.根據(jù)權(quán)利要求5所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述S3的獲取步驟S301:基于所述能量耦合因子與所述諧波影響傾向參數(shù)，對(duì)兩者進(jìn)行逐點(diǎn)乘積運(yùn)算，將每個(gè)樣本點(diǎn)的能量耦合因子數(shù)值與對(duì)應(yīng)所述諧波影響傾向參數(shù)數(shù)值相乘，將乘積運(yùn)算結(jié)果按樣本編號(hào)升序排列為64維向量，建立參數(shù)融合數(shù)據(jù)組，得到諧波融合系數(shù)；S302:調(diào)用K均值聚類(lèi)算法，設(shè)定初始質(zhì)心數(shù)量為3,以所述諧波融合系數(shù)向量為輸入數(shù)據(jù)，初始化質(zhì)心為隨機(jī)選擇的3個(gè)樣本向量，計(jì)算多樣本向量與質(zhì)心向量的歐氏距離，將樣本歸類(lèi)至距離最近的質(zhì)心簇，更新質(zhì)心為簇內(nèi)樣本均值，迭代至質(zhì)心偏移量小于0.1%,輸出樣本分類(lèi)編號(hào)，生成擾動(dòng)聚類(lèi)標(biāo)簽；S303:調(diào)用K均值聚類(lèi)算法，設(shè)定初始質(zhì)心數(shù)量為3,以所述諧波融合系數(shù)向量為輸入數(shù)據(jù)，初始化質(zhì)心為隨機(jī)選擇的3個(gè)樣本向量，計(jì)算多樣本向量與質(zhì)心向量的歐氏距離，將樣本歸類(lèi)至距離最近的質(zhì)心簇，更新質(zhì)心為簇內(nèi)樣本均值，迭代至質(zhì)心偏移量小于0.1%,輸出樣本分類(lèi)編號(hào)，生成擾動(dòng)聚類(lèi)標(biāo)簽。7.根據(jù)權(quán)利要求6所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述方法還包括：S4:通過(guò)所述主邊帶頻率軌跡向量序列，將相鄰周期頻率增量差值構(gòu)造為頻率遞增幅度序列，對(duì)原始序列進(jìn)行一階累加生成新序列，建立累加序列與時(shí)間變量的灰色微分方程求解趨勢(shì)偏移度，基于預(yù)測(cè)結(jié)果與所述通帶壓制行為決策指令，對(duì)阻帶邊界參數(shù)與通帶安全帶寬執(zhí)行預(yù)加載操作。8.根據(jù)權(quán)利要求7所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述頻率遞增幅度序列具體為正向增量、負(fù)向增量、零增量，所述阻帶邊界參數(shù)具體為上限截止頻率、下限截所述灰色微分方程的發(fā)展系數(shù)與灰色作用量通過(guò)最小二乘法求解，殘差比例閾值設(shè)定為5%。9.根據(jù)權(quán)利要求8所述的一種射頻微波信號(hào)調(diào)制方法，其特征在于，所述S4的獲取步驟S401:基于所述主邊帶頻率軌跡向量序列，提取相鄰周期內(nèi)第n+1周期與第n周期的頻率增量，計(jì)算兩者之間的差值絕對(duì)值，按時(shí)間順序排列差值序列，將所述差值序列中每個(gè)元素除以序列最大值完成歸一化處理，生成頻率遞增量差；S402:調(diào)用灰色預(yù)測(cè)模型，以所述頻率遞增量差序列為輸入數(shù)據(jù)，對(duì)原始序列進(jìn)行一階累加生成新序列，建立累加序列與時(shí)間變量的灰色微分方程，求解方程得到發(fā)展系數(shù)與灰色作用量，計(jì)算預(yù)測(cè)值與原始序列的殘差比例，輸出未來(lái)周期內(nèi)差值變化方向，得到趨勢(shì)偏S403:根據(jù)所述趨勢(shì)偏移度的正負(fù)符號(hào)與幅度值，結(jié)合通帶壓制指令矩陣中的頻段壓4制強(qiáng)度參數(shù)，計(jì)算阻帶邊界左移或右移的步長(zhǎng)系數(shù)，將步長(zhǎng)系數(shù)與通帶安全帶寬基準(zhǔn)值相乘，更新阻帶邊界坐標(biāo)與通帶帶寬數(shù)值，生成預(yù)加載參數(shù)組。10.一種射頻微波信號(hào)調(diào)制系統(tǒng)，其特征在于，根據(jù)權(quán)利要求1-9任一項(xiàng)所述的一種射頻微波信號(hào)調(diào)制方法，所述系統(tǒng)包括：時(shí)頻耦合分析模塊，用于通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與所述時(shí)域峰值間距的差分化運(yùn)算結(jié)果，生成能量耦合因子，將所述能量耦合因子傳遞至擾動(dòng)分級(jí)控制模塊；諧波軌跡建模模塊，用于通過(guò)計(jì)算當(dāng)前調(diào)制周期內(nèi)的主頻點(diǎn)偏移量，結(jié)合能量?jī)A斜向量的歸一化處理結(jié)果，調(diào)用協(xié)方差矩陣對(duì)諧波漂移量與能量?jī)A斜度進(jìn)行加權(quán)融合，輸出諧波影響傾向參數(shù)，同時(shí)基于連續(xù)周期內(nèi)主邊帶中心頻率的浮動(dòng)位置生成主邊帶頻率軌跡向量序列，將所述諧波影響傾向參數(shù)傳遞至擾動(dòng)分級(jí)控制模塊，將所述主邊帶頻率軌跡向量序列傳遞至頻帶預(yù)調(diào)模塊；擾動(dòng)分級(jí)控制模塊，用于通過(guò)所述能量耦合因子與所述諧波影響傾向參數(shù)，采用K均值聚類(lèi)算法對(duì)兩者的乘性融合結(jié)果進(jìn)行分類(lèi)，劃分諧波擾動(dòng)狀態(tài)等級(jí)，生成矩陣控制項(xiàng)并觸發(fā)對(duì)應(yīng)的通帶壓制行為決策指令，將所述矩陣控制項(xiàng)與所述通帶壓制行為決策指令傳遞至頻帶預(yù)調(diào)模塊；頻帶預(yù)調(diào)模塊，用于通過(guò)所述主邊帶頻率軌跡向量序列，將相鄰周期頻率增量差值構(gòu)造為頻率遞增幅度序列，輸入至灰色預(yù)測(cè)模型進(jìn)行趨勢(shì)判定，基于預(yù)測(cè)結(jié)果與所述通帶壓制行為決策指令，對(duì)阻帶邊界參數(shù)與通帶安全帶寬執(zhí)行預(yù)加載操作。5一種射頻微波信號(hào)調(diào)制方法及系統(tǒng)技術(shù)領(lǐng)域[0001]本發(fā)明涉及頻率變換技術(shù)領(lǐng)域，尤其涉及一種射頻微波信號(hào)調(diào)制方法及系統(tǒng)。背景技術(shù)[0002]頻率變換技術(shù)領(lǐng)域包含對(duì)射頻信號(hào)或微波信號(hào)在頻譜位置上的變換處理。該技術(shù)領(lǐng)域的核心內(nèi)容包括信號(hào)的上變頻與下變頻操作，通過(guò)本振信號(hào)與混頻器的協(xié)同作用實(shí)現(xiàn)信號(hào)頻率的轉(zhuǎn)換。在雷達(dá)、通信、電子對(duì)抗及無(wú)線電測(cè)量等應(yīng)用中，頻率變換是確保系統(tǒng)信號(hào)在指定頻段傳輸、接收與處理的基礎(chǔ)手段。整體上，該技術(shù)領(lǐng)域還涉及本振頻率控制、混頻路徑匹配、噪聲抑制及頻譜展寬等多個(gè)環(huán)節(jié)，構(gòu)成了一個(gè)涵蓋信號(hào)生成、頻率映射與多級(jí)調(diào)配的綜合系統(tǒng)。[0003]其中，射頻微波信號(hào)調(diào)制方法是指通過(guò)控制射頻或微波頻段信號(hào)的幅度、相位或頻率特性，實(shí)現(xiàn)對(duì)信號(hào)在時(shí)間或頻率域的有序調(diào)整的技術(shù)方案。主要針對(duì)在特定通信鏈路或射頻處理模塊中需要對(duì)原始基帶信號(hào)進(jìn)行調(diào)制以適應(yīng)特定頻段傳輸?shù)膯?wèn)題，具體采用了基于鎖相環(huán)結(jié)構(gòu)實(shí)現(xiàn)頻率合成的方式，結(jié)合高階模擬混頻器輸出調(diào)制信號(hào)，通過(guò)精密控制本振頻率與中頻信號(hào)的疊加順序完成調(diào)制過(guò)程。此外，該方法還結(jié)合溫度補(bǔ)償電路調(diào)整信號(hào)幅度的線性響應(yīng)，并通過(guò)諧波抑制濾波器對(duì)調(diào)制輸出頻譜進(jìn)行限定處理，從而完成整個(gè)調(diào)制流程的頻率變換要求。[0004]現(xiàn)有技術(shù)依賴(lài)鎖相環(huán)與固定濾波器組合抑制干擾，未建立時(shí)頻域參數(shù)的動(dòng)態(tài)關(guān)聯(lián)機(jī)制。例如，當(dāng)信號(hào)時(shí)域峰值間隔與頻域主邊帶能量分布出現(xiàn)異步變化時(shí)，現(xiàn)有方案無(wú)法識(shí)別時(shí)頻域特征的協(xié)同偏差，導(dǎo)致調(diào)制信號(hào)相位失配。現(xiàn)有諧波抑制濾波器參數(shù)固定，未量化主頻點(diǎn)偏移量與諧波能量分布的耦合關(guān)系，在信號(hào)頻率快速跳變場(chǎng)景下，殘留諧波能量超出濾波器阻帶范圍?，F(xiàn)有技術(shù)未構(gòu)建主邊帶頻率軌跡向量序列，無(wú)法捕捉相鄰周期頻率偏移的連續(xù)性特征，當(dāng)環(huán)境噪聲引發(fā)主邊帶頻率累積漂移時(shí)，阻帶邊界調(diào)整滯后于實(shí)際干擾變化，造成通帶信號(hào)能量泄露?，F(xiàn)有溫度補(bǔ)償電路僅修正幅度線性響應(yīng)，未同步關(guān)聯(lián)頻率偏移趨勢(shì)與能量?jī)A斜度參數(shù)，導(dǎo)致干擾抑制措施缺乏動(dòng)態(tài)適配性。發(fā)明內(nèi)容[0005]本發(fā)明的目的是解決現(xiàn)有技術(shù)中存在的缺點(diǎn)，而提出的一種射頻微波信號(hào)調(diào)制方法及系統(tǒng)。[0006]為了實(shí)現(xiàn)上述目的，本發(fā)明采用了如下技術(shù)方案：一種射頻微波信號(hào)調(diào)制方法，包括以下步驟：S1:通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與所述時(shí)域峰值間距的差分化運(yùn)算結(jié)果，通過(guò)差異化窗口偏移量下的皮爾遜相關(guān)系數(shù)提取同S2:通過(guò)計(jì)算當(dāng)前調(diào)制周期內(nèi)的主頻點(diǎn)偏移量，基于離散傅里葉變換的基頻分量6相位差計(jì)算諧波漂移量，結(jié)合能量?jī)A斜向量的歸一化處理結(jié)果，調(diào)用協(xié)方差矩陣對(duì)諧波漂移量與能量?jī)A斜度進(jìn)行加權(quán)融合，輸出諧波影響傾向參數(shù)，同時(shí)基于連續(xù)周期內(nèi)主邊帶中心頻率的浮動(dòng)位置生成主邊帶頻率軌跡向量序列；S3:基于所述能量耦合因子與所述諧波影響傾向參數(shù)，將兩者進(jìn)行逐點(diǎn)乘積運(yùn)算并按樣本編號(hào)升序排列為64維向量，采用K均值聚類(lèi)算法對(duì)兩者的乘性融合結(jié)果進(jìn)行分類(lèi)，劃分諧波擾動(dòng)狀態(tài)等級(jí)，生成矩陣控制項(xiàng)并觸發(fā)對(duì)應(yīng)的通帶壓制行為決策指令。[0007]作為本發(fā)明的進(jìn)一步方案，所述能量耦合因子具體為梯度差值、峰值間隔差值、同步因子，所述諧波影響傾向參數(shù)包括加權(quán)偏移量、能量?jī)A斜系數(shù)、協(xié)方差權(quán)重，所述主邊帶頻率軌跡向量序列具體為中心頻率坐標(biāo)、遷移速率、周期索引，所述諧波擾動(dòng)狀態(tài)等級(jí)具體為高干擾交疊狀態(tài)、低干擾孤立狀態(tài)、中干擾過(guò)渡狀態(tài)，所述矩陣控制項(xiàng)包括通帶壓縮指[0008]作為本發(fā)明的進(jìn)一步方案，所述能量耦合因子具體為梯度差值、峰值間隔差值、同步因子，所述諧波影響傾向參數(shù)包括加權(quán)偏移量、能量?jī)A斜系數(shù)、協(xié)方差權(quán)重，所述主邊帶頻率軌跡向量序列具體為中心頻率坐標(biāo)、遷移速率、周期索引，所述諧波擾動(dòng)狀態(tài)等級(jí)具體為高干擾交疊狀態(tài)、低干擾孤立狀態(tài)、中干擾過(guò)渡狀態(tài)，所述矩陣控制項(xiàng)包括通帶壓縮指[0009]作為本發(fā)明的進(jìn)一步方案，所述能量耦合因子具體為梯度差值、峰值間隔差值、同步因子，所述諧波影響傾向參數(shù)包括加權(quán)偏移量、能量?jī)A斜系數(shù)、協(xié)方差權(quán)重，所述主邊帶頻率軌跡向量序列具體為中心頻率坐標(biāo)、遷移速率、周期索引，所述諧波擾動(dòng)狀態(tài)等級(jí)具體為高干擾交疊狀態(tài)、低干擾孤立狀態(tài)、中干擾過(guò)渡狀態(tài)，所述矩陣控制項(xiàng)包括通帶壓縮指[0010]作為本發(fā)明的進(jìn)一步方案，所述主邊帶頻率軌跡向量序列的獲取步驟具體為：S201:基于當(dāng)前調(diào)制周期的主頻點(diǎn)數(shù)據(jù)，采用離散傅里葉變換計(jì)算基頻分量與基準(zhǔn)頻率的相位差，提取基頻幅值波動(dòng)量，通過(guò)線性差值計(jì)算得到主頻偏移量，對(duì)能量?jī)A斜向量的幅值分量進(jìn)行極差歸一化處理，將多分量壓縮至0-1區(qū)間，生成歸一化能量?jī)A斜向量；S202:將所述主頻偏移量映射至諧波漂移量維度，基于協(xié)方差矩陣的對(duì)角線元素權(quán)重，將諧波漂移量與所述歸一化能量?jī)A斜向量的幅值乘積作為融合因子，采用標(biāo)準(zhǔn)差動(dòng)態(tài)修正融合因子權(quán)重系數(shù)，構(gòu)建諧波漂移量與能量?jī)A斜度的線性組合，輸出諧波影響傾向S203:提取連續(xù)10個(gè)調(diào)制周期的主邊帶中心頻率均值，將所述諧波影響傾向參數(shù)作為修正因子，對(duì)相鄰周期均值差值進(jìn)行增益補(bǔ)償，基于三次樣條插值函數(shù)對(duì)補(bǔ)償后的離散點(diǎn)序列進(jìn)行軌跡擬合，并設(shè)定三次樣條插值函數(shù)在相鄰節(jié)點(diǎn)處一階導(dǎo)數(shù)連續(xù)，生成主邊帶頻率軌跡向量序列。[0011]作為本發(fā)明的進(jìn)一步方案，所述主邊帶頻率軌跡向量序列的獲取步驟具體為：S201:基于當(dāng)前調(diào)制周期的主頻點(diǎn)數(shù)據(jù)，采用離散傅里葉變換計(jì)算基頻分量與基準(zhǔn)頻率的相位差，提取基頻幅值波動(dòng)量，通過(guò)線性差值計(jì)算得到主頻偏移量，對(duì)能量?jī)A斜向量的幅值分量進(jìn)行極差歸一化處理，將多分量壓縮至0-1區(qū)間，生成歸一化能量?jī)A斜向量；S202:將所述主頻偏移量映射至諧波漂移量維度，基于協(xié)方差矩陣的對(duì)角線元素權(quán)重，將諧波漂移量與所述歸一化能量?jī)A斜向量的幅值乘積作為融合因子，采用標(biāo)準(zhǔn)差動(dòng)7態(tài)修正融合因子權(quán)重系數(shù)，構(gòu)建諧波漂移量與能量?jī)A斜度的線性組合，輸出諧波影響傾向S203:提取連續(xù)10個(gè)調(diào)制周期的主邊帶中心頻率均值，將所述諧波影響傾向參數(shù)作為修正因子，對(duì)相鄰周期均值差值進(jìn)行增益補(bǔ)償，基于三次樣條插值函數(shù)對(duì)補(bǔ)償后的離散點(diǎn)序列進(jìn)行軌跡擬合，并設(shè)定三次樣條插值函數(shù)在相鄰節(jié)點(diǎn)處一階導(dǎo)數(shù)連續(xù)，生成主邊帶頻率軌跡向量序列。S4:通過(guò)所述主邊帶頻率軌跡向量序列，將相鄰周期頻率增量差值構(gòu)造為頻率遞增幅度序列，對(duì)原始序列進(jìn)行一階累加生成新序列，建立累加序列與時(shí)間變量的灰色微分方程求解趨勢(shì)偏移度，基于預(yù)測(cè)結(jié)果與所述通帶壓制行為決策指令，對(duì)阻帶邊界參數(shù)與通帶安全帶寬執(zhí)行預(yù)加載操作。[0013]作為本發(fā)明的進(jìn)一步方案，所述頻率遞增幅度序列具體為正向增量、負(fù)向增量、零增量，所述阻帶邊界參數(shù)具體為上限截止頻率、下限截止頻率，所述通帶安全帶寬包括高頻所述灰色微分方程的發(fā)展系數(shù)與灰色作用量通過(guò)最小二乘法求解，殘差比例閾值設(shè)定為5%。[0014]作為本發(fā)明的進(jìn)一步方案，所述S4的獲取步驟具體為：S401:基于所述主邊帶頻率軌跡向量序列，提取相鄰周期內(nèi)第n+1周期與第n周期的頻率增量，計(jì)算兩者之間的差值絕對(duì)值，按時(shí)間順序排列差值序列，將所述差值序列中每個(gè)元素除以序列最大值完成歸一化處理，生成頻率遞增量差；S402:調(diào)用灰色預(yù)測(cè)模型，以所述頻率遞增量差序列為輸入數(shù)據(jù)，對(duì)原始序列進(jìn)行一階累加生成新序列，建立累加序列與時(shí)間變量的灰色微分方程，求解方程得到發(fā)展系數(shù)與灰色作用量，計(jì)算預(yù)測(cè)值與原始序列的殘差比例，輸出未來(lái)周期內(nèi)差值變化方向，得到趨勢(shì)偏移度；S403:根據(jù)所述趨勢(shì)偏移度的正負(fù)符號(hào)與幅度值，結(jié)合通帶壓制指令矩陣中的頻段壓制強(qiáng)度參數(shù)，計(jì)算阻帶邊界左移或右移的步長(zhǎng)系數(shù)，將步長(zhǎng)系數(shù)與通帶安全帶寬基準(zhǔn)值相乘，更新阻帶邊界坐標(biāo)與通帶帶寬數(shù)值，生成預(yù)加載參數(shù)組。[0015]一種射頻微波信號(hào)調(diào)制系統(tǒng)，所述一種射頻微波信號(hào)調(diào)制系統(tǒng)用于執(zhí)行上述一種射頻微波信號(hào)調(diào)制方法，所述系統(tǒng)包括：時(shí)頻耦合分析模塊，用于通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與所述時(shí)域峰值間距的差分化運(yùn)算結(jié)果，生成能量耦合因子，將所述能量耦合因子傳遞至擾動(dòng)分級(jí)控制模塊；諧波軌跡建模模塊，用于通過(guò)計(jì)算當(dāng)前調(diào)制周期內(nèi)的主頻點(diǎn)偏移量，結(jié)合能量?jī)A斜向量的歸一化處理結(jié)果，調(diào)用協(xié)方差矩陣對(duì)諧波漂移量與能量?jī)A斜度進(jìn)行加權(quán)融合，輸出諧波影響傾向參數(shù)，同時(shí)基于連續(xù)周期內(nèi)主邊帶中心頻率的浮動(dòng)位置生成主邊帶頻率軌跡向量序列，將所述諧波影響傾向參數(shù)傳遞至擾動(dòng)分級(jí)控制模塊，將所述主邊帶頻率軌跡向量序列傳遞至頻帶預(yù)調(diào)模塊；擾動(dòng)分級(jí)控制模塊，用于通過(guò)所述能量耦合因子與所述諧波影響傾向參數(shù)，采用K8均值聚類(lèi)算法對(duì)兩者的乘性融合結(jié)果進(jìn)行分類(lèi)，劃分諧波擾動(dòng)狀態(tài)等級(jí)，生成矩陣控制項(xiàng)并觸發(fā)對(duì)應(yīng)的通帶壓制行為決策指令，將所述矩陣控制項(xiàng)與所述通帶壓制行為決策指令傳遞至頻帶預(yù)調(diào)模塊；頻帶預(yù)調(diào)模塊，用于通過(guò)所述主邊帶頻率軌跡向量序列，將相鄰周期頻率增量差值構(gòu)造為頻率遞增幅度序列，輸入至灰色預(yù)測(cè)模型進(jìn)行趨勢(shì)判定，基于預(yù)測(cè)結(jié)果與所述通帶壓制行為決策指令，對(duì)阻帶邊界參數(shù)與通帶安全帶寬執(zhí)行預(yù)加載操作。本發(fā)明中，通過(guò)時(shí)域峰值間距與頻域主邊帶能量比值的互相關(guān)分析，提取周期性特征并生成能量耦合因子，將時(shí)頻域參數(shù)動(dòng)態(tài)關(guān)聯(lián)，解決信號(hào)同步性偏差問(wèn)題。主頻點(diǎn)偏移量與能量?jī)A斜向量的協(xié)方差矩陣融合生成諧波影響傾向參數(shù)，量化諧波漂移對(duì)主邊帶能量分布的干擾程度，結(jié)合主邊帶頻率軌跡向量序列精準(zhǔn)定位干擾源演化路徑?；诰垲?lèi)算法對(duì)能量耦合因子與諧波傾向參數(shù)分類(lèi)劃分?jǐn)_動(dòng)等級(jí)，實(shí)現(xiàn)干擾抑制策略的自適應(yīng)匹配。利用灰色預(yù)測(cè)模型分析主邊帶頻率增量差值序列的趨勢(shì)特征，預(yù)判頻譜沖突方向并預(yù)加載阻帶參數(shù)，優(yōu)化通帶壓制響應(yīng)速度。該方法構(gòu)建時(shí)頻聯(lián)合分析、諧波動(dòng)態(tài)建模、干擾狀態(tài)分級(jí)及頻率趨勢(shì)預(yù)判的閉環(huán)調(diào)控體系，提升調(diào)制信號(hào)在復(fù)雜干擾環(huán)境下的抗諧波漂移能力，降低頻譜資源沖突風(fēng)險(xiǎn)。附圖說(shuō)明[0017]圖1為本發(fā)明的工作流程示意圖；圖2為本發(fā)明能量耦合因子的獲取步驟流程圖；圖3為本發(fā)明主邊帶頻率軌跡向量序列的獲取步驟流程圖；圖4為本發(fā)明S3的獲取步驟流程圖；圖5為本發(fā)明S4的獲取步驟流程圖。具體實(shí)施方式[0018]為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合附圖及實(shí)施例，對(duì)本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說(shuō)明。應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅僅用以解釋本發(fā)明，并不用于限定本發(fā)明。示的方位或位置關(guān)系，僅是為了便于描述本發(fā)明和簡(jiǎn)化描述，而不是指示或暗示所指的裝置或元件必須具有特定的方位、以特定的方位構(gòu)造和操作，因此不能理解為對(duì)本發(fā)明的限[0020]實(shí)施例一：請(qǐng)參閱圖1,本發(fā)明提供一種技術(shù)方案：一種射頻微波信號(hào)調(diào)制方法，包括以下步驟：S1:通過(guò)檢測(cè)射頻信號(hào)的時(shí)域峰值間距與頻域主邊帶能量比值，將兩者輸入至互相關(guān)算法進(jìn)行同步性分析，提取時(shí)域峰值序列的周期性特征，基于頻域能量密度梯度值與時(shí)域峰值間距的差分化運(yùn)算結(jié)果，通過(guò)差異化窗口偏移量下的皮爾遜相關(guān)系數(shù)提取同步性9S2:通過(guò)計(jì)算當(dāng)前調(diào)制周期內(nèi)的主頻點(diǎn)偏移量，基于離散傅里葉變換的基頻分量相位差計(jì)算諧波漂移量，結(jié)合能量?jī)A斜向量的歸一化處理結(jié)果，調(diào)用協(xié)方差矩陣對(duì)諧波漂移量與能量?jī)A斜度進(jìn)行加權(quán)融合，輸出諧波影響傾向參數(shù)，同時(shí)基于連續(xù)周期內(nèi)主邊帶中心頻率的浮動(dòng)位置生成主邊帶頻率軌跡向量序列；S3:基于能量耦合因子與諧波影響傾向參數(shù)，將兩者進(jìn)行逐點(diǎn)乘積運(yùn)算并按樣本編號(hào)升序排列為64維向量，采用K均值聚類(lèi)算法對(duì)兩者的乘性融合結(jié)果進(jìn)行分類(lèi)，劃分諧波擾動(dòng)狀態(tài)等級(jí)，生成矩陣控制項(xiàng)并觸發(fā)對(duì)應(yīng)的通帶壓制行為決策指令；S4:通過(guò)主邊帶頻率軌跡向量序列，將相鄰周期頻率增量差值構(gòu)造為頻率遞增幅度序列，對(duì)原始序列進(jìn)行一階累加生成新序列，建立累加序列與時(shí)間變量的灰色微分方程求解趨勢(shì)偏移度，基于預(yù)測(cè)結(jié)果與通帶壓制行為決策指令，對(duì)阻帶邊界參數(shù)與通帶安全帶寬執(zhí)行預(yù)加載操作。[0021]能量耦合因子具體為梯度差值、峰值間隔差值、同步因子，諧波影響傾向參數(shù)包括加權(quán)偏移量、能量?jī)A斜系數(shù)、協(xié)方差權(quán)重，主邊帶頻率軌跡向量序列具體為中心頻率坐標(biāo)、遷移速率、周期索引，諧波擾動(dòng)狀態(tài)等級(jí)具體為高干擾交疊狀態(tài)、低干擾孤立狀態(tài)、中干擾過(guò)渡狀態(tài)，矩陣控制項(xiàng)包括通帶壓縮指令、阻帶遷移指令、峰值壓制指令，頻率遞增幅度序列具體為正向增量、負(fù)向增量、零增量，阻帶邊界參數(shù)具體為上限截止頻率、下通帶安全帶寬包括高頻擴(kuò)展余量、低頻緩沖余量。[0022]差異化窗口偏移量為滑動(dòng)信號(hào)相對(duì)于基準(zhǔn)信號(hào)的時(shí)延步進(jìn)量，步進(jìn)量范圍為當(dāng)前調(diào)制周期的1/10至1/5;諧波漂移量為基頻分量與基準(zhǔn)頻率的相位差對(duì)應(yīng)的頻率偏移量，通過(guò)離散傅里葉變換提取基頻幅值波動(dòng)量，并以線性差值計(jì)算得到；灰色微分方程的發(fā)展系數(shù)與灰色作用量通過(guò)最小二乘法求解，殘差比例閾值設(shè)定為5%。[0023]請(qǐng)參閱圖2,能量耦合因子的獲取步驟具體為：S101:檢測(cè)射頻信號(hào)時(shí)域波形的局部振幅極值點(diǎn)，采用滑動(dòng)窗口極值檢測(cè)法遍歷波形數(shù)據(jù)，標(biāo)記振幅高于前后相鄰采樣點(diǎn)的位置，記錄極值點(diǎn)對(duì)應(yīng)時(shí)間戳，計(jì)算相鄰極值點(diǎn)的時(shí)間間隔，構(gòu)建時(shí)域峰值間距序列并歸一化至相同時(shí)間軸，對(duì)頻域主邊帶范圍進(jìn)行頻點(diǎn)能量累加，并計(jì)算與次邊帶頻點(diǎn)能量總和的比值，生成主邊帶能量比值；檢測(cè)射頻信號(hào)時(shí)域波形的局部振幅極值點(diǎn)，此過(guò)程通過(guò)采用滑動(dòng)窗口極值檢測(cè)方法遍歷已采集的波形數(shù)據(jù)實(shí)現(xiàn)，具體而言，設(shè)定一個(gè)包含三個(gè)連續(xù)采樣點(diǎn)的滑動(dòng)窗口，當(dāng)窗口中心點(diǎn)的振幅值均大于其緊鄰的前一個(gè)采樣點(diǎn)和后一個(gè)采樣點(diǎn)的振幅時(shí)，該中心點(diǎn)被標(biāo)記為一個(gè)局部振幅極大值點(diǎn)，并記錄此極值點(diǎn)對(duì)應(yīng)的時(shí)間戳，例如，在對(duì)一段持續(xù)1000納秒(ns)的射頻信號(hào)以10Gs/s的速率采樣后，獲得了10000個(gè)振幅數(shù)據(jù)點(diǎn)，遍歷這些數(shù)據(jù)點(diǎn)，識(shí)別出所有滿足上述條件的局部振幅極大值點(diǎn)，設(shè)在此1000ns的片段內(nèi)，識(shí)別出5個(gè)極值點(diǎn)，△t?=T?-T?=950.3-780.5=169.8ns構(gòu)建時(shí)域峰值間距序列SAt=[205.6,184.3,240.4,169.8](單位：ns),為確保不同觀測(cè)片段或不同信號(hào)源獲取2.44GHz至2.46GHz,此范圍的設(shè)定依據(jù)是該無(wú)線通信系統(tǒng)的信道帶寬標(biāo)準(zhǔn)，次邊帶范圍定義為主邊帶兩側(cè)各5MHz的頻段，即2.435GHz至2.44GHz以及2.46GHz至2.465GHz,對(duì)在主邊帶范圍內(nèi)的所有離散頻點(diǎn)上的信號(hào)能量(幅度的Rmain/sub=Emain/Esub=(2.8×10-5)/(0.4×1采樣點(diǎn)序號(hào)時(shí)間戳(ns)振幅(mV)是否極值點(diǎn)相鄰極值點(diǎn)時(shí)間間隔(ns)..·..否否否是(T2)否.·.否是(T3)否否是(T4)否否是(T5)否時(shí)間間隔。[0025]S102:基于時(shí)域峰值間距序列與主邊帶能量比值，將序列數(shù)據(jù)歸一化至相同時(shí)間基于S101獲取的歸一化時(shí)域峰值間距序列S'△t=[0.2056,0.1843,0.2404,0.1698]與在連續(xù)觀測(cè)中獲得的一系列主邊帶能量比值，此處為主邊帶能量比值同樣在對(duì)應(yīng)S'△t各元素的時(shí)間段內(nèi)計(jì)算獲得，形成序列Se=[7.0,7.2,6.8,7.3],將這兩個(gè)長(zhǎng)度相同的序列數(shù)據(jù)歸一化至相同時(shí)間軸的操作已在S101中通過(guò)對(duì)時(shí)域峰值間距序列進(jìn)行歸一化完成，主邊帶能量比值序列本身即是在文=(0.2056+0.1843+0.2404+0.1698)/4=0.8001/4=0.200025;y=(7.0+7.2+6.8+7.3)/4=28.3/4=7計(jì)算(x;-)序列：[0.005575,-0.015725,0.040375,-0.030225];計(jì)算(y;-y)序列：[-0.075,0.125,-0.275,0.225];=(0.005575)(-0.075)+(-0.015725)(0.12+(0.040375)(-0.275)+(-0.030225)(0.225);=-0.000418125-0.001965625-0.∑(x?-)2=(0.005575)2+(-0.015725)2+(0.040375)2+(-0=0.000031080625+0.000247275625+0.0016∑(y?-y)2=(-0.075)2+(0.125)2+(-0.275)2=0.005625+0.015625+0.075625+0.05062對(duì)偏移量為-1(S'△右移一位，比較元素減少)和+1S′左移一位)也執(zhí)行類(lèi)列為[-0.852,-0.994375,-0.798],從中提取相關(guān)系數(shù)曲線的全局最大值，這里的“全關(guān)系數(shù)度梯度序列G=[G?,G?,…,GN-1]中的每一項(xiàng)G;=Ei+1-E,例如，序列(單位：10-7W/Hz,為簡(jiǎn)化表示，后續(xù)運(yùn)算中省略單位，但最終結(jié)果會(huì)附加)為D=(0.298+0.998-0.002+0.698+1.298)/5=3.29/5=0.658,設(shè)定一個(gè)201.781=208.361,此處的RMS△t單位是ns,WD·D是一個(gè)調(diào)整后的數(shù)值，最終S201:基于當(dāng)前調(diào)制周期的主頻點(diǎn)數(shù)據(jù)，采用離散傅里葉變換計(jì)算基頻分量與基準(zhǔn)頻率的相位差，提取基頻幅值波動(dòng)量，通過(guò)線性差值計(jì)算得到主頻偏移量，對(duì)能量?jī)A斜向量的幅值分量進(jìn)行極差歸一化處理，將多分量壓縮至0-1區(qū)間，生成歸一化能量?jī)A斜向量；基于當(dāng)前調(diào)制周期的主頻點(diǎn)數(shù)據(jù)，在某雷達(dá)系統(tǒng)中，一個(gè)調(diào)制周期內(nèi)通過(guò)高精度頻譜分析儀測(cè)得主頻點(diǎn)為f?=9410.5MHz,系統(tǒng)的標(biāo)稱(chēng)基準(zhǔn)頻率為fref=9410.0MHz,采用離散傅里葉變換(DFT)對(duì)主頻點(diǎn)附近的窄帶信號(hào)進(jìn)行精細(xì)分析，計(jì)算得到f?處的復(fù)頻譜值為X(f?),其相位為φ(f?),同時(shí)參考基準(zhǔn)頻率fref處的相位(或一個(gè)預(yù)設(shè)的參考相位)波動(dòng)量，例如，連續(xù)監(jiān)測(cè)10個(gè)脈沖的基頻幅值，其標(biāo)準(zhǔn)差為0=0.005V(相對(duì)于平均幅值1V),通過(guò)線性差值計(jì)算得到主頻偏移量△fmain,此處基于相位差進(jìn)行計(jì)算，相位-頻率轉(zhuǎn)換系數(shù)kφ→f由系統(tǒng)設(shè)計(jì)參數(shù)(如VCO調(diào)諧靈敏度)確定，設(shè)kφ→f=30MHz/,則主頻偏移450kHz,接下來(lái)處理從頻譜分析中獲取的能量?jī)A斜向量，該向量表征主邊帶內(nèi)能量分布的不對(duì)稱(chēng)性，其幅值分量為Vtilt=[2.1,3.5,2.8](單位：dB),對(duì)此向量的幅值分量進(jìn)行極差歸一化處理，首先找出最小值Vmin=2.1dB和最大值Vmax=3.5dB,對(duì)每個(gè)分量v;應(yīng)用v′?=(2.1-2.1)/(3.5-2.1)=0v′?=(3.5-2.1)/(3.5-2.1)=1.4v′?=(2.8-2.1)/(3.5-2.1)=0.7/1生成歸一化能量?jī)A斜向量S202:將主頻偏移量映射至諧波漂移量維度，基于協(xié)方差矩陣的對(duì)角線元素權(quán)重，將諧波漂移量與歸一化能量?jī)A斜向量的幅值乘積作為融合因子，采用標(biāo)準(zhǔn)差動(dòng)態(tài)修正融合因子權(quán)重系數(shù)，構(gòu)建諧波漂移量與能量?jī)A斜度的線性組合，輸出諧波影響傾向參數(shù)；將S201中獲得的主頻偏移量△fmain=450kHz映射至諧波漂移量維度，諧波漂移量Hdrift定義為主頻偏移量相對(duì)于一個(gè)參考帶寬的比例，參考帶寬Bref設(shè)定為系統(tǒng)工作帶寬的1%,若系統(tǒng)工作帶寬為50MHz,則Bref=0.5MHz=500kHz,因此，諧波漂移量Harift=△fmain/Bref=450kHz/500kHz=0.9,此為一個(gè)無(wú)單位的相對(duì)值，基于一個(gè)預(yù)先建立的協(xié)方差矩陣的對(duì)角線元素，這些元素反映了歸一化能量?jī)A斜向量各分量歷史波動(dòng)的重要性，設(shè)對(duì)應(yīng)V′的三個(gè)分量的權(quán)重Wdiag=[0.2,0.5,0.3],這些權(quán)重通過(guò)對(duì)初始融合因子,采用為系統(tǒng)較穩(wěn)定，賦予較高權(quán)重，反之則較低，權(quán)重修正系數(shù)W計(jì)算規(guī)則為：,其中kcorr是一個(gè)調(diào)節(jié)因子，設(shè)為2.0,則Wc=1+(0.1-0.05)×2.0=1+0.05×2.0=1+0.1=1.1,修正后的融合因子Ffused=W.×Finitial=1.1×0.585=0.6435,構(gòu)建諧波漂移量與能量?jī)A斜度的優(yōu)化得到的系數(shù)，確保兩者對(duì)最終參數(shù)的貢獻(xiàn)均衡，設(shè)定α=0.6和β=0.4,(α+β=1),Ptendency=0.6×0.9+0.4×0.6435=0.54+0.2574=0.7974周期均值之間的差值進(jìn)行增益補(bǔ)償，計(jì)算相鄰周期的頻率均值差值△f=fi+1-Fj,△f?=9410.52-9410.50=0.02MHMHz…(計(jì)算所有9個(gè)差值)△f=9410.58-9410.59=-0.01MHz補(bǔ)償后的差值△F′?=-0.01×(1+0.1×0.7974)≈-0.0107974MHz,依此計(jì)算得到補(bǔ)償后的差值序列△f′,并由此重構(gòu)補(bǔ)償后的頻率均值序列f',其中f?=9410.5215948-0.0107974=9410.5107974(依此類(lèi)推計(jì)算至f′10)之間構(gòu)造三次多項(xiàng)式S;(t),并設(shè)定三次樣條插值函數(shù)在各內(nèi)部節(jié)點(diǎn)t;處具有連續(xù)的一階和基于S103獲取的能量耦合因子Kcoupling=208.361和S202獲取的諧波影響傾向參數(shù)Ptendency=0.7974,對(duì)這兩個(gè)參數(shù)進(jìn)行逐點(diǎn)乘積運(yùn)算，在實(shí)際獲取的數(shù)據(jù)，從而得到64個(gè)能量耦合因子值Kcoupling;和64個(gè)諧波影響傾向參數(shù)值Ptendencyj,其中j=1,…,64.表2能量耦合因子與諧波影響傾向參數(shù)示例(前5個(gè)樣本):樣本編號(hào)能量耦合因子諧波影響傾向參數(shù)乘積12345M?=208.361×0.7974≈166.139,對(duì)所有64個(gè)乘積值M?,M?,…,M64,將這些乘積值按樣本編號(hào)升序排列，形成一個(gè)64維的向量(即一個(gè)包含64個(gè)元素的一維數(shù)組)M=[166.139,170.642,161.694,171.131,164.138,…,M?4],的參數(shù)融合數(shù)據(jù)組，也稱(chēng)為諧波融合系數(shù)。[0032]S302:調(diào)用K均值聚類(lèi)算法，設(shè)定初始質(zhì)心數(shù)量為3,以諧波融合系數(shù)向量為輸入數(shù)據(jù)，初始化質(zhì)心為隨機(jī)選擇的3個(gè)樣本向量，計(jì)算多樣本向量與質(zhì)心向量的歐氏距離，將樣本歸類(lèi)至距離最近的質(zhì)心簇，更新質(zhì)心為簇內(nèi)樣本均值，迭代至質(zhì)心偏移量小于0.1%,輸出樣本分類(lèi)編號(hào)，生成擾動(dòng)聚類(lèi)標(biāo)簽；調(diào)用K均值聚類(lèi)算法，設(shè)定初始質(zhì)心數(shù)量為3,以S301步驟中生成的諧波融合系數(shù)序列M=[M?,M?,…,M?4]作為輸入數(shù)據(jù)(每個(gè)M視為一個(gè)一維數(shù)據(jù)點(diǎn)),初始化質(zhì)心的方法是隨機(jī)從M中選擇3個(gè)不同的值作為初始質(zhì)心C?,C?,C?,設(shè)隨機(jī)選擇的初始質(zhì)心為C?=M?0=160.5,C?=M?5=175.8,c?=M?0=150.2,對(duì)于數(shù)據(jù)集中的每一個(gè)數(shù)據(jù)點(diǎn)M;,計(jì)算其與3個(gè)質(zhì)心的歐氏距離(在一維情況下即為差的絕對(duì)值),并將M;歸類(lèi)至與其距離最近的質(zhì)心所在的簇，例如，對(duì)于M?=166.139,其與三個(gè)質(zhì)心的距離分別為d(M?,c?)=|166.139-150.2|=15.939,由于M?與c?的距離最小，將M?歸入簇1,完成所有64個(gè)數(shù)據(jù)點(diǎn)的歸類(lèi)后，更新每個(gè)簇的質(zhì)心，新質(zhì)心為該簇內(nèi)所有數(shù)據(jù)點(diǎn)的算術(shù)平均值，此分配和更新過(guò)程不斷迭代，迭代停止的條件是所有質(zhì)心的位置在連續(xù)兩次迭代之間的偏移量均小于質(zhì)心值的0.1%,即對(duì)于每個(gè)質(zhì)心Ck,,當(dāng)滿足此條件時(shí)，算法收斂，輸出每個(gè)樣本所屬的分類(lèi)編號(hào)(1,2,或3),這些分類(lèi)編號(hào)序列構(gòu)成了生成的擾動(dòng)聚類(lèi)標(biāo)簽。[0033]S303:調(diào)用K均值聚類(lèi)算法，設(shè)定初始質(zhì)心數(shù)量為3,以諧波融合系數(shù)向量為輸入數(shù)據(jù)，初始化質(zhì)心為隨機(jī)選擇的3個(gè)樣本向量，計(jì)算多樣本向量與質(zhì)心向量的歐氏距離，將樣本歸類(lèi)至距離最近的質(zhì)心簇，更新質(zhì)心為簇內(nèi)樣本均值，迭代至質(zhì)心偏移量小于0.1%,輸出樣本分類(lèi)編號(hào)，生成擾動(dòng)聚類(lèi)標(biāo)簽。[0034]此步驟的描述與S302完全一致，其目的是執(zhí)行K均值聚類(lèi)以生成擾動(dòng)聚類(lèi)標(biāo)簽，因此，其具體執(zhí)行過(guò)程、數(shù)據(jù)處理方式及算例參照S302中的詳細(xì)說(shuō)明，同樣以S301生成的諧波融合系數(shù)序列M為輸入，設(shè)定聚類(lèi)數(shù)量為3,通過(guò)隨機(jī)選擇初始質(zhì)心，并迭代計(jì)算距離、分配樣本、更新質(zhì)心，直至質(zhì)心偏移量小于0.1%的收斂標(biāo)準(zhǔn)達(dá)成，最終輸出各樣本的聚類(lèi)標(biāo)簽。[0035]請(qǐng)參閱圖5,S4的獲取步驟具體為：S401:基于主邊帶頻率軌跡向量序列，提取相鄰周期內(nèi)第n+1周期與第n周期的頻率增量，計(jì)算兩者之間的差值絕對(duì)值，按時(shí)間順序排列差值序列，將差值序列中每個(gè)元素除以序列最大值完成歸一化處理，生成頻率遞增量差；此步驟的描述與S302完全一致，其目的是執(zhí)行K均值聚類(lèi)以生成擾動(dòng)聚類(lèi)標(biāo)簽，因調(diào)用GM(1,1)灰色預(yù)測(cè)模型，以S401生成的頻率遞增量差序列x(◎)(1)=0.81,x()(2)=0.875,x)(3)=0.375,x?(4)=1.0,對(duì)x(◎)進(jìn)行一階x(1)(3)=1.685+0.375x(◎)(k)+az1)(k)=b,其中z(①(k)=0.5(x①)(k)+x(k-1));(0.81-0.544/(-0.109))e-(-0.109)(4)+0.544(O)(5)=(1)(5)-(1)(4)=3.98此預(yù)測(cè)值0)(5)=0.9211即為得到的趨勢(shì)偏移度。[0037]S403:根據(jù)趨勢(shì)偏移度的正負(fù)符號(hào)與幅度值，結(jié)合通帶壓制指令矩陣中的頻段壓制強(qiáng)度參數(shù)，計(jì)算阻帶邊界左移或右移的步長(zhǎng)系數(shù)，將步長(zhǎng)系數(shù)與通帶安全帶寬基準(zhǔn)值相乘，更新阻帶邊界坐標(biāo)與通帶帶寬數(shù)值，生成預(yù)加載參數(shù)組。[0038]根據(jù)S402得到的趨勢(shì)偏移度Toffset=0.9211,其為正值，幅度為0.9211,結(jié)合一個(gè)通帶壓制指令矩陣(如下表3所示)中的頻段壓制強(qiáng)度參數(shù)，表3通帶壓制指令參照表：趨勢(shì)偏移度ToffsetToffset區(qū)間頻段壓制強(qiáng)度參數(shù)SpSpToffset<0Toffset<0(負(fù)向趨勢(shì))如表3所示，該表根據(jù)趨勢(shì)偏移度的不同區(qū)間設(shè)定了相應(yīng)的頻段壓制強(qiáng)度參數(shù)，由于Toffset=0.