位相角 0 ～ $6\pi$を使用して,らせんの $$x$$座標と $$y$$座標を定義します。らせんをプロットします。

t = linspace(0, 6 *π,201);x = t /π。* cos (t);y = t /π。* sin (t);情节(x, y)

関数量化を使用して, $$x$$座標と $$y$$座標から,らせんの位相角を求めます。関数量化は, $-\pi$～ $\pi$の閉区間にある角度値を返します。

P =量化(y、x);情节(t, P)

このプロットには不連続点のあることがわかります。打开を使用して不連続性を排除します。打开は,Pの連続する要素間の位相差が跳びのしきい値である $\pi$ラジアン以上の場合に, $\pm 2\pi$の倍数を加算します。シフト後の位相角问は,区間0 ~ $6\pi$にあります。

Q =打开(P);情节(t, Q)

周波数応答の位相角をシフトします。この位相曲線には2つの跳びがあります。最初の跳びはW = 3とW = 3.4の間の3.4250ラジアンで2番目の跳びはW = 5とW = 5.4の間の6.3420ラジアンです。位相曲線をプロットします。

清晰;关闭所有；W = [0:0.4:3, 4:0.4:5, 4:0.4:7];P = [-1.5723 -1.5747 -1.5790 -1.5852 -1.5922 -1.6044 -1.6269 -1.6998 1.7252 1.5989 1.5916 1.5708 1.5582 -4.7838 -4.8143 -4.8456 -4.8764 -4.9002];情节(W、P,“bo - - - - - -”）

打开で既定の跳びのしきい値 $$\pi$$ラジアンを使用して,位相角をシフトします。シフト後の位相曲線をプロットします。両方の跳びは,跳びのしきい値 $$\pi$$ラジアンより大きいため,シフトされます。

情节(W,打开(P),“ro - - - - - -”）

次に,跳びのしきい値に5ラジアンを使用して位相角をシフトします。シフト後の位相曲線をプロットします。最初の跳びは,跳びのしきい値5ラジアンより小さいため,シフトされません。

情节(W,打开(P, 5),“ro - - - - - -”）

位相角を含む2列の行列Pを定義します。

P = [0 7.07;0.19 - 0.98;6.67 - 1.18;0.59 - 1.37;0.78 - 1.56)

P =5×20 7.0700 0.1900 0.9800 6.6700 1.1800 0.5900 1.3700 0.7800 1.5600

位相角P(1、2)= 7.07およびP (3,1) = 6.67には,残りのデータと比較して, $$\pi$$より大きい位相差があります。

まず,要素を列単位で比較して,位相角をアンラップします。引数昏暗的を1に指定します。2番目の引数を［］に指定することで,既定の跳びのしきい値 $$\pi$$を使用します。

昏暗的= 1;P1 =打开(P,[],昏暗的)

P1 =5×20 7.0700 0.1900 7.2632 0.3868 7.4632 0.900 7.6532 0.7800 7.8432

列単位ではなく行単位で位相角をシフトするには,昏暗的を1ではなく2に指定します。

昏暗的= 2;P2 =打开(P1,[],昏暗的)

P2 =5×20 0.7868 0.1900 0.9800 0.3868 1.1800 0.900 1.3700 0.7800 1.5600