الكثافة الطيفية

تُظهر الكثافة الطيفية للضوء الفلوري كدالة للطول الموجي البصري قممًا عند الانتقالات الذرية، كما هو موضح بالأسهم المرقمة.
يُظهر شكل الموجة الصوتية بمرور الوقت (يسار) طيفًا واسعًا لقوة الصوت (يمين).

في معالجة الإشارات ، طيف القدرةSxx(و){\displaystyle S_{xx}(f)}إشارة زمنية مستمرةx(ت){\displaystyle x(t)}يصف توزيع الطاقة إلى مكونات الترددو{\displaystyle f}تكوين تلك الإشارة. [ 1 ] يُظهر تحليل فورييه أنه يمكن تحليل أي إشارة فيزيائية إلى توزيع ترددات على نطاق متصل، حيث قد تتركز بعض الطاقة عند ترددات منفصلة. يُطلق على المتوسط ​​الإحصائي لطاقة أو قدرة أي نوع من الإشارات (بما في ذلك الضوضاء ) عند تحليلها من حيث محتواها الترددي، اسم الكثافة الطيفية .

عندما تتركز طاقة الإشارة حول فترة زمنية محدودة، خاصةً إذا كانت طاقتها الكلية محدودة، يمكن حساب كثافة الطاقة الطيفية . والأكثر شيوعًا هو استخدام كثافة القدرة الطيفية (PSD، أو ببساطة طيف القدرة )، والتي تنطبق على الإشارات الموجودة على مدار الزمن ، أو على مدار فترة زمنية طويلة بما يكفي (خاصةً بالنسبة لمدة القياس) بحيث يمكن اعتبارها فترة زمنية غير محدودة. تشير كثافة القدرة الطيفية حينها إلى توزيع القدرة الطيفية الذي يمكن إيجاده، نظرًا لأن الطاقة الكلية لمثل هذه الإشارة على مدار الزمن ستكون عادةً غير محدودة. ينتج عن جمع أو تكامل المكونات الطيفية القدرة الكلية (في عملية فيزيائية) أو التباين (في عملية إحصائية)، وهو ما يُطابق ما يمكن الحصول عليه من خلال تكاملx2(ت){\displaystyle x^{2}(t)}على مدى المجال الزمني، كما تمليه نظرية بارسيفال . [ 1 ]

طيف العملية الفيزيائيةx(ت){\displaystyle x(t)}غالباً ما تحتوي على معلومات أساسية حول طبيعةx{\displaystyle x}على سبيل المثال، يمكن تحديد درجة الصوت ونبرته لآلة موسيقية من خلال التحليل الطيفي. ويتحدد لون مصدر الضوء بواسطة طيف المجال الكهربائي للموجة الكهرومغناطيسية .هـ(ت){\displaystyle E(t)}نظرًا لتذبذبها بتردد عالٍ للغاية، يتطلب الحصول على طيف من بيانات السلاسل الزمنية كهذه استخدام تحويل فورييه ، وتعميمات مبنية على تحليل فورييه. في كثير من الحالات، لا يُسجّل المجال الزمني مباشرةً في التطبيق العملي، كما هو الحال عند استخدام موشور تشتيت للحصول على طيف الضوء في مطياف ، أو عند إدراك الصوت من خلال تأثيره على مستقبلات السمع في الأذن الداخلية، حيث يكون كل مستقبل حساسًا لتردد معين.

مع ذلك، تركز هذه المقالة على الحالات التي تكون فيها السلسلة الزمنية معروفة (على الأقل إحصائيًا) أو تُقاس مباشرةً (مثلًا بواسطة ميكروفون يُسجل بواسطة حاسوب). يُعد طيف القدرة مهمًا في معالجة الإشارات الإحصائية وفي الدراسة الإحصائية للعمليات العشوائية ، فضلًا عن العديد من فروع الفيزياء والهندسة الأخرى . عادةً ما تكون العملية دالةً للزمن، ولكن يمكن مناقشة البيانات في المجال المكاني بالمثل عند تحليلها بدلالة التردد المكاني . [ 1 ]

الوحدات

في الفيزياء ، قد تكون الإشارة موجة، مثل موجة كهرومغناطيسية أو موجة صوتية أو اهتزاز آلية. يصف طيف كثافة القدرة (PSD) للإشارة كثافة قدرة الإشارة كدالة للتردد. ويُعبر عن طيف كثافة القدرة عادةً بوحدة واط لكل هرتز (W/Hz) في النظام الدولي للوحدات . [ 2 ]

عندما تُعرَّف الإشارة بدلالة جهد متغير مع الزمن فقط، على سبيل المثال، لا توجد قدرة محددة مرتبطة بهذا الجهد. في هذه الحالة، تُحسب "القدرة" ببساطة بدلالة مربع الإشارة، حيث تتناسب هذه القيمة دائمًا مع القدرة الفعلية التي تُوصلها الإشارة إلى مقاومة معينة . لذا، يمكن استخدام وحدة V² ⋅Hz⁻¹ لكثافة القدرة الطيفية ( PSD). أما كثافة الطاقة الطيفية (ESD) فتُقاس بوحدة V² ⋅s⋅Hz⁻¹ ، لأن الطاقة هي القدرة مضروبة في الزمن (مثلًا، واط-ساعة ). [ 3 ]

في الحالة العامة، تكون وحدة كثافة القدرة الطيفية (PSD) هي نسبة وحدة التباين لكل وحدة تردد؛ فعلى سبيل المثال، سلسلة من قيم الإزاحة (بالمتر) مع مرور الوقت (بالثواني) يكون لها كثافة قدرة طيفية بوحدة م² / هرتز. في تحليل الاهتزازات العشوائية ، يمكن استخدام وحدة g₀² هرتز⁻¹ لكثافة القدرة الطيفية للتسارع ، حيث g₀ ترمز إلى الجاذبية الأرضية القياسية . [ 4 ]

رياضياً، ليس من الضروري تحديد أبعاد فيزيائية للإشارة أو للمتغير المستقل. في المناقشة التالية، سيبقى معنى x ( t ) غير محدد، ولكن سيُفترض أن المتغير المستقل هو الزمن.

من جانب واحد مقابل من جانبين

