पुनरावर्ती न्यूनतम वर्ग फ़िल्टर

पुनरावर्ती न्यूनतम वर्ग (आरएलएस) एक अनुकूली फ़िल्टर एल्गोरिथ्म है जो पुनरावर्ती रूप से उन गुणांकों को ढूंढता है जो इनपुट संकेत से संबंधित भारित रैखिक न्यूनतम वर्ग लागत फ़ंक्शन को कम करते हैं। यह दृष्टिकोण अन्य एल्गोरिदम जैसे कि न्यूनतम माध्य वर्ग (एलएमएस) के विपरीत है जिसका लक्ष्य माध्य वर्ग त्रुटि को कम करना है। आरएलएस की व्युत्पत्ति में, इनपुट संकेतों को नियतात्मक माना जाता है, जबकि एलएमएस और इसी तरह के एल्गोरिदम के लिए उन्हें स्टोकेस्टिक माना जाता है। अपने अधिकांश प्रतिस्पर्धियों की तुलना में, आरएलएस अत्यंत तीव्र अभिसरण प्रदर्शित करता है। हालाँकि, यह लाभ उच्च कम्प्यूटेशनल सम्मिश्रता की कीमत पर मिलता है।

प्रेरणा

आरएलएस की खोज गॉस ने की थी, लेकिन 1950 तक अप्रयुक्त या नजरअंदाज कर दिया गया था, जब प्लैकेट ने 1821 में गॉस के मूल कार्य को फिर से खोजा था। सामान्य तौर पर, आरएलएस का उपयोग किसी भी समस्या को हल करने के लिए किया जा सकता है जिसे अनुकूली फिल्टर द्वारा हल किया जा सकता है। उदाहरण के लिए, मान लीजिए कि एक संकेत $d(n)$ एक प्रतिध्वनि वाले, ध्वनि वाले चैनल पर प्रसारित होता है जिसके कारण इसे प्राप्त किया जाता है

$x(n)=\sum _{k=0}^{q}b_{n}(k)d(n-k)+v(n)$

जहां $v(n)$ योगात्मक शोर का प्रतिनिधित्व करता है। आरएलएस फ़िल्टर का उद्देश्य $p+1$-टैप FIR फ़िल्टर के उपयोग से वांछित संकेत $d(n)$ को पुनर्प्राप्त करना है, $\mathbf {w}$:

$d(n)\approx \sum _{k=0}^{p}w(k)x(n-k)=\mathbf {w} ^{\mathit {T}}\mathbf {x} _{n}$

जहां $\mathbf {x} _{n}=[x(n)\quad x(n-1)\quad \ldots \quad x(n-p)]^{T}$ कॉलम सदिश है जिसमें $x(n)$ के सबसे हाल के नमूने $p+1$ शामिल हैं। प्राप्त वांछित संकेत का अनुमान है

${\hat {d}}(n)=\sum _{k=0}^{p}w_{n}(k)x(n-k)=\mathbf {w} _{n}^{\mathit {T}}\mathbf {x} _{n}$

