الطور والتردد اللحظي

في عام 1922، ووفقًا لناهين، عرّف جون رينشو كارسون التردد اللحظي للإشارة بأنه " المشتق الزمني لزاوية طور الإشارة". في تعديل التردد ، يصف التردد اللحظي التردد الذي يتغير أعلى وأسفل تردد الموجة الحاملة، عند تردد النغمة الصوتية. [ 1 ]

يُعدّ كلٌّ من الطور اللحظي والتردد من المفاهيم المهمة في معالجة الإشارات ، ويظهران في سياق تمثيل وتحليل الدوال المتغيرة مع الزمن. [ 2 ] الطور اللحظي (المعروف أيضًا بالطور المحلي أو ببساطة الطور ) لدالة ذات قيم مركبة s ( t ) هو الدالة ذات القيم الحقيقية :

φ(ت)=arg{s(ت)}،{\displaystyle \varphi (t)=\arg\{s(t)\},}

حيث يمثل arg دالة الوسيط المركب . التردد اللحظي هو معدل التغير الزمني للطور اللحظي.

وبالنسبة للدالة الحقيقية s ( t )، يتم تحديدها من التمثيل التحليلي للدالة ، s a ( t ): [ 3 ]

φ(ت)=arg{sأ(ت)}=arg{s(ت)+جs^(ت)}،{\displaystyle {\begin{aligned}\varphi (t)&=\arg\{s_{\mathrm {a} }(t)\}\\[4pt]&=\arg\{s(t)+j{\hat {s}}(t)\},\end{aligned}}}

أينs^(ت){\displaystyle {\hat {s}}(t)}يمثل تحويل هيلبرت لـ s ( t ).

عندما تُقيّد الدالة φ ( t ) بقيمتها الأساسية ، إما في الفترة [−π , π ] أو [ 0 , 2π ] ، تُسمى طورًا مُلتفًا . وإلا تُسمى طورًا غير مُلتف ، وهي دالة متصلة للمتغير t ، بافتراض أن s <sub>a</sub> ( t ) دالة متصلة لـ t . ما لم يُذكر خلاف ذلك، يُستنتج الشكل المتصل.

الطور اللحظي مقابل الزمن. تحتوي الدالة على انقطاعين حقيقيين بزاوية 180 درجة عند الزمنين 21 و59، مما يدل على عبور السعة للصفر. أما "الانقطاعات" بزاوية 360 درجة عند الأزمنة 19 و37 و91 فهي ناتجة عن التفاف الطور.
الطور اللحظي لموجة معدلة التردد: MSK (مفتاح الإزاحة الأدنى). يتم ببساطة تكرار رسم بياني "ملفوف" بزاوية 360 درجة عموديًا مرتين إضافيتين، مما يخلق وهم رسم بياني غير ملفوف، ولكن باستخدام 3 × 360 درجة فقط من المحور الرأسي.

أمثلة

المثال 1

s(ت)=أكوس(ωت+θ)،{\displaystyle s(t)=A\cos(\omega t+\theta ),}

حيث ω > 0.

sأ(ت)=أهـج(ωت+θ)،φ(ت)=ωت+θ.{\displaystyle {\begin{align}s_{\mathrm {a} }(t)&=Ae^{j(\omega t+\theta )},\\\varphi (t)&=\omega t+\theta .\end{محاذاة}}}

في هذا المثال البسيط للموجة الجيبية، يُشار إلى الثابت θ عادةً باسم الطور أو إزاحة الطور . φ ( t ) دالة للزمن، بينما θ ليست كذلك. في المثال التالي، نلاحظ أيضًا أن إزاحة الطور لموجة جيبية حقيقية القيمة تكون غامضة ما لم يتم تحديد مرجع (الجيب أو جيب التمام). أما φ ( t ) فهي مُعرَّفة بشكل واضح لا لبس فيه.

المثال 2

s(ت)=أالخطيئة(ωت)=أكوس(ωت-π2)،{\displaystyle s(t)=A\sin(\omega t)=A\cos \left(\omega t-{\frac {\pi }{2}}\right),}

حيث ω > 0.

sأ(ت)=أهـج(ωت-π2)،φ(ت)=ωت-π2.{\displaystyle {\begin{aligned}s_{\mathrm {a} }(t)&=Ae^{j\left(\omega t-{\frac {\pi }{2}}\right)},\\\varphi (t)&=\omega t-{\frac {\pi }{2}}.\end{aligned}}}

في كلا المثالين ، تتوافق القيم العظمى المحلية لـ s ( t ) مع φ ( t ) =  2πN لقيم N الصحيحة . وهذا له تطبيقات في مجال رؤية الحاسوب . 

التركيبات

يُعرَّف التردد الزاوي اللحظي على النحو التالي:

