MATLAB関数の角度、アンラップ



Matlab Function Angle



一つ、角度

位相角



  1. 文法
    P =角度(Z)
  2. 説明
    P = angle(Z)は、複素数配列Zの各要素の位相角(ラジアン)を返します。角度は±πの間です。
  3. 複素数Zの場合、振幅Rと位相角シータは次の式で与えられます。
    R =絶対値(Z)
    θ=角度(Z)
  4. に基づくことができます
    Z =R。* exp(i * theta)は、元の複素数Zに変換されます。

  5. 行列要素の位相角
    複素数値の行列を作成し、各要素の位相角を計算します。
    Z = [1-1 i 2 + 1i 3-1i 4 + 1i
    1 + 2i 2-2i 3 + 2i 4-2i
    1-3i 2 + 3i 3-3i 4 + 3i
    1 + 4i 2-4i 3 + 4i 4-4i];
  6. P =角度(Z)
    P = 4×4

-0.7854 0.4636 -0.3218 0.2450
1.1071 -0.7854 0.5880 -0.4636
-1.2490 0.9828 -0.7854 0.6435
1.3258 -1.1071 0.9273 -0.7854

  1. アルゴリズム
    角度関数は、複素数z = x + iyを使用し、atan2(y、x)を計算して、正のx軸上のxy平面上で、原点から点(x、y)までの光線を見つけます。この位相角は、複素対数の虚数部でもあります。

総括する



1.アンラップとは

システムの位相周波数特性を計算するには、アークタンジェント関数を使用します。コンピューターのアークタンジェント関数は、第1象限と第2象限の角度が0〜piであり、第3象限と第4象限の角度が0〜-piであることを指定します。角度が0から2piに変化しても、実際の結果が0〜piで、次に-pi〜0の場合、w = piでジャンプが発生し、ジャンプ振幅は2piで、位相巻線と呼ばれます。アンラップ(w)は、位相がpiでジャンプしないようにアンラップすることであり、それによって実際の位相変化を反映します。

2、アンラップ機能



アンラップ関数は、データフェーズのジャンプをチェックアウトし、ジャンプを修正します。実際、データがジャンプするかどうかをチェックするときは、標準のunwrap(pha、tol)があります。この通行料が標準です。デフォルトでは、標準はpiです。つまり、unwrapは、データの前後の2つのポイントの差がtolを超えていることを確認すると、ジャンプがあると見なします。次に、後続のデータに2piを加えたものまたは2piを引いたものが連続データになるようにデータが処理されます。ほとんどの場合、標準のpiを使用する方が適切であるため、2番目のパラメーターをデフォルトにすることができます。特別な場合には、tolを設定することで識別基準を調整できます。
3、アンラップ使用

unwrap関数は、数値のシーケンスだけでなく、行列に対しても機能して、行列の各行または列の位相補正を実装できるため、完全な形式はunwrap(pha、tol、 番号 )。最後のパラメーターdimは、行列の行または列をアンラップするかどうかを示します。各列を操作する場合、3番目のパラメーターをデフォルトにして、各行を操作する場合は、3番目のパラメーターに2を入力します。たとえば、unwrap(pha、[]、2)は、次の位相補正を実装することを意味します。デフォルトの検出ジャンプ標準を使用した、行列phiの位相データの各行。

2、アンラップ

  1. アンラップ
    位相角を修正して、より滑らかな位相画像を生成します
  2. 構文
    Q =アンラップ(P)
    Q =アンラップ(P、tol)
    Q =アンラップ(P、[]、dim)
    Q =アンラップ(P、tol、dim)
  3. 説明
    Pの連続する要素間の絶対ジャンプが、デフォルトのジャンプ許容値であるπラジアン以上の場合、
    Q = unwrap(P)は、±2πの倍数を追加することにより、ベクトルPのラジアン位相角を補正します。
    Pが行列の場合、展開すると、操作は列で実行されます。
    Pが多次元配列の場合、展開は最初の非1次元で実行されます。
  4. Q = unwrap(P、tol)は、デフォルト値πの代わりにジャンプ許容値tolを使用します。
    Q = unwrap(P、[]、dim)デフォルトの許容値を使用して、ケラレのコーナーに沿ってアンラップします。
    Q = unwrap(P、tol、dim)は、tolのジャンプ許容値を使用します。

次の位相データは、3次伝達関数の周波数応答から取得されます。位相曲線は、w = 3.0とw = 3.5の間で3.5873ラジアンで-1.8621から1.7252に遷移します。

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

アンラップを使用して位相角を修正すると、得られるジャンプは2.6959であり、デフォルトのジャンプ許容値πよりも小さくなります。グラフは、元の曲線に新しい曲線をプロットします。

例2

配列P要素での不連続性を除いて、滑らかに増加する位相角を特徴とします(3,1)および(1,2)

P = [ 0 7.0686 1.5708 2.3562 0.1963 0.9817 1.7671 2.5525 6.6759 1.1781 1.9635 2.7489 0.5890 1.3744 2.1598 2.9452 ]

関数Q = unwrap(P)これらの不連続性を排除します。

Q = 0 7.0686 1.5708 2.3562 0.1963 7.2649 1.7671 2.5525 0.3927 7.4613 1.9635 2.7489 0.5890 7.6576 2.1598 2.9452