फ़िल्टर $\mathbf {w}$ के मापदंडों का अनुमान लगाना है, और प्रत्येक समय $n$ पर हम वर्तमान अनुमान को $\mathbf {w} _{n}$और अनुकूलित न्यूनतम-वर्ग अनुमान को $\mathbf {w} _{n+1}$के रूप में संदर्भित करते हैं। जैसा कि नीचे दिखाया गया है, $\mathbf {w} _{n}$भी एक कॉलम सदिश है, और ट्रांसपोज़, $\mathbf {w} _{n}^{\mathit {T}}$ एक पंक्ति सदिश है। मैट्रिक्स गुणनफल $\mathbf {w} _{n}^{\mathit {T}}\mathbf {x} _{n}$ (जो कि डॉट गुणनफल है $\mathbf {w} _{n}$और $\mathbf {x} _{n}$${\hat {d}}(n)$, एक अदिश राशि है। अनुमान "अच्छा" है यदि ${\hat {d}}(n)-d(n)$ कुछ न्यूनतम वर्ग अर्थों में परिमाण में छोटा है।

जैसे-जैसे समय बढ़ता है, $\mathbf {w} _{n+1}$के लिए नए अनुमान को खोजने के लिए न्यूनतम वर्ग एल्गोरिथ्म को पूरी तरह से दोबारा करने से बचना चाहिए, $\mathbf {w} _{n}$ के संदर्भ में।

आरएलएस एल्गोरिदम का लाभ यह है कि मैट्रिक्स को उलटने की कोई आवश्यकता नहीं है, जिससे कम्प्यूटेशनल लागत बचती है। एक अन्य लाभ यह है कि यह कलमन फ़िल्टर जैसे परिणामों के पीछे अंतर्ज्ञान प्रदान करता है।

चर्चा

आरएलएस फ़िल्टर के पीछे का विचार फ़िल्टर गुणांक $\mathbf {w} _{n}$ का उचित चयन करके और नए डेटा आने पर फ़िल्टर को अपडेट करके लागत फ़ंक्शन $C$ को कम करना है। त्रुटि संकेत $e(n)$ और वांछित सिग्नल $d(n)$ को नीचे ऋणात्मक प्रतिक्रिया आरेख में परिभाषित किया गया है:

त्रुटि अनुमान के माध्यम से फ़िल्टर गुणांक ${\hat {d}}(n)$ पर परोक्ष रूप से निर्भर करती है:

$e(n)=d(n)-{\hat {d}}(n)$

भारित न्यूनतम वर्ग त्रुटि फ़ंक्शन $C$- जिस लागत फ़ंक्शन को हम कम करना चाहते हैं - वह एक फ़ंक्शन $e(n)$ है इसलिए फ़िल्टर गुणांक पर भी निर्भर है:

$C(\mathbf {w} _{n})=\sum _{i=0}^{n}\lambda ^{n-i}e^{2}(i)$

जहाँ $0<\lambda \leq 1$ "विस्मृति कारक" है जो पुराने त्रुटि नमूनों को तेजी से कम महत्व देता है।

सभी प्रविष्टियों के लिए आंशिक व्युत्पन्न लेकर लागत फ़ंक्शन को कम किया जाता है $k$ गुणांक सदिश का $\mathbf {w} _{n}$ और परिणामों को शून्य पर सेट करना

${\frac {\partial C(\mathbf {w} _{n})}{\partial w_{n}(k)}}=\sum _{i=0}^{n}2\lambda ^{n-i}e(i)\cdot {\frac {\partial e(i)}{\partial w_{n}(k)}}=-\sum _{i=0}^{n}2\lambda ^{n-i}e(i)\,x(i-k)=0\qquad k=0,1,\ldots ,p$

अगला, बदलें $e(n)$ त्रुटि संकेत की परिभाषा के साथ

$\sum _{i=0}^{n}\lambda ^{n-i}\left[d(i)-\sum _{\ell =0}^{p}w_{n}(\ell )x(i-\ell )\right]x(i-k)=0\qquad k=0,1,\ldots ,p$

समीकरण को पुनर्व्यवस्थित करने से परिणाम प्राप्त होते हैं

$\sum _{\ell =0}^{p}w_{n}(\ell )\left[\sum _{i=0}^{n}\lambda ^{n-i}\,x(i-l)x(i-k)\right]=\sum _{i=0}^{n}\lambda ^{n-i}d(i)x(i-k)\qquad k=0,1,\ldots ,p$

इस रूप को मैट्रिक्स के संदर्भ में व्यक्त किया जा सकता है

$\mathbf {R} _{x}(n)\,\mathbf {w} _{n}=\mathbf {r} _{dx}(n)$

जहाँ $\mathbf {R} _{x}(n)$ के लिए भारित नमूना माध्य और नमूना सहप्रसरण मैट्रिक्स है $x(n)$, और $\mathbf {r} _{dx}(n)$ के बीच क्रॉस-सहप्रसरण के लिए समतुल्य अनुमान है $d(n)$ और $x(n)$. इस अभिव्यक्ति के आधार पर हम ऐसे गुणांक पाते हैं जो लागत फ़ंक्शन को कम करते हैं

$\mathbf {w} _{n}=\mathbf {R} _{x}^{-1}(n)\,\mathbf {r} _{dx}(n)$

यह चर्चा का मुख्य परिणाम है.

=== λ=== चुनना छोटे $\lambda$ है, सहप्रसरण मैट्रिक्स में पिछले नमूनों का योगदान उतना ही छोटा है। यह फ़िल्टर को हाल के नमूनों के प्रति अधिक संवेदनशील बनाता है, जिसका अर्थ है फ़िल्टर गुणांक में अधिक उतार-चढ़ाव। $\lambda =1$ h> केस को ग्रोइंग विंडो आरएलएस एल्गोरिदम के रूप में जाना जाता है। व्यवहार में, $\lambda$ आमतौर पर 0.98 और 1 के बीच चुना जाता है। टाइप- II अधिकतम संभावना अनुमान का उपयोग करके इष्टतम $\lambda$ डेटा के एक सेट से अनुमान लगाया जा सकता है।

पुनरावर्ती एल्गोरिथ्म

चर्चा के परिणामस्वरूप गुणांक सदिश निर्धारित करने के लिए एक एकल समीकरण तैयार हुआ जो लागत फ़ंक्शन को न्यूनतम करता है। इस अनुभाग में हम प्रपत्र का पुनरावर्ती समाधान प्राप्त करना चाहते हैं

$\mathbf {w} _{n}=\mathbf {w} _{n-1}+\Delta \mathbf {w} _{n-1}$

जहाँ $\Delta \mathbf {w} _{n-1}$ समय पर एक सुधार कारक है ${n-1}$. हम क्रॉस सहप्रसरण को व्यक्त करके पुनरावर्ती एल्गोरिदम की व्युत्पत्ति शुरू करते हैं $\mathbf {r} _{dx}(n)$ के अनुसार $\mathbf {r} _{dx}(n-1)$

$\mathbf {r} _{dx}(n)$	$=\sum _{i=0}^{n}\lambda ^{n-i}d(i)\mathbf {x} (i)$
	$=\sum _{i=0}^{n-1}\lambda ^{n-i}d(i)\mathbf {x} (i)+\lambda ^{0}d(n)\mathbf {x} (n)$
	$=\lambda \mathbf {r} _{dx}(n-1)+d(n)\mathbf {x} (n)$

जहाँ $\mathbf {x} (i)$ है ${p+1}$ आयामी डेटा सदिश

$\mathbf {x} (i)=[x(i),x(i-1),\dots ,x(i-p)]^{T}$

वैसे ही हम व्यक्त करते हैं $\mathbf {R} _{x}(n)$ के अनुसार $\mathbf {R} _{x}(n-1)$ द्वारा

$\mathbf {R} _{x}(n)$	$=\sum _{i=0}^{n}\lambda ^{n-i}\mathbf {x} (i)\mathbf {x} ^{T}(i)$
	$=\lambda \mathbf {R} _{x}(n-1)+\mathbf {x} (n)\mathbf {x} ^{T}(n)$

गुणांक सदिश उत्पन्न करने के लिए हम नियतात्मक ऑटो-सहप्रसरण मैट्रिक्स के व्युत्क्रम में रुचि रखते हैं। उस कार्य के लिए वुडबरी मैट्रिक्स पहचान काम आती है। साथ

$A$	$=\lambda \mathbf {R} _{x}(n-1)$ is $(p+1)$-by-$(p+1)$
$U$	$=\mathbf {x} (n)$ is $(p+1)$-by-1 (column vector)
$V$	$=\mathbf {x} ^{T}(n)$ is 1-by-$(p+1)$ (row vector)
$C$	$=\mathbf {I} _{1}$ is the 1-by-1 identity matrix

वुडबरी मैट्रिक्स पहचान इस प्रकार है

$\mathbf {R} _{x}^{-1}(n)$	$=$	$\left[\lambda \mathbf {R} _{x}(n-1)+\mathbf {x} (n)\mathbf {x} ^{T}(n)\right]^{-1}$
	$=$	$\lambda ^{-1}\mathbf {R} _{x}^{-1}(n-1)$
		$-\lambda ^{-1}\mathbf {R} _{x}^{-1}(n-1)\mathbf {x} (n)$
		$\left\{1+\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {R} _{x}^{-1}(n-1)\mathbf {x} (n)\right\}^{-1}\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {R} _{x}^{-1}(n-1)$

मानक साहित्य के अनुरूप आने के लिए, हम परिभाषित करते हैं

$\mathbf {P} (n)$	$=\mathbf {R} _{x}^{-1}(n)$
	$=\lambda ^{-1}\mathbf {P} (n-1)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)$

