पीरियडोग्राम: Difference between revisions
No edit summary |
No edit summary |
||
Line 1: | Line 1: | ||
{{Short description|Estimate of the spectral density of a signal}} | {{Short description|Estimate of the spectral density of a signal}} | ||
पीरियडोग्राम सिग्नल के [[वर्णक्रमीय घनत्व]] का अनुमान है। यह शब्द 1898 में [[आर्थर शूस्टर]] द्वारा गढ़ा गया था।<ref name="Schuster"/> आज, पीरियडोग्राम अधिक परिष्कृत विधियों का घटक है ([[वर्णक्रमीय अनुमान]] देखें)। [[फिल्टर के लिए]] और [[खिड़की के कार्य]] के आयाम बनाम आवृत्ति विशेषताओं की जांच करने के लिए यह सबसे सामान्य उपकरण है। [[ स्पेकट्रूम विशेष्यग्य ]] को पीरियडोग्राम के समय-अनुक्रम के रूप में भी प्रयुक्त किया जाता है। | |||
== परिभाषा == | == परिभाषा == | ||
आज कम से कम दो अलग-अलग परिभाषाएँ उपयोग में हैं।<ref name="Comparison"/> उनमें से एक में समय-औसत सम्मिलित है,<ref name="Wolfram"/> और एक नहीं।<ref name="Matlab"/> समय-औसत भी अन्य लेखों (बार्टलेट की विधि और वेल्च की विधि) का दायरा है। यह लेख समय-औसत के बारे में नहीं है। यहाँ ब्याज की परिभाषा यह है कि | आज कम से कम दो अलग-अलग परिभाषाएँ उपयोग में हैं।<ref name="Comparison"/> उनमें से एक में समय-औसत सम्मिलित है,<ref name="Wolfram"/> और एक नहीं।<ref name="Matlab"/> समय-औसत भी अन्य लेखों (बार्टलेट की विधि और वेल्च की विधि) का दायरा है। यह लेख समय-औसत के बारे में नहीं है। यहाँ ब्याज की परिभाषा यह है कि सतत कार्य की शक्ति वर्णक्रमीय घनत्व, अपने ऑटो-सहसंबंध फलन का [[फूरियर रूपांतरण]] है (फूरियर ट्रांसफॉर्म#क्रॉस-सहसंबंध प्रमेय| क्रॉस-सहसंबंध प्रमेय, स्पेक्ट्रल घनत्व# पावर स्पेक्ट्रल घनत्व, और वीनर-खिनचिन प्रमेय देखें): | ||
:<math>\mathcal{F}\{x(t)\circledast x^*(-t)\} = X(f)\cdot X^*(f) = \left| X(f) \right|^2.</math> | :<math>\mathcal{F}\{x(t)\circledast x^*(-t)\} = X(f)\cdot X^*(f) = \left| X(f) \right|^2.</math> | ||
Line 11: | Line 11: | ||
[[File:Periodogram.png|thumb|400px|पीरियोडोग्राम विधि द्वारा परिकलित दो साइनसोइडल आधार कार्यों का एक पावर स्पेक्ट्रम (परिमाण-वर्ग)।]] | [[File:Periodogram.png|thumb|400px|पीरियोडोग्राम विधि द्वारा परिकलित दो साइनसोइडल आधार कार्यों का एक पावर स्पेक्ट्रम (परिमाण-वर्ग)।]] | ||
[[File:Periodogram windows.png|thumb|400px|दो पावर स्पेक्ट्रा (परिमाण-वर्ग) (आयताकार और हैमिंग [[विंडो फंक्शन|विंडो फलन]] प्लस पृष्ठभूमि शोर), पीरियोडोग्राम विधि द्वारा गणना की गई।]]पैरामीटर {{mvar|T,}} के पर्याप्त छोटे मूल्यों के लिए के लिए क्षेत्र {{math|''X''(''f'')}} | [[File:Periodogram windows.png|thumb|400px|दो पावर स्पेक्ट्रा (परिमाण-वर्ग) (आयताकार और हैमिंग [[विंडो फंक्शन|विंडो फलन]] प्लस पृष्ठभूमि शोर), पीरियोडोग्राम विधि द्वारा गणना की गई।]]पैरामीटर {{mvar|T,}} के पर्याप्त छोटे मूल्यों के लिए के लिए क्षेत्र {{math|''X''(''f'')}} मनमाना ढंग- सही सन्निकटन क्षेत्र में देखा जा सकता है<math>-\tfrac{1}{2T} < f < \tfrac{1}{2T}</math> फलन का: | ||
:<math>X_{1/T}(f)\ \triangleq \sum_{k=-\infty}^{\infty} X\left(f - k/T\right),</math> | :<math>X_{1/T}(f)\ \triangleq \sum_{k=-\infty}^{\infty} X\left(f - k/T\right),</math> | ||
Line 19: | Line 19: | ||
:<math>X_{1/T}\left(\tfrac{k}{NT}\right) = \sum_{n=-\infty}^\infty \underbrace{T\cdot x(nT)}_{x[n]}\cdot e^{-i 2\pi \frac{kn}{N}},</math> | :<math>X_{1/T}\left(\tfrac{k}{NT}\right) = \sum_{n=-\infty}^\infty \underbrace{T\cdot x(nT)}_{x[n]}\cdot e^{-i 2\pi \frac{kn}{N}},</math> | ||
जहाँ {{mvar|k}} | जहाँ {{mvar|k}} पूर्णांक है। की आवधिकता <math>e^{-i 2\pi \frac{kn}{N}}</math>[[असतत फूरियर रूपांतरण]] के संदर्भ में इसे बहुत सरलता से लिखा जा सकता है: | ||
:<math>X_{1/T}\left(\tfrac{k}{NT}\right) = \underbrace{\sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}},}_{DFT}\quad \scriptstyle{\text{(sum over any }n\text{-sequence of length }N)},</math> | :<math>X_{1/T}\left(\tfrac{k}{NT}\right) = \underbrace{\sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}},}_{DFT}\quad \scriptstyle{\text{(sum over any }n\text{-sequence of length }N)},</math> | ||
जहाँ <math>x_{_N}</math> | जहाँ <math>x_{_N}</math> आवधिक योग है:<math>x_{_N}[n]\ \triangleq \sum_{m=-\infty}^{\infty} x[n-mN].</math> | ||
जब सभी पूर्णांकों के लिए मूल्यांकन किया जाता है, {{mvar|k}},और के बीच 0 {{mvar|N}}-1, सरणी: | जब सभी पूर्णांकों के लिए मूल्यांकन किया जाता है, {{mvar|k}},और के बीच 0 {{mvar|N}}-1, सरणी: | ||
:<math>S\left(\tfrac{k}{NT}\right) = \left| \sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}} \right|^2</math> | :<math>S\left(\tfrac{k}{NT}\right) = \left| \sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}} \right|^2</math> | ||
पीरियोग्राम है।<ref name="Matlab"/><ref name=Oppenheim/><ref name=Rabiner/> | |||
== अनुप्रयोग == | == अनुप्रयोग == | ||
[[File:Periodogram for Proxima Centauri b.jpg|thumb|[[प्रॉक्सिमा सेंटौरी बी]] के लिए पीरियडोग्राम नीचे दिखाया गया है।<ref>{{cite web|title=Do-it-yourself Science — is Proxima c hiding in this graph?|url=https://www.eso.org/public/images/potw1737a/|website=www.eso.org|access-date=11 September 2017}}</ref>]]जब | [[File:Periodogram for Proxima Centauri b.jpg|thumb|[[प्रॉक्सिमा सेंटौरी बी]] के लिए पीरियडोग्राम नीचे दिखाया गया है।<ref>{{cite web|title=Do-it-yourself Science — is Proxima c hiding in this graph?|url=https://www.eso.org/public/images/potw1737a/|website=www.eso.org|access-date=11 September 2017}}</ref>]]जब प्राथमिकी फिल्टर या विंडो फलन की विस्तृत विशेषताओं की जांच करने के लिए पीरियडोग्राम का उपयोग किया जाता है, तो पैरामीटर {{mvar|N}} को {{math|''x''[''n'']}} गैर-शून्य अवधि के कई गुणकों के रूप में चुना जाता है अनुक्रम, की गैर-शून्य अवधि के कई गुणकों के रूप में चुना जाता है, जिसे शून्य-पैडिंग कहा जाता है (देखें {{slink|Discrete-time Fourier transform|डीटीएफटी का नमूना लेना|nopage=y}}).{{ | ||
efn-ua|{{mvar|N}} is designated {{mvar|NFFT}} in the Matlab and Octave applications. | efn-ua|{{mvar|N}} is designated {{mvar|NFFT}} in the Matlab and Octave applications. | ||
}} जब [[फ़िल्टर बैंक]] को प्रयुक्त करने के लिए तो {{math|''x''[''n'']}} अनुक्रम की गैर-शून्य अवधि के कई उप-गुणक होते हैं (देखें {{slink|Discrete-time Fourier transform|डीटीएफटी का नमूना लेना|nopage=y}}). | }} जब [[फ़िल्टर बैंक]] को प्रयुक्त करने के लिए तो {{math|''x''[''n'']}} अनुक्रम की गैर-शून्य अवधि के कई उप-गुणक होते हैं (देखें {{slink|Discrete-time Fourier transform|डीटीएफटी का नमूना लेना|nopage=y}}). | ||
पीरियडोग्राम की कमियों में से एक यह है कि गणना में उपयोग किए जाने वाले नमूनों की संख्या बढ़ने पर दी गई [[आवृत्ति]] पर विचरण कम नहीं होता है। यह कम सिग्नल-टू-शोर अनुपात पर नोइज़लाइक सिग्नल या यहां तक कि साइनसोइड्स का विश्लेषण करने के लिए आवश्यक औसत प्रदान नहीं करता है। विंडो फलन और फ़िल्टर आवेग प्रतिक्रियाएँ नीरव हैं, लेकिन कई अन्य संकेतों के लिए वर्णक्रमीय अनुमान के अधिक परिष्कृत तरीकों की आवश्यकता होती है। प्रक्रिया के भागों के रूप में विकल्पों में से दो पीरियडोग्राम का उपयोग करते हैं: | पीरियडोग्राम की कमियों में से एक यह है कि गणना में उपयोग किए जाने वाले नमूनों की संख्या बढ़ने पर दी गई [[आवृत्ति]] पर विचरण कम नहीं होता है। यह कम सिग्नल-टू-शोर अनुपात पर नोइज़लाइक सिग्नल या यहां तक कि साइनसोइड्स का विश्लेषण करने के लिए आवश्यक औसत प्रदान नहीं करता है। विंडो फलन और फ़िल्टर आवेग प्रतिक्रियाएँ नीरव हैं, लेकिन कई अन्य संकेतों के लिए वर्णक्रमीय अनुमान के अधिक परिष्कृत तरीकों की आवश्यकता होती है। प्रक्रिया के भागों के रूप में विकल्पों में से दो पीरियडोग्राम का उपयोग करते हैं: | ||
*औसत पीरियडोग्राम की विधि,<ref name="Engelberg"/> अधिक सामान्यतः वेल्च की विधि के रूप में जाना जाता है,<ref name="Welch"/><ref name="Welch2"/> | *औसत पीरियडोग्राम की विधि,<ref name="Engelberg"/> अधिक सामान्यतः वेल्च की विधि के रूप में जाना जाता है,<ref name="Welch"/><ref name="Welch2"/> लंबे x [n] अनुक्रम को कई छोटे, और संभवतः अतिव्यापी, बाद में विभाजित करता है। यह प्रत्येक के विंडोड पीरियडोग्राम की गणना करता है, और सरणी औसत की गणना करता है, यद्दपि सरणी जहां प्रत्येक तत्व सभी पीरियडोग्राम के संबंधित तत्वों का औसत होता है। [[स्थिर प्रक्रिया]]ओं के लिए, यह प्रत्येक तत्व के शोर भिन्नता को पीरियोडोग्राम की संख्या के व्युत्क्रम के बराबर लगभग कारक से कम कर देता है। | ||
*[[ चौरसाई | स्मूथिंग]] समय के अतिरिक्त आवृत्ति में | *[[ चौरसाई | स्मूथिंग]] समय के अतिरिक्त आवृत्ति में औसत तकनीक है। स्मूथेड पीरियडोग्राम को कभी-कभी स्पेक्ट्रल प्लॉट के रूप में जाना जाता है।<ref name="smoothing"/><ref name="dataplot"/> | ||
पेरियोडोग्राम-आधारित तकनीकें छोटे पूर्वाग्रहों का परिचय देती हैं जो कुछ अनुप्रयोगों में अस्वीकार्य हैं। अन्य तकनीकें जो पीरियडोग्राम्स पर विश्वास नहीं करती हैं, [[वर्णक्रमीय घनत्व अनुमान]] लेख में प्रस्तुत की जाती हैं। | पेरियोडोग्राम-आधारित तकनीकें छोटे पूर्वाग्रहों का परिचय देती हैं जो कुछ अनुप्रयोगों में अस्वीकार्य हैं। अन्य तकनीकें जो पीरियडोग्राम्स पर विश्वास नहीं करती हैं, [[वर्णक्रमीय घनत्व अनुमान]] लेख में प्रस्तुत की जाती हैं। | ||
Line 49: | Line 49: | ||
* असतत-समय फूरियर रूपांतरण | * असतत-समय फूरियर रूपांतरण | ||
*[[कम से कम वर्ग वर्णक्रमीय विश्लेषण]], डेटा में पीरियडोग्राम की गणना के लिए जो समान दूरी पर नहीं है | *[[कम से कम वर्ग वर्णक्रमीय विश्लेषण]], डेटा में पीरियडोग्राम की गणना के लिए जो समान दूरी पर नहीं है | ||
*MUSIC_(एल्गोरिदम) (MUSIC), | *MUSIC_(एल्गोरिदम) (MUSIC), लोकप्रिय पैरामीट्रिक [[सुपर-रिज़ॉल्यूशन इमेजिंग]] विधि | ||
* [[एसएएमवी (एल्गोरिदम)]] | * [[एसएएमवी (एल्गोरिदम)]] | ||
Revision as of 11:25, 3 April 2023
पीरियडोग्राम सिग्नल के वर्णक्रमीय घनत्व का अनुमान है। यह शब्द 1898 में आर्थर शूस्टर द्वारा गढ़ा गया था।[1] आज, पीरियडोग्राम अधिक परिष्कृत विधियों का घटक है (वर्णक्रमीय अनुमान देखें)। फिल्टर के लिए और खिड़की के कार्य के आयाम बनाम आवृत्ति विशेषताओं की जांच करने के लिए यह सबसे सामान्य उपकरण है। स्पेकट्रूम विशेष्यग्य को पीरियडोग्राम के समय-अनुक्रम के रूप में भी प्रयुक्त किया जाता है।
परिभाषा
आज कम से कम दो अलग-अलग परिभाषाएँ उपयोग में हैं।[2] उनमें से एक में समय-औसत सम्मिलित है,[3] और एक नहीं।[4] समय-औसत भी अन्य लेखों (बार्टलेट की विधि और वेल्च की विधि) का दायरा है। यह लेख समय-औसत के बारे में नहीं है। यहाँ ब्याज की परिभाषा यह है कि सतत कार्य की शक्ति वर्णक्रमीय घनत्व, अपने ऑटो-सहसंबंध फलन का फूरियर रूपांतरण है (फूरियर ट्रांसफॉर्म#क्रॉस-सहसंबंध प्रमेय| क्रॉस-सहसंबंध प्रमेय, स्पेक्ट्रल घनत्व# पावर स्पेक्ट्रल घनत्व, और वीनर-खिनचिन प्रमेय देखें):
संगणना
पैरामीटर T, के पर्याप्त छोटे मूल्यों के लिए के लिए क्षेत्र X(f) मनमाना ढंग- सही सन्निकटन क्षेत्र में देखा जा सकता है फलन का:
जो नमूनों x(nT) द्वारा सही रूप से निर्धारित किया जाता है जो x(t) गैर-शून्य अवधि का विस्तार करता है (असतत-समय फूरियर रूपांतरण देखें)।
और पैरामीटर N, के पर्याप्त बड़े मूल्यों के लिए प्रपत्र के योग द्वारा मनमाने ढंग से निकट आवृत्ति पर मूल्यांकन किया जा सकता है:
जहाँ k पूर्णांक है। की आवधिकता असतत फूरियर रूपांतरण के संदर्भ में इसे बहुत सरलता से लिखा जा सकता है:
जहाँ आवधिक योग है:
जब सभी पूर्णांकों के लिए मूल्यांकन किया जाता है, k,और के बीच 0 N-1, सरणी:
अनुप्रयोग
जब प्राथमिकी फिल्टर या विंडो फलन की विस्तृत विशेषताओं की जांच करने के लिए पीरियडोग्राम का उपयोग किया जाता है, तो पैरामीटर N को x[n] गैर-शून्य अवधि के कई गुणकों के रूप में चुना जाता है अनुक्रम, की गैर-शून्य अवधि के कई गुणकों के रूप में चुना जाता है, जिसे शून्य-पैडिंग कहा जाता है (देखें § डीटीएफटी का नमूना लेना).[upper-alpha 1] जब फ़िल्टर बैंक को प्रयुक्त करने के लिए तो x[n] अनुक्रम की गैर-शून्य अवधि के कई उप-गुणक होते हैं (देखें § डीटीएफटी का नमूना लेना).
पीरियडोग्राम की कमियों में से एक यह है कि गणना में उपयोग किए जाने वाले नमूनों की संख्या बढ़ने पर दी गई आवृत्ति पर विचरण कम नहीं होता है। यह कम सिग्नल-टू-शोर अनुपात पर नोइज़लाइक सिग्नल या यहां तक कि साइनसोइड्स का विश्लेषण करने के लिए आवश्यक औसत प्रदान नहीं करता है। विंडो फलन और फ़िल्टर आवेग प्रतिक्रियाएँ नीरव हैं, लेकिन कई अन्य संकेतों के लिए वर्णक्रमीय अनुमान के अधिक परिष्कृत तरीकों की आवश्यकता होती है। प्रक्रिया के भागों के रूप में विकल्पों में से दो पीरियडोग्राम का उपयोग करते हैं:
- औसत पीरियडोग्राम की विधि,[8] अधिक सामान्यतः वेल्च की विधि के रूप में जाना जाता है,[9][10] लंबे x [n] अनुक्रम को कई छोटे, और संभवतः अतिव्यापी, बाद में विभाजित करता है। यह प्रत्येक के विंडोड पीरियडोग्राम की गणना करता है, और सरणी औसत की गणना करता है, यद्दपि सरणी जहां प्रत्येक तत्व सभी पीरियडोग्राम के संबंधित तत्वों का औसत होता है। स्थिर प्रक्रियाओं के लिए, यह प्रत्येक तत्व के शोर भिन्नता को पीरियोडोग्राम की संख्या के व्युत्क्रम के बराबर लगभग कारक से कम कर देता है।
- स्मूथिंग समय के अतिरिक्त आवृत्ति में औसत तकनीक है। स्मूथेड पीरियडोग्राम को कभी-कभी स्पेक्ट्रल प्लॉट के रूप में जाना जाता है।[11][12]
पेरियोडोग्राम-आधारित तकनीकें छोटे पूर्वाग्रहों का परिचय देती हैं जो कुछ अनुप्रयोगों में अस्वीकार्य हैं। अन्य तकनीकें जो पीरियडोग्राम्स पर विश्वास नहीं करती हैं, वर्णक्रमीय घनत्व अनुमान लेख में प्रस्तुत की जाती हैं।
यह भी देखें
- मिलान फ़िल्टर
- रेडॉन ट्रांसफ़ॉर्म (रेडॉन ट्रांसफ़ॉर्म)
- वेल्च की विधि
- बार्टलेट की विधि
- असतत-समय फूरियर रूपांतरण
- कम से कम वर्ग वर्णक्रमीय विश्लेषण, डेटा में पीरियडोग्राम की गणना के लिए जो समान दूरी पर नहीं है
- MUSIC_(एल्गोरिदम) (MUSIC), लोकप्रिय पैरामीट्रिक सुपर-रिज़ॉल्यूशन इमेजिंग विधि
- एसएएमवी (एल्गोरिदम)
टिप्पणियाँ
- ↑ N is designated NFFT in the Matlab and Octave applications.
संदर्भ
- ↑ Schuster, Arthur (January 1898). "On the investigation of hidden periodicities with application to a supposed 26 day period of meteorological phenomena" (PDF). Terrestrial Magnetism. 3 (1): 13–41. Bibcode:1898TeMag...3...13S. doi:10.1029/TM003i001p00013.
It is convenient to have a word for some representation of a variable quantity which shall correspond to the 'spectrum' of a luminous radiation. I propose the word periodogram, and define it more particularly in the following way.
- ↑ McSweeney, Laura A. (2004-05-14). "Comparison of periodogram tests". Journal of Statistical Computation and Simulation. online ($50). 76 (4): 357–369. doi:10.1080/10629360500107618. S2CID 120439605.
- ↑ "Periodogram—Wolfram Language Documentation".
- ↑ 4.0 4.1 "Periodogram power spectral density estimate - MATLAB periodogram".
- ↑ Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). Discrete-time signal processing (2nd ed.). Upper Saddle River, N.J.: Prentice Hall. p. 732 (10.55). ISBN 0-13-754920-2.
- ↑ Rabiner, Lawrence R.; Gold, Bernard (1975). "6.18". Theory and application of digital signal processing. Englewood Cliffs, N.J.: Prentice-Hall. pp. 415. ISBN 0-13-914101-4.
- ↑ "Do-it-yourself Science — is Proxima c hiding in this graph?". www.eso.org. Retrieved 11 September 2017.
- ↑ Engelberg, S. (2008), Digital Signal Processing: An Experimental Approach, Springer, Chap. 7 p. 56
- ↑ Welch, Peter D. (June 1967). "The use of Fast Fourier Transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms". IEEE Transactions on Audio and Electroacoustics. AU-15 (2): 70–73. Bibcode:1967ITAE...15...70W. doi:10.1109/TAU.1967.1161901.
- ↑ "Welch's power spectral density estimate - MATLAB pwelch".
- ↑ Spectral Plot, from the NIST Engineering Statistics Handbook.
- ↑
"DATAPLOT Reference Manual" (PDF). NIST.gov. National Institute of Standards and Technology (NIST). 1997-03-11. Retrieved 2019-06-14.
The spectral plot is essentially a "smoothed" periodogram where the smoothing is done in the frequency domain.
अग्रिम पठन
- Box, George E. P.; Jenkins, Gwilym M. (1976). Time series analysis: Forecasting and control. San Francisco: Holden-Day.
- Scargle, J.D. (December 15, 1982). "Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data". Astrophysical Journal, Part 1. 263: 835–853. Bibcode:1982ApJ...263..835S. doi:10.1086/160554.
- Vaughan, Simon; Uttley, Philip (2006). "Detecting X-ray QPOs in active galaxies". Advances in Space Research. 38 (7): 1405–1408. arXiv:astro-ph/0506456. Bibcode:2006AdSpR..38.1405V. doi:10.1016/j.asr.2005.02.064. S2CID 21054467.