يمكن أن تكون كثافة القدرة الطيفية (PSD) دالة أحادية الجانب للترددات الموجبة فقط، أو دالة ثنائية الجانب للترددات الموجبة والسالبة ولكن بنصف السعة فقط. وتكون كثافة القدرة الطيفية للضوضاء أحادية الجانب عمومًا في الهندسة وثنائية الجانب في الفيزياء. [ 5 ]

تعريف

كثافة الطيف الطاقي

في معالجة الإشارات ، طاقة الإشارةx(ت){\displaystyle x(t)}يُعطى بواسطة هـ-|x(ت)|2 دت.{\displaystyle E\triangleq \int _{-\infty }^{\infty }\left|x(t)\right|^{2}\ dt.} بافتراض أن الطاقة الكلية محدودة (أيx(ت){\displaystyle x(t)}(دالة قابلة للتكامل التربيعي ) تسمح بتطبيق نظرية بارسيفال (أو نظرية بلانشيريل ). [ 6 ] أي، -|x(ت)|2دت=-|x^(و)|2دو،{\displaystyle \int _{-\infty }^{\infty }|x(t)|^{2}\,dt=\int _{-\infty }^{\infty }\left|{\hat {x}}(f)\right|^{2}\,df,} أين x^(و)=-هـ-أنا2πوتx(ت) دت،{\displaystyle {\hat {x}}(f)=\int _{-\infty }^{\infty }e^{-i2\pi ft}x(t)\ dt,} هو تحويل فورييه لـx(ت){\displaystyle x(t)}بترددو{\displaystyle f}( بالهرتز ). [ 7 ] تنطبق النظرية أيضًا على حالات الزمن المتقطع. بما أن التكامل على الجانب الأيسر هو طاقة الإشارة، فإن قيمة|x^(و)|2دو{\displaystyle \left|{\hat {x}}(f)\right|^{2}df}يمكن تفسيرها على أنها دالة كثافة مضروبة في فاصل تردد صغير للغاية، يصف الطاقة الموجودة في الإشارة عند الترددو{\displaystyle f}في نطاق الترددو+دو{\displaystyle f+df}.

وبالتالي، فإن كثافة الطاقة الطيفية لـx(ت){\displaystyle x(t)}يُعرَّف على أنه [ 8 ]

الوظيفةS¯xx(و){\displaystyle {\bar {S}}_{xx}(f)}والارتباط الذاتي لـx(ت){\displaystyle x(t)}تشكل زوج تحويل فورييه، وهي نتيجة تُعرف أيضًا باسم نظرية وينر-خينشين (انظر أيضًا مخطط الفترة ).

كمثال مادي لكيفية قياس كثافة الطاقة الطيفية لإشارة ما، لنفترضV(ت){\displaystyle V(t)}يمثل الجهد ( بالفولت ) لنبضة كهربائية تنتشر على طول خط نقل ذي مقاومةZ{\displaystyle Z}ولنفترض أن الخط ينتهي بمقاومة مطابقة (بحيث يتم توصيل كل طاقة النبضة إلى المقاومة ولا ينعكس أي منها). وفقًا لقانون أوم ، فإن القدرة المُسلَّمة إلى المقاومة في الزمن tت{\displaystyle t}يساويV(ت)2/Z{\displaystyle V(t)^{2}/Z}لذا، يتم إيجاد الطاقة الكلية عن طريق التكامل.V(ت)2/Z{\displaystyle V(t)^{2}/Z}بالنسبة للزمن خلال مدة النبضة. لإيجاد قيمة كثافة الطاقة الطيفيةS¯xx(و){\displaystyle {\bar {S}}_{xx}(f)}بترددو{\displaystyle f}يمكن إدخال مرشح تمرير نطاقي بين خط النقل والمقاوم، بحيث يمرر هذا المرشح نطاقًا ضيقًا من الترددات فقط.Δو{\displaystyle \Delta f}على سبيل المثال، بالقرب من التردد المطلوب، ثم قم بقياس الطاقة الكليةهـ(و){\displaystyle E(f)}تبددت عبر المقاوم. قيمة كثافة الطاقة الطيفية عندو{\displaystyle f}ثم يُقدّر أنهـ(و)/Δو{\displaystyle E(f)/\Delta f}في هذا المثال، بما أن القوةV(ت)2/Z{\displaystyle V(t)^{2}/Z}لها الوحدة V 2 ⋅Ω −1 ، الطاقةهـ(و){\displaystyle E(f)}لها الوحدة V 2 ⋅s⋅Ω −1  = J ، ومن ثم التقديرهـ(و)/Δو{\displaystyle E(f)/\Delta f}تبلغ وحدة كثافة الطاقة الطيفية J⋅Hz −1 . في كثير من الحالات، من الشائع حذف خطوة القسمة علىZ{\displaystyle Z}وبالتالي فإن كثافة الطيف الطاقي لها الوحدة V 2 ⋅s·Hz −1 .

يمكن تعميم هذا التعريف بطريقة مباشرة ليشمل إشارة منفصلة ذات عدد لا نهائي قابل للعد من القيمxن{\displaystyle x_{n}}مثل إشارة يتم أخذ عينات منها في أوقات منفصلةتن=ت0+(نΔت){\displaystyle t_{n}=t_{0}+(n\,\Delta t)}: S¯xx(و)=ليمشمال(Δت)2|ن=-شمالشمالxنهـ-أنا2πونΔت|2|x^د(و)|2،{\displaystyle {\bar {S}}_{xx}(f)=\lim _{N\to \infty }(\Delta t)^{2}\underbrace {\left|\sum _{n=-N}^{N}x_{n}e^{-i2\pi fn\,\Delta t}\right|^{2}} _{\left|{\hat {x}}_{d}(f)\right|^{2}},} أينx^د(و){\displaystyle {\hat {x}}_{d}(f)}هو تحويل فورييه المتقطع لـxن.{\displaystyle x_{n}.} فترة أخذ العيناتΔت{\displaystyle \Delta t}يلزم ذلك للحفاظ على الوحدة الفيزيائية الصحيحة ولضمان استعادة الحالة المستمرة في الحدΔت0{\displaystyle \Delta t\to 0}لكن في العلوم الرياضية، غالباً ما تُحدد الفترة بـ 1، مما يُبسط النتائج على حساب عموميتها. (انظر أيضاً: التردد المُعَيَّر (الوحدة) )

كثافة القدرة الطيفية

طيف القدرة لتباين درجة حرارة إشعاع الخلفية الكونية الميكروي المقاس بدلالة المقياس الزاوي. الخط المتصل هو نموذج نظري للمقارنة.

يُعدّ التعريف المذكور أعلاه لكثافة الطيف الطاقي مناسبًا للإشارات العابرة (الإشارات النبضية) التي تتركز طاقتها حول نافذة زمنية واحدة؛ وفي هذه الحالة، تكون تحويلات فورييه للإشارات موجودة عمومًا. أما بالنسبة للإشارات المستمرة على مدار الزمن، فيجب تعريف كثافة الطيف القدرة (PSD) الموجودة للعمليات المستقرة ؛ وهذا يصف كيفية توزيع قدرة الإشارة أو السلسلة الزمنية على التردد، كما في المثال البسيط المذكور سابقًا. هنا، يمكن أن تكون القدرة هي القدرة الفيزيائية الفعلية، أو غالبًا، لتسهيل التعامل مع الإشارات المجردة، تُعرّف ببساطة على أنها مربع قيمة الإشارة. على سبيل المثال، يدرس الإحصائيون تباين دالة ما مع مرور الوقت.x(ت){\displaystyle x(t)}(أو على متغير مستقل آخر)، وباستخدام تشبيه بالإشارات الكهربائية (من بين عمليات فيزيائية أخرى)، من المعتاد الإشارة إليه باسم طيف القدرة حتى عندما لا تكون هناك قدرة فيزيائية متضمنة. إذا قام المرء بإنشاء مصدر جهد فيزيائي يتبعx(ت){\displaystyle x(t)}وإذا طبقنا ذلك على طرفي مقاومة قيمتها أوم واحد ، فإن القدرة اللحظية المبددة في تلك المقاومة ستُعطى بالعلاقة التالية:x2(ت){\displaystyle x^{2}(t)}واط .

متوسط ​​القدرةP{\displaystyle P}إشارةx(ت){\displaystyle x(t)}وبالتالي، يُعطى المتوسط ​​الزمني الإجمالي بالمعادلة التالية، حيث تمثل الفترةتي{\displaystyle T}يتمحور حول وقت عشوائيت=ت0{\displaystyle t=t_{0}}: P=ليمتي1تيت0-تي/2ت0+تي/2|x(ت)|2دت{\displaystyle P=\lim _{T\to \infty }{\frac {1}{T}}\int _{t_{0}-T/2}^{t_{0}+T/2}\left|x(t)\right|^{2}\,dt}

عندما يكون من الأنسب التعامل مع الحدود الزمنية في الإشارة نفسها بدلاً من الحدود الزمنية في حدود التكامل، يمكن كتابة متوسط ​​القدرة أيضًا على النحو التالي: P=ليمتي1تي-|xتي(ت)|2دت،{\displaystyle P=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }\left|x_{T}(t)\right|^{2}\,dt,} أينxتي(ت)=x(ت)wتي(ت){\displaystyle x_{T}(t)=x(t)w_{T}(t)}وwتي(ت){\displaystyle w_{T}(t)}الوحدة داخل الفترة التعسفية وصفر في أي مكان آخر.