जहां लाभ सदिश $g(n)$ है

$\mathbf {g} (n)$	$=\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)\left\{1+\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)\right\}^{-1}$
	$=\mathbf {P} (n-1)\mathbf {x} (n)\left\{\lambda +\mathbf {x} ^{T}(n)\mathbf {P} (n-1)\mathbf {x} (n)\right\}^{-1}$

आगे बढ़ने से पहले ये लाना जरूरी है $\mathbf {g} (n)$ दूसरे रूप में

$\mathbf {g} (n)\left\{1+\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)\right\}$	$=\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)$
$\mathbf {g} (n)+\mathbf {g} (n)\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)$	$=\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)$

बायीं ओर का दूसरा पद घटाने पर प्राप्त होता है

$\mathbf {g} (n)$	$=\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)\mathbf {x} (n)$
	$=\lambda ^{-1}\left[\mathbf {P} (n-1)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\mathbf {P} (n-1)\right]\mathbf {x} (n)$

की पुनरावर्ती परिभाषा के साथ $\mathbf {P} (n)$ वांछित प्रपत्र इस प्रकार है

$\mathbf {g} (n)=\mathbf {P} (n)\mathbf {x} (n)$

अब हम रिकर्सन पूरा करने के लिए तैयार हैं। चर्चा के अनुसार

