8的概率分布；
[0022] 并計算sfb中的比例因子0至8的概率分布向量，記為％，其中，％的維數為 IX9, ^中的第1個元素的值為Sfb中的比例因子0的概率分布，仏中的第2個元素的值 為sfb中的比例因子1的概率分布，依次類推，中的第9個元素的值為sfb中的比例因子 8的概率分布；
[0023] ⑥_6、計算々與々的差值向量，記為AP，處AP代表比例因子〇至 8的使用概率增量向量；然后將AP中的9個元素的值按序作為訓練樣本集合中的第k個 子樣本的第22個特征值至第30個特征值；
[0024] ⑥-7、計算々與$6的商值向量，記為dP，獲，dP代表比例因子0至8的 使用概率增長倍數向量；然后將dP中的9個元素的值按序作為訓練樣本集合中的第k個子 樣本的第31個特征值至第39個特征值；
[0025] ⑥_8、計算比例因子0的使用概率增長倍數與比例因子1至8各自的使用概率增 長倍數的差值，將得到的8個差值按序作為訓練樣本集合中的第k個子樣本的第40個特征 值至第47個特征值；
[0026] ⑥-9、將訓練樣本集合中的第k個子樣本的47個特征值構成行向量Fk。
[0027] 所述的步驟⑧中一個待檢測MP3音頻的47個特征值的提取過程為：
[0028] ⑧-1、將待檢測MP3音頻中利用長窗編碼方式進行處理的幀定義為長窗編碼幀； 然后利用步驟③-1中所使用的MP3解碼器對待檢測MP3音頻中的每幀進行解碼處理，得到 待檢測MP3音頻解碼后的WAV音頻，在解碼處理過程中，提取出待檢測MP3音頻中的每幀長 窗編碼幀的位置；接著根據待檢測MP3音頻中的每幀長窗編碼幀的位置，獲取待檢測MP3音 頻中的所有長窗編碼幀的位置的比例因子矩陣，記為sfa'，其中，如果待檢測MP3音頻為單 聲道音頻，則sfa'的維數為2w'X21，如果待檢測MP3音頻為雙聲道音頻，則sfa'的維數為 4w'X21，w'表示待檢測MP3音頻中的長窗編碼幀的總幀數；
[0029] ⑧-2、利用步驟②中所使用的MP3編碼器對步驟⑧-1中得到的解碼后的WAV音 頻進行編碼處理，得到新的MP3音頻，在編碼處理過程中，根據步驟⑧-1中提取出的待檢測 MP3音頻中的每幀長窗編碼幀的位置，獲取步驟⑧-1中得到的解碼后的WAV音頻在編碼處 理過程中的所有長窗編碼幀的位置的比例因子矩陣，記為sfb'，其中，本步驟中的MP3編碼 器的編碼比特率與待檢測MP3音頻的獲取過程中最后一次壓縮編碼所使用的MP3編碼器的 編碼比特率相同，如果待檢測MP3音頻為單聲道音頻，則sfb'的維數為2w'X21，如果待檢 測MP3音頻為雙聲道音頻，則sfb'的維數為4w'X21 ;
[0030] ⑧-3、計算sfa'中的比例因子1至5轉移為sfb'中的比例因子0至3的轉移概 率矩陣，記為P'，其中，P'的維數為5X4,P'中的第1行第1列元素的值為Sfa'中的比例 因子1轉移為sfb'中的比例因子0的轉移概率，P'中的第1行第2列元素的值為Sfa'中 的比例因子1轉移為sfb'中的比例因子1的轉移概率，P'中的第1行第3列元素的值為 sfa'中的比例因子1轉移為sfb'中的比例因子2的轉移概率，P'中的第1行第4列元素的 值為sfa'中的比例因子1轉移為sfb'中的比例因子3的轉移概率，依次類推，P'中的第5 行第4列元素的值為sfa'中的比例因子5轉移為sfb'中的比例因子3的轉移概率；然后將 P'中的20個元素的值按序作為待檢測MP3音頻的前20個特征值；
[0031] ⑧_4、計算sfa'與sfb'的差值矩陣，記為Asf'，Asf' =sfa'_sfb' ；然后計算 Asf'中的所有元素的值的均值，將該均值作為待檢測MP3音頻的第21個特征值；
[0032] ⑧-5、計算sfa'中的比例因子0至8的概率分布向量，記為％，其中，?的維數 為IX9, ％中的第1個元素的值為sfa'中的比例因子0的概率分布，％中的第2個元素 的值為sfa'中的比例因子1的概率分布，依次類推，中的第9個元素的值為sfa'中的比 例因子8的概率分布；
[0033] 并計算Sfb'中的比例因子0至8的概率分布向量，記為％，其中，&的維數為 1X9, 中的第1個元素的值為sfb'中的比例因子0的概率分布，中的第2個元素的 值為sfb'中的比例因子1的概率分布，依次類推，％中的第9個元素的值為sfb'中的比例 因子8的概率分布；
[0034] ⑧-6、計算與.?，的差值向量，記為AP'，處'=匕.-AP'代表比例因子 〇至8的使用概率增量向量；然后將AP'中的9個元素的值按序作為待檢測MP3音頻的第 22個特征值至第30個特征值；
[0035] ⑧_7、計算％,與％的商值向量，記為dP，，代表比例因子0至 8的使用概率增長倍數向量；然后將dP'中的9個元素的值按序作為待檢測MP3音頻的第 31個特征值至第39個特征值；
[0036] ⑧-8、計算比例因子0的使用概率增長倍數與比例因子1至8各自的使用概率增 長倍數的差值，將得到的8個差值按序作為待檢測MP3音頻的第40個特征值至第47個特 征值。
[0037] 與現有技術相比，本發(fā)明的優(yōu)點在于：
[0038] 1)通過獲取未壓縮的WAV音頻樣本對應的一次壓縮MP3音頻樣本、二次壓縮MP3 音頻樣本和三次壓縮MP3音頻樣本，且限制獲取二次壓縮MP3音頻樣本的第一次壓縮和獲 取三次壓縮MP3音頻樣本的第一次壓縮所使用的MP3編碼器的編碼比特率相同，獲取一次 壓縮MP3音頻樣本、獲取二次壓縮MP3音頻樣本的第二次壓縮、獲取三次壓縮MP3音頻樣本 的第二次壓縮和第三次壓縮所使用的MP3編碼器的編碼比特率相同，利用這些MP3音頻樣 本各自的特征值訓練得到的訓練模板能夠很好的判斷待檢測的MP3音頻樣本為一次壓縮 的MP3音頻樣本或為二次壓縮的MP3音頻樣本還是為三次壓縮的MP3音頻樣本，實現了MP3 音頻壓縮歷史檢測。
[0039] 2)本發(fā)明方法通過分析MP3編碼基本原理，根據長窗比例因子在MP3多次壓縮過 程中的變化規(guī)律，提取基于長窗比例因子的特征值作為MP3音頻壓縮歷史檢測的依據。通 過利用長窗比例因子的特征可有效降低檢測方法的復雜度。
[0040] 3)本發(fā)明通過提取兩次壓縮前后的長窗比例因子的差值均值，以及利用概率統(tǒng)計 方法與單步馬爾科夫轉移概率分別提取長窗比例因子的概率分布和轉移概率作為該檢測 方法的有效特征，這些特征的提取和利用可有效提高檢測方法的準確率。
【附圖說明】
[0041] 圖1為本發(fā)明方法的總體實現框圖。
【具體實施方式】
[0042] 以下結合附圖實施例對本發(fā)明作進一步詳細描述。
[0043] 本發(fā)明提出的一種MP3音頻壓縮歷史檢測方法，其總體實現框圖如圖1所示，其包 括以下步驟：
[0044] ①選取N個風格不同且未壓縮的WAV音頻樣本，其中，N彡10,如取N= 300。在 此，未壓縮的WAV音頻樣本可以選擇blues、pop、classical、country和folk等風格的音 頻。
[0045] ②獲取與每個未壓縮的WAV音頻樣本對應的一次壓縮MP3音頻樣本，具體過程為： 利用現有的MP3編碼器對每個未壓縮的WAV音頻樣本進行壓縮編碼，得到每個未壓縮的WAV 音頻樣本經過一次壓縮后形成的一次壓縮MP3音頻樣本。在此，現有的MP3編碼器如可采用 lame3. 99. 5MP3 編碼器，編碼比特率為brlpbrh可以為 64kbps、80kbps、96kbps、112kbps、 128kbps、160kbps以及192kbps中的一種編碼比特率。
[0046] ③獲取與每個未壓縮的WAV音頻樣本對應的二次壓縮MP3音頻樣本，具體過程為： ③ -1、利用MP3編碼器對每個未壓縮的WAV音頻樣本進行壓縮編碼，得到每個未壓縮的WAV 音頻樣本經過一次壓縮后形成的一次壓縮MP3音頻樣本，再利用現有的MP3解碼器將得到 的每個一次壓縮MP3音頻樣本解壓成WAV音頻樣本；③-2、利用MP3編碼器對步驟③-1中 解壓成的每個WAV音頻樣本進行壓縮編碼，得到與每個未壓縮的WAV音頻樣本對應的二次 壓縮MP3音頻樣本。
[0047] 其中，所述的步驟③-1和所述的步驟③-2中使用的MP3編碼器均與所述的步驟 ②中使用的MP3編碼器相同，且所述的步驟③-1中的MP3編碼器的編碼比特率與所述的步 驟②中的MP3編碼器的編碼比特率相同或不相同，所述的步驟③-2中的MP3編碼器的編碼 比特率與所述的步驟②中的MP3編碼器的編碼比特率相同。假如步驟②中使用的MP3編碼 器的編碼比特率為brliS128kbps，則步驟③-1中的MP3編碼器的編碼比特率br2i可以 為 64kbps、80kbps、96kbps、112kbps、128kbps、160kbps以及 192kbps中的一種編碼比特率， 步驟③-2中的MP3編碼器的編碼比特率br22S128kbps。
[0048] 在此，現有的MP3解碼器如可采用lame3. 99. 5MP3解碼器。
[0049] ④獲取與每個未壓縮的WAV音頻樣本對應的三次壓縮MP3音頻樣本，具體過程為： ④ -1、利用MP3編碼器對每個未壓縮的WAV音頻樣本進行壓縮編碼，得到每個未壓縮的WAV 音頻樣本經過一次壓縮后形成的一次壓縮MP3音頻樣本，再利用MP3解碼器將得到的每個 一次壓縮MP3音頻樣本解壓成WAV音頻樣本；④-2、利用MP3編碼器對步驟④-1中解壓成 的每個WAV音頻樣本進行壓縮編碼，得到與每個未壓縮的WAV音頻樣本對應的二次壓縮MP3 音頻樣本，再利用MP3解碼器將得到的每個二次壓縮MP3音頻樣本解壓成WAV音頻樣本； ④-3、利用MP3編碼器對步驟④-2中解壓成的每個WAV音頻樣本進行壓縮編碼，得到與每 個未壓縮的WAV音頻樣本對應的三次壓縮MP3音頻樣本。
[0050] 其中，所述的步驟④-1、所述的步驟④_2和