matlab之unwrap函数
要計(jì)算一個(gè)系統(tǒng)相頻特性,就要用到反正切函數(shù),計(jì)算機(jī)中反正切函數(shù)規(guī)定,在一、二象限中的角度為0~pi,三四象限的角度為0~-pi。
若一個(gè)角度從0變到2pi,但實(shí)際得到的結(jié)果是0~pi,再由-pi~0,在w=pi處發(fā)生跳變,跳變幅度為2pi,這就叫相位的卷繞。
unwrap(w)就是解卷繞,使相位在pi處不發(fā)生跳變,從而反應(yīng)出真實(shí)的相位變化 "
查看 unwrap 的幫助文檔可以發(fā)現(xiàn) unwrap 還可以輸入一個(gè)參數(shù) tol,默認(rèn)tol = pi。也可以根據(jù)情況修改下 tol。
以下純屬個(gè)人見解,若有錯(cuò)誤,謝謝指正:
unwrap也就是如果前后兩個(gè)數(shù)據(jù)之差大于默認(rèn)值(一般為π,也可自己在調(diào)用函數(shù)時(shí)設(shè)定),那么就把第二個(gè)值加上或者減去2π,反正最終結(jié)果是它倆之差小于默認(rèn)值即可。【特別注意】這個(gè)函數(shù)并不是把[-π,π]范圍的數(shù)據(jù)變成[0,2π]的數(shù)據(jù),感覺有點(diǎn)類似于讓數(shù)據(jù)具有一定的單調(diào)性。
看看matlab里面的實(shí)例:
?w = [0:.2:3,3.5:1:10];?
p = [ ? ?0
? ? ?-1.5728
? ? ?-1.5747
? ? ?-1.5772
? ? ?-1.5790
? ? ?-1.5816
? ? ?-1.5852
? ? ?-1.5877
? ? ?-1.5922
? ? ?-1.5976
? ? ?-1.6044
? ? ?-1.6129
? ? ?-1.6269
? ? ?-1.6512
? ? ?-1.6998
? ? ?-1.8621
? ? ? 1.7252
? ? ? 1.6124
? ? ? 1.5930
? ? ? 1.5916
? ? ? 1.5708
? ? ? 1.5708
? ? ? 1.5708 ];
semilogx(w,p,'b*-'), hold
當(dāng)我進(jìn)行unwrap以后:
可以對(duì)比這兩張圖,第一張明顯有跳變,而第二張就比較平滑,類似于單調(diào)性比較好。
此外unwrap的調(diào)用方法有三種,可以自己設(shè)置跳變的閾值,默認(rèn)是π。
總結(jié)
以上是生活随笔為你收集整理的matlab之unwrap函数的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 安卓 日志 路径
- 下一篇: linux里的run-level,lin