متىP{\displaystyle P}إذا كانت قيمة التكامل غير صفرية، فيجب أن ينمو التكامل إلى ما لا نهاية على الأقل بنفس سرعةتي{\displaystyle T}نعم. هذا هو السبب في أننا لا نستطيع استخدام طاقة الإشارة، وهي ذلك التكامل المتباعد.

في تحليل محتوى التردد للإشارةx(ت){\displaystyle x(t)}قد يرغب المرء في حساب تحويل فورييه العاديx^(و){\displaystyle {\hat {x}}(f)}مع ذلك، بالنسبة للعديد من الإشارات المهمة، لا يوجد تحويل فورييه عادي بشكل رسمي. [ ملاحظة 1 ] ومع ذلك، في ظل شروط مناسبة، لا تزال بعض تعميمات تحويل فورييه (مثل تحويل فورييه-ستيلتيس ) متوافقة مع نظرية بارسيفال . على هذا النحو، P=ليمتي1تي-|x^تي(و)|2دو،{\displaystyle P=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }|{\hat {x}}_{T}(f)|^{2}\,df,} حيث يحدد التكامل كثافة الطيف القدرة : [ 9 ] [ 10 ]

ثم تسمح نظرية الالتفاف بالنظر إلى|x^تي(و)|2{\displaystyle |{\hat {x}}_{T}(f)|^{2}}باعتبارها تحويل فورييه للالتفاف الزمني لـxتي*(-ت){\displaystyle x_{T}^{*}(-t)}وxتي(ت){\displaystyle x_{T}(t)}، حيث يمثل * المرافق المركب.