$\mathbf {w} _{n}$	$=\mathbf {P} (n)\,\mathbf {r} _{dx}(n)$
	$=\lambda \mathbf {P} (n)\,\mathbf {r} _{dx}(n-1)+d(n)\mathbf {P} (n)\,\mathbf {x} (n)$

दूसरा चरण की पुनरावर्ती परिभाषा से अनुसरण करता है $\mathbf {r} _{dx}(n)$. आगे हम की पुनरावर्ती परिभाषा को शामिल करते हैं $\mathbf {P} (n)$ के वैकल्पिक रूप के साथ $\mathbf {g} (n)$ और पाओ

$\mathbf {w} _{n}$	$=\lambda \left[\lambda ^{-1}\mathbf {P} (n-1)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)\right]\mathbf {r} _{dx}(n-1)+d(n)\mathbf {g} (n)$
	$=\mathbf {P} (n-1)\mathbf {r} _{dx}(n-1)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\mathbf {P} (n-1)\mathbf {r} _{dx}(n-1)+d(n)\mathbf {g} (n)$
	$=\mathbf {P} (n-1)\mathbf {r} _{dx}(n-1)+\mathbf {g} (n)\left[d(n)-\mathbf {x} ^{T}(n)\mathbf {P} (n-1)\mathbf {r} _{dx}(n-1)\right]$

साथ $\mathbf {w} _{n-1}=\mathbf {P} (n-1)\mathbf {r} _{dx}(n-1)$ हम अद्यतन समीकरण पर पहुँचते हैं

$\mathbf {w} _{n}$	$=\mathbf {w} _{n-1}+\mathbf {g} (n)\left[d(n)-\mathbf {x} ^{T}(n)\mathbf {w} _{n-1}\right]$
	$=\mathbf {w} _{n-1}+\mathbf {g} (n)\alpha (n)$

जहाँ $\alpha (n)=d(n)-\mathbf {x} ^{T}(n)\mathbf {w} _{n-1}$ एक प्राथमिकता और एक पश्चवर्ती त्रुटि है। इसकी तुलना पिछली त्रुटि से करें; फ़िल्टर अद्यतन होने के बाद गणना की गई त्रुटि:

$e(n)=d(n)-\mathbf {x} ^{T}(n)\mathbf {w} _{n}$

इसका मतलब है कि हमें सुधार कारक मिल गया है

$\Delta \mathbf {w} _{n-1}=\mathbf {g} (n)\alpha (n)$

यह सहज रूप से संतोषजनक परिणाम इंगित करता है कि सुधार कारक त्रुटि और लाभ सदिश दोनों के लिए सीधे आनुपातिक है, जो भार कारक के माध्यम से नियंत्रित करता है कि कितनी संवेदनशीलता वांछित है, $\lambda$.