ω(ت)=دφ(ت)دت،{\displaystyle \omega (t)={\frac {d\varphi (t)}{dt}},}

ويُعرَّف التردد اللحظي (العادي) على النحو التالي:

و(ت)=12πω(ت)=12πدφ(ت)دت{\displaystyle f(t)={\frac {1}{2\pi }}\omega (t)={\frac {1}{2\pi }}{\frac {d\varphi (t)}{dt}}}

حيث يجب أن تكون φ ( t ) هي المرحلة غير الملفوفة ؛ وإلا، إذا كانت φ ( t ) ملفوفة، فإن الانقطاعات في φ ( t ) ستؤدي إلى نبضات ديراك دلتا في f ( t ).

العملية العكسية، التي تكشف دائماً عن الطور، هي:

φ(ت)=-تω(τ)دτ=2π-تو(τ)دτ=-0ω(τ)دτ+0تω(τ)دτ=φ(0)+0تω(τ)دτ.// )\,d\tau +\int _{0}^{t}\omega (\tau )\,d\tau \\[5pt]&=\varphi (0)+\int _{0}^{t}\omega (\tau )\,d\tau .\end{محاذاة}}}

يمكن اشتقاق هذا التردد اللحظي، ω ( t )، مباشرة من الأجزاء الحقيقية والخيالية لـ sa ( t ) ، بدلاً من arg المركب دون القلق بشأن فك الطور.