لإثبات الادعاء الوارد في المعادلة 2 أدناه، سنجد تعبيرًا لـ[x^تي(و)]*{\displaystyle [{\hat {x}}_{T}(f)]^{*}}سيكون ذلك مفيدًا لهذا الغرض. في الواقع، سنوضح ذلك.[x^تي(و)]*=F{xتي*(-ت)}{\displaystyle [{\hat {x}}_{T}(f)]^{*}={\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}}ابدأ بملاحظة أن F{xتي*(-ت)}=-xتي*(-ت)هـ-أنا2πوتدت{\displaystyle {\begin{aligned}{\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}&=\int _{-\infty }^{\infty }x_{T}^{*}(-t)e^{-i2\pi ft}dt\end{aligned}}} ودعz=-ت{\displaystyle z=-t}، لهذا السببz-{\displaystyle z\rightarrow -\infty }متىت{\displaystyle t\rightarrow \infty }والعكس صحيح. -xتي*(-ت)هـ-أنا2πوتدت=-xتي*(z)هـأنا2πوz(-دz)=-xتي*(z)هـأنا2πوzدz=-xتي*(ت)هـأنا2πوتدت{\displaystyle {\begin{aligned}\int _{-\infty }^{\infty }x_{T}^{*}(-t)e^{-i2\pi ft}dt&=\int _{\infty }^{-\infty }x_{T}^{*}(z)e^{i2\pi fz}\left(-dz\right)\\&=\int _{-\infty }^{\infty }x_{T}^{*}(z)e^{i2\pi fz}dz\\&=\int _{-\infty }^{\infty }x_{T}^{*}(t)e^{i2\pi ft}dt\end{aligned}}} حيث تم استخدام في السطر الأخيرz{\displaystyle z}وت{\displaystyle t}كونها متغيرات وهمية. لذلك، لدينا F{xتي*(-ت)}=-xتي*(-ت)هـ-أنا2πوتدت=-xتي*(ت)هـأنا2πوتدت=-xتي*(ت)[هـ-أنا2πوت]*دت=[-xتي(ت)هـ-أنا2πوتدت]*=[F{xتي(ت)}]*=[x^تي(و)]*{\displaystyle {\begin{aligned}{\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}&=\int _{-\infty }^{\infty }x_{T}^{*}(-t)e^{-i2\pi ft}dt\\&=\int _{-\infty }^{\infty }x_{T}^{*}(t)e^{i2\pi ft}dt\\&=\int _{-\infty }^{\infty }x_{T}^{*}(t)[e^{-i2\pi ft}]^{*}dt\\&=\left[\int _{-\infty }^{\infty }x_{T}(t)e^{-i2\pi ft}dt\right]^{*}\\&=\left[{\mathcal {F}}\left\{x_{T}(t)\right\}\right]^{*}\\&=\left[{\hat {x}}_{T}(f)\right]^{*}\end{aligned}}} qed

والآن، دعونا نبرهن على صحة الادعاء الوارد في المعادلة 2 باستخدام المتطابقة الموضحة. بالإضافة إلى ذلك، سنقوم بالاستبدال.u(ت)=xتي*(-ت){\displaystyle u(t)=x_{T}^{*}(-t)}وبهذه الطريقة، لدينا: |x^تي(و)|2=[x^تي(و)]*x^تي(و)=F{xتي*(-ت)}F{xتي(ت)}=F{u(ت)}F{xتي(ت)}=F{u(ت)*xتي(ت)}=-[-u(τ-ت)xتي(ت)دت]هـ-أنا2πوτدτ=-[-xتي*(ت-τ)xتي(ت)دت]هـ-أنا2πوτ دτ،{\displaystyle {\begin{aligned}\left|{\hat {x}}_{T}(f)\right|^{2}&=[{\hat {x}}_{T}(f)]^{*}\cdot {\hat {x}}_{T}(f)\\&={\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}\cdot {\mathcal {F}}\left\{x_{T}(t)\right\}\\&={\mathcal {F}}\left\{u(t)\right\}\cdot {\mathcal {F}}\left\{x_{T}(t)\right\}\\&={\mathcal {F}}\left\{u(t)\mathbin {\mathbf {*} } x_{T}(t)\right\}\\&=\int _{-\infty }^{\infty }\left[\int _{-\infty }^{\infty }u(\tau -t)x_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau \\&=\int _{-\infty }^{\infty }\left[\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }\ d\tau ,\end{aligned}}} حيث تم استخدام نظرية الالتفاف عند الانتقال من السطر الثالث إلى السطر الرابع.

الآن، إذا قسمنا عملية الالتفاف الزمني المذكورة أعلاه على الدورةتي{\displaystyle T}واعتبر الحد كما يليتي{\displaystyle T\rightarrow \infty }، فتصبح دالة الارتباط الذاتي للإشارة غير المُقاسة.x(ت){\displaystyle x(t)}، والذي يُشار إليه بـRxx(τ){\displaystyle R_{xx}(\tau )}بشرط أنx(ت){\displaystyle x(t)}هي إرجودية ، وهذا صحيح في معظم الحالات العملية، ولكن ليس في جميعها. [ ملاحظة 2 ]ليمتي1تي|x^تي(و)|2=-[ليمتي1تي-xتي*(ت-τ)xتي(ت)دت]هـ-أنا2πوτ دτ=-Rxx(τ)هـ-أنا2πوτدτ{\displaystyle \lim _{T\to \infty }{\frac {1}{T}}\left|{\hat {x}}_{T}(f)\right|^{2}=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }\ d\tau =\int _{-\infty }^{\infty }R_{xx}(\tau )e^{-i2\pi f\tau }d\tau }

بافتراض خاصية الإرجودية لـx(ت){\displaystyle x(t)}ويمكن إيجاد كثافة القدرة الطيفية مرة أخرى كتحويل فورييه لدالة الارتباط الذاتيRxx{\displaystyle R_{xx}}، وهي خاصية تُعرف باسم نظرية وينر-خينشين . [ 11 ]

يستخدم العديد من المؤلفين هذه العلاقة لتعريف كثافة القدرة الطيفية بدلالة دالة الارتباط الذاتي بدلاً من تحويل فورييه للإشارة كما فعلنا. [ 12 ]

قوة الإشارة في نطاق تردد معين[و1،و2]{\displaystyle [f_{1},f_{2}]}، أين0<و1<و2{\displaystyle 0<f_{1}<f_{2}}يمكن حسابها عن طريق التكامل على التردد. بما أنSxx(-و)=Sxx(و){\displaystyle S_{xx}(-f)=S_{xx}(f)}، ويمكن إسناد قدر متساوٍ من الطاقة إلى نطاقات التردد الموجبة والسالبة، وهو ما يفسر العامل 2 بالشكل التالي (تعتمد هذه العوامل البسيطة على الاصطلاحات المستخدمة): Pمحدود النطاق=2و1و2Sxx(و)دو{\displaystyle P_{\textsf {band-limited}}=2\int _{f_{1}}^{f_{2}}S_{xx}(f)\,df} وبشكل أعم، يمكن استخدام تقنيات مماثلة لتقدير الكثافة الطيفية المتغيرة مع الزمن. في هذه الحالة، الفترة الزمنيةتي{\displaystyle T}تكون محدودة وليست قريبة من اللانهاية. ينتج عن ذلك انخفاض في التغطية الطيفية والدقة، حيث أن الترددات الأقل من1/تي{\displaystyle 1/T}لا يتم أخذ عينات منها، والنتائج عند ترددات ليست مضاعفات صحيحة لـ1/تي{\displaystyle 1/T}ليست مستقلة. باستخدام سلسلة زمنية واحدة فقط، سيكون طيف القدرة المُقدَّر "مشوشًا" للغاية؛ ومع ذلك، يمكن التخفيف من ذلك إذا كان من الممكن تقييم القيمة المتوقعة (في المعادلة أعلاه) باستخدام عدد كبير (أو لا نهائي) من الأطياف قصيرة المدى التي تتوافق مع مجموعات إحصائية من تحقيقاتx(ت){\displaystyle x(t)}تم تقييمها خلال الفترة الزمنية المحددة.