आरएलएस एल्गोरिदम सारांश

पी-वें ऑर्डर आरएलएस फ़िल्टर के लिए आरएलएस एल्गोरिदम को संक्षेप में प्रस्तुत किया जा सकता है

Parameters:	$p=$ filter order
	$\lambda =$ forgetting factor
	$\delta =$ value to initialize $\mathbf {P} (0)$
Initialization:	$\mathbf {w} (0)=0$,
	$x(k)=0,k=-p,\dots ,-1$,
	$d(k)=0,k=-p,\dots ,-1$
	$\mathbf {P} (0)=\delta I$ where $I$ is the identity matrix of rank $p+1$
Computation:	For $n=1,2,\dots$
	$\mathbf {x} (n)=\left[{\begin{matrix}x(n)\\x(n-1)\\\vdots \\x(n-p)\end{matrix}}\right]$
	$\alpha (n)=d(n)-\mathbf {x} ^{T}(n)\mathbf {w} (n-1)$
	$\mathbf {g} (n)=\mathbf {P} (n-1)\mathbf {x} (n)\left\{\lambda +\mathbf {x} ^{T}(n)\mathbf {P} (n-1)\mathbf {x} (n)\right\}^{-1}$
	$\mathbf {P} (n)=\lambda ^{-1}\mathbf {P} (n-1)-\mathbf {g} (n)\mathbf {x} ^{T}(n)\lambda ^{-1}\mathbf {P} (n-1)$
	$\mathbf {w} (n)=\mathbf {w} (n-1)+\,\alpha (n)\mathbf {g} (n)$.

के लिए प्रत्यावर्तन $P$ बीजगणितीय रिकाटी समीकरण का अनुसरण करता है और इस प्रकार कलमन फिल्टर के समानांतर खींचता है।

जाली पुनरावर्ती न्यूनतम वर्ग फ़िल्टर (एलआरएलएस)

जाली पुनरावर्ती न्यूनतम वर्ग अनुकूली फ़िल्टर मानक आरएलएस से संबंधित है, सिवाय इसके कि इसके लिए कम अंकगणितीय संचालन (आदेश एन) की आवश्यकता होती है। यह पारंपरिक एलएमएस एल्गोरिदम पर अतिरिक्त लाभ प्रदान करता है जैसे कि तेज अभिसरण दर, मॉड्यूलर संरचना, और इनपुट सहसंबंध मैट्रिक्स के ईजेनवैल्यू प्रसार में भिन्नता के प्रति असंवेदनशीलता। वर्णित एलआरएलएस एल्गोरिदम पिछली त्रुटियों पर आधारित है और इसमें सामान्यीकृत फॉर्म शामिल है। व्युत्पत्ति मानक आरएलएस एल्गोरिथ्म के समान है और की परिभाषा पर आधारित है $d(k)\,\!$. आगे की भविष्यवाणी के मामले में, हमारे पास है $d(k)=x(k)\,\!$ इनपुट संकेत के साथ $x(k-1)\,\!$ सबसे अद्यतित नमूने के रूप में। पिछड़ी भविष्यवाणी का मामला है $d(k)=x(k-i-1)\,\!$, जहां i अतीत में नमूने का सूचकांक है जिसकी हम भविष्यवाणी करना चाहते हैं, और इनपुट संकेत $x(k)\,\!$ सबसे ताज़ा नमूना है.

पैरामीटर सारांश

$\kappa _{f}(k,i)\,\!$ अग्र परावर्तन गुणांक है

$\kappa _{b}(k,i)\,\!$ पश्चगामी परावर्तन गुणांक है

$e_{f}(k,i)\,\!$ तात्कालिक पूर्ववर्ती अग्रगामी भविष्यवाणी त्रुटि का प्रतिनिधित्व करता है

$e_{b}(k,i)\,\!$ तात्कालिक पश्चगामी पूर्वानुमान त्रुटि का प्रतिनिधित्व करता है

$\xi _{b_{\min }}^{d}(k,i)\,\!$ न्यूनतम न्यूनतम-वर्ग पिछड़ा पूर्वानुमान त्रुटि है

