• 沒有找到結果。

第四章 改良式本質模態分解法之架構與分析結果

4.2 模擬訊號分析

4.2.2 非平穩訊號

(1) 振幅呈指數衰減之正弦波調幅訊號 )

3(t

S 為振幅呈指數衰減之正弦波調幅訊號,表示式如下:

) 6 sin(

15 )

( 0.4

3 t e t

S t  (4.3)

圖 4.10 訊號S3(t)原始波形

參數設定為預設值。圖 4.11 為利用立方雲線聯結中值點之分解結果,可 以明顯的看出S3(t)訊號之 IMF1 可視為近似為原訊號之分量,這與 2.1.2 節 描述一致:IMF 為一調幅(AM)訊號,因此 RIMD 也會視S3(t)訊號為一IMF 分量。計算相同條件下,EMD 與 RIMD 立方雲線拆解出之 IMF1 與原訊號 間之mse,如表 4.4 所示。

圖 4.11 S3(t)經RIMD 之分解結果

表 4.4 S3(t)分量mse 之比較表

分量 EMD RIMD 立方雲線

IMF1:15e0.4tsin(6t) 0.011 0.004

(2) 振幅呈指數衰減與平穩之正弦波訊號的疊加

S3(t)加入一頻率為7Hz 之正弦波疊加成訊號S4(t),表示式如下:

) 6 sin(

15 ) 14 sin(

11 )

( 0.4

4 t t e t

S    t  (4.4)

圖 4.12 訊號S4(t)原始波形

參數設定為預設值。利用立方雲線聯結中值點,分別以斜率比值與平均 值求得邊界極值之分解結果作比較,如圖 4.13 與圖 4.14 所示。可以明顯看 出S4(t)訊號前兩個 IMFs 可視為近似 7Hz 與振幅呈指數衰減之 3Hz 之正弦 波,比較兩者IMFs 之波形,可以發現前者整體較為對稱。比較相同條件下,

三者間(EMD、RIMD 斜率比值以及 RIMD 平均值)拆解出之前兩個 IMFs 與 原分量間之mse,結果如表 4.5 所示,可知利用 EMD 以及 RIMD 斜率比值 整體效果較佳。

圖 4.13 S4(t)經RIMD 之分解結果(斜率比值求邊界極值)

圖 4.14 S4(t)經RIMD 之分解結果(平均值求邊界極值)

表 4.5 S4(t)分量mse 之比較表

RIMD 立方雲線

分量 EMD

斜率比值 平均值

IMF1:11sin(14t) 0.414 0.463 0.398 IMF2:15e0.4tsin(6t) 0.431 0.411 3.879

(3) 兩振幅呈指數衰減之正弦波訊號的疊加

訊號S5(t)為兩振幅呈指數衰減之正弦波訊號的疊加,表示式如下:

) 6 sin(

13 ) 14 sin(

15 )

( 0.3 0.5

5 t e t e t

S t   t  (4.5)

圖 4.15 訊號S5(t)原始波形

參數設定為預設值。利用立方雲線聯結中值點,分別以斜率比值與平均 值求得邊界極值之分解結果作比較,如圖 4.16 與圖 4.17 所示。可以看出

)

5(t

S 訊號前兩個 IMFs 可視為近似振幅呈指數衰減之 7Hz 與 3Hz 正弦波,

且前者波形平穩性與對稱性表現較好。比較相同條件下,三者間(EMD、RIMD 斜率比值以及RIMD 平均值)拆解出之前兩個 IMFs 與原分量間之 mse,結果 如表 4.6 所示,可知 RIMD 斜率比值表現最佳。

圖 4.16 S5(t)經RIMD 之分解結果(斜率比值求邊界極值)

圖 4.17 S5(t)經RIMD 之分解結果(平均值求邊界極值)

表 4.6 S5(t)分量mse 之比較表

RIMD 立方雲線

分量 EMD

斜率比值 平均值

IMF1:15e0.3tsin(14t) 0.847 0.508 0.516 IMF2:13e0.5tsin(6t) 0.698 0.398 2.666

(4) 頻率遞增之正弦波調頻訊號 )