كما هو الحال مع كثافة الطيف الطاقي، يمكن تعميم تعريف كثافة الطيف القدرة ليشمل متغيرات الزمن المنفصلةxن{\displaystyle x_{n}}كما في السابق، يمكننا اعتبار نافذة من-شمالنشمال{\displaystyle -N\leq n\leq N}مع أخذ عينات من الإشارة في أوقات منفصلةتن=ت0+(نΔت){\displaystyle t_{n}=t_{0}+(n\,\Delta t)}لفترة قياس إجماليةتي=(2شمال+1)Δت{\displaystyle T=(2N+1)\,\Delta t}. Sxx(و)=ليمشمال(Δت)2تي|ن=-شمالشمالxنهـ-أنا2πونΔت|2{\displaystyle S_{xx}(f)=\lim _{N\to \infty }{\frac {(\Delta t)^{2}}{T}}\left|\sum _{n=-N}^{N}x_{n}e^{-i2\pi fn\,\Delta t}\right|^{2}} لاحظ أنه يمكن الحصول على تقدير واحد لكثافة القدرة الطيفية من خلال عدد محدود من عمليات أخذ العينات. وكما كان الحال سابقًا، يتم الحصول على كثافة القدرة الطيفية الفعلية عندماشمال{\displaystyle N}(وبالتاليتي{\displaystyle T}عندما يقترب عدد القياسات من اللانهاية، يتم تطبيق القيمة المتوقعة رسميًا. في التطبيقات العملية، يُجرى عادةً حساب متوسط ​​كثافة القدرة الطيفية (PSD) لعدد محدود من القياسات عبر العديد من التجارب للحصول على تقدير أدق لكثافة القدرة الطيفية النظرية للعملية الفيزيائية الكامنة وراء القياسات الفردية. تُسمى كثافة القدرة الطيفية المحسوبة هذه أحيانًا مخطط الفترة . يتقارب مخطط الفترة هذا مع كثافة القدرة الطيفية الحقيقية مع ازدياد عدد التقديرات وفترة حساب المتوسط.تي{\displaystyle T}يقترب من اللانهاية. [ 13 ]

إذا كان لكل من الإشارتين كثافة طيفية للقدرة، فيمكن حساب الكثافة الطيفية المتقاطعة بشكل مماثل؛ وبما أن الكثافة الطيفية للقدرة مرتبطة بالارتباط الذاتي، فإن الكثافة الطيفية المتقاطعة مرتبطة بالارتباط المتبادل .

خصائص كثافة القدرة الطيفية

تتضمن بعض خصائص PSD ما يلي: [ 14 ]

  • يكون طيف القدرة دائمًا حقيقيًا وغير سالب، كما أن طيف العملية ذات القيم الحقيقية هو أيضًا دالة زوجية للتردد:Sxx(-و)=Sxx(و){\displaystyle S_{xx}(-f)=S_{xx}(f)}.
  • بالنسبة لعملية عشوائية مستمرة x(t)، يمكن إعادة بناء دالة الارتباط الذاتي R xx ( t ) من طيف قدرتها S xx (f) باستخدام تحويل فورييه العكسي
  • باستخدام نظرية بارسيفال ، يمكن حساب العزم الثاني (متوسط ​​القدرة) لعملية ما عن طريق تكامل طيف القدرة على جميع الترددات: P=هـ(x2)=-Sxx(و)دو{\displaystyle P=\operatorname {E} (x^{2})=\int _{-\infty }^{\infty }\!S_{xx}(f)\,df}
  • بالنسبة لعملية حقيقية x ( t ) ذات كثافة طيفية للقدرة Sxx(و){\displaystyle S_{xx}(f)}يمكن للمرء حساب الطيف المتكامل أو توزيع القدرة الطيفيةF(و){\displaystyle F(f)}، والتي تحدد متوسط ​​القدرة المحدودة النطاق الموجودة في الترددات من التيار المستمر إلى f باستخدام: [ 15 ]F(و)=20وSxx(و)دو.{\displaystyle F(f)=2\int _{0}^{f}S_{xx}(f')\,df'.} لاحظ أن التعبير السابق للطاقة الكلية (تباين الإشارة) هو حالة خاصة حيث f → ∞ . 

كثافة القدرة الطيفية المتقاطعة

بافتراض إشارتينx(ت){\displaystyle x(t)}وy(ت){\displaystyle y(t)}، وكل منها يمتلك كثافة طيفية للقدرةSxx(و){\displaystyle S_{xx}(f)}وSyy(و){\displaystyle S_{yy}(f)}من الممكن تعريف كثافة القدرة الطيفية المتقاطعة ( CPSD ) أو الكثافة الطيفية المتقاطعة ( CSD ). ولنبدأ بالنظر إلى متوسط ​​قدرة هذه الإشارة المركبة. P=ليمتي1تي-[xتي(ت)+yتي(ت)]*[xتي(ت)+yتي(ت)]دت=ليمتي1تي-|xتي(ت)|2+xتي*(ت)yتي(ت)+yتي*(ت)xتي(ت)+|yتي(ت)|2دت{\displaystyle {\begin{aligned}P&=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }\left[x_{T}(t)+y_{T}(t)\right]^{*}\left[x_{T}(t)+y_{T}(t)\right]dt\\&=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }|x_{T}(t)|^{2}+x_{T}^{*}(t)y_{T}(t)+y_{T}^{*}(t)x_{T}(t)+|y_{T}(t)|^{2}dt\\\end{aligned}}}

