一種基于同步壓縮變換重構(gòu)的幅值校正方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于信號(hào)處理技術(shù)領(lǐng)域，具體設(shè)及一種基于同步壓縮變換重構(gòu)的幅值校正 方法。
【背景技術(shù)】
[0002] 時(shí)頻分析方法是處理非平穩(wěn)信號(hào)的有力工具，目前常用的有經(jīng)典的時(shí)頻方法，如 短時(shí)傅里葉變換和小波變換，近幾年有學(xué)者通過(guò)運(yùn)些經(jīng)典的方法，得出了一種新的時(shí)頻分 析方法：同步壓縮變換。該方法類(lèi)似經(jīng)驗(yàn)?zāi)Ｊ椒纸?EMD)，可W進(jìn)行信號(hào)分解，得到單成分信 號(hào)的時(shí)域圖，即對(duì)信號(hào)進(jìn)行重構(gòu)。
[0003] 根據(jù)文獻(xiàn)【latsenko D,Mcclintock P V E,Stefanovska A.Linear and synchrosqueezed time-frequency representations revisited:Overview,standards of use.resolution.reconstruction , concentration,and algorithms[J].Digital Signal Processing，2015，42.】中可W總結(jié)出一般的同步壓縮變換，但是一般的同步壓縮 變換中使用了瞬時(shí)頻率的近似值，該近似值大大降低了時(shí)域重構(gòu)信號(hào)的精度，所得到的信 號(hào)的幅值相比原始信號(hào)具有很大的差別。
【發(fā)明內(nèi)容】
[0004] 為了克服上述現(xiàn)有技術(shù)的缺點(diǎn)，本發(fā)明的目的在于提供一種基于同步壓縮變換重 構(gòu)的幅值校正方法，校正所得到信號(hào)的幅值，提高時(shí)域重構(gòu)信號(hào)的精度。
[000引為了達(dá)到上述目的，本發(fā)明采取的技術(shù)方案為：
[0006] -種基于同步壓縮變換重構(gòu)的幅值校正方法，包括W下步驟：
[0007] 步驟一，對(duì)線性調(diào)頻信號(hào)s(t)進(jìn)行短時(shí)傅里葉變換（STFT)，得到短時(shí)傅里葉變換 系數(shù)Vs，其中信號(hào)的形式為s(t)=Asin(2JT(i) (t))，A為信號(hào)的幅值，t為時(shí)間參數(shù)，Φ (t)為 信號(hào)的相位函數(shù)；
[000引步驟二，用得到的短時(shí)傅里葉變換系數(shù)Vs估算瞬時(shí)頻率，即
得 到估計(jì)瞬時(shí)頻率氣，在該估計(jì)瞬時(shí)頻率的基礎(chǔ)上進(jìn)行改進(jìn)，用ω =? +0 0旅一了)來(lái)計(jì) 算準(zhǔn)確瞬時(shí)頻率S，其中Φ"(τ)為相位函數(shù)的二階導(dǎo)數(shù)，τ為過(guò)渡參數(shù)，不影響結(jié)果；
[0009] 步驟Ξ，將信號(hào)的所需要分析的頻率范圍設(shè)定為[Wmin，Wmax]，將該范圍劃分為多 個(gè)Wwk為中屯、頻率，ΚΔ ω為間隔的頻率段，通過(guò)c〇k=化-k〇)A ω來(lái)計(jì)算中屯、頻率，其中
,Δ ω在不加大計(jì)算量的前提下盡可能的 取小值；
[0010] 步驟四，整合前Ξ步計(jì)算出來(lái)的參數(shù)V S、竊、ω k和Δ ω，通過(guò)公式
完成同步壓縮變換，即將所有對(duì)應(yīng)同一 頻率段的準(zhǔn)確瞬時(shí)頻率的短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為同步壓縮系數(shù)，如果某幾個(gè)準(zhǔn)確瞬時(shí) 頻率對(duì)應(yīng)第k個(gè)頻率段，那么就將相應(yīng)的所有短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為和該頻率段相對(duì) 應(yīng)的第k個(gè)同步壓縮系數(shù)，從而得到同步壓縮變換系數(shù)，公式中j是計(jì)數(shù)參數(shù)，在該公式中只 有計(jì)數(shù)作用，不影響結(jié)果，％為通過(guò)步驟二計(jì)算出來(lái)的第j個(gè)準(zhǔn)確瞬時(shí)頻率，Vs(t，W^為在 時(shí)間為t和頻率為ω川寸的短時(shí)傅里葉變換系數(shù)，而ω J是屬于所需分析的頻率范圍，為已知 量；
[0011] 步驟五，提取同步壓縮變換的系數(shù)脊線，并將脊線附近的系數(shù)進(jìn)行整合，得到相應(yīng) 的信號(hào)成分。
[0012] 所述的步驟一中線性調(diào)頻信號(hào)s(t)為單成分信號(hào)或多成分信號(hào)。
[001引所述的步驟ミ中Δ ω為0.巧化.2之間。
[0014] 本發(fā)明的有益效果為：由于本發(fā)明在步驟二中對(duì)瞬時(shí)頻率的計(jì)算進(jìn)行了改進(jìn)，能 夠得到準(zhǔn)確的瞬時(shí)頻率，因此提高了時(shí)域重構(gòu)信號(hào)的精度。
【附圖說(shuō)明】
[0015] 圖1為使用本發(fā)明方法對(duì)單成分信號(hào)進(jìn)行信號(hào)重構(gòu)，其中（a)、（b)和（C)圖分別為 仿真出的四階相位信號(hào)、Ξ階相位信號(hào)和二階相位信號(hào)。
[0016] 圖2為使用一般的同步壓縮變換對(duì)單成分信號(hào)進(jìn)行信號(hào)重構(gòu)，其中（a)、（b)和(C) 圖分別為仿真出的四階相位信號(hào)、Ξ階相位信號(hào)和二階相位信號(hào)。
[0017] 圖3為使用本發(fā)明方法對(duì)多成分信號(hào)進(jìn)行分解得到的各個(gè)成分信號(hào)與原始信號(hào)的 對(duì)比圖。
[0018] 圖4為使用一般的同步壓縮變換對(duì)多成分信號(hào)進(jìn)行分解得到的各個(gè)成分信號(hào)與原 始信號(hào)的對(duì)比圖。
【具體實(shí)施方式】
[0019] 下面結(jié)合附圖和實(shí)施例對(duì)本發(fā)明作進(jìn)一步詳細(xì)闡述。
[0020] 實(shí)施例一
[0021] -種基于同步壓縮變換重構(gòu)的幅值校正方法，包括W下步驟：
[0022] 步驟一，分別仿真不同階相位的單成分信號(hào)，其中單成分信號(hào)分別為四階相位信 號(hào) s(t) = sin(23i(250t+100t4))，s階相位信號(hào) s(t) = sin(23i(250t+100t3))和二階相位信 號(hào)s(t) = sin(化(250t+l(K)t2))，對(duì)線性調(diào)頻信號(hào)s(t)進(jìn)行短時(shí)傅里葉變換(STFT)，得到短 時(shí)傅里葉變換系數(shù)Vs;
[0023] 步驟二，用得到的短時(shí)傅里葉變換系數(shù)Vs估算瞬時(shí)頻率，即
得 到估計(jì)瞬時(shí)頻率聳，在該估計(jì)瞬時(shí)頻率的基礎(chǔ)上進(jìn)行改進(jìn)，用
農(nóng)計(jì) 算準(zhǔn)確瞬時(shí)頻率諒，其中Φ"(τ)為相位函數(shù)的二階導(dǎo)數(shù)，τ為過(guò)渡參數(shù)，不影響結(jié)果；
[0024] 步驟Ξ，將信號(hào)的所需要分析的頻率范圍設(shè)定為[Wmin，Wmax]，將該范圍劃分為多 個(gè)Wwk為中屯、頻率，ΚΔ ω為間隔的頻率段，通過(guò)c〇k=化-ko)A ω來(lái)計(jì)算中屯、頻率，其中
，.Δ ω -般應(yīng)該在不加大計(jì)算量的前提下盡 可能的取小值，本實(shí)施例中為0.135;
[002引步驟四，整合前Ξ步計(jì)算出來(lái)的參數(shù)Vs、?、Wk和Δ ω，通過(guò)公式
完成同步壓縮變換，即將所有對(duì)應(yīng)同一 頻率段的準(zhǔn)確瞬時(shí)頻率的短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為同步壓縮系數(shù)，即如果某幾個(gè)準(zhǔn)確瞬 時(shí)頻率對(duì)應(yīng)第k個(gè)頻率段，那么就將相應(yīng)的所有短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為和該頻率段相 對(duì)應(yīng)的第k個(gè)同步壓縮系數(shù)，從而得到同步壓縮變換系數(shù)，如圖1中（a)、（b)、（c)中的SWFT 圖；其中j是計(jì)數(shù)參數(shù)，在該公式中只有計(jì)數(shù)作用，不影響結(jié)果，鳴為通過(guò)步驟二計(jì)算出來(lái)的 第j個(gè)準(zhǔn)確瞬時(shí)頻率，Vs(t，Wj)為在時(shí)間為t和頻率為Wj時(shí)的短時(shí)傅里葉變換系數(shù)，而COj 是屬于所需分析的頻率范圍，為已知量；
[0026] 步驟五，提取同步壓縮變換的系數(shù)脊線，并將脊線附近的系數(shù)進(jìn)行整合，得到相應(yīng) 的信號(hào)成分，如圖1中(a)、(b)、（c)中的重構(gòu)與原始對(duì)比圖。
[0027] 由一般的同步壓縮變換所得到的結(jié)果如圖2所示，對(duì)比圖1和圖2,可W看出本發(fā)明 方法相對(duì)一般的方法得到的重構(gòu)信號(hào)與原始信號(hào)更加一致，幅值精度得到了極大提高。
[002引實(shí)施例二
[0029] -種基于同步壓縮變換重構(gòu)的幅值校正方法，包括W下步驟：
[0030] 步驟一，仿真一個(gè)雙成分信號(hào) s = sl+s2,其中 sl = sin(化(250t+100t3))，s2 = sin (2JT(l〇〇t+l(K)t2))，其中對(duì)信號(hào)s(t)進(jìn)行短時(shí)傅里葉變換（STFT)，得到短時(shí)傅里葉變換系 數(shù)Vs;
[0031] 步驟二，用得到的短時(shí)傅里葉變換系數(shù)Vs估算瞬時(shí)頻率，即
得到 估計(jì)瞬時(shí)頻率4，在該估計(jì)瞬時(shí)頻率的基礎(chǔ)上進(jìn)行改進(jìn)，用
來(lái)計(jì)算 準(zhǔn)確瞬時(shí)頻率ffl，:其中φ"(τ)為相位函數(shù)的二階導(dǎo)數(shù)，τ為過(guò)渡參數(shù)，不影響結(jié)果；
[0032] 步驟Ξ，將信號(hào)的所需要分析的頻率范圍設(shè)定為[Wmin，Wmax]，將該范圍劃分為多 個(gè)Wwk為中屯、頻率，ΚΔ ω為間隔的頻率段，通過(guò)c〇k=化-k〇)A ω來(lái)計(jì)算中屯、頻率，其中
，Δ ω-般應(yīng)該在不加大計(jì)算量的前提下盡 可能的取小值，本實(shí)施例中取為0.135;
[0033] 步驟四，整合前Ξ步計(jì)算出來(lái)的參數(shù)Vs、孩、Wk和Δ ω，通過(guò)公式
完成同步壓縮變換，即將所有對(duì)應(yīng)同一 頻率段的準(zhǔn)確瞬時(shí)頻率的短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為同步壓縮系數(shù)，即如果某幾個(gè)準(zhǔn)確瞬 時(shí)頻率對(duì)應(yīng)第k個(gè)頻率段，那么就將相應(yīng)的所有短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為和該頻率段相 對(duì)應(yīng)的第k個(gè)同步壓縮系數(shù)，從而得到同步壓縮變換系數(shù);其中j是計(jì)數(shù)參數(shù)，在該公式中只 有計(jì)數(shù)作用，不影響結(jié)果，兩;為通過(guò)步驟二計(jì)算出來(lái)的第j個(gè)準(zhǔn)確瞬時(shí)頻率，Vs(t，W^為在 時(shí)間為t和頻率為ω川寸的短時(shí)傅里葉變換系數(shù)，而ω J是屬于所需分析的頻率范圍，為已知 量；
[0034] 步驟五，提取同步壓縮變換的系數(shù)脊線，并將脊線附近的系數(shù)進(jìn)行整合，得到相應(yīng) 的信號(hào)成分，如圖3所示。
[0035] 由一般的同步壓縮變換所得到的結(jié)果如圖4所示，對(duì)比圖3和圖4,可W看出本發(fā)明 方法相對(duì)一般的方法得到的重構(gòu)信號(hào)與原始信號(hào)更加一致，幅值精度得到了極大提高。
【主權(quán)項(xiàng)】
1. 一種基于同步壓縮變換重構(gòu)的幅值校正方法，其特征在于，包括以下步驟： 步驟一，對(duì)線性調(diào)頻信號(hào)S(t)進(jìn)行短時(shí)傅里葉變換（STFT)，得到短時(shí)傅里葉變換系數(shù) Vs，其中信號(hào)的形式為s(t)=A Sin(23i(Ht))，A為信號(hào)的幅值，t為時(shí)間參數(shù)，Φ(〇為信號(hào) 的相位函數(shù)； 步驟二，用得到的短時(shí)傅里葉變換系數(shù)1估算瞬時(shí)頻率，即，得到估計(jì) 瞬時(shí)頻率?,,在該估計(jì)瞬時(shí)頻率的基礎(chǔ)上進(jìn)行改進(jìn)，用丨來(lái)計(jì)算準(zhǔn)確瞬 時(shí)頻率其中φ〃（τ)為相位函數(shù)的二階導(dǎo)數(shù)，τ為過(guò)渡參數(shù)，不影響結(jié)果； 步驟三，將信號(hào)的所需要分析的頻率范圍設(shè)定為[?min，comax]，將該范圍劃分為多個(gè)以 wk為中心頻率，以Δ ω為間隔的頻率段，通過(guò)cok= (k-k〇) Δ ω來(lái)計(jì)算中心頻率，其中，△ω在不加大計(jì)算量的前提下盡可能的 取小值； 步驟四，整合前三步計(jì)算出來(lái)的參數(shù)Vs、孩:、wk和Δ ω，通過(guò)公式完成同步壓縮變換，即將所有對(duì)應(yīng)同 一頻率段的準(zhǔn)確瞬時(shí)頻率的短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為同步壓縮系數(shù)，如果某幾個(gè)準(zhǔn)確瞬 時(shí)頻率對(duì)應(yīng)第k個(gè)頻率段，那么就將相應(yīng)的所有短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為和該頻率段相 對(duì)應(yīng)的第k個(gè)同步壓縮系數(shù)，從而得到同步壓縮變換系數(shù)，公式中j是計(jì)數(shù)參數(shù)，在該公式中 只有計(jì)數(shù)作用，不影響結(jié)果，％為通過(guò)步驟二計(jì)算出來(lái)的第j個(gè)準(zhǔn)確瞬時(shí)頻率，V s(t，〇^)為 在時(shí)間為t和頻率為時(shí)的短時(shí)傅里葉變換系數(shù)，而是屬于所需分析的頻率范圍，為已 知量； 步驟五，提取同步壓縮變換的系數(shù)脊線，并將脊線附近的系數(shù)進(jìn)行整合，得到相應(yīng)的信 號(hào)成分。2. 根據(jù)權(quán)利要求1所述的一種基于同步壓縮變換重構(gòu)的幅值校正方法，其特征在于:所 述的步驟一中線性調(diào)頻信號(hào)s(t)為單成分信號(hào)或多成分信號(hào)。3. 根據(jù)權(quán)利要求1所述的一種基于同步壓縮變換重構(gòu)的幅值校正方法，其特征在于:所 述的步驟三中Δ ω為0.1到0.2之間。
【專利摘要】一種基于同步壓縮變換重構(gòu)的幅值校正方法，對(duì)線性調(diào)頻信號(hào)進(jìn)行短時(shí)傅里葉變換，得到短時(shí)傅里葉變換系數(shù)，然后估算瞬時(shí)頻率，在該估計(jì)瞬時(shí)頻率的基礎(chǔ)上進(jìn)行改進(jìn)，計(jì)算準(zhǔn)確瞬時(shí)頻率，再計(jì)算中心頻率，然后完成同步壓縮變換，即將所有對(duì)應(yīng)同一頻率段的準(zhǔn)確瞬時(shí)頻率的短時(shí)傅里葉變換系數(shù)轉(zhuǎn)化為同步壓縮系數(shù)，最后提取同步壓縮變換的系數(shù)脊線，并將脊線附近的系數(shù)進(jìn)行整合，得到相應(yīng)的信號(hào)成分，由于本發(fā)明對(duì)瞬時(shí)頻率的計(jì)算進(jìn)行了改進(jìn)，能夠得到準(zhǔn)確的瞬時(shí)頻率，因此提高了時(shí)域重構(gòu)信號(hào)的精度。
【IPC分類(lèi)】G06F17/14
【公開(kāi)號(hào)】CN105446940
【申請(qǐng)?zhí)枴緾N201510821373
【發(fā)明人】廖與禾, 李帥, 陳鵬, 孫鵬, 林京, 雷亞國(guó)
【申請(qǐng)人】西安交通大學(xué)
【公開(kāi)日】2016年3月30日
【申請(qǐng)日】2015年11月23日