6(t

S 為頻率遞增之正弦波調頻訊號,表示式如下:

) ) 5 14 sin((

11 )

6(t t t

S    (4.6)

圖 4.18 訊號S6(t)原始波形

參數設定為預設值。圖 4.19 為利用立方雲線聯結中值點之分解結果,

可以明顯的看出S6(t)訊號之IMF1 可視為近似為原訊號之分量;這也與 2.1.2 節描述一致:IMF 為一調頻(FM)訊號,因此 RIMD 也會視S6(t)訊號為一IMF 分量。計算相同條件下,EMD 與 RIMD 立方雲線拆解出之 IMF1 與原訊號 間之mse,如表 4.7 所示。

表 4.7 S6(t)分量mse 之比較表

分量 EMD RIMD 立方雲線

IMF1:11sin((145t)t) 0.027 0.028

圖 4.19 S6(t)經RIMD 之分解結果

(5) 兩頻率遞增之正弦波調頻訊號的疊加

訊號S7(t)為兩頻率遞增之正弦波的疊加,表示式如下:

) ) 2 6 sin((

5 ) ) 5 14 sin((

11 )

7(t t t t t

S       (4.7)

圖 4.20 訊號S7(t)原始波形

參數設定為預設值。利用立方雲線聯結中值點,分別以斜率比值與平均 值求得邊界極值之分解結果作比較,如圖 4.21 與圖 4.22 所示,可以看出

)

7(t

S 訊號前兩個 IMFs 可視為兩近似頻率遞增之正弦波調頻訊號。前者

IMF1、IMF2 之波形在後方邊界處皆有小振盪;至於後者則是在 IMF1 有較 佳(對稱性)之表現。比較相同條件下,三者間(EMD、RIMD 斜率比值以及 RIMD 平均值)拆解出之前兩個 IMFs 與原分量間之 mse,結果如表 4.8 所示,

可以發現EMD 整體表現較佳。

圖 4.21 S7(t)經RIMD 之分解結果(斜率比值求邊界極值)

圖 4.22 S7(t)經RIMD 之分解結果(平均比值求邊界極值)

表 4.8 S7(t)分量mse 之比較表

RIMD 立方雲線

分量 EMD

斜率比值 平均值

IMF1:11sin((145t)t) 0.120 0.225 0.120 IMF2:5sin((62t)t) 0.117 0.218 0.756

(6) 振幅呈指數衰減與頻率遞增之正弦波訊號的疊加

訊號S8(t)為振幅呈指數衰減與頻率遞增之正弦波的疊加,表示式如下:

) ) 5 14 sin((

11 ) 6 sin(

13 )

( 0.5

8 t e t t t

S t     (4.8)

圖 4.23 訊號S8(t)原始波形

參數設定為預設值。利用立方雲線聯結中值點,分別以斜率比值與平均 值求得邊界極值之分解結果作比較,如圖 4.24 與圖 4.25 所示。可以看出

)

8(t

S 訊號前兩個IMFs 可視為頻率遞增與振幅呈指數衰減之正弦波訊號,前 者波形在整體之平穩性與對稱性較好;後者依然在IMF1 有較佳之表現。比 較相同條件下,三者間(EMD、RIMD 斜率比值以及 RIMD 平均值)拆解出之 前兩個IMFs 與原分量間之 mse,結果如表 4.9 所示,可以發現 RIMD 斜率 比值表現最佳。

圖 4.24 S8(t)經RIMD 之分解結果(斜率比值求邊界極值)

圖 4.25 S8(t)經RIMD 之分解結果(平均值求邊界極值)

表 4.9 S8(t)分量mse 之比較表

RIMD 立方雲線

分量 EMD

斜率比值 平均值

IMF1:11sin((145t)t) 0.359 0.288 0.221 IMF2:13e0.5tsin(6t) 0.655 0.545 2.410

相關文件