باستخدام نفس الرموز والأساليب المستخدمة في اشتقاق كثافة القدرة الطيفية، نستغل نظرية بارسيفال ونحصل على Sxy(و)=ليمتي1تي[x^تي*(و)y^تي(و)]Syx(و)=ليمتي1تي[y^تي*(و)x^تي(و)]{\displaystyle {\begin{aligned}S_{xy}(f)&=\lim _{T\to \infty }{\frac {1}{T}}\left[{\hat {x}}_{T}^{*}(f){\hat {y}}_{T}(f)\right]&S_{yx}(f)&=\lim _{T\to \infty }{\frac {1}{T}}\left[{\hat {y}}_{T}^{*}(f){\hat {x}}_{T}(f)\right]\end{aligned}}} حيث، مرة أخرى، مساهماتSxx(و){\displaystyle S_{xx}(f)}وSyy(و){\displaystyle S_{yy}(f)}هذا مفهوم بالفعل. لاحظ أنSxy*(و)=Syx(و){\displaystyle S_{xy}^{*}(f)=S_{yx}(f)}لذا، فإن المساهمة الكاملة في القوة المتقاطعة، بشكل عام، تأتي من ضعف الجزء الحقيقي لأي من توزيعات كثافة القدرة المتقاطعة الفردية . وكما في السابق، نعيد صياغة هذه النواتج من هنا على أنها تحويل فورييه لعملية التفاف زمني، والتي عند قسمتها على الدورة وأخذها إلى النهايةتي{\displaystyle T\to \infty }يصبح تحويل فورييه لدالة الارتباط المتبادل . [ 16 ]Sxy(و)=-[ليمتي1تي-xتي*(ت-τ)yتي(ت)دت]هـ-أنا2πوτدτ=-Rxy(τ)هـ-أنا2πوτدτSyx(و)=-[ليمتي1تي-yتي*(ت-τ)xتي(ت)دت]هـ-أنا2πوτدτ=-Ryx(τ)هـ-أنا2πوτدτ،{\displaystyle {\begin{aligned}S_{xy}(f)&=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )y_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau =\int _{-\infty }^{\infty }R_{xy}(\tau )e^{-i2\pi f\tau }d\tau \\S_{yx}(f)&=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }y_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau =\int _{-\infty }^{\infty }R_{yx}(\tau )e^{-i2\pi f\tau }d\tau ,\end{aligned}}} أينRxy(τ){\displaystyle R_{xy}(\tau )}هو الارتباط المتبادل لـx(ت){\displaystyle x(t)}معy(ت){\displaystyle y(t)}وRyx(τ){\displaystyle R_{yx}(\tau )}هو الارتباط المتبادل لـy(ت){\displaystyle y(t)}معx(ت){\displaystyle x(t)}في ضوء ذلك، يُنظر إلى قسم خدمات الحماية الشخصية على أنه حالة خاصة من قسم خدمات الحماية المدنية.x(ت)=y(ت){\displaystyle x(t)=y(t)}. لوx(ت){\displaystyle x(t)}وy(ت){\displaystyle y(t)}هي إشارات حقيقية (مثل الجهد أو التيار)، وتحويلات فورييه الخاصة بهاx^(و){\displaystyle {\hat {x}}(f)}وy^(و){\displaystyle {\hat {y}}(f)}عادةً ما تقتصر الترددات على الترددات الموجبة وفقًا للعرف. لذلك، في معالجة الإشارات النموذجية، فإنّ كثافة طيف القدرة الكاملة (CPSD) هي مجرد واحدة من كثافات طيف القدرة الكاملة (CPSD) مضروبة في اثنين. CPSDممتلىء=2Sxy(و)=2Syx(و){\displaystyle \operatorname {CPSD} _{\text{Full}}=2S_{xy}(f)=2S_{yx}(f)}

بالنسبة للإشارات المنفصلة x n و y n ، فإن العلاقة بين الكثافة الطيفية المتقاطعة والتباين المتقاطع هي Sxy(و)=ن=-Rxy(τن)هـ-أنا2πوτنΔτ{\displaystyle S_{xy}(f)=\sum _{n=-\infty }^{\infty }R_{xy}(\tau _{n})e^{-i2\pi f\tau _{n}}\,\Delta \tau }

تقدير

يهدف تقدير الكثافة الطيفية إلى تقدير الكثافة الطيفية لإشارة عشوائية من سلسلة من العينات الزمنية. وبحسب المعلومات المتوفرة عن الإشارة، قد تتضمن تقنيات التقدير أساليب بارامترية أو غير بارامترية ، وقد تعتمد على تحليل المجال الزمني أو مجال التردد. على سبيل المثال، تتضمن إحدى التقنيات البارامترية الشائعة مطابقة المشاهدات مع نموذج الانحدار الذاتي . أما التقنية غير البارامترية الشائعة فهي مخطط الفترة .

