वर्णक्रमीय घनत्व अनुमान: Difference between revisions
(Created page with "{{Short description|Signal processing technique}} {{for|the statistical method|Probability density estimation}} {{broader|Spectral density}} सांख्यिकीय...") |
No edit summary |
||
| Line 2: | Line 2: | ||
{{for|the statistical method|Probability density estimation}} | {{for|the statistical method|Probability density estimation}} | ||
{{broader|Spectral density}} | {{broader|Spectral density}} | ||
[[सांख्यिकीय सिग्नल प्रोसेसिंग]] में, [[वर्णक्रमीय घनत्व]] अनुमान (एसडीई) या केवल वर्णक्रमीय अनुमान का लक्ष्य सिग्नल के समय नमूनों के अनुक्रम से सिग्नल के वर्णक्रमीय घनत्व (जिसे [[पावर स्पेक्ट्रम]] के रूप में भी जाना जाता है) का अनुमान लगाना है।<ref>[[Peter Stoica|P Stoica]] and R Moses, Spectral Analysis of Signals, Prentice Hall, 2005.</ref> सहज रूप से कहें तो, वर्णक्रमीय घनत्व सिग्नल की [[आवृत्ति]] सामग्री को दर्शाता है। वर्णक्रमीय घनत्व का अनुमान लगाने का | [[सांख्यिकीय सिग्नल प्रोसेसिंग]] में, [[वर्णक्रमीय घनत्व]] अनुमान (एसडीई) या केवल वर्णक्रमीय अनुमान का लक्ष्य सिग्नल के समय नमूनों के अनुक्रम से सिग्नल के वर्णक्रमीय घनत्व (जिसे [[पावर स्पेक्ट्रम]] के रूप में भी जाना जाता है) का अनुमान लगाना है।<ref>[[Peter Stoica|P Stoica]] and R Moses, Spectral Analysis of Signals, Prentice Hall, 2005.</ref> सहज रूप से कहें तो, वर्णक्रमीय घनत्व सिग्नल की [[आवृत्ति]] सामग्री को दर्शाता है। वर्णक्रमीय घनत्व का अनुमान लगाने का उद्देश्य इन आवधिकों के अनुरूप आवृत्तियों पर चोटियों को देखकर, डेटा में किसी भी आवधिक फ़ंक्शन का पता लगाना है। | ||
कुछ एसडीई तकनीकें मानती हैं कि | कुछ एसडीई तकनीकें मानती हैं कि सिग्नल सीमित (आमतौर पर छोटी) संख्या में उत्पन्न आवृत्तियों और शोर से बना होता है और उत्पन्न आवृत्तियों के स्थान और तीव्रता का पता लगाने की कोशिश करता है। अन्य लोग घटकों की संख्या पर कोई धारणा नहीं बनाते हैं और संपूर्ण उत्पादन स्पेक्ट्रम का अनुमान लगाना चाहते हैं। | ||
== सिंहावलोकन == | == सिंहावलोकन == | ||
| Line 10: | Line 10: | ||
[[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>) इस उप-विभाजन के लिए विशेष रूप से उपयुक्त हैं। गैर-आवधिक कार्यों के विश्लेषण के लिए सामान्य गणितीय तकनीकें [[फूरियर विश्लेषण]] की श्रेणी में आती हैं। | ||
किसी फ़ंक्शन का [[फूरियर रूपांतरण]] | किसी फ़ंक्शन का [[फूरियर रूपांतरण]] आवृत्ति स्पेक्ट्रम उत्पन्न करता है जिसमें मूल सिग्नल के बारे में सारी जानकारी होती है, लेकिन अलग रूप में। इसका मतलब यह है कि मूल फ़ंक्शन को व्युत्क्रम फूरियर रूपांतरण द्वारा पूरी तरह से पुनर्निर्मित (संश्लेषित) किया जा सकता है। सही पुनर्निर्माण के लिए, स्पेक्ट्रम विश्लेषक को प्रत्येक आवृत्ति घटक के [[आयाम]] और [[चरण]] (तरंगों) दोनों को संरक्षित करना होगा। जानकारी के इन दो टुकड़ों को 2-आयामी वेक्टर के रूप में, [[जटिल संख्या]] के रूप में, या ध्रुवीय निर्देशांक में परिमाण (आयाम) और चरण के रूप में (यानी, चरण के रूप में) दर्शाया जा सकता [[नमूनाकरण (सिग्नल प्रोसेसिंग)]] में सामान्य तकनीक वर्ग आयाम, या [[शक्ति (भौतिकी)]] पर विचार करना है; इस मामले में परिणामी प्लॉट को पावर स्पेक्ट्रम के रूप में जाना जाता है। | ||
उत्क्रमणीयता के कारण, फूरियर रूपांतरण को समय के बजाय आवृत्ति के संदर्भ में फ़ंक्शन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह | उत्क्रमणीयता के कारण, फूरियर रूपांतरण को समय के बजाय आवृत्ति के संदर्भ में फ़ंक्शन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह आवृत्ति डोमेन प्रतिनिधित्व है। समय डोमेन में निष्पादित किए जा सकने वाले रैखिक परिचालनों में ऐसे समकक्ष होते हैं जिन्हें अक्सर आवृत्ति डोमेन में अधिक आसानी से निष्पादित किया जा सकता है। फ़्रिक्वेंसी विश्लेषण रैखिक और गैर-रेखीय दोनों, विभिन्न समय-डोमेन संचालन के प्रभावों की समझ और व्याख्या को भी सरल बनाता है। उदाहरण के लिए, केवल गैर-रैखिकता|गैर-रैखिक या समय-संस्करण प्रणाली|समय-संस्करण संचालन ही आवृत्ति स्पेक्ट्रम में नई आवृत्तियाँ बना सकते हैं। | ||
व्यवहार में, लगभग सभी सॉफ्टवेयर और इलेक्ट्रॉनिक उपकरण जो आवृत्ति स्पेक्ट्रा उत्पन्न करते हैं, | व्यवहार में, लगभग सभी सॉफ्टवेयर और इलेक्ट्रॉनिक उपकरण जो आवृत्ति स्पेक्ट्रा उत्पन्न करते हैं, [[उलटा फूरियर रूपांतरण]] (डीएफटी) का उपयोग करते हैं, जो सिग्नल के नमूने (सिग्नल प्रोसेसिंग) पर काम करता है, और जो पूर्ण अभिन्न समाधान के लिए गणितीय अनुमान प्रदान करता है। डीएफटी लगभग हमेशा कुशल एल्गोरिदम द्वारा कार्यान्वित किया जाता है जिसे [[[[असतत फूरियर रूपांतरण]]]] (एफएफटी) कहा जाता है। डीएफटी के वर्ग-परिमाण घटकों की सरणी प्रकार का पावर स्पेक्ट्रम है जिसे [[ periodogram ]] कहा जाता है, जिसका व्यापक रूप से [[आवेग प्रतिक्रिया]] और [[विंडो फ़ंक्शन]] जैसे शोर-मुक्त कार्यों की आवृत्ति विशेषताओं की जांच के लिए उपयोग किया जाता है। लेकिन कम सिग्नल-टू-शोर अनुपात पर शोर जैसे संकेतों या यहां तक कि साइनसोइड्स पर लागू होने पर पीरियोडोग्राम प्रसंस्करण-लाभ प्रदान नहीं करता है। दूसरे शब्दों में, किसी दी गई आवृत्ति पर इसके वर्णक्रमीय अनुमान का विचरण कम नहीं होता है क्योंकि गणना में उपयोग किए गए नमूनों की संख्या बढ़ जाती है। इसे समय के साथ औसत करके (वेल्च की विधि) कम किया जा सकता है<ref>{{Citation|last=Welch|first=P. D.|title=The use of Fast Fourier Transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms|journal=IEEE Transactions on Audio and Electroacoustics|volume=AU-15|issue=2|year=1967|pages=70–73|doi=10.1109/TAU.1967.1161901|bibcode=1967ITAE...15...70W}}</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 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 41: | Line 41: | ||
* [[मूविंग-एवरेज मॉडल]] (एमए) अनुमान, जो मानता है कि एनवां नमूना पिछले पी नमूनों में शोर शर्तों के साथ सहसंबद्ध है। | * [[मूविंग-एवरेज मॉडल]] (एमए) अनुमान, जो मानता है कि एनवां नमूना पिछले पी नमूनों में शोर शर्तों के साथ सहसंबद्ध है। | ||
* [[ऑटोरेग्रेसिव मूविंग एवरेज]] (एआरएमए) अनुमान, जो एआर और एमए मॉडल का सामान्यीकरण करता है। | * [[ऑटोरेग्रेसिव मूविंग एवरेज]] (एआरएमए) अनुमान, जो एआर और एमए मॉडल का सामान्यीकरण करता है। | ||
* [[संगीत (एल्गोरिदम)]] (संगीत) | * [[संगीत (एल्गोरिदम)]] (संगीत) लोकप्रिय [[सुपर-रिज़ॉल्यूशन इमेजिंग]] विधि है। | ||
* अधिकतम एन्ट्रापी वर्णक्रमीय आकलन | * अधिकतम एन्ट्रापी वर्णक्रमीय आकलन पूर्ण-ध्रुव विधि है जो एसडीई के लिए उपयोगी है जब एकल वर्णक्रमीय विशेषताएं, जैसे तेज चोटियां, अपेक्षित होती हैं। | ||
और अंत में अर्ध-पैरामीट्रिक तकनीकों के कुछ उदाहरण: | और अंत में अर्ध-पैरामीट्रिक तकनीकों के कुछ उदाहरण: | ||
* विरल पुनरावृत्तीय सहप्रसरण-आधारित अनुमान (स्पाइस) अनुमान,<ref name=":0" />और अधिक सामान्यीकृत <math>(r,q)</math>-मसाला।<ref>{{Cite journal|last1=Sward|first1=Johan|last2=Adalbjornsson|first2=Stefan Ingi|last3=Jakobsson|first3=Andreas|date=March 2017|title=विरल पुनरावृत्तीय सहप्रसरण-आधारित अनुमानक का एक सामान्यीकरण|url=http://dx.doi.org/10.1109/icassp.2017.7952898|journal=2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)|pages=3954–3958 |publisher=IEEE|doi=10.1109/icassp.2017.7952898|isbn=978-1-5090-4117-6 |s2cid=5640068 }}</ref> *पुनरावृत्तीय अनुकूली दृष्टिकोण (आईएए) अनुमान।<ref>{{Cite journal|last1=Yardibi|first1=Tarik|last2=Li|first2=Jian|last3=Stoica|first3=Petre|last4=Xue|first4=Ming|last5=Baggeroer|first5=Arthur B.|date=January 2010|title=Source Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares|url=https://ieeexplore.ieee.org/document/5417172|journal=IEEE Transactions on Aerospace and Electronic Systems|volume=46|issue=1|pages=425–443|doi=10.1109/TAES.2010.5417172|bibcode=2010ITAES..46..425Y |hdl=1721.1/59588 |s2cid=18834345 |issn=0018-9251|hdl-access=free}}</ref> *[[लैस्सो (सांख्यिकी)]], कम से कम वर्ग वर्णक्रमीय विश्लेषण के समान लेकिन | * विरल पुनरावृत्तीय सहप्रसरण-आधारित अनुमान (स्पाइस) अनुमान,<ref name=":0" />और अधिक सामान्यीकृत <math>(r,q)</math>-मसाला।<ref>{{Cite journal|last1=Sward|first1=Johan|last2=Adalbjornsson|first2=Stefan Ingi|last3=Jakobsson|first3=Andreas|date=March 2017|title=विरल पुनरावृत्तीय सहप्रसरण-आधारित अनुमानक का एक सामान्यीकरण|url=http://dx.doi.org/10.1109/icassp.2017.7952898|journal=2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)|pages=3954–3958 |publisher=IEEE|doi=10.1109/icassp.2017.7952898|isbn=978-1-5090-4117-6 |s2cid=5640068 }}</ref> *पुनरावृत्तीय अनुकूली दृष्टिकोण (आईएए) अनुमान।<ref>{{Cite journal|last1=Yardibi|first1=Tarik|last2=Li|first2=Jian|last3=Stoica|first3=Petre|last4=Xue|first4=Ming|last5=Baggeroer|first5=Arthur B.|date=January 2010|title=Source Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares|url=https://ieeexplore.ieee.org/document/5417172|journal=IEEE Transactions on Aerospace and Electronic Systems|volume=46|issue=1|pages=425–443|doi=10.1109/TAES.2010.5417172|bibcode=2010ITAES..46..425Y |hdl=1721.1/59588 |s2cid=18834345 |issn=0018-9251|hdl-access=free}}</ref> *[[लैस्सो (सांख्यिकी)]], कम से कम वर्ग वर्णक्रमीय विश्लेषण के समान लेकिन विरल दंड लागू करने के साथ।<ref>{{Cite journal|last1=Panahi|first1=Ashkan|last2=Viberg|first2=Mats|date=February 2011|title=LASSO-आधारित DOA आकलन पद्धति के रिज़ॉल्यूशन पर|url=http://dx.doi.org/10.1109/wsa.2011.5741938|journal=2011 International ITG Workshop on Smart Antennas|pages=1–5 |publisher=IEEE|doi=10.1109/wsa.2011.5741938|isbn=978-1-61284-075-8 |s2cid=7013162 }}</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> | ||
जहां <math>\phi_1,\ldots,\phi_p</math> निश्चित गुणांक हैं और <math>\epsilon_t</math> शून्य माध्य और नवीनता विचरण वाली | जहां <math>\phi_1,\ldots,\phi_p</math> निश्चित गुणांक हैं और <math>\epsilon_t</math> शून्य माध्य और नवीनता विचरण वाली श्वेत रव प्रक्रिया है <math>\sigma^2_p</math>. इस प्रक्रिया के लिए एसडीएफ है | ||
:<math> | :<math> | ||
| Line 66: | Line 66: | ||
* ऑटोरेग्रेसिव मॉडल#यूल-वॉकर समीकरण|यूल-वॉकर अनुमानक यूल-वॉकर समीकरणों को पुनरावर्ती रूप से हल करके पाए जाते हैं <math>\text{AR}(p)</math> प्रक्रिया | * ऑटोरेग्रेसिव मॉडल#यूल-वॉकर समीकरण|यूल-वॉकर अनुमानक यूल-वॉकर समीकरणों को पुनरावर्ती रूप से हल करके पाए जाते हैं <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> प्रतिगमन समस्या के रूप में प्रक्रिया करें और आगे-पीछे विधि का उपयोग करके उस समस्या को हल करें। वे बर्ग अनुमानकर्ताओं के साथ प्रतिस्पर्धी हैं। | ||
* अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें | * अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें अरेखीय अनुकूलन शामिल है और यह पहले तीन की तुलना में अधिक जटिल है। | ||
वैकल्पिक पैरामीट्रिक तरीकों में [[ चलती औसत मॉडल ]] (एमए) और पूर्ण ऑटोरेग्रेसिव मूविंग एवरेज मॉडल (एआरएमए) में फिट होना शामिल है। | वैकल्पिक पैरामीट्रिक तरीकों में [[ चलती औसत मॉडल ]] (एमए) और पूर्ण ऑटोरेग्रेसिव मूविंग एवरेज मॉडल (एआरएमए) में फिट होना शामिल है। | ||
| Line 82: | Line 82: | ||
=== एकाधिक स्वर === | === एकाधिक स्वर === | ||
सिग्नल के लिए | सिग्नल के लिए विशिष्ट मॉडल <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> शोर के कारण वर्णक्रमीय घनत्व फ़ंक्शन के अलावा आवेग कार्य भी होता है। | ||
आवृत्ति अनुमान के लिए सबसे आम तरीकों में इन घटकों को निकालने के लिए शोर [[रैखिक उपस्थान]] की पहचान करना शामिल है। ये विधियाँ | आवृत्ति अनुमान के लिए सबसे आम तरीकों में इन घटकों को निकालने के लिए शोर [[रैखिक उपस्थान]] की पहचान करना शामिल है। ये विधियाँ सिग्नल उप-स्थान और शोर उप-स्थान में ऑटोसहसंबंध मैट्रिक्स के [[Eigendecomposition]] पर आधारित हैं। इन उप-स्थानों की पहचान होने के बाद, शोर उप-स्थान से घटक आवृत्तियों को खोजने के लिए आवृत्ति अनुमान फ़ंक्शन का उपयोग किया जाता है। शोर उप-स्थान आधारित आवृत्ति अनुमान की सबसे लोकप्रिय विधियाँ हैं पिसारेंको हार्मोनिक अपघटन|पिसारेंको की विधि, मल्टीपल सिग्नल वर्गीकरण (संगीत) विधि, ईजेनवेक्टर विधि और न्यूनतम मानक विधि। | ||
; पिसारेंको हार्मोनिक अपघटन|पिसारेंको की विधि: <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 96: | Line 96: | ||
== उदाहरण गणना == | == उदाहरण गणना == | ||
कल्पना करना <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 118: | Line 118: | ||
:<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>x(t)</math> आवृत्ति के साथ घटक से आ रहा है <math>\nu_k</math> है <math>\tfrac{1}{2}A_k^2.</math> ये सभी योगदान औसत शक्ति में जुड़ जाते हैं <math>x(t).</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>S</math> चरणीय फ़ंक्शन है, जो नीरस रूप से घटता नहीं है। इसकी छलांग अवधि (रिंग) घटकों की आवृत्तियों पर होती है <math>x</math>, और प्रत्येक छलांग का मूल्य उस घटक की शक्ति या भिन्नता है। | ||
विचरण स्वयं के साथ डेटा का सहप्रसरण है। यदि हम अब उसी डेटा पर विचार करें लेकिन थोड़े अंतराल के साथ <math>\tau</math>, हम इसका सहप्रसरण ले सकते हैं <math>x(t)</math> साथ <math>x(t + \tau)</math>, और इसे स्वतःसहसंबंध फ़ंक्शन के रूप में परिभाषित करें <math>c</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>x</math>: | ऐसा दिखाया जा सकता है <math>c</math> समान अवधियों के साथ आवधिक घटकों में विघटित किया जा सकता है <math>x</math>: | ||
| Line 134: | Line 134: | ||
यह वास्तव में का वर्णक्रमीय अपघटन है <math>c</math> विभिन्न आवृत्तियों पर, और शक्ति के वितरण से संबंधित है <math>x</math> आवृत्तियों पर: आवृत्ति घटक का आयाम <math>c</math> सिग्नल की औसत शक्ति में इसका योगदान है। | यह वास्तव में का वर्णक्रमीय अपघटन है <math>c</math> विभिन्न आवृत्तियों पर, और शक्ति के वितरण से संबंधित है <math>x</math> आवृत्तियों पर: आवृत्ति घटक का आयाम <math>c</math> सिग्नल की औसत शक्ति में इसका योगदान है। | ||
इस उदाहरण का पावर स्पेक्ट्रम निरंतर नहीं है, और इसलिए इसका कोई व्युत्पन्न नहीं है, और इसलिए इस सिग्नल में पावर स्पेक्ट्रल घनत्व फ़ंक्शन नहीं है। सामान्य तौर पर, पावर स्पेक्ट्रम आम तौर पर दो भागों का योग होगा: | इस उदाहरण का पावर स्पेक्ट्रम निरंतर नहीं है, और इसलिए इसका कोई व्युत्पन्न नहीं है, और इसलिए इस सिग्नल में पावर स्पेक्ट्रल घनत्व फ़ंक्शन नहीं है। सामान्य तौर पर, पावर स्पेक्ट्रम आम तौर पर दो भागों का योग होगा: लाइन स्पेक्ट्रम जैसे कि इस उदाहरण में, जो निरंतर नहीं है और इसमें घनत्व फ़ंक्शन नहीं है, और अवशेष, जो बिल्कुल निरंतर है और इसमें घनत्व फ़ंक्शन होता है . | ||
== यह भी देखें == | == यह भी देखें == | ||
Revision as of 14:40, 8 July 2023
सांख्यिकीय सिग्नल प्रोसेसिंग में, वर्णक्रमीय घनत्व अनुमान (एसडीई) या केवल वर्णक्रमीय अनुमान का लक्ष्य सिग्नल के समय नमूनों के अनुक्रम से सिग्नल के वर्णक्रमीय घनत्व (जिसे पावर स्पेक्ट्रम के रूप में भी जाना जाता है) का अनुमान लगाना है।[1] सहज रूप से कहें तो, वर्णक्रमीय घनत्व सिग्नल की आवृत्ति सामग्री को दर्शाता है। वर्णक्रमीय घनत्व का अनुमान लगाने का उद्देश्य इन आवधिकों के अनुरूप आवृत्तियों पर चोटियों को देखकर, डेटा में किसी भी आवधिक फ़ंक्शन का पता लगाना है।
कुछ एसडीई तकनीकें मानती हैं कि सिग्नल सीमित (आमतौर पर छोटी) संख्या में उत्पन्न आवृत्तियों और शोर से बना होता है और उत्पन्न आवृत्तियों के स्थान और तीव्रता का पता लगाने की कोशिश करता है। अन्य लोग घटकों की संख्या पर कोई धारणा नहीं बनाते हैं और संपूर्ण उत्पादन स्पेक्ट्रम का अनुमान लगाना चाहते हैं।
सिंहावलोकन
स्पेक्ट्रम विश्लेषण, जिसे आवृत्ति डोमेन विश्लेषण या वर्णक्रमीय घनत्व अनुमान भी कहा जाता है, जटिल सिग्नल को सरल भागों में विघटित करने की तकनीकी प्रक्रिया है। जैसा कि ऊपर वर्णित है, कई भौतिक प्रक्रियाओं को कई व्यक्तिगत आवृत्ति घटकों के योग के रूप में सबसे अच्छा वर्णित किया गया है। कोई भी प्रक्रिया जो विभिन्न मात्राओं (जैसे आयाम, शक्तियाँ, तीव्रता) बनाम आवृत्ति (या चरण (तरंगें)) की मात्रा निर्धारित करती है, उसे स्पेक्ट्रम विश्लेषण कहा जा सकता है।
स्पेक्ट्रम विश्लेषण पूरे सिग्नल पर किया जा सकता है। वैकल्पिक रूप से, सिग्नल को छोटे खंडों (कभी-कभी फ़्रेम कहा जाता है) में तोड़ा जा सकता है, और स्पेक्ट्रम विश्लेषण को इन व्यक्तिगत खंडों पर लागू किया जा सकता है। आवधिक कार्य (जैसे ) इस उप-विभाजन के लिए विशेष रूप से उपयुक्त हैं। गैर-आवधिक कार्यों के विश्लेषण के लिए सामान्य गणितीय तकनीकें फूरियर विश्लेषण की श्रेणी में आती हैं।
किसी फ़ंक्शन का फूरियर रूपांतरण आवृत्ति स्पेक्ट्रम उत्पन्न करता है जिसमें मूल सिग्नल के बारे में सारी जानकारी होती है, लेकिन अलग रूप में। इसका मतलब यह है कि मूल फ़ंक्शन को व्युत्क्रम फूरियर रूपांतरण द्वारा पूरी तरह से पुनर्निर्मित (संश्लेषित) किया जा सकता है। सही पुनर्निर्माण के लिए, स्पेक्ट्रम विश्लेषक को प्रत्येक आवृत्ति घटक के आयाम और चरण (तरंगों) दोनों को संरक्षित करना होगा। जानकारी के इन दो टुकड़ों को 2-आयामी वेक्टर के रूप में, जटिल संख्या के रूप में, या ध्रुवीय निर्देशांक में परिमाण (आयाम) और चरण के रूप में (यानी, चरण के रूप में) दर्शाया जा सकता नमूनाकरण (सिग्नल प्रोसेसिंग) में सामान्य तकनीक वर्ग आयाम, या शक्ति (भौतिकी) पर विचार करना है; इस मामले में परिणामी प्लॉट को पावर स्पेक्ट्रम के रूप में जाना जाता है।
उत्क्रमणीयता के कारण, फूरियर रूपांतरण को समय के बजाय आवृत्ति के संदर्भ में फ़ंक्शन का प्रतिनिधित्व कहा जाता है; इस प्रकार, यह आवृत्ति डोमेन प्रतिनिधित्व है। समय डोमेन में निष्पादित किए जा सकने वाले रैखिक परिचालनों में ऐसे समकक्ष होते हैं जिन्हें अक्सर आवृत्ति डोमेन में अधिक आसानी से निष्पादित किया जा सकता है। फ़्रिक्वेंसी विश्लेषण रैखिक और गैर-रेखीय दोनों, विभिन्न समय-डोमेन संचालन के प्रभावों की समझ और व्याख्या को भी सरल बनाता है। उदाहरण के लिए, केवल गैर-रैखिकता|गैर-रैखिक या समय-संस्करण प्रणाली|समय-संस्करण संचालन ही आवृत्ति स्पेक्ट्रम में नई आवृत्तियाँ बना सकते हैं।
व्यवहार में, लगभग सभी सॉफ्टवेयर और इलेक्ट्रॉनिक उपकरण जो आवृत्ति स्पेक्ट्रा उत्पन्न करते हैं, उलटा फूरियर रूपांतरण (डीएफटी) का उपयोग करते हैं, जो सिग्नल के नमूने (सिग्नल प्रोसेसिंग) पर काम करता है, और जो पूर्ण अभिन्न समाधान के लिए गणितीय अनुमान प्रदान करता है। डीएफटी लगभग हमेशा कुशल एल्गोरिदम द्वारा कार्यान्वित किया जाता है जिसे [[असतत फूरियर रूपांतरण]] (एफएफटी) कहा जाता है। डीएफटी के वर्ग-परिमाण घटकों की सरणी प्रकार का पावर स्पेक्ट्रम है जिसे periodogram कहा जाता है, जिसका व्यापक रूप से आवेग प्रतिक्रिया और विंडो फ़ंक्शन जैसे शोर-मुक्त कार्यों की आवृत्ति विशेषताओं की जांच के लिए उपयोग किया जाता है। लेकिन कम सिग्नल-टू-शोर अनुपात पर शोर जैसे संकेतों या यहां तक कि साइनसोइड्स पर लागू होने पर पीरियोडोग्राम प्रसंस्करण-लाभ प्रदान नहीं करता है। दूसरे शब्दों में, किसी दी गई आवृत्ति पर इसके वर्णक्रमीय अनुमान का विचरण कम नहीं होता है क्योंकि गणना में उपयोग किए गए नमूनों की संख्या बढ़ जाती है। इसे समय के साथ औसत करके (वेल्च की विधि) कम किया जा सकता है[2]) या अधिक आवृत्ति (चौरसाई )। वर्णक्रमीय घनत्व आकलन (एसडीई) के लिए वेल्च की विधि का व्यापक रूप से उपयोग किया जाता है। हालाँकि, पीरियोडोग्राम-आधारित तकनीकें छोटे पूर्वाग्रह पेश करती हैं जो कुछ अनुप्रयोगों में अस्वीकार्य हैं। इसलिए अन्य विकल्प अगले भाग में प्रस्तुत किए गए हैं।
तकनीक
बुनियादी आवर्त सारणी की कमियों को कम करने के लिए वर्णक्रमीय आकलन की कई अन्य तकनीकें विकसित की गई हैं। इन तकनीकों को आम तौर पर गैर-पैरामीट्रिक सांख्यिकी|गैर-पैरामीट्रिक, पैरामीट्रिक अनुमान, और हाल ही में सेमीपैरामीट्रिक मॉडल|अर्ध-पैरामीट्रिक (जिसे विरल भी कहा जाता है) विधियों में विभाजित किया जा सकता है।[3] गैर-पैरामीट्रिक दृष्टिकोण स्पष्ट रूप से सहप्रसरण या प्रक्रिया के स्पेक्ट्रम का अनुमान लगाते हैं, बिना यह माने कि प्रक्रिया में कोई विशेष संरचना है। बुनियादी अनुप्रयोगों (उदाहरण के लिए वेल्च की विधि) के लिए उपयोग में आने वाले कुछ सबसे आम अनुमानक गैर-पैरामीट्रिक अनुमानक हैं जो पीरियोडोग्राम से निकटता से संबंधित हैं। इसके विपरीत, पैरामीट्रिक दृष्टिकोण यह मानते हैं कि अंतर्निहित स्थिर प्रक्रिया में निश्चित संरचना होती है जिसे कम संख्या में मापदंडों का उपयोग करके वर्णित किया जा सकता है (उदाहरण के लिए, ऑटोरेग्रेसिव मूविंग एवरेज मॉडल | ऑटो-रिग्रेसिव या मूविंग एवरेज मॉडल का उपयोग करके)। इन दृष्टिकोणों में, कार्य उस मॉडल के मापदंडों का अनुमान लगाना है जो स्टोकेस्टिक प्रक्रिया का वर्णन करता है। अर्ध-पैरामीट्रिक विधियों का उपयोग करते समय, अंतर्निहित प्रक्रिया को गैर-पैरामीट्रिक ढांचे का उपयोग करके मॉडलिंग किया जाता है, अतिरिक्त धारणा के साथ कि मॉडल के गैर-शून्य घटकों की संख्या छोटी है (यानी, मॉडल विरल है)। गुम डेटा पुनर्प्राप्ति के लिए भी इसी तरह के तरीकों का उपयोग किया जा सकता है [4] साथ ही संपीड़ित संवेदन।
गैर-पैरामीट्रिक वर्णक्रमीय घनत्व आकलन तकनीकों की आंशिक सूची निम्नलिखित है:
- पीरियोडोग्राम, असतत फूरियर रूपांतरण का मापांक वर्ग
- लोम्ब-स्कार्गल पीरियोडोग्राम, जिसके लिए डेटा को समान रूप से स्थान देने की आवश्यकता नहीं है
- बार्टलेट की विधि वर्णक्रमीय घनत्व अनुमान के विचरण को कम करने के लिए सिग्नल के कई खंडों से लिए गए पीरियडोग्राम का औसत है
- वेल्च की विधि बार्टलेट की विधि का विंडो संस्करण है जो ओवरलैपिंग सेगमेंट का उपयोग करती है
- मल्टीटेपर पीरियडोग्राम-आधारित विधि है जो वर्णक्रमीय घनत्व अनुमान के विचरण को कम करने के लिए वर्णक्रमीय घनत्व का स्वतंत्र अनुमान बनाने के लिए कई टेपर या विंडो का उपयोग करती है।
- न्यूनतम-वर्ग वर्णक्रमीय विश्लेषण, ज्ञात आवृत्तियों के अनुरूप न्यूनतम वर्गों पर आधारित
- गैर-समान असतत फूरियर ट्रांसफॉर्म का उपयोग तब किया जाता है जब सिग्नल नमूने असमान रूप से समय श्रृंखला में होते हैं
- एकवचन स्पेक्ट्रम विश्लेषण गैरपैरामीट्रिक विधि है जो वर्णक्रमीय घनत्व का अनुमान लगाने के लिए सहप्रसरण मैट्रिक्स के एकवचन मूल्य अपघटन का उपयोग करता है
- अल्पकालीन फूरियर रूपांतरण
- सूचना क्षेत्र सिद्धांत#क्रिटिकल फिल्टर सूचना क्षेत्र सिद्धांत पर आधारित गैर-पैरामीट्रिक विधि है जो शोर, अपूर्ण डेटा और वाद्य प्रतिक्रिया कार्यों से निपट सकती है।
नीचे पैरामीट्रिक तकनीकों की आंशिक सूची दी गई है:
- ऑटोरेग्रेसिव मॉडल (एआर) अनुमान, जो मानता है कि एनवां नमूना पिछले पी नमूनों के साथ सहसंबद्ध है।
- मूविंग-एवरेज मॉडल (एमए) अनुमान, जो मानता है कि एनवां नमूना पिछले पी नमूनों में शोर शर्तों के साथ सहसंबद्ध है।
- ऑटोरेग्रेसिव मूविंग एवरेज (एआरएमए) अनुमान, जो एआर और एमए मॉडल का सामान्यीकरण करता है।
- संगीत (एल्गोरिदम) (संगीत) लोकप्रिय सुपर-रिज़ॉल्यूशन इमेजिंग विधि है।
- अधिकतम एन्ट्रापी वर्णक्रमीय आकलन पूर्ण-ध्रुव विधि है जो एसडीई के लिए उपयोगी है जब एकल वर्णक्रमीय विशेषताएं, जैसे तेज चोटियां, अपेक्षित होती हैं।
और अंत में अर्ध-पैरामीट्रिक तकनीकों के कुछ उदाहरण:
- विरल पुनरावृत्तीय सहप्रसरण-आधारित अनुमान (स्पाइस) अनुमान,[3]और अधिक सामान्यीकृत -मसाला।[5] *पुनरावृत्तीय अनुकूली दृष्टिकोण (आईएए) अनुमान।[6] *लैस्सो (सांख्यिकी), कम से कम वर्ग वर्णक्रमीय विश्लेषण के समान लेकिन विरल दंड लागू करने के साथ।[7]
पैरामीट्रिक अनुमान
पैरामीट्रिक वर्णक्रमीय अनुमान में, कोई यह मानता है कि सिग्नल स्थिर प्रक्रिया द्वारा तैयार किया गया है जिसमें वर्णक्रमीय घनत्व फ़ंक्शन (एसडीएफ) है यह आवृत्ति का कार्य है और पैरामीटर .[8] फिर अनुमान की समस्या इन मापदंडों का अनुमान लगाने में से बन जाती है।
पैरामीट्रिक एसडीएफ अनुमान का सबसे सामान्य रूप मॉडल के रूप में ऑटोरेग्रेसिव मॉडल का उपयोग करता है आदेश की .[8]: 392 संकेत अनुक्रम शून्य माध्य का पालन करना प्रक्रिया समीकरण को संतुष्ट करती है
जहां निश्चित गुणांक हैं और शून्य माध्य और नवीनता विचरण वाली श्वेत रव प्रक्रिया है . इस प्रक्रिया के लिए एसडीएफ है
साथ नमूनाकरण समय अंतराल और नाइक्विस्ट आवृत्ति.
मापदंडों का अनुमान लगाने के लिए कई दृष्टिकोण हैं की प्रक्रिया और इस प्रकार वर्णक्रमीय घनत्व:[8]: 452-453
- ऑटोरेग्रेसिव मॉडल#यूल-वॉकर समीकरण|यूल-वॉकर अनुमानक यूल-वॉकर समीकरणों को पुनरावर्ती रूप से हल करके पाए जाते हैं प्रक्रिया
- बर्ग अनुमानक यूल-वॉकर समीकरणों को सामान्य न्यूनतम वर्ग समस्या के रूप में मानकर पाए जाते हैं। बर्ग अनुमानकों को आम तौर पर यूल-वॉकर अनुमानकों से बेहतर माना जाता है।[8]: 452 बर्ग ने इन्हें अधिकतम एन्ट्रॉपी वर्णक्रमीय अनुमान के साथ जोड़ा।[9]
- आगे-पीछे न्यूनतम-वर्ग अनुमानक का व्यवहार करते हैं प्रतिगमन समस्या के रूप में प्रक्रिया करें और आगे-पीछे विधि का उपयोग करके उस समस्या को हल करें। वे बर्ग अनुमानकर्ताओं के साथ प्रतिस्पर्धी हैं।
- अधिकतम संभावना अनुमानक अधिकतम संभावना दृष्टिकोण का उपयोग करके मापदंडों का अनुमान लगाते हैं। इसमें अरेखीय अनुकूलन शामिल है और यह पहले तीन की तुलना में अधिक जटिल है।
वैकल्पिक पैरामीट्रिक तरीकों में चलती औसत मॉडल (एमए) और पूर्ण ऑटोरेग्रेसिव मूविंग एवरेज मॉडल (एआरएमए) में फिट होना शामिल है।
आवृत्ति अनुमान
फ़्रिक्वेंसी अनुमान अनुमान सिद्धांत की प्रक्रिया है जो घटकों की संख्या के बारे में दी गई धारणाओं के शोर की उपस्थिति में अंकीय संकेत प्रक्रिया की आवृत्ति, आयाम और चरण-शिफ्ट है।[10] यह उपरोक्त सामान्य तरीकों के विपरीत है, जो घटकों के बारे में पूर्व धारणा नहीं बनाते हैं।
एकल स्वर
यदि कोई केवल सबसे ऊंची आवृत्ति का अनुमान लगाना चाहता है, तो वह पिच का पता लगाने का एल्गोरिदम का उपयोग कर सकता है। यदि प्रमुख आवृत्ति समय के साथ बदलती है, तो समस्या तात्कालिक आवृत्ति के अनुमान की हो जाती है जैसा कि समय-आवृत्ति प्रतिनिधित्व में परिभाषित किया गया है। तात्कालिक आवृत्ति अनुमान के तरीकों में विग्नर-विले वितरण और उच्च क्रम अस्पष्टता कार्यों पर आधारित तरीके शामिल हैं।[11] यदि कोई प्राप्त सिग्नल के सभी (संभवतः जटिल) आवृत्ति घटकों (संचरित सिग्नल और शोर सहित) को जानना चाहता है, तो वह मल्टी-टोन दृष्टिकोण का उपयोग करता है।
एकाधिक स्वर
सिग्नल के लिए विशिष्ट मॉडल का योग होता है सफ़ेद शोर की उपस्थिति में जटिल घातांक,
- .
की शक्ति वर्णक्रमीय घनत्व से बना है शोर के कारण वर्णक्रमीय घनत्व फ़ंक्शन के अलावा आवेग कार्य भी होता है।
आवृत्ति अनुमान के लिए सबसे आम तरीकों में इन घटकों को निकालने के लिए शोर रैखिक उपस्थान की पहचान करना शामिल है। ये विधियाँ सिग्नल उप-स्थान और शोर उप-स्थान में ऑटोसहसंबंध मैट्रिक्स के Eigendecomposition पर आधारित हैं। इन उप-स्थानों की पहचान होने के बाद, शोर उप-स्थान से घटक आवृत्तियों को खोजने के लिए आवृत्ति अनुमान फ़ंक्शन का उपयोग किया जाता है। शोर उप-स्थान आधारित आवृत्ति अनुमान की सबसे लोकप्रिय विधियाँ हैं पिसारेंको हार्मोनिक अपघटन|पिसारेंको की विधि, मल्टीपल सिग्नल वर्गीकरण (संगीत) विधि, ईजेनवेक्टर विधि और न्यूनतम मानक विधि।
- पिसारेंको हार्मोनिक अपघटन|पिसारेंको की विधि
- एकाधिक सिग्नल वर्गीकरण
- ,
- आइजेनवेक्टर विधि
- न्यूनतम मानक विधि
उदाहरण गणना
कल्पना करना , से को शून्य माध्य वाली समय श्रृंखला (अलग समय) है। मान लीजिए कि यह आवधिक घटकों की सीमित संख्या का योग है (सभी आवृत्तियाँ सकारात्मक हैं):
का विचरण जैसा कि ऊपर दिया गया है, शून्य-माध्य फ़ंक्शन के लिए है
यदि ये डेटा विद्युत सिग्नल से लिए गए नमूने थे, तो यह इसकी औसत शक्ति होगी (शक्ति प्रति यूनिट समय ऊर्जा है, इसलिए यदि ऊर्जा आयाम वर्ग के अनुरूप है तो यह विचरण के अनुरूप है)।
अब, सरलता के लिए, मान लीजिए कि संकेत समय में अनंत रूप से फैलता है, इसलिए हम सीमा को पार कर जाते हैं यदि औसत शक्ति सीमित है, जो वास्तविकता में लगभग हमेशा मामला होता है, तो निम्न सीमा मौजूद होती है और डेटा का भिन्नता होती है।
फिर से, सरलता के लिए, हम निरंतर समय पर जाएंगे, और मान लेंगे कि संकेत दोनों दिशाओं में समय में अनंत रूप से फैलता है। तब ये दो सूत्र बन जाते हैं
और
मूल माध्य का वर्ग है , तो का विचरण है इसलिए, की औसत शक्ति में योगदान आवृत्ति के साथ घटक से आ रहा है है ये सभी योगदान औसत शक्ति में जुड़ जाते हैं फिर आवृत्ति के फलन के रूप में शक्ति है और इसका सांख्यिकीय संचयी वितरण कार्य होगा
चरणीय फ़ंक्शन है, जो नीरस रूप से घटता नहीं है। इसकी छलांग अवधि (रिंग) घटकों की आवृत्तियों पर होती है , और प्रत्येक छलांग का मूल्य उस घटक की शक्ति या भिन्नता है।
विचरण स्वयं के साथ डेटा का सहप्रसरण है। यदि हम अब उसी डेटा पर विचार करें लेकिन थोड़े अंतराल के साथ , हम इसका सहप्रसरण ले सकते हैं साथ , और इसे स्वतःसहसंबंध फ़ंक्शन के रूप में परिभाषित करें सिग्नल (या डेटा) का :
यदि यह अस्तित्व में है, तो यह सम कार्य है यदि औसत शक्ति परिबद्ध है, तो सर्वत्र विद्यमान है, परिमित है और सीमाबद्ध है जो डेटा की औसत शक्ति या विचरण है।
ऐसा दिखाया जा सकता है समान अवधियों के साथ आवधिक घटकों में विघटित किया जा सकता है :
यह वास्तव में का वर्णक्रमीय अपघटन है विभिन्न आवृत्तियों पर, और शक्ति के वितरण से संबंधित है आवृत्तियों पर: आवृत्ति घटक का आयाम सिग्नल की औसत शक्ति में इसका योगदान है।
इस उदाहरण का पावर स्पेक्ट्रम निरंतर नहीं है, और इसलिए इसका कोई व्युत्पन्न नहीं है, और इसलिए इस सिग्नल में पावर स्पेक्ट्रल घनत्व फ़ंक्शन नहीं है। सामान्य तौर पर, पावर स्पेक्ट्रम आम तौर पर दो भागों का योग होगा: लाइन स्पेक्ट्रम जैसे कि इस उदाहरण में, जो निरंतर नहीं है और इसमें घनत्व फ़ंक्शन नहीं है, और अवशेष, जो बिल्कुल निरंतर है और इसमें घनत्व फ़ंक्शन होता है .
यह भी देखें
- बहुआयामी वर्णक्रमीय अनुमान
- पीरियोडोग्राम
- सिगस्पेक
- spectrogram
- समय-आवृत्ति विश्लेषण
- समय-आवृत्ति प्रतिनिधित्व
- कम संभावना
- वर्णक्रमीय विद्युत वितरण
संदर्भ
- ↑ P Stoica and R Moses, Spectral Analysis of Signals, Prentice Hall, 2005.
- ↑ 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
- ↑ 3.0 3.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.
- ↑ 8.0 8.1 8.2 8.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.