$\xi _{f_{\min }}^{d}(k,i)\,\!$ न्यूनतम न्यूनतम-वर्ग अग्रेषित पूर्वानुमान त्रुटि है

$\gamma (k,i)\,\!$ प्राथमिक और पश्चवर्ती त्रुटियों के बीच एक रूपांतरण कारक है

$v_{i}(k)\,\!$ फीडफॉरवर्ड गुणक गुणांक हैं।

$\varepsilon \,\!$ एक छोटा धनात्मक स्थिरांक है जो 0.01 हो सकता है

एलआरएलएस एल्गोरिदम सारांश

एलआरएलएस फ़िल्टर के लिए एल्गोरिदम को संक्षेप में प्रस्तुत किया जा सकता है

Initialization:
	For ${\textstyle i=0,1,\ldots ,N}$
	$\delta (-1,i)=\delta _{D}(-1,i)=0\,\!$ (if ${\textstyle x(k)=0}$ for ${\textstyle k<0}$)
	$\xi _{b_{\min }}^{d}(-1,i)=\xi _{f_{\min }}^{d}(-1,i)=\varepsilon$
	$\gamma (-1,i)=1\,\!$
	$e_{b}(-1,i)=0\,\!$
	End
Computation:
	For ${\textstyle k\geq 0}$
	$\gamma (k,0)=1\,\!$
	$e_{b}(k,0)=e_{f}(k,0)=x(k)\,\!$
	$\xi _{b_{\min }}^{d}(k,0)=\xi _{f_{\min }}^{d}(k,0)=x^{2}(k)+\lambda \xi _{f_{\min }}^{d}(k-1,0)\,\!$
	$e(k,0)=d(k)\,\!$
	For ${\textstyle i=0,1,\ldots ,N}$
	$\delta (k,i)=\lambda \delta (k-1,i)+{\frac {e_{b}(k-1,i)e_{f}(k,i)}{\gamma (k-1,i)}}$
	$\gamma (k,i+1)=\gamma (k,i)-{\frac {e_{b}^{2}(k,i)}{\xi _{b_{\min }}^{d}(k,i)}}$
	$\kappa _{b}(k,i)={\frac {\delta (k,i)}{\xi _{f_{\min }}^{d}(k,i)}}$
	$\kappa _{f}(k,i)={\frac {\delta (k,i)}{\xi _{b_{\min }}^{d}(k-1,i)}}$
	$e_{b}(k,i+1)=e_{b}(k-1,i)-\kappa _{b}(k,i)e_{f}(k,i)\,\!$
	$e_{f}(k,i+1)=e_{f}(k,i)-\kappa _{f}(k,i)e_{b}(k-1,i)\,\!$
	$\xi _{b_{\min }}^{d}(k,i+1)=\xi _{b_{\min }}^{d}(k-1,i)-\delta (k,i)\kappa _{b}(k,i)$
	$\xi _{f_{\min }}^{d}(k,i+1)=\xi _{f_{\min }}^{d}(k,i)-\delta (k,i)\kappa _{f}(k,i)$
	Feedforward filtering
	$\delta _{D}(k,i)=\lambda \delta _{D}(k-1,i)+{\frac {e(k,i)e_{b}(k,i)}{\gamma (k,i)}}$
	$v_{i}(k)={\frac {\delta _{D}(k,i)}{\xi _{b_{\min }}^{d}(k,i)}}$
	$e(k,i+1)=e(k,i)-v_{i}(k)e_{b}(k,i)\,\!$
	End
	End

सामान्यीकृत जाली पुनरावर्ती न्यूनतम वर्ग फ़िल्टर (एनएलआरएलएस)

एलआरएलएस के सामान्यीकृत रूप में कम पुनरावृत्ति और चर होते हैं। इसकी गणना एल्गोरिदम के आंतरिक चर में सामान्यीकरण लागू करके की जा सकती है जो उनके परिमाण को एक से सीमित रखेगा। इसका उपयोग आम तौर पर वास्तविक समय के अनुप्रयोगों में नहीं किया जाता है क्योंकि विभाजन और वर्ग-रूट संचालन की संख्या उच्च कम्प्यूटेशनल लोड के साथ आती है।

एनएलआरएलएस एल्गोरिदम सारांश