عادةً ما تُقدَّر الكثافة الطيفية باستخدام طرق تحويل فورييه (مثل طريقة ويلش )، ولكن يمكن أيضًا استخدام تقنيات أخرى مثل طريقة الإنتروبيا القصوى . في عام 2026، تم إثبات التقدير غير المتحيز للكثافات الطيفية وتعميماتها ذات الرتب الأعلى في سياق الأطياف المتعددة استنادًا إلى تحويلات فورييه المُنَظَّمة وإحصائية K. [ 17 ]

  • المركز الطيفي للإشارة هو نقطة المنتصف لدالة الكثافة الطيفية الخاصة بها، أي التردد الذي يقسم التوزيع إلى جزأين متساويين.
  • يمثل تردد الحافة الطيفية (SEF)، والذي يُعبر عنه عادةً بـ "SEF x"، التردد الذي تقع تحته نسبة x بالمئة من إجمالي طاقة إشارة معينة؛ وعادةً ما تتراوح قيمة x بين 75 و95. وهو مقياس شائع الاستخدام في مراقبة تخطيط كهربية الدماغ (EEG) ، حيث يُستخدم SEF لتقدير عمق التخدير ومراحل النوم . [ 18 ] [ 19 ]
  • الغلاف الطيفي هو منحنى غلاف كثافة الطيف. وهو يصف نقطة زمنية واحدة (نافذة واحدة، على وجه الدقة). على سبيل المثال، في الاستشعار عن بعد باستخدام مطياف ، يمثل الغلاف الطيفي لمعلم ما حدود خصائصه الطيفية ، كما هو محدد بنطاق مستويات السطوع في كل نطاق من النطاقات الطيفية ذات الأهمية.
  • الكثافة الطيفية دالة للتردد، وليست دالة للزمن. مع ذلك، يمكن حساب الكثافة الطيفية لنافذة صغيرة من إشارة أطول، ورسمها بيانيًا مقابل الزمن المرتبط بتلك النافذة. يُسمى هذا الرسم البياني مخططًا طيفيًا . وهو أساس عدد من تقنيات التحليل الطيفي، مثل تحويل فورييه قصير المدى والمويجات .
  • يشير مصطلح "الطيف" عمومًا إلى كثافة القدرة الطيفية، كما ذُكر سابقًا، والتي تُصوّر توزيع محتوى الإشارة على التردد. بالنسبة لدوال النقل (مثل مخطط بود ، ونبضة التردد المتغير )، يمكن تمثيل استجابة التردد الكاملة بيانيًا في جزأين: القدرة مقابل التردد، والطور مقابل التردد - كثافة طيف الطور ، أو طيف الطور ، أو الطور الطيفي . في حالات أقل شيوعًا، قد يكون الجزآن هما الجزء الحقيقي والجزء التخيلي لدالة النقل. يجب عدم الخلط بين هذا وبين استجابة التردد لدالة النقل، والتي تتضمن أيضًا الطور (أو ما يُكافئه، جزءًا حقيقيًا وجزءًا تخيليًا) كدالة للتردد. استجابة النبضة في المجال الزمنيح(ت){\displaystyle h(t)}لا يمكن استخلاصها بشكل فريد من كثافة القدرة الطيفية وحدها دون جزء الطور. على الرغم من أن هذه أيضًا أزواج تحويل فورييه، إلا أنه لا يوجد تناظر (كما هو الحال في الارتباط الذاتي ) يُجبر تحويل فورييه على أن يكون ذا قيمة حقيقية. انظر: نبضة فائقة القصر#الطور الطيفي ، ضوضاء الطور ، تأخير المجموعة .
  • أحيانًا يُصادف المرء كثافة طيفية للسعة ( ASD )، وهي الجذر التربيعي لكثافة القدرة الطيفية (PSD). تُقاس كثافة الطيف للسعة لإشارة الجهد بوحدة V⋅Hz −1/2 . [ 20 ] يُعدّ هذا مفيدًا عندما يكون شكل الطيف ثابتًا نسبيًا، إذ تتناسب تغيرات كثافة الطيف للسعة مع تغيرات مستوى جهد الإشارة نفسه. ولكن من الناحية الرياضية، يُفضّل استخدام كثافة القدرة الطيفية، لأنه في هذه الحالة فقط تكون المساحة تحت المنحنى ذات دلالة من حيث القدرة الفعلية على جميع الترددات أو على نطاق ترددي مُحدد.

التطبيقات