φ(ت)=arg{sأ(ت)}=atan2(أنام[sأ(ت)]،Rهـ[sأ(ت)])+2م1π=دالة الظل العكسي(أنام[sأ(ت)]Rهـ[sأ(ت)])+م2π{\displaystyle {\begin{align}\varphi (t)&=\arg\{s_{\mathrm {a} }(t)\}\\[4pt]&=\operatorname {atan2} ({\mathcal {Im}[s_{\mathrm {a} }(t)]،{\mathcal {Re}}[s_{\mathrm {a} }(t)])+2m_{1}\pi \\[4pt]&=\arctan \left({\frac {{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{{\mathcal {Re}}[s_{\mathrm {a} }(t)]}}\right)+m_{2}\pi \end{محاذاة}}}

يمثل كل من 2m1π و m2π مضاعفات صحيحة للعدد π اللازمة لفك الطور. عند قيم الزمن t ، حيث لا يوجد تغيير في العدد الصحيح m2 ، تكون مشتقة φ ( t ) هي

ω(ت)=دφ(ت)دت=ددتدالة الظل العكسي(أنام[sأ(ت)]Rهـ[sأ(ت)])=11+(أنام[sأ(ت)]Rهـ[sأ(ت)])2ددت(أنام[sأ(ت)]Rهـ[sأ(ت)])=Rهـ[sأ(ت)]دأنام[sأ(ت)]دت-أنام[sأ(ت)]دRهـ[sأ(ت)]دت(Rهـ[sأ(ت)])2+(أنام[sأ(ت)])2=1|sأ(ت)|2(Rهـ[sأ(ت)]دأنام[sأ(ت)]دت-أنام[sأ(ت)]دRهـ[sأ(ت)]دت)=1(s(ت))2+(s^(ت))2(s(ت)دs^(ت)دت-s^(ت)دs(ت)دت){\displaystyle {\begin{aligned}\omega (t)={\frac {d\varphi (t)}{dt}}&={\frac {d}{dt}}\arctan \left({\frac {{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{{\mathcal {Re}}[s_{\mathrm {a} }(t)]}}\right)\\[3pt]&={\frac {1}{1+\left({\frac {{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{{\mathcal {Re}}[s_{\mathrm {a} }(t)]}}\right)^{2}}}{\frac {d}{dt}}\left({\frac {{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{{\mathcal {Re}}[s_{\mathrm {a} }(t)]}}\right)\\[3pt]&={\frac {{\mathcal {Re}}[s_{\mathrm {a} }(t)]{\frac {d{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{dt}}-{\mathcal {Im}}[s_{\mathrm {a} }(t)]{\frac {d{\mathcal {Re}}[s_{\mathrm {a} }(t)]}{dt}}}{({\mathcal {Re}}[s_{\mathrm {a} }(t)])^{2}+({\mathcal {Im}}[s_{\mathrm {a} }(t)])^{2}}}\\[3pt]&={\frac {1}{|s_{\mathrm {a} }(t)|^{2}}}\left({\mathcal {Re}}[s_{\mathrm {a} }(t)]{\frac {d{\mathcal {Im}}[s_{\mathrm {a} }(t)]}{dt}}-{\mathcal {Im}}[s_{\mathrm {a} }(t)]{\frac {d{\mathcal {Re}}[s_{\mathrm {a} }(t)]}{dt}}\right)\\[3pt]&={\frac {1}{(s(t))^{2}+\left({\hat {s}}(t)\right)^{2}}}\left(s(t){\frac {d{\hat {s}}(t)}{dt}}-{\hat {s}}(t){\frac {ds(t)}{dt}}\right)\end{aligned}}}

بالنسبة للدوال ذات الزمن المتقطع، يمكن كتابة ذلك على شكل استدعاء ذاتي:

φ[ن]=φ[ن-1]+ω[ن]=φ[ن-1]+arg{sأ[ن]}-arg{sأ[ن-1]}Δφ[ن]=φ[ن-1]+arg{sأ[ن]sأ[ن-1]}{\displaystyle {\begin{aligned}\varphi [n]&=\varphi [n-1]+\omega [n]\\&=\varphi [n-1]+\underbrace {\arg\{s_{\mathrm {a} }[n]\}-\arg\{s_{\mathrm {a} }[n-1]\}} _{\Delta \varphi [n]}\\&=\varphi [n-1]+\arg \left\{{\frac {s_{\mathrm {a} }[n]}{s_{\mathrm {a} }[n-1]}}\right\}\\\end{aligned}}}

يمكن إزالة الانقطاعات بإضافة 2π عندما يكون Δφ [ n ] ≤ −π ، وطرح 2π عندما يكون Δφ [ n ] >  π . يسمح ذلك لـ φ [ n ] بالتراكم بلا حدود، مما ينتج عنه طور لحظي غير ملتف. الصيغة المكافئة التي تستبدل عملية باقي القسمة على بعملية ضرب عددية مركبة هي: 

φ[ن]=φ[ن-1]+arg{sأ[ن]sأ*[ن-1]}،{\displaystyle \varphi [n]=\varphi [n-1]+\arg\{s_{\mathrm {a} }[n]\,s_{\mathrm {a} }^{*}[n-1]\},}

حيث تشير العلامة النجمية إلى المرافق المركب . التردد اللحظي المتقطع (بوحدات راديان لكل عينة) هو ببساطة تقدم الطور لتلك العينة.

ω[ن]=arg{sأ[ن]sأ*[ن-1]}.{\displaystyle \omega [n]=\arg\{s_{\mathrm {a} }[n]\,s_{\mathrm {a} }^{*}[n-1]\}.}

تمثيل معقد

في بعض التطبيقات، مثل حساب متوسط ​​قيم الطور في عدة لحظات زمنية، قد يكون من المفيد تحويل كل قيمة إلى عدد مركب، أو تمثيل متجهي: [ 4 ]

هـأناφ(ت)=sأ(ت)|sأ(ت)|=كوس(φ(ت))+أناالخطيئة(φ(ت)).{\displaystyle e^{i\varphi (t)}={\frac {s_{\mathrm {a} }(t)}{|s_{\mathrm {a} }(t)|}}=\cos(\varphi (t))+i\sin(\varphi (t)).}

يشبه هذا التمثيل تمثيل الطور الملتف من حيث أنه لا يميز بين مضاعفات 2π في الطور، ولكنه يشبه تمثيل الطور غير الملتف لكونه متصلاً. ويمكن الحصول على متوسط ​​الطور المتجهي كقيمة وسيطة لمجموع الأعداد المركبة دون مراعاة الالتفاف.

انظر أيضاً

مراجع

  1. ناهين، بول (2024). الراديو الرياضي: داخل سحر AM وFM والنطاق الجانبي الأحادي . برينستون: مطبعة جامعة برينستون. ص 210-213 . ISBN  9780691235318.
  2. سيديتش، إي.؛ ديوروفيتش، آي.؛ ستانكوفيتش، إل. (أغسطس 2008). "تحليل الأداء الكمي لـ Scalogram كمُقدِّر للتردد اللحظي". معاملات IEEE في معالجة الإشارات . 56 (8): 3837-3845 . Bibcode : 2008ITSP...56.3837S . doi : 10.1109/TSP.2008.924856 . ISSN 1053-587X . S2CID 16396084 .  
  3. بلاكليدج، جوناثان م. (2006). معالجة الإشارات الرقمية: الأساليب الرياضية والحسابية، وتطوير البرمجيات والتطبيقات ( الطبعة الثانية). دار وودهيد للنشر. ص 134. ISBN   1904275265.
  4. وانغ، س. (2014). "طريقة محسّنة لفكّ الطور الموجّه بالجودة وتطبيقاتها في التصوير بالرنين المغناطيسي" . التقدم في بحوث الكهرومغناطيسية . 145 : 273-286 . doi : 10.2528/PIER14021005 .

للمزيد من القراءة

  • كوهين، ليون (1995). تحليل التردد الزمني . برنتيس هول.
  • غرانلوند؛ كنوتسون (1995). معالجة الإشارات لرؤية الحاسوب . دار نشر كلوير الأكاديمية.