एनएलआरएलएस फ़िल्टर के लिए एल्गोरिदम को संक्षेप में प्रस्तुत किया जा सकता है

Initialization:
	For ${\textstyle i=0,1,\ldots ,N.}$
	${\overline {\delta }}(-1,i)=0\,\!$ (if ${\textstyle x(k)=d(k)=0}$ for ${\textstyle k<0}$)
	${\overline {\delta }}_{D}(-1,i)=0\,\!$
	${\overline {e}}_{b}(-1,i)=0\,\!$
	End
	$\sigma _{x}^{2}(-1)=\lambda \sigma _{d}^{2}(-1)=\varepsilon \,\!$
Computation:
	For ${\textstyle k\geq 0}$
	$\sigma _{x}^{2}(k)=\lambda \sigma _{x}^{2}(k-1)+x^{2}(k)\,\!$ (Input signal energy)
	$\sigma _{d}^{2}(k)=\lambda \sigma _{d}^{2}(k-1)+d^{2}(k)\,\!$ (Reference signal energy)
	${\overline {e}}_{b}(k,0)={\overline {e}}_{f}(k,0)={\frac {x(k)}{\sigma _{x}(k)}}\,\!$
	${\overline {e}}(k,0)={\frac {d(k)}{\sigma _{d}(k)}}\,\!$
	For ${\textstyle i=0,1,\ldots ,N}$
	${\overline {\delta }}(k,i)=\delta (k-1,i){\sqrt {(1-{\overline {e}}_{b}^{2}(k-1,i))(1-{\overline {e}}_{f}^{2}(k,i))}}+{\overline {e}}_{b}(k-1,i){\overline {e}}_{f}(k,i)$
	${\overline {e}}_{b}(k,i+1)={\frac {{\overline {e}}_{b}(k-1,i)-{\overline {\delta }}(k,i){\overline {e}}_{f}(k,i)}{\sqrt {(1-{\overline {\delta }}^{2}(k,i))(1-{\overline {e}}_{f}^{2}(k,i))}}}$
	${\overline {e}}_{f}(k,i+1)={\frac {{\overline {e}}_{f}(k,i)-{\overline {\delta }}(k,i){\overline {e}}_{b}(k-1,i)}{\sqrt {(1-{\overline {\delta }}^{2}(k,i))(1-{\overline {e}}_{b}^{2}(k-1,i))}}}$
	Feedforward filter
	${\overline {\delta }}_{D}(k,i)={\overline {\delta }}_{D}(k-1,i){\sqrt {(1-{\overline {e}}_{b}^{2}(k,i))(1-{\overline {e}}^{2}(k,i))}}+{\overline {e}}(k,i){\overline {e}}_{b}(k,i)$
	${\overline {e}}(k,i+1)={\frac {1}{\sqrt {(1-{\overline {e}}_{b}^{2}(k,i))(1-{\overline {\delta }}_{D}^{2}(k,i))}}}[{\overline {e}}(k,i)-{\overline {\delta }}_{D}(k,i){\overline {e}}_{b}(k,i)]$
	End
	End

यह भी देखें

• अनुकूली फ़िल्टर
• कर्नेल अनुकूली फ़िल्टर
• न्यूनतम माध्य वर्ग फ़िल्टर
• शून्य-बल तुल्यकारक

संदर्भ

• Hayes, Monson H. (1996). "9.4: Recursive Least Squares". Statistical Digital Signal Processing and Modeling. Wiley. p. 541. ISBN 0-471-59431-8.
• Simon Haykin, Adaptive Filter Theory, Prentice Hall, 2002, ISBN 0-13-048434-2
• M.H.A Davis, R.B. Vinter, Stochastic Modelling and Control, Springer, 1985, ISBN 0-412-16200-8
• Weifeng Liu, Jose Principe and Simon Haykin, Kernel Adaptive Filtering: A Comprehensive Introduction, John Wiley, 2010, ISBN 0-470-44753-2
• R.L.Plackett, Some Theorems in Least Squares, Biometrika, 1950, 37, 149–157, ISSN 0006-3444
• C.F.Gauss, Theoria combinationis observationum erroribus minimis obnoxiae, 1821, Werke, 4. Gottinge

टिप्पणियाँ

[Category:Statistical signal processi