أي إشارة يمكن تمثيلها كمتغير يتغير مع الزمن لها طيف ترددي خاص بها. يشمل ذلك إشارات مألوفة مثل الضوء المرئي (الذي يُدرك كلون )، والنوتات الموسيقية (التي تُدرك كدرجة صوتيةوالراديو/التلفزيون (المحدد بتردده، أو أحيانًا بطول موجته )، وحتى دوران الأرض المنتظم. عند عرض هذه الإشارات في شكل طيف ترددي، تتضح جوانب معينة من الإشارات المستقبلة أو العمليات الأساسية التي تُنتجها. في بعض الحالات، قد يتضمن الطيف الترددي ذروة واضحة تُقابل مكون الموجة الجيبية . بالإضافة إلى ذلك، قد توجد ذروات تُقابل توافقيات الذروة الأساسية، مما يشير إلى إشارة دورية ليست جيبية بسيطة. أو قد يُظهر الطيف المتصل فترات ترددية ضيقة مُعززة بشدة تُقابل الرنين، أو فترات ترددية ذات طاقة شبه معدومة كما هو الحال في مرشح الحجب .

الهندسة الكهربائية

مخطط طيفي لإشارة راديو FM، حيث يمثل المحور الأفقي التردد، بينما يمثل المحور الرأسي الزمن المتزايد لأعلى.

يُعدّ مفهوم طيف القدرة للإشارة واستخدامه أساسيًا في الهندسة الكهربائية ، لا سيما في أنظمة الاتصالات الإلكترونية ، بما في ذلك الاتصالات اللاسلكية والرادارات والأنظمة ذات الصلة، بالإضافة إلى تقنية الاستشعار عن بُعد السلبي . وتُستخدم أجهزة إلكترونية تُسمى محللات الطيف لمراقبة وقياس أطياف قدرة الإشارات.

يقيس محلل الطيف مقدار تحويل فورييه قصير المدى (STFT) لإشارة دخل. إذا أمكن اعتبار الإشارة المراد تحليلها عملية ثابتة، فإن تحويل فورييه قصير المدى يُعد تقديرًا جيدًا ومُنعّمًا لكثافة طيف قدرتها.

علم الكونيات

يتم تحديد التقلبات البدائية ، وهي اختلافات الكثافة في الكون المبكر، من خلال طيف القدرة الذي يعطي قوة الاختلافات كدالة للمقياس المكاني.

انظر أيضاً

ملحوظات

  1. لا يزال بعض المؤلفين، على سبيل المثال ( Risken & Frank 1996 ، ص 30) ، يستخدمون تحويل فورييه غير المعياري بطريقة رسمية لصياغة تعريف لكثافة القدرة الطيفية  x^(ω)x^*(ω)=2πو(ω)دلتا(ω-ω)،{\displaystyle \langle {\hat {x}}(\omega ){\hat {x}}^{\ast }(\omega ')\rangle =2\pi f(\omega )\delta (\omega -\omega '),} أيندلتا(ω-ω){\displaystyle \delta (\omega -\omega ')}هي دالة ديراك دلتا . قد تكون هذه العبارات الرسمية مفيدة أحيانًا لتوجيه الحدس، ولكن يجب استخدامها دائمًا بأقصى درجات الحذر.
  2. تُضفي نظرية وينر -خينشين معنىً على هذه الصيغة لأي عملية ثابتة بالمعنى الواسع في ظل فرضيات أضعف:Rxx{\displaystyle R_{xx}} لا يشترط أن تكون قابلة للتكامل المطلق، يكفي أن تكون موجودة. لكن لا يمكن تفسير التكامل بالطريقة المعتادة. تصبح الصيغة منطقية أيضًا إذا فُسِّرت على أنها تتضمن توزيعات (بمعنى لوران شوارتز ، وليس بمعنى دالة التوزيع التراكمي الإحصائي ) بدلًا من الدوال.Rxx{\displaystyle R_{xx}}إذا كانت متصلة، فيمكن استخدام نظرية بوخنر لإثبات أن تحويل فورييه الخاص بها موجود كمقياس موجب ، ودالة التوزيع الخاصة به هي F (ولكن ليس بالضرورة كدالة وليس بالضرورة أن تمتلك كثافة احتمالية).
  1. 1 2 3 P Stoica & R Moses (2005). "التحليل الطيفي للإشارات" (PDF) .
  2. مارال 2004 .
  3. نورتون وكارتشوب 2003 .
  4. بيروليني 2007 ، ص 83.
  5. باشوتا، روديجر (5 أبريل 2005). "كثافة القدرة الطيفية" . rp-photonics.com . مؤرشف من الأصل بتاريخ 15 أبريل 2024. تم الاطلاع عليه بتاريخ 26 يونيو 2024 .
  6. ^ أوبنهايم وفيرغيز 2016 ، ص. 60.
  7. شتاين 2000 ، ص 108، 115.
  8. ^ أوبنهايم وفيرغيز 2016 ، ص. 14.
  9. ^ أوبنهايم وفيرغيز 2016 ، ص 422-423.
  10. ميلر وتشيلدرز 2012 ، ص 429-431.
  11. ميلر وتشيلدرز 2012 ، ص 433.
  12. دينيس وارد ريكر (2003). معالجة إشارات الصدى . سبرينغر. ISBN 978-1-4020-7395-3.
  13. براون وهوانغ 1997 .
  14. ميلر وتشيلدرز 2012 ، ص 431.
  15. دافنبورت وروت 1987 .
  16. ويليام د. بيني (2009). "دورة معالجة الإشارات، الفصل 7" .
  17. سيف، ماركوس؛ قربانيتيماد، أرمين؛ فاغنر، فابيان؛ هاغيله، دانيال (2026). "التقدير الصحيح للأطياف ذات الرتبة الأعلى: من التحديات النظرية إلى التنفيذ العملي متعدد القنوات في SignalSnap" . معالجة الإشارات الرقمية . 173 105893. doi : 10.1016/j.dsp.2026.105893 .
  18. إيرانمانش ورودريغيز -فيليغاس 2017 .
  19. امتياز ورودريغيز -فيليغاس 2014 .
  20. مايكل سيرنا وأودري ف. هارفي (2000). "أساسيات تحليل وقياس الإشارات باستخدام تحويل فورييه السريع" (ملف PDF) . مؤرشف من الأصل في 15 سبتمبر 2012.

مراجع

  • بيروليني، أليساندرو (2007). هندسة الموثوقية . برلين؛ نيويورك: سبرينغر ساينس آند بيزنس ميديا. ISBN 978-3-540-49388-4.
  • براون، روبرت غروفر؛ هوانغ، باتريك واي سي (1997). مقدمة في الإشارات العشوائية وتطبيق ترشيح كالمان مع تمارين وحلول ماتلاب . نيويورك: وايلي-ليس. ISBN 978-0-471-12839-7.
  • دافنبورت، ويلبر ب. (الابن)؛ روت، ويليام ل. (1987). مقدمة في نظرية الإشارات العشوائية والضوضاء . نيويورك: مطبعة وايلي-IEEE. ISBN 978-0-87942-235-6.
  • إمتياز، سيد أنس؛ رودريغيز-فيليغاس، إستر (2014). "خوارزمية منخفضة التكلفة الحسابية للكشف عن نوم حركة العين السريعة باستخدام تخطيط كهربية الدماغ أحادي القناة" . حوليات الهندسة الطبية الحيوية . 42 (11): 2344-2359 . doi : 10.1007/s10439-014-1085-6 . PMC 4204008. PMID 25113231 .  
  • إيرانمانش، سام؛ رودريغيز-فيليغاس، إستر (2017). "نظام كشف مغزل النوم منخفض الطاقة للغاية على شريحة". معاملات IEEE في الدوائر والأنظمة الطبية الحيوية . 11 (4): 858-866 . Bibcode : 2017ITBC...11..858I . doi : 10.1109/TBCAS.2017.2690908 . hdl : 10044/1/46059 . PMID: 28541914. S2CID : 206608057 .  
  • مارال، جيرارد (2004). شبكات VSAT . غرب ساسكس، إنجلترا؛ هوبوكين، نيوجيرسي: وايلي. ISBN 978-0-470-86684-9.
  • ميلر، سكوت؛ تشيلدرز، دونالد (2012). الاحتمالات والعمليات العشوائية . بوسطن، ماساتشوستس: أكاديميك برس. ISBN 978-0-12-386981-4. OCLC 696092052 . 
  • نورتون، إم بي؛ كارتشوب، دي جي (2003). أساسيات تحليل الضوضاء والاهتزازات للمهندسين . كامبريدج: مطبعة جامعة كامبريدج. ISBN 978-0-521-49913-2.
  • أوبنهايم، آلان ف.؛ فيرغيز، جورج س. (2016). الإشارات والأنظمة والاستدلال . بوسطن: بيرسون. ISBN 978-0-13-394328-3.
  • ريكن، هانس؛ فرانك، حتى (1996). معادلة فوكر-بلانك . نيويورك: سبرينغر للعلوم والإعلام التجاري. رقم ISBN 978-3-540-61530-9.
  • شتاين، جوناثان Y. (2000). معالجة الإشارات الرقمية . نيويورك واينهايم: وايلي إنترساينس. رقم ISBN 978-0-471-29546-4.