स्पेक्ट्रल घनत्व अनुमान: Difference between revisions
No edit summary |
No edit summary |
||
Line 5: | Line 5: | ||
== अवलोकन == | == अवलोकन == | ||
[[Image:Voice waveform and spectrum.png|thumb|आवाज तरंग और इसकी आवृत्ति | [[Image:Voice waveform and spectrum.png|thumb|आवाज तरंग और इसकी आवृत्ति वर्णक्रम का उदाहरण]] | ||
[[Image:triangle-td and fd.png|thumb|एक आवधिक तरंग (त्रिकोण तरंग) और इसकी आवृत्ति | [[Image:triangle-td and fd.png|thumb|एक आवधिक तरंग (त्रिकोण तरंग) और इसकी आवृत्ति वर्णक्रम, 220 हर्ट्ज पर एक मौलिक आवृत्ति दिखा रही है जिसके बाद 220 हर्ट्ज के गुणक (हार्मोनिक्स) हैं।]] | ||
[[File:Comparison of periodogram and Welch methods of spectral density estimation.png|thumb|तुलना के लिए, संगीत के एक खंड की शक्ति वर्णक्रमीय घनत्व दो अलग-अलग तरीकों से अनुमानित है।]]वर्णक्रम विश्लेषण, जिसे [[वर्णक्रमीय अनुमान]] विश्लेषण या वर्णक्रमीय घनत्व अनुमान के रूप में भी जाना जाता है, एक जटिल संकेत को सरल भागों में विघटित करने की तकनीकी प्रक्रिया है। जैसा कि ऊपर बताया गया है, कई भौतिक प्रक्रियाओं को कई अलग-अलग आवृत्ति घटकों के योग के रूप में वर्णित किया गया है। कोई भी प्रक्रिया जो विभिन्न राशियों (जैसे आयाम, शक्तियाँ, तीव्रता) बनाम आवृत्ति (या चरण (तरंगें)) की मात्रा निर्धारित करती है, उसे वर्णक्रम विश्लेषण कहा जा सकता है। | [[File:Comparison of periodogram and Welch methods of spectral density estimation.png|thumb|तुलना के लिए, संगीत के एक खंड की शक्ति वर्णक्रमीय घनत्व दो अलग-अलग तरीकों से अनुमानित है।]]वर्णक्रम विश्लेषण, जिसे [[वर्णक्रमीय अनुमान]] विश्लेषण या वर्णक्रमीय घनत्व अनुमान के रूप में भी जाना जाता है, एक जटिल संकेत को सरल भागों में विघटित करने की तकनीकी प्रक्रिया है। जैसा कि ऊपर बताया गया है, कई भौतिक प्रक्रियाओं को कई अलग-अलग आवृत्ति घटकों के योग के रूप में वर्णित किया गया है। कोई भी प्रक्रिया जो विभिन्न राशियों (जैसे आयाम, शक्तियाँ, तीव्रता) बनाम आवृत्ति (या चरण (तरंगें)) की मात्रा निर्धारित करती है, उसे वर्णक्रम विश्लेषण कहा जा सकता है। | ||
वर्णक्रम विश्लेषण पूरे संकेत पर किया जा सकता है। वैकल्पिक रूप से, एक संकेत को छोटे खंडों में तोड़ा जा सकता है (कभी-कभी ''फ्रेम'' कहा जाता है), और इन अलग-अलग खंडों पर वर्णक्रम विश्लेषण लागू किया जा सकता है। आवधिक कार्य (जैसे <math>\sin (t)</math>) इस उप-विभाजन के लिए विशेष रूप से उपयुक्त हैं। गैर-आवधिक कार्यों के विश्लेषण के लिए सामान्य गणितीय विधियाँ [[फूरियर विश्लेषण|सांध्वनिक विश्लेषण]] की श्रेणी में आती हैं। | वर्णक्रम विश्लेषण पूरे संकेत पर किया जा सकता है। वैकल्पिक रूप से, एक संकेत को छोटे खंडों में तोड़ा जा सकता है (कभी-कभी ''फ्रेम'' कहा जाता है), और इन अलग-अलग खंडों पर वर्णक्रम विश्लेषण लागू किया जा सकता है। आवधिक कार्य (जैसे <math>\sin (t)</math>) इस उप-विभाजन के लिए विशेष रूप से उपयुक्त हैं। गैर-आवधिक कार्यों के विश्लेषण के लिए सामान्य गणितीय विधियाँ [[फूरियर विश्लेषण|सांध्वनिक विश्लेषण]] की श्रेणी में आती हैं। | ||
फलन का [[फूरियर रूपांतरण|सांध्वनिक]] [[फूरियर रूपांतरण|रूपांतरण]] एक आवृत्ति वर्णक्रम उत्पन्न करता है, जिसमें मूल संकेत के विषय में सभी जानकारी होती है, परन्तु एक अलग रूप में। इसका तात्पर्य यह है, कि [[उलटा फूरियर रूपांतरण|प्रतिलोम संध्वनिक रूपांतरण]] द्वारा मूल कार्य को पूरी तरह से पुनर्निर्मित (संश्लेषित) किया जा सकता है। सही पुनर्निर्माण के लिए, | फलन का [[फूरियर रूपांतरण|सांध्वनिक]] [[फूरियर रूपांतरण|रूपांतरण]] एक आवृत्ति वर्णक्रम उत्पन्न करता है, जिसमें मूल संकेत के विषय में सभी जानकारी होती है, परन्तु एक अलग रूप में। इसका तात्पर्य यह है, कि [[उलटा फूरियर रूपांतरण|प्रतिलोम संध्वनिक रूपांतरण]] द्वारा मूल कार्य को पूरी तरह से पुनर्निर्मित (संश्लेषित) किया जा सकता है। सही पुनर्निर्माण के लिए, वर्णक्रम विश्लेषक को प्रत्येक आवृत्ति घटक के [[आयाम]] और [[चरण]] (तरंगों) दोनों को संरक्षित करना चाहिए। जानकारी के इन दो टुकड़ों को 2-आयामी वेक्टर के रूप में, एक [[जटिल संख्या]] के रूप में, या परिमाण (आयाम) और चरण के रूप में ध्रुवीय निर्देशांक (अर्थात, एक चरण के रूप में) के रूप में दर्शाया जा सकता है। संकेत प्रसंस्करण में एक सामान्य विधि वर्ग आयाम, या [[शक्ति (भौतिकी)]] पर विचार करना है; इस विषय में परिणामी भूखंड को शक्ति वर्णक्रम कहा जाता है। | ||
प्रतिवर्तीता के कारण, सान्ध्वनिक रूपांतरण को समय के अतिरिक्त आवृत्ति के संदर्भ में फलन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह एक आवृत्ति कार्यक्षेत्र प्रतिनिधित्व है। रैखिक संचालन जो समय कार्यक्षेत्र में किए जा सकते हैं, उनके समकक्ष हैं, जो प्रायः आवृत्ति कार्यक्षेत्र में अधिक सरलता से किए जा सकते हैं। बारंबारता विश्लेषण भी रैखिक और गैर-रैखिक दोनों तरह के विभिन्न समय -कार्यक्षेत्र संचालनों के प्रभावों की समझ और व्याख्या को सरल करता है। उदाहरण के लिए, केवल गैर-रैखिक या [[समय-भिन्न प्रणाली]] आवृत्ति वर्णक्रम में नई आवृत्तियां निर्मित कर सकते हैं। | प्रतिवर्तीता के कारण, सान्ध्वनिक रूपांतरण को समय के अतिरिक्त आवृत्ति के संदर्भ में फलन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह एक आवृत्ति कार्यक्षेत्र प्रतिनिधित्व है। रैखिक संचालन जो समय कार्यक्षेत्र में किए जा सकते हैं, उनके समकक्ष हैं, जो प्रायः आवृत्ति कार्यक्षेत्र में अधिक सरलता से किए जा सकते हैं। बारंबारता विश्लेषण भी रैखिक और गैर-रैखिक दोनों तरह के विभिन्न समय -कार्यक्षेत्र संचालनों के प्रभावों की समझ और व्याख्या को सरल करता है। उदाहरण के लिए, केवल गैर-रैखिक या [[समय-भिन्न प्रणाली]] आवृत्ति वर्णक्रम में नई आवृत्तियां निर्मित कर सकते हैं। | ||
Line 19: | Line 19: | ||
== तकनीक == | == तकनीक == | ||
आधारभूत पीरियडोग्राम के नुकसान को कम करने के लिए स्पेक्ट्रल आकलन के लिए कई अन्य तकनीकों का विकास किया गया है। इन तकनीकों को आम तौर पर गैर-आवर्तिता वक्र आंकड़ों में विभाजित किया जा सकता है।<ref name=":0">{{Cite journal|last1=Stoica|first1=Petre|last2=Babu|first2=Prabhu|last3=Li|first3=Jian|date=January 2011|title=अलग-अलग मॉडल में विरल पैरामीटर अनुमान की नई विधि और अनियमित रूप से सैंपल किए गए डेटा के स्पेक्ट्रल विश्लेषण के लिए इसका उपयोग|url=https://ieeexplore.ieee.org/document/5599897|journal=IEEE Transactions on Signal Processing|volume=59|issue=1|pages=35–47|doi=10.1109/TSP.2010.2086452|bibcode=2011ITSP...59...35S |s2cid=15936187 |issn=1053-587X}}</ref> गैर- प्राचलिक दृष्टिकोण स्पष्ट रूप से [[सहप्रसरण]] या प्रक्रिया के | आधारभूत पीरियडोग्राम के नुकसान को कम करने के लिए स्पेक्ट्रल आकलन के लिए कई अन्य तकनीकों का विकास किया गया है। इन तकनीकों को आम तौर पर गैर-आवर्तिता वक्र आंकड़ों में विभाजित किया जा सकता है।<ref name=":0">{{Cite journal|last1=Stoica|first1=Petre|last2=Babu|first2=Prabhu|last3=Li|first3=Jian|date=January 2011|title=अलग-अलग मॉडल में विरल पैरामीटर अनुमान की नई विधि और अनियमित रूप से सैंपल किए गए डेटा के स्पेक्ट्रल विश्लेषण के लिए इसका उपयोग|url=https://ieeexplore.ieee.org/document/5599897|journal=IEEE Transactions on Signal Processing|volume=59|issue=1|pages=35–47|doi=10.1109/TSP.2010.2086452|bibcode=2011ITSP...59...35S |s2cid=15936187 |issn=1053-587X}}</ref> गैर- प्राचलिक दृष्टिकोण स्पष्ट रूप से [[सहप्रसरण]] या प्रक्रिया के वर्णक्रम का अनुमान लगाए बिना यह मानते हैं कि प्रक्रिया की कोई विशेष संरचना है। बुनियादी अनुप्रयोगों (जैसे वेल्च की विधि) के लिए उपयोग किए जाने वाले कुछ सबसे आम अनुमानक गैर- प्राचलिक अनुमानक हैं जो पीरियोडोग्राम से निकटता से संबंधित हैं। इसके विपरीत, प्राचलिक दृष्टिकोण मानते हैं कि अंतर्निहित [[स्थिर प्रक्रिया]] में एक निश्चित संरचना होती है जिसे पैरामीटर की एक छोटी संख्या का उपयोग करके वर्णित किया जा सकता है (उदाहरण के लिए, [[ऑटोरेग्रेसिव मूविंग एवरेज मॉडल|ऑटोरेग्रेसिव मूविंग एवरेज प्रारूप]] | ऑटो-रिग्रेसिव या मूविंग एवरेज प्रारूप का उपयोग करके)। इन दृष्टिकोणों में, कार्य प्रारूप के मापदंडों का अनुमान लगाना है जो स्टोकेस्टिक प्रक्रिया का वर्णन करता है। अर्ध- प्राचलिक विधियों का उपयोग करते समय, अंतर्निहित प्रक्रिया को एक गैर- प्राचलिक ढांचे का उपयोग करके तैयार किया जाता है, अतिरिक्त धारणा के साथ कि प्रारूप के गैर-शून्य घटकों की संख्या छोटी है (यानी, प्रारूप विरल है)। लापता डेटा रिकवरी के लिए भी इसी तरह के तरीकों का इस्तेमाल किया जा सकता है <ref>{{Cite journal|last1=Stoica|first1=Petre|last2=Li|first2=Jian|last3=Ling|first3=Jun|last4=Cheng|first4=Yubo|date=April 2009|title=एक गैर पैरामीट्रिक पुनरावृत्त अनुकूली दृष्टिकोण के माध्यम से लापता डेटा रिकवरी|url=http://dx.doi.org/10.1109/icassp.2009.4960347|journal=2009 IEEE International Conference on Acoustics, Speech and Signal Processing|pages=3369–3372 |publisher=IEEE|doi=10.1109/icassp.2009.4960347|isbn=978-1-4244-2353-8 }}</ref> साथ ही संपीड़ित संवेदन। | ||
निम्नलिखित गैर- प्राचलिक वर्णक्रमीय घनत्व अनुमान तकनीकों की आंशिक सूची है: | निम्नलिखित गैर- प्राचलिक वर्णक्रमीय घनत्व अनुमान तकनीकों की आंशिक सूची है: | ||
Line 29: | Line 29: | ||
* [[कम से कम वर्ग वर्णक्रमीय विश्लेषण]], ज्ञात आवृत्तियों के लिए [[कम से कम वर्गों]] के आधार पर | * [[कम से कम वर्ग वर्णक्रमीय विश्लेषण]], ज्ञात आवृत्तियों के लिए [[कम से कम वर्गों]] के आधार पर | ||
* [[गैर-समान असतत फूरियर रूपांतरण]] का उपयोग तब किया जाता है जब सिग्नल नमूने असमान रूप से समय श्रृंखला में होते हैं | * [[गैर-समान असतत फूरियर रूपांतरण]] का उपयोग तब किया जाता है जब सिग्नल नमूने असमान रूप से समय श्रृंखला में होते हैं | ||
* [[एकवचन स्पेक्ट्रम विश्लेषण]] एक गैर प्राचलिक विधि है जो वर्णक्रमीय घनत्व का अनुमान लगाने के लिए सहप्रसरण मैट्रिक्स के एकवचन मूल्य अपघटन का उपयोग करती है | * [[एकवचन स्पेक्ट्रम विश्लेषण|एकवचन वर्णक्रम विश्लेषण]] एक गैर प्राचलिक विधि है जो वर्णक्रमीय घनत्व का अनुमान लगाने के लिए सहप्रसरण मैट्रिक्स के एकवचन मूल्य अपघटन का उपयोग करती है | ||
* [[शॉर्ट-टाइम फूरियर रूपांतरण]] | * [[शॉर्ट-टाइम फूरियर रूपांतरण]] | ||
*[[सूचना क्षेत्र सिद्धांत]]#महत्वपूर्ण फ़िल्टर सूचना क्षेत्र सिद्धांत पर आधारित एक गैर- प्राचलिक विधि है जो शोर, अपूर्ण डेटा और वाद्य प्रतिक्रिया कार्यों से निपट सकता है | *[[सूचना क्षेत्र सिद्धांत]]#महत्वपूर्ण फ़िल्टर सूचना क्षेत्र सिद्धांत पर आधारित एक गैर- प्राचलिक विधि है जो शोर, अपूर्ण डेटा और वाद्य प्रतिक्रिया कार्यों से निपट सकता है | ||
Line 35: | Line 35: | ||
नीचे प्राचलिक तकनीकों की आंशिक सूची दी गई है: | नीचे प्राचलिक तकनीकों की आंशिक सूची दी गई है: | ||
* [[Autoregressive model]] (AR) आकलन, जो मानता है कि nth नमूना पिछले p नमूनों के साथ सहसंबद्ध है। | * [[Autoregressive model]] (AR) आकलन, जो मानता है कि nth नमूना पिछले p नमूनों के साथ सहसंबद्ध है। | ||
* [[मूविंग-एवरेज मॉडल]] (MA) का अनुमान, जो मानता है कि nth नमूना पिछले p नमूनों में शोर की शर्तों के साथ सहसंबद्ध है। | * [[मूविंग-एवरेज मॉडल|मूविंग-एवरेज प्रारूप]] (MA) का अनुमान, जो मानता है कि nth नमूना पिछले p नमूनों में शोर की शर्तों के साथ सहसंबद्ध है। | ||
* Autoregressive मूविंग एवरेज (ARMA) अनुमान, जो AR और MA | * Autoregressive मूविंग एवरेज (ARMA) अनुमान, जो AR और MA प्रारूप का सामान्यीकरण करता है। | ||
* [[संगीत (एल्गोरिदम)]] (संगीत) एक लोकप्रिय [[सुपर-रिज़ॉल्यूशन इमेजिंग]] विधि है। | * [[संगीत (एल्गोरिदम)]] (संगीत) एक लोकप्रिय [[सुपर-रिज़ॉल्यूशन इमेजिंग]] विधि है। | ||
* अधिकतम एन्ट्रॉपी वर्णक्रमीय अनुमान एसडीई के लिए उपयोगी एक सर्व-ध्रुवीय विधि है जब एकवचन वर्णक्रमीय विशेषताओं, जैसे तेज चोटियों की अपेक्षा की जाती है। | * अधिकतम एन्ट्रॉपी वर्णक्रमीय अनुमान एसडीई के लिए उपयोगी एक सर्व-ध्रुवीय विधि है जब एकवचन वर्णक्रमीय विशेषताओं, जैसे तेज चोटियों की अपेक्षा की जाती है। | ||
Line 48: | Line 48: | ||
प्राचलिक स्पेक्ट्रल अनुमान में, कोई मानता है कि सिग्नल एक स्थिर प्रक्रिया द्वारा तैयार किया गया है जिसमें स्पेक्ट्रल घनत्व फ़ंक्शन (एसडीएफ) है। <math>S(f; a_1, \ldots, a_p)</math> यह आवृत्ति का एक कार्य है <math>f</math> और <math>p</math> पैरामीटर <math>a_1, \ldots, a_p</math>.<ref name=Percival1993>{{cite book |last1=Percival|first1=Donald B.|last2=Walden|first2=Andrew T.|title=भौतिक अनुप्रयोगों के लिए वर्णक्रमीय विश्लेषण|date=1992|publisher=Cambridge University Press|isbn=9780521435413}}</ref> अनुमान समस्या तब इन पैरामीटरों का आकलन करने में से एक बन जाती है। | प्राचलिक स्पेक्ट्रल अनुमान में, कोई मानता है कि सिग्नल एक स्थिर प्रक्रिया द्वारा तैयार किया गया है जिसमें स्पेक्ट्रल घनत्व फ़ंक्शन (एसडीएफ) है। <math>S(f; a_1, \ldots, a_p)</math> यह आवृत्ति का एक कार्य है <math>f</math> और <math>p</math> पैरामीटर <math>a_1, \ldots, a_p</math>.<ref name=Percival1993>{{cite book |last1=Percival|first1=Donald B.|last2=Walden|first2=Andrew T.|title=भौतिक अनुप्रयोगों के लिए वर्णक्रमीय विश्लेषण|date=1992|publisher=Cambridge University Press|isbn=9780521435413}}</ref> अनुमान समस्या तब इन पैरामीटरों का आकलन करने में से एक बन जाती है। | ||
प्राचलिक एसडीएफ अनुमान का सबसे सामान्य रूप एक | प्राचलिक एसडीएफ अनुमान का सबसे सामान्य रूप एक प्रारूप के रूप में एक ऑटोरेग्रेसिव प्रारूप का उपयोग करता है <math>\text{AR}(p)</math> आदेश की <math>p</math>.{{r|Percival1993|page1=392}} एक संकेत अनुक्रम <math>\{Y_t\}</math> शून्य माध्य का पालन करना <math>\text{AR}(p)</math> प्रक्रिया समीकरण को संतुष्ट करती है | ||
:<math>Y_t = \phi_1Y_{t-1} + \phi_2Y_{t-2} + \cdots + \phi_pY_{t-p} + \epsilon_t,</math> | :<math>Y_t = \phi_1Y_{t-1} + \phi_2Y_{t-2} + \cdots + \phi_pY_{t-p} + \epsilon_t,</math> | ||
Line 60: | Line 60: | ||
मापदंडों का अनुमान लगाने के लिए कई दृष्टिकोण हैं <math>\phi_1, \ldots, \phi_p,\sigma^2_p</math> की <math>\text{AR}(p)</math> प्रक्रिया और इस प्रकार वर्णक्रमीय घनत्व:{{r|Percival1993|page1=452-453}} | मापदंडों का अनुमान लगाने के लिए कई दृष्टिकोण हैं <math>\phi_1, \ldots, \phi_p,\sigma^2_p</math> की <math>\text{AR}(p)</math> प्रक्रिया और इस प्रकार वर्णक्रमीय घनत्व:{{r|Percival1993|page1=452-453}} | ||
* ऑटोरिग्रेसिव | * ऑटोरिग्रेसिव प्रारूप#यूल-वॉकर समीकरण|यूल-वॉकर अनुमानक एक के लिए यूल-वॉकर समीकरणों को पुनरावर्ती रूप से हल करके पाए जाते हैं। <math>\text{AR}(p)</math> प्रक्रिया | ||
* बर्ग अनुमानक यूल-वॉकर समीकरणों को सामान्य न्यूनतम वर्ग समस्या के रूप में मानते हुए पाए जाते हैं। बर्ग अनुमानकों को आमतौर पर यूल-वाकर अनुमानकों से बेहतर माना जाता है।{{r|Percival1993|page1=452}} बर्ग ने इन्हें अधिकतम एन्ट्रॉपी स्पेक्ट्रल अनुमान के साथ संबद्ध किया।<ref name=Burg>Burg, J.P. (1967) "Maximum Entropy Spectral Analysis", ''Proceedings of the 37th Meeting of the Society of Exploration Geophysicists'', Oklahoma City, Oklahoma.</ref> | * बर्ग अनुमानक यूल-वॉकर समीकरणों को सामान्य न्यूनतम वर्ग समस्या के रूप में मानते हुए पाए जाते हैं। बर्ग अनुमानकों को आमतौर पर यूल-वाकर अनुमानकों से बेहतर माना जाता है।{{r|Percival1993|page1=452}} बर्ग ने इन्हें अधिकतम एन्ट्रॉपी स्पेक्ट्रल अनुमान के साथ संबद्ध किया।<ref name=Burg>Burg, J.P. (1967) "Maximum Entropy Spectral Analysis", ''Proceedings of the 37th Meeting of the Society of Exploration Geophysicists'', Oklahoma City, Oklahoma.</ref> | ||
* फॉरवर्ड-बैकवर्ड न्यूनतम-वर्ग अनुमानक इसका इलाज करते हैं <math>\text{AR}(p)</math> प्रतिगमन समस्या के रूप में प्रक्रिया करें और आगे-पीछे विधि का उपयोग करके उस समस्या को हल करें। वे बर्ग अनुमानकों के साथ प्रतिस्पर्धी हैं। | * फॉरवर्ड-बैकवर्ड न्यूनतम-वर्ग अनुमानक इसका इलाज करते हैं <math>\text{AR}(p)</math> प्रतिगमन समस्या के रूप में प्रक्रिया करें और आगे-पीछे विधि का उपयोग करके उस समस्या को हल करें। वे बर्ग अनुमानकों के साथ प्रतिस्पर्धी हैं। | ||
* अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें एक गैर-रैखिक अनुकूलन | * अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें एक गैर-रैखिक अनुकूलन सम्मिलित है और पहले तीन की तुलना में अधिक जटिल है। | ||
वैकल्पिक प्राचलिक विधियों में [[ चलती औसत मॉडल ]] (एमए) और फुल ऑटोरेग्रेसिव मूविंग एवरेज | वैकल्पिक प्राचलिक विधियों में [[ चलती औसत मॉडल | चलती औसत प्रारूप]] (एमए) और फुल ऑटोरेग्रेसिव मूविंग एवरेज प्रारूप (एआरएमए) को फिट करना सम्मिलित है। | ||
== आवृत्ति अनुमान == | == आवृत्ति अनुमान == | ||
आवृत्ति अनुमान, अवयव की संख्या के विषय में दी गई धारणाओं के [[शोर|ध्वनि]] की उपस्थिति में [[ अंकीय संकेत प्रक्रिया |अंकीय संकेत प्रक्रिया]] की आवृत्ति, विपुलता और चरण-परिवर्तन के आकलन सिद्धांत की प्रक्रिया है।<ref>Hayes, Monson H., ''Statistical Digital Signal Processing and Modeling'', John Wiley & Sons, Inc., 1996. {{ISBN|0-471-59431-8}}.</ref> यह उपरोक्त सामान्य विधियों के विपरीत है, जो घटकों के बारे में पूर्व धारणा नहीं बनाते हैं। | |||
=== सिंगल टोन === | === सिंगल टोन === | ||
यदि कोई केवल एक ही सबसे तीव्र आवृत्ति का अनुमान लगाना चाहता है, तो वह [[ पिच डिटेक्शन एल्गोरिदम |ध्वनि पहचान कलन विधि]] का उपयोग कर सकता है। यदि प्रमुख आवृत्ति समय के साथ परिवर्तित है, तो समस्या [[तात्कालिक आवृत्ति]] का अनुमान बन जाती है, जैसा कि समय-आवृत्ति प्रतिनिधित्व में परिभाषित किया गया है। तात्क्षणिक आवृत्ति अनुमान के विधियों में विग्नर-विल वितरण और उच्च क्रम अस्पष्टता कार्यों पर आधारित सम्मिलित हैं।<ref name=Lerga>{{cite web|last=Lerga|first=Jonatan|title=सिग्नल तात्कालिक आवृत्ति अनुमान विधियों का अवलोकन|url=http://www.fer.unizg.hr/_download/repository/Jonatan_Lerga_-_kvalifikacijski_rad.pdf|publisher=University of Rijeka|access-date=22 March 2014}}</ref> | |||
यदि कोई केवल एक ही सबसे | |||
यदि कोई प्राप्त | यदि कोई प्राप्त संकेत (संचारित सिग्नल और ध्वनि सहित) के सभी (संभवतः जटिल) आवृत्ति घटकों को जानना चाहता है, तो एक बहु-टोन दृष्टिकोण का उपयोग करता है। | ||
=== एकाधिक स्वर === | === एकाधिक स्वर === | ||
एक संकेत के लिए एक विशिष्ट | एक संकेत के लिए एक विशिष्ट प्रारूप <math>x(n)</math> का योग होता है <math>p</math> सफेद ध्वनि की उपस्थिति में जटिल घातांक, <math>w(n)</math> | ||
:<math>x(n) = \sum_{i=1}^p A_i e^{j n \omega_i} + w(n)</math>. | :<math>x(n) = \sum_{i=1}^p A_i e^{j n \omega_i} + w(n)</math>. | ||
की शक्ति वर्णक्रमीय घनत्व <math>x(n)</math> से बना है <math>p</math> शोर के कारण वर्णक्रमीय घनत्व समारोह के अलावा आवेग कार्य करता है। | की शक्ति वर्णक्रमीय घनत्व <math>x(n)</math> से बना है <math>p</math> शोर के कारण वर्णक्रमीय घनत्व समारोह के अलावा आवेग कार्य करता है। | ||
आवृत्ति अनुमान के लिए सबसे | आवृत्ति अनुमान के लिए सबसे सरलतम विधियों में इन घटकों को निकालने के लिए ध्वनि रैखिक उप-स्थान की पहचान करना सम्मिलित है। ये विधियाँ एक संकेत उप-स्थान और एक ध्वनि उप-क्षेत्र में स्वतःसंबंध मैट्रिक्स के [[एजकम्पोजीशन]] पर आधारित हैं। इन उप-स्थानों की पहचान के पश्चात, ध्वनि उप-स्थान से घटक आवृत्तियों को खोजने के लिए एक आवृत्ति अनुमान फलन का उपयोग किया जाता है। ध्वनि उप-स्थान आधारित आवृत्ति आकलन के सबसे लोकप्रिय विधि हैं, पिसारेंको हार्मोनिक, पिसारेंको की विधि, एकाधिक संकेत वर्गीकरण (संगीत) विधि, आइजन्वेक्टर विधि और न्यूनतम मानदंड विधि। | ||
; पिसारेंको हार्मोनिक अपघटन | पिसारेंको की विधि: <math>\hat{P}_\text{PHD}\left(e^{j \omega}\right) = \frac{1}{\left|\mathbf{e}^H \mathbf{v}_\text{min}\right|^2}</math> | ; पिसारेंको हार्मोनिक अपघटन | पिसारेंको की विधि: <math>\hat{P}_\text{PHD}\left(e^{j \omega}\right) = \frac{1}{\left|\mathbf{e}^H \mathbf{v}_\text{min}\right|^2}</math> | ||
Line 92: | Line 92: | ||
== उदाहरण गणना == | == उदाहरण गणना == | ||
कल्पना करना <math>x_n</math>, से <math>n=0</math> को <math>N-1</math> शून्य माध्य के साथ एक समय श्रृंखला (असतत समय) है। मान लीजिए कि यह आवधिक घटकों की एक परिमित संख्या का योग है (सभी आवृत्तियाँ धनात्मक हैं): | कल्पना करना <math>x_n</math>, से <math>n=0</math> को <math>N-1</math> शून्य माध्य के साथ एक समय श्रृंखला (असतत समय) है। मान लीजिए कि यह आवधिक घटकों की एक परिमित संख्या का योग है, (सभी आवृत्तियाँ धनात्मक हैं): | ||
:<math>\begin{align} | :<math>\begin{align} | ||
Line 102: | Line 102: | ||
:<math>\frac{1}{N} \sum_{n=0}^{N-1} x_n^2.</math> | :<math>\frac{1}{N} \sum_{n=0}^{N-1} x_n^2.</math> | ||
यदि ये | यदि ये आँकड़े एक विद्युत संकेत से लिए गए नमूने थे, तो यह इसकी औसत शक्ति होगी (शक्ति प्रति इकाई समय में ऊर्जा है, इसलिए यह विचरण के अनुरूप है यदि ऊर्जा आयाम वर्ग के अनुरूप है)। | ||
अब, | अब, सरलता के लिए, मान लें कि सिग्नल समय में अनंतता तक फैलता है, इसलिए हम इसे सीमा के रूप में ले जाते हैं <math>N\to \infty.</math> यदि औसत शक्ति परिबद्ध है, जो वास्तविकता में लगभग सदैव ही होती है, तो निम्नलिखित सीमा उपस्थित होती है, और यह आँकड़े का विचरण है। | ||
:<math>\lim_{N \to \infty} \frac{1}{N} \sum_{n=0}^{N-1} x_n^2.</math> | :<math>\lim_{N \to \infty} \frac{1}{N} \sum_{n=0}^{N-1} x_n^2.</math> | ||
तत्पश्चात, सरलता के लिए, हम निरंतर समय को पास करेंगे, और यह मानेंगे कि संकेत दोनों दिशाओं में समय में असीमित रूप से विस्तारित होता है। तब ये दो सूत्र बन जाते हैं | |||
:<math>x(t) = \sum_k A_k \sin(2\pi\nu_k t + \phi_k)</math> | :<math>x(t) = \sum_k A_k \sin(2\pi\nu_k t + \phi_k)</math> | ||
Line 113: | Line 113: | ||
:<math>\lim_{T\to\infty} \frac{1}{2T} \int_{-T}^T x(t)^2 dt.</math> | :<math>\lim_{T\to\infty} \frac{1}{2T} \int_{-T}^T x(t)^2 dt.</math> | ||
का मूल माध्य वर्ग <math>\sin</math> है <math>1/\sqrt{2}</math>, इसलिए का विचरण <math>A_k \sin(2\pi\nu_k t + \phi_k)</math> है <math>\tfrac{1}{2} A_k^2.</math> | का मूल माध्य वर्ग <math>\sin</math> है <math>1/\sqrt{2}</math>, इसलिए का विचरण <math>A_k \sin(2\pi\nu_k t + \phi_k)</math> है <math>\tfrac{1}{2} A_k^2.</math>, की औसत शक्ति में योगदान <math>x(t)</math> आवृत्ति के साथ घटक से आ रहा है <math>\nu_k</math> है <math>\tfrac{1}{2}A_k^2.</math> ये सभी योगदान की औसत शक्ति को जोड़ते हैं <math>x(t).</math> | ||
इसके पश्चात आवृत्ति के कार्य के रूप में शक्ति है <math>\tfrac{1}{2}A_k^2,</math> और इसका सांख्यिकीय [[संचयी वितरण कार्य]] <math>S(\nu)</math> होगा | |||
:<math>S(\nu) = \sum _ {k : \nu_k < \nu} \frac{1}{2} A_k^2.</math> | :<math>S(\nu) = \sum _ {k : \nu_k < \nu} \frac{1}{2} A_k^2.</math> | ||
<math>S</math> एक [[ समारोह की ओर कदम बढ़ाएं ]] है, नीरस रूप से गैर-घटता है। इसकी छलांग अवधि (रिंग) के घटकों की आवृत्तियों पर होती है <math>x</math>, और प्रत्येक छलांग का मान उस घटक की शक्ति या भिन्नता है। | <math>S</math> एक [[ समारोह की ओर कदम बढ़ाएं | समारोह की ओर कदम बढ़ाएं]] है, नीरस रूप से गैर-घटता है। इसकी छलांग अवधि (रिंग) के घटकों की आवृत्तियों पर होती है <math>x</math>, और प्रत्येक छलांग का मान उस घटक की शक्ति या भिन्नता है। | ||
प्रसरण स्वयं के साथ | प्रसरण स्वयं के साथ आँकड़े का सहप्रसरण है। यदि हम अब समान आँकड़े पर विचार करें, परन्तु एक अंतराल के साथ <math>\tau</math>,को हम सहप्रसरण ले सकते हैं, <math>x(t)</math> साथ <math>x(t + \tau)</math>, और इसे स्वतः सहसंबंध फलन के रूप में परिभाषित करें <math>c</math> संकेत (या डेटा) के <math>x</math>: | ||
:<math>c(\tau) = \lim_{T\to\infty} \frac{1}{2T} \int_{-T}^T x(t) x(t + \tau) dt.</math> | :<math>c(\tau) = \lim_{T\to\infty} \frac{1}{2T} \int_{-T}^T x(t) x(t + \tau) dt.</math> | ||
यदि यह अस्तित्व में है, तो यह | यदि यह अस्तित्व में है, तो यह <math>\tau.</math> का एक कार्य भी है, यदि औसत शक्ति परिबद्ध है, तब <math>c</math> हर जगह उपस्थित है एवं <math>c(0)</math> परिमित है, और इससे घिरा हुआ है जो आँकड़े की औसत शक्ति या भिन्नता है। | ||
यह दिखाया जा सकता है <math>c</math> के रूप में एक ही अवधि के साथ आवधिक घटकों में विघटित किया जा सकता है | यह दिखाया जा सकता है,कि <math>c</math> के रूप में <math>x</math>: एक ही अवधि के साथ आवधिक घटकों में विघटित किया जा सकता है: | ||
:<math>c(\tau) = \sum_k \frac{1}{2} A_k^2 \cos(2\pi\nu_k\tau).</math> | :<math>c(\tau) = \sum_k \frac{1}{2} A_k^2 \cos(2\pi\nu_k\tau).</math> | ||
यह वास्तव में वर्णक्रमीय अपघटन है, <math>c</math> विभिन्न आवृत्तियों पर | यह वास्तव में वर्णक्रमीय अपघटन है, और <math>c</math> विभिन्न आवृत्तियों पर शक्ति के वितरण से संबंधित है, <math>x</math> आवृत्तियों पर: के एक आवृत्ति घटक का आयाम <math>c</math> सिग्नल की औसत शक्ति में इसका योगदान होता है। | ||
इस उदाहरण | इस उदाहरण की शक्ति वर्णक्रम निरंतर नहीं है, इसलिए इसका अवकलन नहीं होता है, और इसलिए इस संकेत का शक्ति स्पेक्ट्रल घनत्व फलन नहीं होता है। सामान्य रूप से, शक्ति वर्णक्रम दो हिस्सों का योग होता है: एक ऐसी रेखा वर्णक्रम जैसा इस उदाहरण में है, जो निरंतर नहीं है, और जिसका घनत्व फलन नहीं होता है, और एक शेष, जो पूरी तरह से निरंतर होता है,एवं जिसका घनत्व फलन होता है। | ||
== यह भी देखें == | == यह भी देखें == |
Revision as of 22:10, 28 June 2023
सांख्यिकीय संकेत प्रसंस्करण में, वर्णक्रमीय घनत्व अनुमान (एसडीई) या केवल वर्णक्रमीय अनुमान का लक्ष्य सिग्नल के समय नमूनों के अनुक्रम से सिग्नल के वर्णक्रमीय घनत्व (जिसे पावर स्पेक्ट्रल घनत्व के रूप में भी जाना जाता है) का अनुमान लगाना है।[1] सहज रूप से कहें तो, वर्णक्रमीय घनत्व सिग्नल की आवृत्ति सामग्री को दर्शाता है। वर्णक्रमीय घनत्व का अनुमान लगाने का एक उद्देश्य इन आवधिकों के अनुरूप आवृत्तियों पर चोटियों को देखकर, आँकड़े में किसी भी आवधिकता का पता लगाना है।
कुछ एसडीई तकनीकों का मानना है, कि एक संकेत उत्पन्न आवृत्तियों की एक सीमित (सामान्यतः छोटी) संख्या और ध्वनि से बना होता है, और उत्पन्न आवृत्तियों के स्थान और तीव्रता को खोजने का प्रयास करता है। अन्य घटकों की संख्या पर कोई धारणा नहीं बनाते हैं, और संपूर्ण उत्पादक उत्पन्न करने का अनुमान लगाना चाहते हैं।
अवलोकन
वर्णक्रम विश्लेषण, जिसे वर्णक्रमीय अनुमान विश्लेषण या वर्णक्रमीय घनत्व अनुमान के रूप में भी जाना जाता है, एक जटिल संकेत को सरल भागों में विघटित करने की तकनीकी प्रक्रिया है। जैसा कि ऊपर बताया गया है, कई भौतिक प्रक्रियाओं को कई अलग-अलग आवृत्ति घटकों के योग के रूप में वर्णित किया गया है। कोई भी प्रक्रिया जो विभिन्न राशियों (जैसे आयाम, शक्तियाँ, तीव्रता) बनाम आवृत्ति (या चरण (तरंगें)) की मात्रा निर्धारित करती है, उसे वर्णक्रम विश्लेषण कहा जा सकता है।
वर्णक्रम विश्लेषण पूरे संकेत पर किया जा सकता है। वैकल्पिक रूप से, एक संकेत को छोटे खंडों में तोड़ा जा सकता है (कभी-कभी फ्रेम कहा जाता है), और इन अलग-अलग खंडों पर वर्णक्रम विश्लेषण लागू किया जा सकता है। आवधिक कार्य (जैसे ) इस उप-विभाजन के लिए विशेष रूप से उपयुक्त हैं। गैर-आवधिक कार्यों के विश्लेषण के लिए सामान्य गणितीय विधियाँ सांध्वनिक विश्लेषण की श्रेणी में आती हैं।
फलन का सांध्वनिक रूपांतरण एक आवृत्ति वर्णक्रम उत्पन्न करता है, जिसमें मूल संकेत के विषय में सभी जानकारी होती है, परन्तु एक अलग रूप में। इसका तात्पर्य यह है, कि प्रतिलोम संध्वनिक रूपांतरण द्वारा मूल कार्य को पूरी तरह से पुनर्निर्मित (संश्लेषित) किया जा सकता है। सही पुनर्निर्माण के लिए, वर्णक्रम विश्लेषक को प्रत्येक आवृत्ति घटक के आयाम और चरण (तरंगों) दोनों को संरक्षित करना चाहिए। जानकारी के इन दो टुकड़ों को 2-आयामी वेक्टर के रूप में, एक जटिल संख्या के रूप में, या परिमाण (आयाम) और चरण के रूप में ध्रुवीय निर्देशांक (अर्थात, एक चरण के रूप में) के रूप में दर्शाया जा सकता है। संकेत प्रसंस्करण में एक सामान्य विधि वर्ग आयाम, या शक्ति (भौतिकी) पर विचार करना है; इस विषय में परिणामी भूखंड को शक्ति वर्णक्रम कहा जाता है।
प्रतिवर्तीता के कारण, सान्ध्वनिक रूपांतरण को समय के अतिरिक्त आवृत्ति के संदर्भ में फलन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह एक आवृत्ति कार्यक्षेत्र प्रतिनिधित्व है। रैखिक संचालन जो समय कार्यक्षेत्र में किए जा सकते हैं, उनके समकक्ष हैं, जो प्रायः आवृत्ति कार्यक्षेत्र में अधिक सरलता से किए जा सकते हैं। बारंबारता विश्लेषण भी रैखिक और गैर-रैखिक दोनों तरह के विभिन्न समय -कार्यक्षेत्र संचालनों के प्रभावों की समझ और व्याख्या को सरल करता है। उदाहरण के लिए, केवल गैर-रैखिक या समय-भिन्न प्रणाली आवृत्ति वर्णक्रम में नई आवृत्तियां निर्मित कर सकते हैं।
व्यवहार में, लगभग सभी सॉफ्टवेयर और वैद्युत्कीय उपकरण जो आवृत्ति स्पेक्ट्रा उत्पन्न करते हैं, एक विशेष फूरियर परिवर्तन (डीएफटी) का उपयोग करते हैं, जो संकेत के नमूनों पर कार्य करता है, और पूर्ण समाकलित समाधान के लिए एक गणितीय अनुमान प्रदान करता है। जो संकेत के नमूनाकरण पर संचालित होता है, और जो पूर्ण अभिन्न समाधान के लिए गणितीय समीपता प्रदान करता है। डीएफटी लगभग अनिवार्य रूप से [[असतत फूरियर रूपांतरण]] (एफएफटी) नामक एक कुशल कलन विधि द्वारा कार्यान्वित किया जाता है। एक डीएफटी के वर्ग-परिमाण घटकों की सारणी एक प्रकार का शक्ति वर्णक्रम है, जिसे पीरियोग्राम कहा जाता है, जिसका व्यापक रूप से आवेग प्रतिक्रिया और विंडो फलन जैसे ध्वनि-मुक्त कार्यों की आवृत्ति विशेषताओं की जांच के लिए उपयोग किया जाता है। परन्तु कम सिग्नल-से-ध्वनि अनुपात पर ध्वनि जैसा संकेत या यहां तक कि ज्यावक्र पर लागू होने पर पीरियडोग्राम प्रसंस्करण-लाभ प्रदान नहीं करता है। दूसरे शब्दों में, किसी दी गई आवृत्ति पर इसके वर्णक्रमीय अनुमान का प्रसरण कम नहीं होता है, क्योंकि संगणना में उपयोग किए जाने वाले नमूनों की संख्या बढ़ जाती है। इसे समय के साथ औसत करके कम किया जा सकता है, (वेल्च की विधि[1]) या अधिक आवृत्ति (चौरसाई)। वेल्च की विधि व्यापक रूप से वर्णक्रमीय घनत्व अनुमान (एसडीई) के लिए उपयोग की जाती है। यद्यपि, पीरियोग्राम-आधारित तकनीकें छोटे पूर्वाग्रहों का परिचय देती हैं जो कुछ अनुप्रयोगों में अस्वीकार्य हैं। इसलिए अन्य विकल्प अगले भाग में प्रस्तुत किए गए हैं।
तकनीक
आधारभूत पीरियडोग्राम के नुकसान को कम करने के लिए स्पेक्ट्रल आकलन के लिए कई अन्य तकनीकों का विकास किया गया है। इन तकनीकों को आम तौर पर गैर-आवर्तिता वक्र आंकड़ों में विभाजित किया जा सकता है।[2] गैर- प्राचलिक दृष्टिकोण स्पष्ट रूप से सहप्रसरण या प्रक्रिया के वर्णक्रम का अनुमान लगाए बिना यह मानते हैं कि प्रक्रिया की कोई विशेष संरचना है। बुनियादी अनुप्रयोगों (जैसे वेल्च की विधि) के लिए उपयोग किए जाने वाले कुछ सबसे आम अनुमानक गैर- प्राचलिक अनुमानक हैं जो पीरियोडोग्राम से निकटता से संबंधित हैं। इसके विपरीत, प्राचलिक दृष्टिकोण मानते हैं कि अंतर्निहित स्थिर प्रक्रिया में एक निश्चित संरचना होती है जिसे पैरामीटर की एक छोटी संख्या का उपयोग करके वर्णित किया जा सकता है (उदाहरण के लिए, ऑटोरेग्रेसिव मूविंग एवरेज प्रारूप | ऑटो-रिग्रेसिव या मूविंग एवरेज प्रारूप का उपयोग करके)। इन दृष्टिकोणों में, कार्य प्रारूप के मापदंडों का अनुमान लगाना है जो स्टोकेस्टिक प्रक्रिया का वर्णन करता है। अर्ध- प्राचलिक विधियों का उपयोग करते समय, अंतर्निहित प्रक्रिया को एक गैर- प्राचलिक ढांचे का उपयोग करके तैयार किया जाता है, अतिरिक्त धारणा के साथ कि प्रारूप के गैर-शून्य घटकों की संख्या छोटी है (यानी, प्रारूप विरल है)। लापता डेटा रिकवरी के लिए भी इसी तरह के तरीकों का इस्तेमाल किया जा सकता है [3] साथ ही संपीड़ित संवेदन।
निम्नलिखित गैर- प्राचलिक वर्णक्रमीय घनत्व अनुमान तकनीकों की आंशिक सूची है:
- पीरियोडोग्राम, असतत फूरियर रूपांतरण का मापांक वर्ग
- लोम्ब-स्कार्ल पीरियोग्राम, जिसके लिए डेटा को समान रूप से स्थान देने की आवश्यकता नहीं है
- बार्टलेट की विधि वर्णक्रमीय घनत्व अनुमान के विचरण को कम करने के लिए सिग्नल के कई खंडों से लिए गए पीरियडोग्राम का औसत है
- वेल्च की विधि बार्टलेट की विधि का एक खिड़की वाला संस्करण है जो ओवरलैपिंग सेगमेंट का उपयोग करता है
- मल्टीटेपर एक पीरियडोग्राम-आधारित विधि है जो वर्णक्रमीय घनत्व अनुमान के विचरण को कम करने के लिए वर्णक्रमीय घनत्व के स्वतंत्र अनुमान बनाने के लिए कई टेपर्स या विंडो का उपयोग करती है।
- कम से कम वर्ग वर्णक्रमीय विश्लेषण, ज्ञात आवृत्तियों के लिए कम से कम वर्गों के आधार पर
- गैर-समान असतत फूरियर रूपांतरण का उपयोग तब किया जाता है जब सिग्नल नमूने असमान रूप से समय श्रृंखला में होते हैं
- एकवचन वर्णक्रम विश्लेषण एक गैर प्राचलिक विधि है जो वर्णक्रमीय घनत्व का अनुमान लगाने के लिए सहप्रसरण मैट्रिक्स के एकवचन मूल्य अपघटन का उपयोग करती है
- शॉर्ट-टाइम फूरियर रूपांतरण
- सूचना क्षेत्र सिद्धांत#महत्वपूर्ण फ़िल्टर सूचना क्षेत्र सिद्धांत पर आधारित एक गैर- प्राचलिक विधि है जो शोर, अपूर्ण डेटा और वाद्य प्रतिक्रिया कार्यों से निपट सकता है
नीचे प्राचलिक तकनीकों की आंशिक सूची दी गई है:
- Autoregressive model (AR) आकलन, जो मानता है कि nth नमूना पिछले p नमूनों के साथ सहसंबद्ध है।
- मूविंग-एवरेज प्रारूप (MA) का अनुमान, जो मानता है कि nth नमूना पिछले p नमूनों में शोर की शर्तों के साथ सहसंबद्ध है।
- Autoregressive मूविंग एवरेज (ARMA) अनुमान, जो AR और MA प्रारूप का सामान्यीकरण करता है।
- संगीत (एल्गोरिदम) (संगीत) एक लोकप्रिय सुपर-रिज़ॉल्यूशन इमेजिंग विधि है।
- अधिकतम एन्ट्रॉपी वर्णक्रमीय अनुमान एसडीई के लिए उपयोगी एक सर्व-ध्रुवीय विधि है जब एकवचन वर्णक्रमीय विशेषताओं, जैसे तेज चोटियों की अपेक्षा की जाती है।
और अंत में अर्ध- प्राचलिक तकनीकों के कुछ उदाहरण:
- स्पार्स इटरेटिव कोवैरियंस-आधारित अनुमान (स्पाइस) अनुमान,[2]और अधिक सामान्यीकृत -मसाला।[4] * पुनरावृत्त अनुकूली दृष्टिकोण (आईएए) अनुमान।[5] *लास्सो (सांख्यिकी), कम से कम वर्ग वर्णक्रमीय विश्लेषण के समान परन्तु एक विरलता लागू करने वाले दंड के साथ।[6]
प्राचलिक अनुमान
प्राचलिक स्पेक्ट्रल अनुमान में, कोई मानता है कि सिग्नल एक स्थिर प्रक्रिया द्वारा तैयार किया गया है जिसमें स्पेक्ट्रल घनत्व फ़ंक्शन (एसडीएफ) है। यह आवृत्ति का एक कार्य है और पैरामीटर .[7] अनुमान समस्या तब इन पैरामीटरों का आकलन करने में से एक बन जाती है।
प्राचलिक एसडीएफ अनुमान का सबसे सामान्य रूप एक प्रारूप के रूप में एक ऑटोरेग्रेसिव प्रारूप का उपयोग करता है आदेश की .[7]: 392 एक संकेत अनुक्रम शून्य माध्य का पालन करना प्रक्रिया समीकरण को संतुष्ट करती है
जहां निश्चित गुणांक हैं और शून्य माध्य और नवीनता विचरण के साथ एक श्वेत शोर प्रक्रिया है . इस प्रक्रिया के लिए एसडीएफ है
साथ नमूना समय अंतराल और Nyquist आवृत्ति।
मापदंडों का अनुमान लगाने के लिए कई दृष्टिकोण हैं की प्रक्रिया और इस प्रकार वर्णक्रमीय घनत्व:[7]: 452-453
- ऑटोरिग्रेसिव प्रारूप#यूल-वॉकर समीकरण|यूल-वॉकर अनुमानक एक के लिए यूल-वॉकर समीकरणों को पुनरावर्ती रूप से हल करके पाए जाते हैं। प्रक्रिया
- बर्ग अनुमानक यूल-वॉकर समीकरणों को सामान्य न्यूनतम वर्ग समस्या के रूप में मानते हुए पाए जाते हैं। बर्ग अनुमानकों को आमतौर पर यूल-वाकर अनुमानकों से बेहतर माना जाता है।[7]: 452 बर्ग ने इन्हें अधिकतम एन्ट्रॉपी स्पेक्ट्रल अनुमान के साथ संबद्ध किया।[8]
- फॉरवर्ड-बैकवर्ड न्यूनतम-वर्ग अनुमानक इसका इलाज करते हैं प्रतिगमन समस्या के रूप में प्रक्रिया करें और आगे-पीछे विधि का उपयोग करके उस समस्या को हल करें। वे बर्ग अनुमानकों के साथ प्रतिस्पर्धी हैं।
- अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें एक गैर-रैखिक अनुकूलन सम्मिलित है और पहले तीन की तुलना में अधिक जटिल है।
वैकल्पिक प्राचलिक विधियों में चलती औसत प्रारूप (एमए) और फुल ऑटोरेग्रेसिव मूविंग एवरेज प्रारूप (एआरएमए) को फिट करना सम्मिलित है।
आवृत्ति अनुमान
आवृत्ति अनुमान, अवयव की संख्या के विषय में दी गई धारणाओं के ध्वनि की उपस्थिति में अंकीय संकेत प्रक्रिया की आवृत्ति, विपुलता और चरण-परिवर्तन के आकलन सिद्धांत की प्रक्रिया है।[9] यह उपरोक्त सामान्य विधियों के विपरीत है, जो घटकों के बारे में पूर्व धारणा नहीं बनाते हैं।
सिंगल टोन
यदि कोई केवल एक ही सबसे तीव्र आवृत्ति का अनुमान लगाना चाहता है, तो वह ध्वनि पहचान कलन विधि का उपयोग कर सकता है। यदि प्रमुख आवृत्ति समय के साथ परिवर्तित है, तो समस्या तात्कालिक आवृत्ति का अनुमान बन जाती है, जैसा कि समय-आवृत्ति प्रतिनिधित्व में परिभाषित किया गया है। तात्क्षणिक आवृत्ति अनुमान के विधियों में विग्नर-विल वितरण और उच्च क्रम अस्पष्टता कार्यों पर आधारित सम्मिलित हैं।[10]
यदि कोई प्राप्त संकेत (संचारित सिग्नल और ध्वनि सहित) के सभी (संभवतः जटिल) आवृत्ति घटकों को जानना चाहता है, तो एक बहु-टोन दृष्टिकोण का उपयोग करता है।
एकाधिक स्वर
एक संकेत के लिए एक विशिष्ट प्रारूप का योग होता है सफेद ध्वनि की उपस्थिति में जटिल घातांक,
- .
की शक्ति वर्णक्रमीय घनत्व से बना है शोर के कारण वर्णक्रमीय घनत्व समारोह के अलावा आवेग कार्य करता है।
आवृत्ति अनुमान के लिए सबसे सरलतम विधियों में इन घटकों को निकालने के लिए ध्वनि रैखिक उप-स्थान की पहचान करना सम्मिलित है। ये विधियाँ एक संकेत उप-स्थान और एक ध्वनि उप-क्षेत्र में स्वतःसंबंध मैट्रिक्स के एजकम्पोजीशन पर आधारित हैं। इन उप-स्थानों की पहचान के पश्चात, ध्वनि उप-स्थान से घटक आवृत्तियों को खोजने के लिए एक आवृत्ति अनुमान फलन का उपयोग किया जाता है। ध्वनि उप-स्थान आधारित आवृत्ति आकलन के सबसे लोकप्रिय विधि हैं, पिसारेंको हार्मोनिक, पिसारेंको की विधि, एकाधिक संकेत वर्गीकरण (संगीत) विधि, आइजन्वेक्टर विधि और न्यूनतम मानदंड विधि।
- पिसारेंको हार्मोनिक अपघटन | पिसारेंको की विधि
- एकाधिक संकेत वर्गीकरण
- ,
- ईजेनवेक्टर विधि
- न्यूनतम मानदंड विधि
उदाहरण गणना
कल्पना करना , से को शून्य माध्य के साथ एक समय श्रृंखला (असतत समय) है। मान लीजिए कि यह आवधिक घटकों की एक परिमित संख्या का योग है, (सभी आवृत्तियाँ धनात्मक हैं):
का विचरण ऊपर दिए गए शून्य-माध्य फ़ंक्शन के लिए है
यदि ये आँकड़े एक विद्युत संकेत से लिए गए नमूने थे, तो यह इसकी औसत शक्ति होगी (शक्ति प्रति इकाई समय में ऊर्जा है, इसलिए यह विचरण के अनुरूप है यदि ऊर्जा आयाम वर्ग के अनुरूप है)।
अब, सरलता के लिए, मान लें कि सिग्नल समय में अनंतता तक फैलता है, इसलिए हम इसे सीमा के रूप में ले जाते हैं यदि औसत शक्ति परिबद्ध है, जो वास्तविकता में लगभग सदैव ही होती है, तो निम्नलिखित सीमा उपस्थित होती है, और यह आँकड़े का विचरण है।
तत्पश्चात, सरलता के लिए, हम निरंतर समय को पास करेंगे, और यह मानेंगे कि संकेत दोनों दिशाओं में समय में असीमित रूप से विस्तारित होता है। तब ये दो सूत्र बन जाते हैं
और
का मूल माध्य वर्ग है , इसलिए का विचरण है , की औसत शक्ति में योगदान आवृत्ति के साथ घटक से आ रहा है है ये सभी योगदान की औसत शक्ति को जोड़ते हैं
इसके पश्चात आवृत्ति के कार्य के रूप में शक्ति है और इसका सांख्यिकीय संचयी वितरण कार्य होगा
एक समारोह की ओर कदम बढ़ाएं है, नीरस रूप से गैर-घटता है। इसकी छलांग अवधि (रिंग) के घटकों की आवृत्तियों पर होती है , और प्रत्येक छलांग का मान उस घटक की शक्ति या भिन्नता है।
प्रसरण स्वयं के साथ आँकड़े का सहप्रसरण है। यदि हम अब समान आँकड़े पर विचार करें, परन्तु एक अंतराल के साथ ,को हम सहप्रसरण ले सकते हैं, साथ , और इसे स्वतः सहसंबंध फलन के रूप में परिभाषित करें संकेत (या डेटा) के :
यदि यह अस्तित्व में है, तो यह का एक कार्य भी है, यदि औसत शक्ति परिबद्ध है, तब हर जगह उपस्थित है एवं परिमित है, और इससे घिरा हुआ है जो आँकड़े की औसत शक्ति या भिन्नता है।
यह दिखाया जा सकता है,कि के रूप में : एक ही अवधि के साथ आवधिक घटकों में विघटित किया जा सकता है:
यह वास्तव में वर्णक्रमीय अपघटन है, और विभिन्न आवृत्तियों पर शक्ति के वितरण से संबंधित है, आवृत्तियों पर: के एक आवृत्ति घटक का आयाम सिग्नल की औसत शक्ति में इसका योगदान होता है।
इस उदाहरण की शक्ति वर्णक्रम निरंतर नहीं है, इसलिए इसका अवकलन नहीं होता है, और इसलिए इस संकेत का शक्ति स्पेक्ट्रल घनत्व फलन नहीं होता है। सामान्य रूप से, शक्ति वर्णक्रम दो हिस्सों का योग होता है: एक ऐसी रेखा वर्णक्रम जैसा इस उदाहरण में है, जो निरंतर नहीं है, और जिसका घनत्व फलन नहीं होता है, और एक शेष, जो पूरी तरह से निरंतर होता है,एवं जिसका घनत्व फलन होता है।
यह भी देखें
- बहुआयामी वर्णक्रमीय अनुमान
- पीरियोडोग्राम
- सिगस्पेक
- spectrogram
- समय-आवृत्ति विश्लेषण
- समय-आवृत्ति प्रतिनिधित्व
- कम संभावना
- वर्णक्रमीय बिजली वितरण
संदर्भ
- ↑ Welch, P. D. (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
- ↑ 2.0 2.1 Stoica, Petre; Babu, Prabhu; Li, Jian (January 2011). "अलग-अलग मॉडल में विरल पैरामीटर अनुमान की नई विधि और अनियमित रूप से सैंपल किए गए डेटा के स्पेक्ट्रल विश्लेषण के लिए इसका उपयोग". IEEE Transactions on Signal Processing. 59 (1): 35–47. Bibcode:2011ITSP...59...35S. doi:10.1109/TSP.2010.2086452. ISSN 1053-587X. S2CID 15936187.
- ↑ Stoica, Petre; Li, Jian; Ling, Jun; Cheng, Yubo (April 2009). "एक गैर पैरामीट्रिक पुनरावृत्त अनुकूली दृष्टिकोण के माध्यम से लापता डेटा रिकवरी". 2009 IEEE International Conference on Acoustics, Speech and Signal Processing. IEEE: 3369–3372. doi:10.1109/icassp.2009.4960347. ISBN 978-1-4244-2353-8.
- ↑ Sward, Johan; Adalbjornsson, Stefan Ingi; Jakobsson, Andreas (March 2017). "विरल पुनरावृत्त सहप्रसरण-आधारित अनुमानक का एक सामान्यीकरण". 2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP). IEEE: 3954–3958. doi:10.1109/icassp.2017.7952898. ISBN 978-1-5090-4117-6. S2CID 5640068.
- ↑ Yardibi, Tarik; Li, Jian; Stoica, Petre; Xue, Ming; Baggeroer, Arthur B. (January 2010). "Source Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares". IEEE Transactions on Aerospace and Electronic Systems. 46 (1): 425–443. Bibcode:2010ITAES..46..425Y. doi:10.1109/TAES.2010.5417172. hdl:1721.1/59588. ISSN 0018-9251. S2CID 18834345.
- ↑ Panahi, Ashkan; Viberg, Mats (February 2011). "LASSO-आधारित DOA आकलन पद्धति के संकल्प पर". 2011 International ITG Workshop on Smart Antennas. IEEE: 1–5. doi:10.1109/wsa.2011.5741938. ISBN 978-1-61284-075-8. S2CID 7013162.
- ↑ 7.0 7.1 7.2 7.3 Percival, Donald B.; Walden, Andrew T. (1992). भौतिक अनुप्रयोगों के लिए वर्णक्रमीय विश्लेषण. Cambridge University Press. ISBN 9780521435413.
- ↑ Burg, J.P. (1967) "Maximum Entropy Spectral Analysis", Proceedings of the 37th Meeting of the Society of Exploration Geophysicists, Oklahoma City, Oklahoma.
- ↑ Hayes, Monson H., Statistical Digital Signal Processing and Modeling, John Wiley & Sons, Inc., 1996. ISBN 0-471-59431-8.
- ↑ Lerga, Jonatan. "सिग्नल तात्कालिक आवृत्ति अनुमान विधियों का अवलोकन" (PDF). University of Rijeka. Retrieved 22 March 2014.
अग्रिम पठन
- Porat, B. (1994). Digital Processing of Random Signals: Theory & Methods. Prentice Hall. ISBN 978-0-13-063751-2.
- Priestley, M.B. (1991). Spectral Analysis and Time Series. Academic Press. ISBN 978-0-12-564922-3.
- Stoica, P.; Moses, R. (2005). Spectral Analysis of Signals. Prentice Hall. ISBN 978-0-13-113956-5.
- Thomson, D. J. (1982). "Spectrum estimation and harmonic analysis". Proceedings of the IEEE. 70 (9): 1055–1096. Bibcode:1982IEEEP..70.1055T. CiteSeerX 10.1.1.471.1278. doi:10.1109/PROC.1982.12433. S2CID 290772.