चिर्प जेड-रूपांतरण: Difference between revisions
(Created page with "{{refimprove|date=January 2013}} चिरप z-परिणत (सीजेडटी) असतत फूरियर रूपांतरण (डीएफटी)...") |
m (9 revisions imported from alpha:चिर्प_जेड-रूपांतरण) |
||
(8 intermediate revisions by 4 users not shown) | |||
Line 1: | Line 1: | ||
'''चिर्प जेड-रूपांतरण''' (सीजेडटी) [[असतत फूरियर रूपांतरण]] (डीएफटी) का एक सामान्यीकरण है। जबकि डीएफटी यूनिट सर्कल के साथ समान रूप से दूरी वाले बिंदुओं पर जेड-रूपांतरण का नमूना लेता है, [[ एस विमान |एस प्लेन]] में सीधी रेखाओं के अनुरूप, जेड-प्लेन में सर्पिल आर्क्स के साथ चहचहाना जेड-रूपांतरण नमूने लेता है।<ref name=Shilling>[https://krex.k-state.edu/dspace/bitstream/handle/2097/7844/LD2668R41972S43.pdf A study of the Chirp Z-transform and its applications] - Shilling, Steve Alan</ref><ref>{{Cite web|url=http://www.mathworks.com/help/signal/ref/czt.html|title=चिरप जेड-ट्रांसफॉर्म - MATLAB czt|website=www.mathworks.com|access-date=2016-09-22}}</ref> इस प्रकार से डीएफटी, वास्तविक डीएफटी और ज़ूम डीएफटी की गणना सीजेडटी के विशेष स्तिथियों के रूप में की जा सकती है। | |||
विशेष रूप से, | विशेष रूप से, चिर्प Z रूपांतरण, Z रूपांतरण की गणना बिंदु z<sub>k</sub> की एक सीमित संख्या पर करता है एक [[लघुगणकीय सर्पिल]] समोच्च के साथ, इसे इस प्रकार परिभाषित किया गया है:<ref name=Shilling/><ref>{{Cite web|url=http://prod.sandia.gov/techlib/access-control.cgi/2005/057084.pdf|title=Chirp Z-Transform Spectral Zoom Optimization with MATLAB®|last=Martin|first=Grant D.|date=November 2005|website=|access-date=}}</ref> | ||
:<math>X_k = \sum_{n=0}^{N-1} x(n) z_{k}^{-n} </math> | :<math>X_k = \sum_{n=0}^{N-1} x(n) z_{k}^{-n} </math> | ||
:<math>z_k = A\cdot W^{-k}, k=0,1,\dots,M-1</math> | :<math>z_k = A\cdot W^{-k}, k=0,1,\dots,M-1</math> | ||
जहां A | जहां A सम्मिश्र प्रारंभिक बिंदु है, W बिंदुओं के बीच सम्मिश्र अनुपात है, और M गणना किए जाने वाले बिंदुओं की संख्या है। | ||
डीएफटी की तरह, चिर्प जेड-रूपांतरण की गणना O(''n'' log ''n'') ऑपरेशंस में की जा सकती है जहां n=\max(M,N) है। | |||
व्युत्क्रम चिर्प जेड-रूपांतरण (आईसीजेडटी) के लिए एक O(''n'' log ''n'') एल्गोरिदम 2003 और 2019 में वर्णित किया गया था।<ref>{{Cite thesis |type=PhD |last=Bostan |first=Alin |date=2003 |title=बुनियादी कंप्यूटर बीजगणित संचालन के लिए कुशल एल्गोरिदम|publisher=Ecole polytechnique |url=https://specfun.inria.fr/bostan/these/These.pdf}}</ref><ref>{{Cite journal|last=Bostan|first=Alin|last2=Schost|first2=Éric|date=2005|title=बिंदुओं के विशेष सेट पर बहुपद मूल्यांकन और प्रक्षेप|journal=Journal of Complexity|language=en|volume=21|issue=4|pages=420–446|doi=10.1016/j.jco.2004.09.009|doi-access=free}}</ref> <ref>[https://scitechdaily.com/engineers-solve-50-year-old-puzzle-in-signal-processing-inverse-chirp-z-transform/ Engineers Solve 50-Year-Old Puzzle in Signal Processing – Inverse Chirp Z-Transform], By IOWA STATE UNIVERSITY OCTOBER 10, 2019</ref> | |||
==ब्लूस्टीन का एल्गोरिदम== | ==ब्लूस्टीन का एल्गोरिदम== | ||
ब्लूस्टीन का एल्गोरिदम<ref name="bluestein">{{Cite journal|last=Bluestein|first=L.|date=1970-12-01|title=असतत फूरियर रूपांतरण की गणना के लिए एक रैखिक फ़िल्टरिंग दृष्टिकोण|journal=IEEE Transactions on Audio and Electroacoustics|volume=18|issue=4|pages=451–455|doi=10.1109/TAU.1970.1162132|issn=0018-9278}}</ref><ref>{{cite web|url=http://www.dsprelated.com/dspbooks/mdft/Bluestein_s_FFT_Algorithm.html |title=ब्लूस्टीन का एफएफटी एल्गोरिदम|publisher=DSPRelated.com}}</ref> सीजेडटी को एक [[कनवल्शन]] के रूप में व्यक्त करता है और तेज़ फूरियर | ब्लूस्टीन का एल्गोरिदम<ref name="bluestein">{{Cite journal|last=Bluestein|first=L.|date=1970-12-01|title=असतत फूरियर रूपांतरण की गणना के लिए एक रैखिक फ़िल्टरिंग दृष्टिकोण|journal=IEEE Transactions on Audio and Electroacoustics|volume=18|issue=4|pages=451–455|doi=10.1109/TAU.1970.1162132|issn=0018-9278}}</ref><ref>{{cite web|url=http://www.dsprelated.com/dspbooks/mdft/Bluestein_s_FFT_Algorithm.html |title=ब्लूस्टीन का एफएफटी एल्गोरिदम|publisher=DSPRelated.com}}</ref> सीजेडटी को एक [[कनवल्शन]] के रूप में व्यक्त करता है और तेज़ फूरियर रूपांतरण/आईएफएफटी का उपयोग करके इसे कुशलतापूर्वक कार्यान्वित करता है। | ||
चूंकि डीएफटी सीजेडटी का | चूंकि डीएफटी सीजेडटी का विशेष स्तिथि है, यह अभाज्य संख्या आकार सहित इच्छानुसार आकारों के [[फास्ट फूरियर ट्रांसफॉर्म|फास्ट फूरियर रूपांतरण]] (डीएफटी) की कुशल गणना की अनुमति देता है। इस प्रकार से (प्राइम साइज के एफएफटी के लिए अन्य एल्गोरिदम, रेडर का एल्गोरिदम, डीएफटी को कनवल्शन के रूप में फिर से लिखकर भी कार्य करता है।) इसकी कल्पना 1968 में [[लियो ब्लूस्टीन]] द्वारा की गई थी।<ref name="bluestein" /> किन्तु ब्लूस्टीन के एल्गोरिदम का उपयोग (एकतरफा) जेड-रूपांतरण (रेबिनर एट अल।, 1969) के आधार पर डीएफटी की तुलना में अधिक सामान्य परिवर्तनों की गणना करने के लिए किया जा सकता है। | ||
याद रखें कि डीएफटी को सूत्र द्वारा परिभाषित किया गया है | याद रखें कि डीएफटी को सूत्र द्वारा परिभाषित किया गया है | ||
Line 21: | Line 19: | ||
:<math> X_k = \sum_{n=0}^{N-1} x_n e^{-\frac{2\pi i}{N} nk } | :<math> X_k = \sum_{n=0}^{N-1} x_n e^{-\frac{2\pi i}{N} nk } | ||
\qquad | \qquad | ||
k = 0,\dots,N-1. </math> | k = 0,\dots,N-1. </math> | ||
यदि हम घातांक में गुणनफल nk को पहचान से प्रतिस्थापित करते हैं | यदि हम घातांक में गुणनफल nk को पहचान से प्रतिस्थापित करते हैं | ||
Line 30: | Line 28: | ||
\qquad | \qquad | ||
k = 0,\dots,N-1. </math> | k = 0,\dots,N-1. </math> | ||
यह | यह योग वास्तव में दो अनुक्रमों a और bn का एक संलयन है जिसे निम्न द्वारा परिभाषित किया गया है: | ||
:<math>a_n = x_n e^{-\frac{\pi i}{N} n^2 }</math> | :<math>a_n = x_n e^{-\frac{\pi i}{N} n^2 }</math> | ||
:<math>b_n = e^{\frac{\pi i}{N} n^2 },</math> | :<math>b_n = e^{\frac{\pi i}{N} n^2 },</math> | ||
कनवल्शन के आउटपुट को N चरण कारकों ''b<sub>k</sub>''<sup>*</sup> से गुणा किया जाता है। वह है: | |||
:<math>X_k = b_k^* \left(\sum_{n=0}^{N-1} a_n b_{k-n}\right) \qquad k = 0,\dots,N-1. </math> | :<math>X_k = b_k^* \left(\sum_{n=0}^{N-1} a_n b_{k-n}\right) \qquad k = 0,\dots,N-1. </math> | ||
यह कनवल्शन, | यह कनवल्शन, परिवर्तन में, [[कनवल्शन प्रमेय]] के माध्यम से एफएफटी की जोड़ी (साथ ही सम्मिश्र [[कलरव]] b<sub>''n''</sub> की पूर्व-गणना की गई एफएफटी)) के साथ किया जा सकता है । इस प्रकार से मुख्य तथ्य यह है कि ये एफएफटी समान लंबाई ''N'' के नहीं हैं: इस तरह के कनवल्शन की गणना एफएफटी से केवल शून्य-पैडिंग द्वारा 2N-1 से अधिक या उसके समान लंबाई तक की जा सकती है। विशेष रूप से, कोई दो या किसी अन्य [[चिकनी संख्या|अत्यधिक समग्र संख्या]] आकार की शक्ति तक पैड कर सकता है, कूली-टुकी एफएफटी एल्गोरिथ्म O(N log N) समय में कूली-टुकी एल्गोरिथ्म जिसके लिए एफएफटी को कुशलतापूर्वक निष्पादित किया जा सकता है। इस प्रकार, ब्लूस्टीन का एल्गोरिदम प्राइम-आकार डीएफटी की गणना करने के लिए एक O(N log N) विधि प्रदान करता है, चूंकि समग्र आकार के लिए कूली-टुकी एल्गोरिदम की तुलना में कई गुना धीमा है। | ||
ब्लूस्टीन के एल्गोरिदम में कनवल्शन के लिए शून्य-पैडिंग का उपयोग कुछ अतिरिक्त टिप्पणी का पात्र है। मान लीजिए कि हम लंबाई M ≥ 2N–1 पर शून्य-पैड करते हैं। इसका | इस प्रकार से ब्लूस्टीन के एल्गोरिदम में कनवल्शन के लिए शून्य-पैडिंग का उपयोग कुछ अतिरिक्त टिप्पणी का पात्र है। मान लीजिए कि हम लंबाई M ≥ 2N–1 पर शून्य-पैड करते हैं। इसका अर्थ यह है कि ''a<sub>n</sub>'' को लंबाई M की एक सरणी ''A<sub>n</sub>'' तक विस्तारित किया गया है, जहां ''A<sub>n</sub>'' = ''a<sub>n</sub>'' 0 ≤ n < N के लिए और ''A<sub>n</sub>'' = 0 अन्यथा - "शून्य-पैडिंग" का सामान्य अर्थ। चूंकि, कनवल्शन में ''b<sub>k</sub>''<sub>–''n''</sub> पद के कारण, ''b<sub>n</sub>'' के लिए n के धनात्मक और ऋणात्मक (ध्यान दें कि ''b''<sub>–''n''</sub> = ''b<sub>n</sub>'') दोनों मान आवश्यक हैं। शून्य-पैडेड सरणी के डीएफटी द्वारा निहित आवधिक सीमाओं का अर्थ है कि -n ''M''–''n'' के समान है। इस प्रकार, ''b<sub>n</sub>'' को लंबाई M की एक सरणी ''B<sub>n</sub>'' तक विस्तारित किया गया है, जहां ''B''<sub>0</sub> = ''b''<sub>0</sub>, ''B<sub>n</sub>'' = ''B<sub>M</sub>''<sub>–''n''</sub> = ''b<sub>n</sub>'' 0 < n < N के लिए, और ''B<sub>n</sub>'' = 0 अन्यथा। सामान्य कनवल्शन प्रमेय के अनुसार, a और b का कनवल्शन प्राप्त करने के लिए A और B को एफएफटी किया जाता है, बिंदुवार गुणा किया जाता है, और विपरीत एफएफटी किया जाता है। | ||
आइए हम इस बारे में और अधिक | आइए हम इस बारे में और अधिक स्पष्ट हों कि डीएफटी के लिए ब्लूस्टीन के एल्गोरिदम में किस प्रकार के कनवल्शन की आवश्यकता है। यदि अनुक्रम b<sub>''n''</sub> अवधि n के साथ N में आवधिक थे, तो यह लंबाई N का चक्रीय घुमाव होगा, और शून्य-पैडिंग केवल कम्प्यूटेशनल सुविधा के लिए होगी। चूंकि, सामान्यतः ऐसा नहीं होता है: | ||
:<math>b_{n+N} = e^{\frac{\pi i}{N} (n+N)^2 } = b_n \left[ e^{\frac{\pi i}{N} (2Nn+N^2) } \right] = (-1)^N b_n .</math> | :<math>b_{n+N} = e^{\frac{\pi i}{N} (n+N)^2 } = b_n \left[ e^{\frac{\pi i}{N} (2Nn+N^2) } \right] = (-1)^N b_n .</math> | ||
इसलिए, N सम और विषम संख्याओं के लिए कनवल्शन चक्रीय है, | इसलिए, N सम और विषम संख्याओं के लिए कनवल्शन चक्रीय है, किन्तु इस स्तिथि में N [[समग्र संख्या]] है और कोई सामान्यतः कूली-टुकी जैसे अधिक कुशल एफएफटी एल्गोरिदम का उपयोग करेगा। चूंकि, N विषम के लिए, फिर b<sub>''n''</sub> एंटीपेरियोडिक फलन है और हमारे पास तकनीकी रूप से लंबाई N का ऋणात्मक चक्रीय घुमाव है। चूंकि, ऊपर बताए अनुसार एक शून्य-पैड ''a<sub>n</sub>'' को कम से कम 2N−1 की लंबाई तक ले जाने पर ऐसे अंतर विलुप्त हो जाते हैं। इसलिए, इसे एक सरल रैखिक कनवल्शन के आउटपुट के उपसमुच्चय के रूप में विचार कदाचित्स सबसे सरल है (अर्थात डेटा का कोई वैचारिक विस्तार, आवधिक या अन्यथा नहीं) है। | ||
== z-परिवर्तन == | == z-परिवर्तन == | ||
ब्लूस्टीन के एल्गोरिदम का उपयोग (एकतरफा) जेड- | ब्लूस्टीन के एल्गोरिदम का उपयोग (एकतरफा) जेड-रूपांतरण (रेबिनर एट अल।, 1969) के आधार पर अधिक सामान्य परिवर्तन की गणना करने के लिए भी किया जा सकता है। विशेष रूप से, यह प्रपत्र के किसी भी परिवर्तन की गणना कर सकता है: | ||
:<math> X_k = \sum_{n=0}^{N-1} x_n z^{nk} | :<math> X_k = \sum_{n=0}^{N-1} x_n z^{nk} | ||
\qquad | \qquad | ||
k = 0,\dots,M-1, </math> | k = 0,\dots,M-1, </math> | ||
एक | एक इच्छानुसार [[जटिल संख्या|सम्मिश्र संख्या]] z के लिए और इनपुट और आउटपुट की भिन्न संख्या N और M के लिए है। ब्लूस्टीन के एल्गोरिदम को देखते हुए, इस तरह के परिवर्तन का उपयोग किया जा सकता है, उदाहरण के लिए, स्पेक्ट्रम के कुछ भाग के अधिक सूक्ष्म अंतर वाले इंटरपोलेशन को प्राप्त करने के लिए (चूंकि ज़ूम एफएफटी के समान आवृत्ति रिज़ॉल्यूशन अभी भी कुल नमूना समय तक सीमित है), इच्छानुसार रूप से से ध्रुवों को बढ़ाएं स्थानांतरण-फ़ंक्शन विश्लेषण आदि। | ||
एल्गोरिदम को | इस प्रकार से एल्गोरिदम को चिर्प ज़ेड-ट्रांसफ़ॉर्म एल्गोरिदम करार दिया गया था, क्योंकि फूरियर-ट्रांसफ़ॉर्म केस (|z| = 1) के लिए, अनुक्रम b<sub>''n''</sub> ऊपर से रैखिक रूप से बढ़ती आवृत्ति का एक सम्मिश्र साइनसॉइड है, जिसे [[राडार]] सिस्टम में (रैखिक) चिर्प कहा जाता है। | ||
==यह भी देखें== | ==यह भी देखें== | ||
Line 62: | Line 60: | ||
==संदर्भ== | ==संदर्भ== | ||
{{reflist}} | {{reflist}} | ||
===सामान्य=== | ===सामान्य=== | ||
* लियो आई. ब्लूस्टीन, असतत फूरियर | * लियो आई. ब्लूस्टीन, असतत फूरियर रूपांतरण की गणना के लिए एक रैखिक फ़िल्टरिंग दृष्टिकोण, पूर्वोत्तर इलेक्ट्रॉनिक्स अनुसंधान और इंजीनियरिंग मीटिंग रिकॉर्ड '10', 218-219 (1968)। | ||
* लॉरेंस आर. रैबिनर, रोनाल्ड डब्ल्यू. शेफ़र, और चार्ल्स एम. रेडर, [https://web.archive.org/web/20130618204835/http://www3.alcatel-lucent.com/bstj/vol48-1969/ आलेख/bstj48-5-1249.pdf | * लॉरेंस आर. रैबिनर, रोनाल्ड डब्ल्यू. शेफ़र, और चार्ल्स एम. रेडर, [https://web.archive.org/web/20130618204835/http://www3.alcatel-lucent.com/bstj/vol48-1969/ आलेख/bstj48-5-1249.pdf चिर्प ज़ेड-रूपांतरण एल्गोरिदम और उसका अनुप्रयोग], बेल सिस्ट। टेक. जे. '48', 1249-1292 (1969)। इसमें भी प्रकाशित: रैबिनर, शेफर, और रेडर, [https://web.archive.org/web/20150220065735/http://cronos.rutgers.edu/~lrr/Reprints/015_czt.pdf द चिर्प ज़ेड-रूपांतरण एल्गोरिथ्म] , आईईईई ट्रांस। ऑडियो इलेक्ट्रोकॉस्टिक्स '17' (2), 86-92 (1969)। | ||
* डी. एच. बेली और पी. एन. स्वार्जट्रॉबर, द फ्रैक्शनल फूरियर | * डी. एच. बेली और पी. एन. स्वार्जट्रॉबर, द फ्रैक्शनल फूरियर रूपांतरण एंड एप्लिकेशन, [[ सियाम समीक्षा |सियाम समीक्षा]] '33', 389-404 (1991)। (ध्यान दें कि z-परिवर्तन के लिए यह शब्दावली गैरमानक है: एक भिन्नात्मक फूरियर परिवर्तन पारंपरिक रूप से एक पूरी तरह से अलग, निरंतर परिवर्तन को संदर्भित करता है।) | ||
* लॉरेंस रैबिनर, द चिर ज़ेड-ट्रांसफ़ॉर्म एल्गोरिदम-सेरेन्डिपिटी में एक पाठ, आईईईई सिग्नल प्रोसेसिंग पत्रिका '21', 118-119 (मार्च 2004)। (ऐतिहासिक टिप्पणी।) | * लॉरेंस रैबिनर, द चिर ज़ेड-ट्रांसफ़ॉर्म एल्गोरिदम-सेरेन्डिपिटी में एक पाठ, आईईईई सिग्नल प्रोसेसिंग पत्रिका '21', 118-119 (मार्च 2004)। (ऐतिहासिक टिप्पणी।) | ||
* व्लादिमीर सुखॉय और अलेक्जेंडर स्टॉयचेव: [https://www.nature.com/articles/s41598-019-50234-9.pdf यूनिट सर्कल से व्युत्क्रम एफएफटी को सामान्य बनाना], (अक्टूबर 2019)। # खुला एक्सेस। | * व्लादिमीर सुखॉय और अलेक्जेंडर स्टॉयचेव: [https://www.nature.com/articles/s41598-019-50234-9.pdf यूनिट सर्कल से व्युत्क्रम एफएफटी को सामान्य बनाना], (अक्टूबर 2019)। # खुला एक्सेस। | ||
Line 83: | Line 78: | ||
[[Category: Machine Translated Page]] | [[Category: Machine Translated Page]] | ||
[[Category:Created On 16/08/2023]] | [[Category:Created On 16/08/2023]] | ||
[[Category:Vigyan Ready]] |
Latest revision as of 09:54, 1 December 2023
चिर्प जेड-रूपांतरण (सीजेडटी) असतत फूरियर रूपांतरण (डीएफटी) का एक सामान्यीकरण है। जबकि डीएफटी यूनिट सर्कल के साथ समान रूप से दूरी वाले बिंदुओं पर जेड-रूपांतरण का नमूना लेता है, एस प्लेन में सीधी रेखाओं के अनुरूप, जेड-प्लेन में सर्पिल आर्क्स के साथ चहचहाना जेड-रूपांतरण नमूने लेता है।[1][2] इस प्रकार से डीएफटी, वास्तविक डीएफटी और ज़ूम डीएफटी की गणना सीजेडटी के विशेष स्तिथियों के रूप में की जा सकती है।
विशेष रूप से, चिर्प Z रूपांतरण, Z रूपांतरण की गणना बिंदु zk की एक सीमित संख्या पर करता है एक लघुगणकीय सर्पिल समोच्च के साथ, इसे इस प्रकार परिभाषित किया गया है:[1][3]
जहां A सम्मिश्र प्रारंभिक बिंदु है, W बिंदुओं के बीच सम्मिश्र अनुपात है, और M गणना किए जाने वाले बिंदुओं की संख्या है।
डीएफटी की तरह, चिर्प जेड-रूपांतरण की गणना O(n log n) ऑपरेशंस में की जा सकती है जहां n=\max(M,N) है।
व्युत्क्रम चिर्प जेड-रूपांतरण (आईसीजेडटी) के लिए एक O(n log n) एल्गोरिदम 2003 और 2019 में वर्णित किया गया था।[4][5] [6]
ब्लूस्टीन का एल्गोरिदम
ब्लूस्टीन का एल्गोरिदम[7][8] सीजेडटी को एक कनवल्शन के रूप में व्यक्त करता है और तेज़ फूरियर रूपांतरण/आईएफएफटी का उपयोग करके इसे कुशलतापूर्वक कार्यान्वित करता है।
चूंकि डीएफटी सीजेडटी का विशेष स्तिथि है, यह अभाज्य संख्या आकार सहित इच्छानुसार आकारों के फास्ट फूरियर रूपांतरण (डीएफटी) की कुशल गणना की अनुमति देता है। इस प्रकार से (प्राइम साइज के एफएफटी के लिए अन्य एल्गोरिदम, रेडर का एल्गोरिदम, डीएफटी को कनवल्शन के रूप में फिर से लिखकर भी कार्य करता है।) इसकी कल्पना 1968 में लियो ब्लूस्टीन द्वारा की गई थी।[7] किन्तु ब्लूस्टीन के एल्गोरिदम का उपयोग (एकतरफा) जेड-रूपांतरण (रेबिनर एट अल।, 1969) के आधार पर डीएफटी की तुलना में अधिक सामान्य परिवर्तनों की गणना करने के लिए किया जा सकता है।
याद रखें कि डीएफटी को सूत्र द्वारा परिभाषित किया गया है
यदि हम घातांक में गुणनफल nk को पहचान से प्रतिस्थापित करते हैं
हम इस प्रकार प्राप्त करते हैं:
यह योग वास्तव में दो अनुक्रमों a और bn का एक संलयन है जिसे निम्न द्वारा परिभाषित किया गया है:
कनवल्शन के आउटपुट को N चरण कारकों bk* से गुणा किया जाता है। वह है:
यह कनवल्शन, परिवर्तन में, कनवल्शन प्रमेय के माध्यम से एफएफटी की जोड़ी (साथ ही सम्मिश्र कलरव bn की पूर्व-गणना की गई एफएफटी)) के साथ किया जा सकता है । इस प्रकार से मुख्य तथ्य यह है कि ये एफएफटी समान लंबाई N के नहीं हैं: इस तरह के कनवल्शन की गणना एफएफटी से केवल शून्य-पैडिंग द्वारा 2N-1 से अधिक या उसके समान लंबाई तक की जा सकती है। विशेष रूप से, कोई दो या किसी अन्य अत्यधिक समग्र संख्या आकार की शक्ति तक पैड कर सकता है, कूली-टुकी एफएफटी एल्गोरिथ्म O(N log N) समय में कूली-टुकी एल्गोरिथ्म जिसके लिए एफएफटी को कुशलतापूर्वक निष्पादित किया जा सकता है। इस प्रकार, ब्लूस्टीन का एल्गोरिदम प्राइम-आकार डीएफटी की गणना करने के लिए एक O(N log N) विधि प्रदान करता है, चूंकि समग्र आकार के लिए कूली-टुकी एल्गोरिदम की तुलना में कई गुना धीमा है।
इस प्रकार से ब्लूस्टीन के एल्गोरिदम में कनवल्शन के लिए शून्य-पैडिंग का उपयोग कुछ अतिरिक्त टिप्पणी का पात्र है। मान लीजिए कि हम लंबाई M ≥ 2N–1 पर शून्य-पैड करते हैं। इसका अर्थ यह है कि an को लंबाई M की एक सरणी An तक विस्तारित किया गया है, जहां An = an 0 ≤ n < N के लिए और An = 0 अन्यथा - "शून्य-पैडिंग" का सामान्य अर्थ। चूंकि, कनवल्शन में bk–n पद के कारण, bn के लिए n के धनात्मक और ऋणात्मक (ध्यान दें कि b–n = bn) दोनों मान आवश्यक हैं। शून्य-पैडेड सरणी के डीएफटी द्वारा निहित आवधिक सीमाओं का अर्थ है कि -n M–n के समान है। इस प्रकार, bn को लंबाई M की एक सरणी Bn तक विस्तारित किया गया है, जहां B0 = b0, Bn = BM–n = bn 0 < n < N के लिए, और Bn = 0 अन्यथा। सामान्य कनवल्शन प्रमेय के अनुसार, a और b का कनवल्शन प्राप्त करने के लिए A और B को एफएफटी किया जाता है, बिंदुवार गुणा किया जाता है, और विपरीत एफएफटी किया जाता है।
आइए हम इस बारे में और अधिक स्पष्ट हों कि डीएफटी के लिए ब्लूस्टीन के एल्गोरिदम में किस प्रकार के कनवल्शन की आवश्यकता है। यदि अनुक्रम bn अवधि n के साथ N में आवधिक थे, तो यह लंबाई N का चक्रीय घुमाव होगा, और शून्य-पैडिंग केवल कम्प्यूटेशनल सुविधा के लिए होगी। चूंकि, सामान्यतः ऐसा नहीं होता है:
इसलिए, N सम और विषम संख्याओं के लिए कनवल्शन चक्रीय है, किन्तु इस स्तिथि में N समग्र संख्या है और कोई सामान्यतः कूली-टुकी जैसे अधिक कुशल एफएफटी एल्गोरिदम का उपयोग करेगा। चूंकि, N विषम के लिए, फिर bn एंटीपेरियोडिक फलन है और हमारे पास तकनीकी रूप से लंबाई N का ऋणात्मक चक्रीय घुमाव है। चूंकि, ऊपर बताए अनुसार एक शून्य-पैड an को कम से कम 2N−1 की लंबाई तक ले जाने पर ऐसे अंतर विलुप्त हो जाते हैं। इसलिए, इसे एक सरल रैखिक कनवल्शन के आउटपुट के उपसमुच्चय के रूप में विचार कदाचित्स सबसे सरल है (अर्थात डेटा का कोई वैचारिक विस्तार, आवधिक या अन्यथा नहीं) है।
z-परिवर्तन
ब्लूस्टीन के एल्गोरिदम का उपयोग (एकतरफा) जेड-रूपांतरण (रेबिनर एट अल।, 1969) के आधार पर अधिक सामान्य परिवर्तन की गणना करने के लिए भी किया जा सकता है। विशेष रूप से, यह प्रपत्र के किसी भी परिवर्तन की गणना कर सकता है:
एक इच्छानुसार सम्मिश्र संख्या z के लिए और इनपुट और आउटपुट की भिन्न संख्या N और M के लिए है। ब्लूस्टीन के एल्गोरिदम को देखते हुए, इस तरह के परिवर्तन का उपयोग किया जा सकता है, उदाहरण के लिए, स्पेक्ट्रम के कुछ भाग के अधिक सूक्ष्म अंतर वाले इंटरपोलेशन को प्राप्त करने के लिए (चूंकि ज़ूम एफएफटी के समान आवृत्ति रिज़ॉल्यूशन अभी भी कुल नमूना समय तक सीमित है), इच्छानुसार रूप से से ध्रुवों को बढ़ाएं स्थानांतरण-फ़ंक्शन विश्लेषण आदि।
इस प्रकार से एल्गोरिदम को चिर्प ज़ेड-ट्रांसफ़ॉर्म एल्गोरिदम करार दिया गया था, क्योंकि फूरियर-ट्रांसफ़ॉर्म केस (|z| = 1) के लिए, अनुक्रम bn ऊपर से रैखिक रूप से बढ़ती आवृत्ति का एक सम्मिश्र साइनसॉइड है, जिसे राडार सिस्टम में (रैखिक) चिर्प कहा जाता है।
यह भी देखें
संदर्भ
- ↑ 1.0 1.1 A study of the Chirp Z-transform and its applications - Shilling, Steve Alan
- ↑ "चिरप जेड-ट्रांसफॉर्म - MATLAB czt". www.mathworks.com. Retrieved 2016-09-22.
- ↑ Martin, Grant D. (November 2005). "Chirp Z-Transform Spectral Zoom Optimization with MATLAB®" (PDF).
- ↑ Bostan, Alin (2003). बुनियादी कंप्यूटर बीजगणित संचालन के लिए कुशल एल्गोरिदम (PDF) (PhD). Ecole polytechnique.
- ↑ Bostan, Alin; Schost, Éric (2005). "बिंदुओं के विशेष सेट पर बहुपद मूल्यांकन और प्रक्षेप". Journal of Complexity (in English). 21 (4): 420–446. doi:10.1016/j.jco.2004.09.009.
- ↑ Engineers Solve 50-Year-Old Puzzle in Signal Processing – Inverse Chirp Z-Transform, By IOWA STATE UNIVERSITY OCTOBER 10, 2019
- ↑ 7.0 7.1 Bluestein, L. (1970-12-01). "असतत फूरियर रूपांतरण की गणना के लिए एक रैखिक फ़िल्टरिंग दृष्टिकोण". IEEE Transactions on Audio and Electroacoustics. 18 (4): 451–455. doi:10.1109/TAU.1970.1162132. ISSN 0018-9278.
- ↑ "ब्लूस्टीन का एफएफटी एल्गोरिदम". DSPRelated.com.
सामान्य
- लियो आई. ब्लूस्टीन, असतत फूरियर रूपांतरण की गणना के लिए एक रैखिक फ़िल्टरिंग दृष्टिकोण, पूर्वोत्तर इलेक्ट्रॉनिक्स अनुसंधान और इंजीनियरिंग मीटिंग रिकॉर्ड '10', 218-219 (1968)।
- लॉरेंस आर. रैबिनर, रोनाल्ड डब्ल्यू. शेफ़र, और चार्ल्स एम. रेडर, आलेख/bstj48-5-1249.pdf चिर्प ज़ेड-रूपांतरण एल्गोरिदम और उसका अनुप्रयोग, बेल सिस्ट। टेक. जे. '48', 1249-1292 (1969)। इसमें भी प्रकाशित: रैबिनर, शेफर, और रेडर, द चिर्प ज़ेड-रूपांतरण एल्गोरिथ्म , आईईईई ट्रांस। ऑडियो इलेक्ट्रोकॉस्टिक्स '17' (2), 86-92 (1969)।
- डी. एच. बेली और पी. एन. स्वार्जट्रॉबर, द फ्रैक्शनल फूरियर रूपांतरण एंड एप्लिकेशन, सियाम समीक्षा '33', 389-404 (1991)। (ध्यान दें कि z-परिवर्तन के लिए यह शब्दावली गैरमानक है: एक भिन्नात्मक फूरियर परिवर्तन पारंपरिक रूप से एक पूरी तरह से अलग, निरंतर परिवर्तन को संदर्भित करता है।)
- लॉरेंस रैबिनर, द चिर ज़ेड-ट्रांसफ़ॉर्म एल्गोरिदम-सेरेन्डिपिटी में एक पाठ, आईईईई सिग्नल प्रोसेसिंग पत्रिका '21', 118-119 (मार्च 2004)। (ऐतिहासिक टिप्पणी।)
- व्लादिमीर सुखॉय और अलेक्जेंडर स्टॉयचेव: यूनिट सर्कल से व्युत्क्रम एफएफटी को सामान्य बनाना, (अक्टूबर 2019)। # खुला एक्सेस।
- व्लादिमीर सुखॉय और अलेक्जेंडर स्टोयचेव: यूनिट सर्कल पर चिर आकृति के लिए आईसीजेडटी एल्गोरिदम का संख्यात्मक त्रुटि विश्लेषण, विज्ञान प्रतिनिधि 10, 4852 (2020)।
बाहरी संबंध
- A DSP algorithm for frequency analysis - the Chirp-Z Transform (CZT)
- Solving a 50-year-old puzzle in signal processing, part two