مرشح إهليلجي

المرشح الإهليلجي (المعروف أيضًا بمرشح كاوير ، نسبةً إلى ويلهلم كاوير ، أو مرشح زولوتاريف ، نسبةً إلى يغور زولوتاريف ) هو مرشح لمعالجة الإشارات يتميز بتوازن التموج في كلٍ من نطاق التمرير ونطاق الإيقاف . يمكن ضبط مقدار التموج في كل نطاق بشكل مستقل، ولا يوجد مرشح آخر من نفس الرتبة يتميز بانتقال أسرع في الكسب بين نطاق التمرير ونطاق الإيقاف ، وذلك لقيم التموج المحددة (سواءً كان التموج متوازنًا أم لا). بدلاً من ذلك، يمكن التخلي عن إمكانية ضبط تموج نطاق التمرير ونطاق الإيقاف بشكل مستقل، وتصميم مرشح يكون أقل حساسية لتغيرات المكونات.

عندما يقترب التموج في نطاق الإيقاف من الصفر، يصبح المرشح مرشح تشيبيشيف من النوع الأول . وعندما يقترب التموج في نطاق التمرير من الصفر، يصبح المرشح مرشح تشيبيشيف من النوع الثاني ، وأخيرًا، عندما تقترب قيمتا التموج من الصفر، يصبح المرشح مرشح باترورث .

يُعطى كسب مرشح القطع الناقص ذي التمرير المنخفض كدالة للتردد الزاوي ω بالصيغة التالية:

جين(ω)=11+ϵ2Rن2(ξ،ω/ω0){\displaystyle G_{n}(\omega )={1 \over {\sqrt {1+\epsilon ^{2}R_{n}^{2}(\xi ,\omega /\omega _{0})}}}}

حيث R n هي الدالة الكسرية الإهليلجية من الرتبة n (المعروفة أحيانًا باسم دالة تشيبيشيف الكسرية) و

ω0{\displaystyle \omega _{0}} تردد القطع
ϵ{\displaystyle \epsilon }هو عامل التموج
ξ{\displaystyle \xi }هو عامل الانتقائية

تحدد قيمة عامل التموج تموج نطاق التمرير، بينما يحدد مزيج عامل التموج وعامل الانتقائية تموج نطاق الإيقاف.

ملكيات

استجابة التردد لمرشح تمرير منخفض إهليلجي من الرتبة الرابعة بقيمتي ε  =  0.5 و ξ  =  1.05. كما هو موضح أيضًا الحد الأدنى للكسب في نطاق التمرير والحد الأقصى للكسب في نطاق الإيقاف، ومنطقة الانتقال بين التردد المعياري 1 و ξ
صورة مقربة لمنطقة الانتقال في الرسم البياني أعلاه.
  • في نطاق التمرير، تتراوح قيمة الدالة الكسرية الإهليلجية بين الصفر والواحد. وبالتالي، يتراوح كسب نطاق التمرير بين 1 و1/1+ϵ2{\displaystyle 1/{\sqrt {1+\epsilon ^{2}}}}.
  • في نطاق التوقف، تتغير الدالة الكسرية الإهليلجية بين اللانهاية وعامل التمييزلن{\displaystyle L_{n}}والذي يُعرَّف على النحو التالي:
لن=Rن(ξ،ξ){\displaystyle L_{n}=R_{n}(\xi ,\xi )\,}
وبالتالي، سيتراوح كسب نطاق الإيقاف بين 0 و1/1+ϵ2لن2{\displaystyle 1/{\sqrt {1+\epsilon ^{2}L_{n}^{2}}}}.
  • في حدودξ{\displaystyle \xi \rightarrow \infty }تتحول الدالة الكسرية الإهليلجية إلى متعددة حدود تشيبيشيف ، وبالتالي يصبح المرشح مرشح تشيبيشيف من النوع الأول ، بمعامل تموج ε
  • بما أن مرشح باترورث هو شكل محدود لمرشح تشيبيشيف، فإنه يترتب على ذلك أنه في حالة الحدξ{\displaystyle \xi \rightarrow \infty }،ω00{\displaystyle \omega _{0}\rightarrow 0}وϵ0{\displaystyle \epsilon \rightarrow 0}بحيثϵRن(ξ،1/ω0)=1{\displaystyle \epsilon \,R_{n}(\xi ,1/\omega _{0})=1}يتحول المرشح إلى مرشح باترورث
  • في حدودξ{\displaystyle \xi \rightarrow \infty }،ϵ0{\displaystyle \epsilon \rightarrow 0}و ω00{\displaystyle \omega _{0}\rightarrow 0}بحيثξω0=1{\displaystyle \xi \omega _{0}=1}وϵلن=α{\displaystyle \epsilon L_{n}=\alpha }يصبح المرشح مرشحًا من نوع تشيبيشيف الثاني مع كسب
جي(ω)=11+1α2تين2(1/ω){\displaystyle G(\omega )={\frac {1}{\sqrt {1+{\frac {1}{\alpha ^{2}T_{n}^{2}(1/\omega )}}}}}}

الأقطاب والأصفار

لوغاريتم القيمة المطلقة لكسب مرشح إهليلجي من الرتبة الثامنة في فضاء التردد المركب (s = σ + jω) حيث ε = 0.5، ξ = 1.05، وω₀ = 1. تمثل النقاط البيضاء الأقطاب، بينما تمثل النقاط السوداء الأصفار. يوجد إجمالاً 16 قطباً و8 أصفار مزدوجة. ما يبدو كقطب واحد وصفر واحد بالقرب من منطقة الانتقال هو في الواقع أربعة أقطاب وصفران مزدوجان كما هو موضح في العرض المكبر أدناه. في هذه الصورة، يشير اللون الأسود إلى كسب 0.0001 أو أقل، بينما يشير اللون الأبيض إلى كسب 10 أو أكثر.
عرض موسع لمنطقة الانتقال في الصورة أعلاه، يوضح الأقطاب الأربعة والصفران المزدوجين.

ستتطابق أصفار كسب المرشح الإهليلجي مع أقطاب الدالة الكسرية الإهليلجية، والتي تم اشتقاقها في المقالة المتعلقة بالدوال الكسرية الإهليلجية .

يمكن اشتقاق أقطاب كسب مرشح إهليلجي بطريقة مشابهة جدًا لاشتقاق أقطاب كسب مرشح تشيبيشيف من النوع الأول . ولتبسيط الأمر، نفترض أن تردد القطع يساوي واحدًا.(ωصم){\displaystyle (\omega _{pm})}ستكون أصفار مقام كسب المرشح الإهليلجي هي جزء من كسب المرشح. باستخدام التردد المركبs=σ+جω{\displaystyle s=\sigma +j\omega }وهذا يعني أن:

1+ϵ2Rن2(-جs،ξ)=0{\displaystyle 1+\epsilon ^{2}R_{n}^{2}(-js,\xi )=0\,}

تعريف-جs=جد(w،1/ξ){\displaystyle -js=\mathrm {cd} (w,1/\xi )}حيث cd() هي دالة جيب التمام الإهليلجية لجاكوبي، وباستخدام تعريف الدوال الكسرية الإهليلجية ينتج ما يلي:

1+ϵ2جد2(نwكنك،1لن)=0{\displaystyle 1+\epsilon ^{2}\mathrm {cd} ^{2}\left({\frac {nwK_{n}}{K}},{\frac {1}{L_{n}}}\right)=0\,}

أينك=ك(1/ξ){\displaystyle K=K(1/\xi )}وكن=ك(1/لن){\displaystyle K_{n}=K(1/L_{n})}حل لإيجاد قيمة w

w=كنكنجد-1(±جϵ،1لن)+مكن{\displaystyle w={\frac {K}{nK_{n}}}\mathrm {cd} ^{-1}\left({\frac {\pm j}{\epsilon }},{\frac {1}{L_{n}}}\right)+{\frac {mK}{n}}}

حيث يتم توضيح القيم المتعددة للدالة العكسية cd() باستخدام الفهرس الصحيح m .

تكون أقطاب دالة الكسب الإهليلجية كما يلي:

sصم=أناجد(w،1/ξ){\displaystyle s_{pm}=i\,\mathrm {cd} (w,1/\xi )\,}

كما هو الحال بالنسبة لكثيرات حدود تشيبيشيف، يمكن التعبير عن ذلك بشكل معقد صريح ( Lutovac & et al. 2001 ، § 12.8).

sصم=أ+جبج{\displaystyle s_{pm}={\frac {a+jb}{c}}}
أ=-ζن1-ζن21-xم21-xم2/ξ2{\displaystyle a=-\zeta _{n}{\sqrt {1-\zeta _{n}^{2}}}{\sqrt {1-x_{m}^{2}}}{\sqrt {1-x_{m}^{2}/\xi ^{2}}}}
ب=xم1-ζن2(1-1/ξ2){\displaystyle b=x_{m}{\sqrt {1-\zeta _{n}^{2}(1-1/\xi ^{2})}}}
ج=1-ζن2+xأنا2ζن2/ξ2{\displaystyle c=1-\zeta _{n}^{2}+x_{i}^{2}\zeta _{n}^{2}/\xi ^{2}}

أينζن{\displaystyle \zeta _{n}}هي وظيفة منن،ϵ{\displaystyle n,\,\epsilon }وξ{\displaystyle \xi }وxم{\displaystyle x_{m}}هي أصفار الدالة الكسرية الإهليلجية.ζن{\displaystyle \zeta _{n}}يمكن التعبير عنها لجميع قيم n بدلالة دوال جاكوبي الإهليلجية، أو جبريًا لبعض الرتب، وخاصة الرتب 1 و2 و3. بالنسبة للرتبتين 1 و2 لدينا

ζ1=11+ϵ2{\displaystyle \zeta _{1}={\frac {1}{\sqrt {1+\epsilon ^{2}}}}}
ζ2=2(1+ت)1+ϵ2+(1-ت)2+ϵ2(1+ت)2{\displaystyle \zeta _{2}={\frac {2}{(1+t){\sqrt {1+\epsilon ^{2}}}+{\sqrt {(1-t)^{2}+\epsilon ^{2}(1+t)^{2}}}}}}

أين

ت=1-1/ξ2{\displaystyle t={\sqrt {1-1/\xi ^{2}}}}

التعبير الجبري لـζ3{\displaystyle \zeta _{3}}الأمر معقد إلى حد ما (انظر Lutovac & et al. (2001 ، § 12.8.1) ).

يمكن استخدام خاصية التداخل للدوال الكسرية الإهليلجية لبناء تعبيرات من رتبة أعلى لـζن{\displaystyle \zeta _{n}}:

ζمن(ξ،ϵ)=ζم(ξ،1ζن2(لم،ϵ)-1){\displaystyle \zeta _{m\cdot n}(\xi ,\epsilon )=\zeta _{m}\left(\xi ,{\sqrt {{\frac {1}{\zeta _{n}^{2}(L_{m},\epsilon )}}-1}}\right)}

أينلم=Rم(ξ،ξ){\displaystyle L_{m}=R_{m}(\xi ,\xi )}.

اعتبارات التصميم

عوامل الجودة المعيارية لأقطاب مرشح إهليلجي من الرتبة الثامنة بمعامل ξ  =  1.1 كدالة لعامل التموج ε . يمثل كل منحنى أربعة أقطاب، حيث أن أزواج الأقطاب المترافقة المركبة وأزواج الأقطاب الموجبة والسالبة لها نفس عامل الجودة. (ينطبق المنحنيان الأزرق والسماوي تقريبًا). يتم تقليل عامل الجودة لجميع الأقطاب في آن واحد عند ε Qmin  =  1  / L n = 0.02323...   

انظر لوتوفاك وآخرون . (2001 ، § 12.11، 13.14) .

تُحدد المرشحات الإهليلجية عمومًا باشتراط قيمة معينة لتموج نطاق التمرير، αp ، وتموج نطاق الإيقاف، αs ، وحدّة القطع. تُحدد هذه القيم الحد الأدنى للرتبة n للدالة الإهليلجية الكسرية التي تصف سلوك المرشح. على وجه التحديد، لنفترض مرشحًا يمرر الترددات الأقل من ωp ويوقف الترددات الأعلى من ωs ، ونُعرّف الدالة المساعدة .u=121-1-(ωصωs)241+1-(ωصωs)24{\displaystyle u={\frac {1}{2}}\cdot {\frac {1-{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}{1+{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}}}لتقييد كسب نطاق التمرير في [1- G p , 1] وكسب نطاق الإيقاف في [0, 1- G s ] ، يجب أن يحقق الترتيب n [ 1 ]نسجل(16جيsجيص)-سجل(u+2u5+15u9+150u13)أين:جيs=10αs/10-1جيص=10αص/10-1{\displaystyle {\begin{aligned}&n\geq {\frac {\log {\left(16{\frac {G_{s}}{G_{p}}}\right)}}{-\log(u+2u^{5}+15u^{9}+150u^{13})}}\\&{\text{Where:}}\\&Gs=10^{\alpha _{s/10}}-1\\&Gp=10^{\alpha _{p/10}}-1\\\end{aligned}}}تنشأ هذه المعادلة من تقريب لدالة التكامل الإهليلجي الكاملة K؛ الصيغة الدقيقة هي [ 2 ]نك(1-τ12)ك(τ2)ك(τ1)ك(1-τ22){\displaystyle n\geq {\frac {K\left({\sqrt {1-\tau _{1}^{2}}}\right)K(\tau _{2})}{K(\tau _{1})K\left({\sqrt {1-\tau _{2}^{2}}}\right)}}}أينτ1=جيصجيsτ2=ωصωs{\displaystyle {\begin{aligned}\tau _{1}&={\sqrt {\frac {G_{p}}{G_{s}}}}\\\tau _{2}&={\frac {\omega _{p}}{\omega _{s}}}\end{aligned}}}من الاعتبارات التصميمية الأخرى حساسية دالة الكسب لقيم المكونات الإلكترونية المستخدمة في بناء المرشح. تتناسب هذه الحساسية عكسيًا مع عامل الجودة ( Q ) لأقطاب دالة نقل المرشح . يُعرَّف عامل الجودة للقطب كما يلي:

سؤال=-|sصم|2Rهـ(sصم)=-12كوس(arg(sصم)){\displaystyle Q=-{\frac {|s_{pm}|}{2\mathrm {Re} (s_{pm})}}=-{\frac {1}{2\cos(\arg(s_{pm}))}}}

وهو مقياس لتأثير القطب على دالة الكسب. بالنسبة للمرشح الإهليلجي، يحدث أنه بالنسبة لرتبة معينة، توجد علاقة بين عامل التموج وعامل الانتقائية تعمل على تقليل عامل الجودة لجميع الأقطاب في دالة النقل في آن واحد:

ϵسؤالمأنان=1لن(ξ){\displaystyle \epsilon _{Qmin}={\frac {1}{\sqrt {L_{n}(\xi )}}}}

ينتج عن ذلك مرشحٌ ذو حساسيةٍ ضئيلةٍ للغاية لتغيرات المكونات، ولكن ستُفقد القدرة على تحديد تموجات نطاقي التمرير والإيقاف بشكلٍ مستقل. بالنسبة لهذه المرشحات، مع ازدياد الرتبة، يقل التموج في كلا النطاقين ويزداد معدل القطع. إذا قررنا استخدام مرشح إهليلجي ذي عامل جودة منخفض (Q) لتحقيق حد أدنى معين من التموج في نطاقات المرشح مع معدل قطع معين، فإن الرتبة المطلوبة ستكون عمومًا أكبر من الرتبة المطلوبة في حال عدم وجود قيد عامل الجودة المنخفض. ستبدو صورة القيمة المطلقة للكسب مشابهةً جدًا للصورة في القسم السابق، باستثناء أن الأقطاب مُرتبة في دائرة بدلًا من قطع ناقص. لن تكون الأقطاب متباعدةً بانتظام، وستكون هناك أصفار على محور ω، على عكس مرشح باترورث ، الذي تُرتّب أقطابه في دائرة متباعدة بانتظام دون أصفار.

مقارنة مع المرشحات الخطية الأخرى

إليكم صورة توضح المرشح الإهليلجي بجانب أنواع أخرى شائعة من المرشحات التي تم الحصول عليها بنفس عدد المعاملات:

كما هو واضح من الصورة، فإن المرشحات الإهليلجية أكثر حدة من جميع المرشحات الأخرى، لكنها تظهر تموجات على كامل عرض النطاق الترددي.

البناء من أصفار نقل تشيبيشيف

تُعتبر نطاقات التوقف في المرشحات الإهليلجية في جوهرها مرشحات تشيبيشيف ذات أصفار نقل، حيث تُرتّب أصفار النقل بطريقة تُنتج نطاق توقف متساوي التموج. وبناءً على ذلك، يُمكن تحويل معادلة خصائص مرشح تشيبيشيف،ك(s){\displaystyle K(s)}يحتوي على أصفار انعكاس تشيبيشيف في البسط ولا يحتوي على أصفار انتقال في المقام، إلى مرشح إهليلجيك(s){\displaystyle K(s)}يحتوي البسط على أصفار انعكاس القطع الناقص وأصفار انتقال القطع الناقص في المقام، وذلك عن طريق إنشاء أصفار انتقال بشكل متكرر من المعكوس المُقاس لأصفار انعكاس تشيبيشيف، ثم إعادة إنشاء نطاق تمرير تشيبيشيف متساوي التموج من أصفار الانتقال، وتكرار ذلك حتى لا تُحدث التكرارات أي تغييرات أخرى ذات دلالة.ك(s){\displaystyle K(s)}[ 3 ] عامل القياس المستخدم ،Ωج{\displaystyle \Omega _{c}}، هي نسبة تردد القطع بين نطاق الإيقاف ونطاق المرور، وتُعرف أيضًا باسم معكوس "عامل الانتقائية". [ 1 ] وبما أن تصميمات الدوائر الإهليلجية تُحدد عادةً بناءً على متطلبات توهين نطاق الإيقاف،Ωج{\displaystyle \Omega _{c}}، ويمكن اشتقاقها من المعادلات التي تحدد الحد الأدنى للرتبة، n، أعلاه.

الωs/ωص{\displaystyle \omega _{s}/\omega _{p}}نسبة،Ωج{\displaystyle \Omega _{c}}يمكن اشتقاق ذلك من خلال حل مسألة الحد الأدنى للترتيب، n ، المذكورة أعلاه بشكل عكسي من n لإيجادΩج{\displaystyle \Omega _{c}}[ 1 ]

ن= عدد الأقطاب (رتبة المرشح)αص و αs=توهين نطاق المرور وتوهين نطاق الإيقاف، على التواليωص و ωs=ترددات نطاق المرور ونطاق الإيقاف، على التواليد=جيsجيص=10αs/10-110αص/10-1q=(16د)(-1/ن)0=-q+u+2u5+15u9+150u13u= الجذر الحقيقي للمعادلة أعلاهك= عامل الانتقائية =1-(1-2u1+2u)4Ωج=ωsωص=1ك=11-(1-2u1+2u)4{\displaystyle {\begin{aligned}n&={\text{ number of poles (order of the filter)}}\\\alpha _{p}{\text{ and }}\alpha _{s}&={\text{pass band and stop band attenuation, respectively}}\\\omega _{p}{\text{ and }}\omega _{s}&={\text{pass band and stop band frequencies, respectively}}\\D&={\frac {G_{s}}{G_{p}}}={\frac {10^{\alpha _{s}/10}-1}{10^{\alpha _{p}/10}-1}}\\q&=(16D)^{(-1/n)}\\0&=-q+u+2u^{5}+15u^{9}+150u^{13}\\u&={\text{ real root of above equation}}\\k&={\text{ the selectivity factor }}={\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}\\\Omega _{c}&={\frac {\omega _{s}}{\omega _{p}}}={\frac {1}{k}}={\frac {1}{\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}}\\\end{aligned}}}

كثيرات الحدود المميزة،ك(s){\displaystyle K(s)}محسوب منΩج{\displaystyle \Omega _{c}}ويمكن بعد ذلك ترجمة متطلبات التوهين إلى كثيرات حدود دالة النقل،جي(s){\displaystyle G(s)}مع الترجمة الكلاسيكية،جي(s)=1/(1+ε2ك(s)ك(-s))|أعمدة LHP{\displaystyle G(s)={\sqrt {1/(1+\varepsilon ^{2}K(s)K(-s))}}{\bigg |}_{\text{LHP poles}}}أينε2=10أص/10.-1.{\displaystyle \varepsilon ^{2}=10^{Ap/10.}-1.}وأص{\displaystyle A_{p}} هو تموج نطاق التمرير. [ 3 ] [ 4 ]

مثال بسيط

صمم مرشحًا إهليلجيًا بنطاق تمرير ذي تموج قدره 1  ديسيبل من 0 إلى 1 راديان/ثانية، ونطاق إيقاف ذي تموج قدره 40  ديسيبل من 1.25 راديان/ثانية على الأقل إلى{\displaystyle \infty }.

بتطبيق الحسابات المذكورة أعلاه لقيمة n قبل تطبيق دالة ceil() ، نجد أن n تساوي 4.83721900، ويتم تقريبها إلى أقرب عدد صحيح، وهو 5، بتطبيق دالة ceil() ، مما يعني أن مرشحًا إهليلجيًا بخمسة أقطاب مطلوب لتلبية متطلبات التصميم المحددة. بتطبيق الحسابات المذكورة أعلاه لـΩج{\displaystyle \Omega _{c}}يلزم تصميم نطاق إيقاف بمقدار 40  ديسيبل بالضبط من التوهين،Ωج{\displaystyle \Omega _{c}}وقد تبين أن قيمته هي 1.2186824.

يمكن تنفيذ دالة الانعكاس المقياسي متعدد الحدود عن طريق نقل كل جذر، s ، إلىΩج/s{\displaystyle \Omega _{c}/s}ويمكن تحقيق ذلك بسهولة عن طريق عكس متعددة الحدود وقياسها بـΩج{\displaystyle \Omega _{c}}كما هو موضح.

أsن+بsن-1...جs2+دs1+هـs0(هـΩجن)sن+(دΩجن-1)sن-1...(جΩج2)s2+(بΩج1)s1+(أΩج0)s0{\displaystyle {\begin{aligned}&as^{n}+bs^{n-1}\dots cs^{2}+ds^{1}+es^{0}\Longrightarrow ({\frac {e}{\Omega _{c}^{n}}})s^{n}+({\frac {d}{\Omega _{c}^{n-1}}})s^{n-1}\dots ({\frac {c}{\Omega _{c}^{2}}})s^{2}+({\frac {b}{\Omega _{c}^{1}}})s^{1}+({\frac {a}{\Omega _{c}^{0}}})s^{0}\\\end{aligned}}}

تكون خطوات تصميم القطع الناقص كما يلي: [ 3 ]

  1. صمم مرشحات تشيبيشيف ذات  تموج نطاق تمرير 1 ديسيبل.
  2. اعكس جميع الانعكاسات الصفرية حولΩج{\displaystyle \Omega _{c}}لإنشاء أصفار الإرسال
  3. قم بإنشاء نطاق تمرير متساوي التموج من أصفار الإرسال باستخدام العملية الموضحة في أصفار إرسال تشيبيشيف.
  4. كرر الخطوتين ٢ و٣ حتى لا يتغير نطاق التمرير ونطاق الإيقاف بشكل ملحوظ. عادةً، ينتج عن ١٥ إلى ٢٥ تكرارًا اختلافات في المعاملات في حدود ١ × ١٠⁻¹⁵.

لتوضيح الخطوات، تبدأ معادلات K(s) أدناه بمعادلة تشيبيشيف K(s) القياسية، ثم تُكرر العملية. تظهر اختلافات واضحة في التكرارات الثلاثة الأولى. عند الوصول إلى 18 تكرارًا، تصبح الاختلافات في K(s) ضئيلة. يمكن إيقاف التكرارات عندما يصبح التغير في معاملات K(s) صغيرًا بما يكفي لتلبية متطلبات دقة التصميم. تم توحيد جميع تكرارات K(s) أدناه بحيث|ك(ج)|=1{\displaystyle |K(j)|=1}ومع ذلك، يمكن تأجيل هذه الخطوة حتى التكرار الأخير، إذا رغب في ذلك.

التكرار 0: ك(s)=16s5+20s3+5s1التكرار 1: ك(s)=9.2965947s5+12.999133s3+4.0025668s0.14167325s4+0.84164496s2+1التكرار الثاني: ك(s)=8.6496472s5+12.270597s3+3.8746611s0.19518773s4+0.94147634s2+1التكرار 17: ك(s)=8.550086786383502s5+12.157269873073034s3+3.854163602012615s0.2043607336740334s4+0.9573802183509494s2+1التكرار 18: ك(s)=8.550086786383422s5+12.157269873072942s3+3.854163602012599s0.2043607336740334s4+0.9573802183509494s2+1{\displaystyle {\begin{aligned}{\text{iteration 0: }}K(s)&={\frac {16s^{5}+20s^{3}+5s}{1}}\\{\text{iteration 1: }}K(s)&={\frac {9.2965947s^{5}+12.999133s^{3}+4.0025668s}{0.14167325s^{4}+0.84164496s^{2}+1}}\\{\text{iteration 2: }}K(s)&={\frac {8.6496472s^{5}+12.270597s^{3}+3.8746611s}{0.19518773s^{4}+0.94147634s^{2}+1}}\\&\vdots \\{\text{iteration 17: }}K(s)&={\frac {8.550086786383502s^{5}+12.157269873073034s^{3}+3.854163602012615s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\{\text{iteration 18: }}K(s)&={\frac {8.550086786383422s^{5}+12.157269873072942s^{3}+3.854163602012599s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\\end{aligned}}}

للعثور علىجي(s){\displaystyle G(s)}دالة التحويل، قم بما يلي. [ 3 ]

ε2=101دب/10.-1.=0.25892541جي(s)=جي(s)جي(-s)|أعمدة LHP=11+ε2ك(s)ك(-s)|أعمدة LHP=ك(s)دهـنك(-s)دهـنك(s)دهـنك(-s)دهـن+ε2ك(s)نuمك(-s)نuم|أعمدة LHP=0.20436073s4+0.95738022s2+1-18.928479s10-53.78661s8-54.942632s6-22.939175s4-1.931467+1|أعمدة LHP{\displaystyle {\begin{aligned}\varepsilon ^{2}&=10^{1dB/10.}-1.=.25892541\\G(s)&={\sqrt {G(s)G(-s)}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {1}{1+\varepsilon ^{2}K(s)K(-s)}}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {K(s)_{den}K(-s)_{den}}{K(s)_{den}K(-s)_{den}+\varepsilon ^{2}K(s)_{num}K(-s)_{num}}}}{\bigg |}_{\text{LHP poles}}\\&={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{{\sqrt {-18.928479s^{10}-53.78661s^{8}-54.942632s^{6}-22.939175s^{4}-1.931467+1}}{\bigg |}_{\text{LHP poles}}}}\\\end{aligned}}}

للحصول علىجي(s){\displaystyle G(s)}من النصف الأيسر من المستوى المركب، حلل البسط والمقام للحصول على الجذور باستخدام خوارزمية إيجاد الجذور . تخلص من جميع الجذور الموجودة في النصف الأيمن من المستوى المركب للمقام، واقسم الجذور المتكررة في البسط على اثنين، ثم أعد بناء المعادلة.جي(s){\displaystyle G(s)}مع الجذور المتبقية. [ 3 ] [ 4 ] عمومًا، قم بالتطبيع|جي(s)|{\displaystyle |G(s)|}إلى 1 فيs=0{\displaystyle s=0}.

جي(s)=0.20436073s4+0.95738022s2+14.3506872s5+4.0174213s4+8.0362343s3+4.9129149s2+3.4288915s+1{\displaystyle {\begin{aligned}&G(s)={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{4.3506872s^{5}+4.0174213s^{4}+8.0362343s^{3}+4.9129149s^{2}+3.4288915s+1}}\\\end{aligned}}}

للتأكد من أن المثالجي(s){\displaystyle G(s)}صحيح، مخططجي(s){\displaystyle G(s)}على امتدادجω{\displaystyle j\omega }يظهر أدناه تموج نطاق تمرير يبلغ 1  ديسيبل، وتردد قطع يبلغ 1 راديان/ثانية، وتوهين نطاق إيقاف يبلغ 40  ديسيبل يبدأ من 1.21868 راديان/ثانية

محاكاة إهليلجية بخمسة أقطاب

حتى تعديلات الطلبات

لا يمكن تنفيذ مرشحات القطع الناقص ذات الرتبة الزوجية، والمُنفذة باستخدام عناصر سلبية، كالمحاثات والمكثفات وخطوط النقل، ذات نهايات متساوية القيمة على كلا الجانبين، باستخدام دالة النقل الإهليلجية التقليدية دون استخدام ملفات مقترنة، وهو أمر قد لا يكون مرغوبًا فيه أو عمليًا. ويعود ذلك إلى عدم القدرة الفيزيائية على استيعاب أصفار انعكاس تشيبيشيف وأصفار النقل ذات الرتبة الزوجية ، والتي ينتج عنها قيم مصفوفة التشتت S12 التي تتجاوز قيمة S12 عندω=0{\displaystyle \omega =0}، وقيم S12 المحدودة الموجودة عندω={\displaystyle \omega =\infty }إذا لم يكن من الممكن تصميم المرشح بزيادة أو تقليل أحد النهايات لاستيعاب نطاق التمرير S12، فيجب تعديل دالة النقل الإهليلجية بحيث يتم نقل صفر الانعكاس ذي الرتبة الزوجية الأدنى إلىω=0{\displaystyle \omega =0}وأعلى رتبة زوجية لنقل البيانات من الصفر إلىω={\displaystyle \omega =\infty }مع الحفاظ على استجابة التموج المتساوي لنطاق التمرير ونطاق الإيقاف. [ 5 ]

يتضمن التعديل المطلوب تعيين كل قطب وصفر لدالة النقل الإهليلجية بطريقة تُعيّن أدنى تردد انعكاس صفر إلى صفر، وأعلى تردد نقل صفر إلى{\displaystyle \infty }والأقطاب والأصفار المتبقية حسب الحاجة للحفاظ على نطاق تمرير ونطاق إيقاف متساوي التموج. ويمكن إيجاد صفر الانعكاس ذي التردد الأدنى عن طريق تحليلك(s){\displaystyle K(s)}يمكن إيجاد البسط، وأعلى تردد إرسال صفري، عن طريق تحليلك(s){\displaystyle K(s)}المقام.

لتحويل أصفار الانعكاس، يتم تطبيق المعادلة التالية على جميع أقطاب وأصفارك(s){\displaystyle K(s)}[ 5 ] بينما من الناحية النظرية ، يمكن إجراء عمليات الترجمة على أي منك(s){\displaystyle K(s)}أوجي(s){\displaystyle G(s)}يجب استخراج أصفار الانعكاس منك(s){\displaystyle K(s)}لذا، من الأفضل عمومًا إجراء عمليات الترجمة علىك(s){\displaystyle K(s)}.

Rأنا=Rأنا2+ωليا21-ωليا2{\displaystyle R_{i}'={\sqrt {\frac {R_{i}^{2}+\omega _{LO}^{2}}{1-\omega _{LO}^{2}}}}}

أين:

Rأنا{\displaystyle R_{i}}هل الدالة الإهليلجية الأصلية صفر أم قطب؟

Rأنا{\displaystyle R_{i}'}يمثل الصفر أو القطب المُعين لدالة التحويل المعدلة ذات الرتبة الزوجية.

ωليا{\displaystyle \omega _{LO}}هو أدنى تردد انعكاس صفري في نطاق المرور.

علامة المكون التخيلي لـRأنا{\displaystyle R_{i}'}يتم تحديده من خلال إشارة الأصلRأنا{\displaystyle R_{i}}.

لتحويل أصفار الإرسال، يتم تطبيق المعادلة التالية على جميع أقطاب وأصفارك(s){\displaystyle K(s)}[ 5 ] بينما من الناحية النظرية ، يمكن إجراء عمليات الترجمة على أي منك(s){\displaystyle K(s)}أوجي(s){\displaystyle G(s)}إذا كان يجب استخراج أصفار الانعكاس منك(s){\displaystyle K(s)}قد يكون من الأجدى إجراء عمليات الترجمة علىك(s){\displaystyle K(s)}.

Rأنا=(ωحأنا2-1)Rأنا2ωحأنا2+Rأنا2{\displaystyle R_{i}'={\sqrt {\frac {(\omega _{HI}^{2}-1)R_{i}^{2}}{\omega _{HI}^{2}+R_{i}^{2}}}}}

أين:

Rأنا{\displaystyle R_{i}}هل الدالة الإهليلجية الأصلية صفر أم قطب؟

Rأنا{\displaystyle R_{i}'}يمثل الصفر أو القطب المُعين لدالة التحويل المعدلة ذات الرتبة الزوجية.

ωحأنا{\displaystyle \omega _{HI}}هو أعلى تردد إرسال صفري في نطاق المرور.

علامة المكون التخيلي لـRأنا{\displaystyle R_{i}'}يتم تحديده من خلال إشارة الأصلRأنا{\displaystyle R_{i}}إذا كان يتم التشغيل علىجي(s){\displaystyle G(s)}إشارة المكون الحقيقي لـRأنا{\displaystyle R_{i}'}يجب أن تكون القيمة سالبة لتتوافق مع متطلبات النصف الأيسر من المستوى.

من المهم ملاحظة أن جميع التطبيقات تتطلب ترجمة نطاقات المرور والإيقاف. على سبيل المثال، لا تتطلب أجهزة الإرسال والاستقبال الشبكية السلبية سوى ترجمة نطاقات الإيقاف ذات الرتبة الزوجية، وتعمل بكفاءة أكبر مع نطاقات المرور ذات الرتبة الزوجية غير المترجمة. [ 5 ]

متىجي(s){\displaystyle G(s)}عند اكتمال العملية، يتم إنشاء دالة نقل متساوية التموج بقيم مصفوفة تشتت لـ S12 تساوي 1ω=0{\displaystyle \omega =0}و 0 عندω={\displaystyle \omega =\infty }، والتي يمكن تنفيذها باستخدام شبكات سلبية متساوية الإنهاء.

يوضح الرسم التوضيحي أدناه مرشحًا إهليلجيًا من الرتبة الثامنة تم تعديله لدعم الشبكات السلبية متساوية الإنهاء من الرتبة الزوجية عن طريق نقل صفر الانعكاس ذي التردد الأدنى من تردد محدود إلى 0 وصفر الإرسال ذي التردد الأعلى إلى{\displaystyle \infty }مع الحفاظ على استجابة ترددية متساوية التموج لنطاق التمرير ونطاق الإيقاف .

رسم توضيحي إهليلجي معدل بترتيب زوجي
رسم توضيحي إهليلجي معدل بترتيب زوجي

الΩج{\displaystyle \Omega _{c}}وحسابات الرتبة المذكورة في فقرة "إنشاء المرشح الإهليلجي" أعلاه خاصة بالمرشحات الإهليلجية غير المعدلة فقط. مع أن تعديلات الرتبة لا تؤثر على توهين نطاق التمرير أو نطاق الإيقاف، إلا أنه من المتوقع وجود أخطاء طفيفة في الرتبة.Ωج{\displaystyle \Omega _{c}}الحسابات. لذلك، من المهم تطبيق تعديلات الترتيب الزوجي في النهاية.ك(s){\displaystyle K(s)}تكتمل التكرارات إذا رغبنا في الحفاظ على توهين نطاقي المرور والإيقاف. إذا تم إنشاء دالة إهليلجية معدلة من رتبة زوجية منΩج{\displaystyle \Omega _{c}}المتطلب الفعليΩج{\displaystyle \Omega _{c}}سيكون أكبر قليلاً من التصميمΩج{\displaystyle \Omega _{c}}وبالمثل، قد ينتج عن عملية حساب الطلب، ن ، قيمة أصغر من الطلب المطلوب فعليًا.

تطبيق الساعة الرملية

مرشح الساعة الرملية هو حالة خاصة من المرشحات حيث تكون أصفار الانعكاس هي مقلوب أصفار الإرسال  حول تردد قطع مُعَيَّر قدره 3.01 ديسيبل عند تردد قطع 1 راديان/ثانية، مما يجعل جميع أقطاب المرشح تقع على دائرة الوحدة . [ 6 ] يتميز تطبيق مرشح الساعة الرملية الإهليلجي عن مرشح تشيبيشيف العكسي بأن نطاق التمرير أكثر استواءً، كما يتميز عن المرشحات الإهليلجية التقليدية بأن تأخير المجموعة يكون أقل حدة عند تردد القطع.

استجابة التردد المتبادلة S11 و S12 على شكل الساعة الرملية
استجابة ترددية متبادلة S11 و S12 لساعة رملية ذات 7 أقطاب

عملية التخليق

أبسط طريقة لتصميم مرشح الساعة الرملية هي تصميم مرشح إهليلجي بتوهين محدد لنطاق التوقف، As ، وتوهين محسوب لنطاق المرور يفي بمتطلبات الشبكة ثنائية المنافذ عديمة الفقد التي تكون معلمات التشتت|S11|2+|S12|2=1{\displaystyle |S_{11}|^{2}+|S_{12}|^{2}=1}[ 7 ] بالإضافة إلى التحويل الحسابي المعروف لقيمة الديسيبل،(Sأناج)دب=20لoز10(|Sأناج|أرأناتح){\displaystyle (S_{ij})_{dB}=20log_{10}(|S_{ij}|_{arith})}، ينتج عن المعالجة الجبرية متطلبات التوهين المحسوبة لنطاق التمرير التالية.

أص=-10سجل10(1.-10(-أs/10)){\displaystyle A_{p}=-10\log _{10}{(1.-10^{(-A_{s}/10)})}}

سيُنتج التردد A p  ، المُعرّف أعلاه، أصفارًا متبادلة للانعكاس والنفاذ حول تردد قطع غير معروف حتى الآن يبلغ 3.01 ديسيبل. لتصميم مرشح إهليلجي بتردد تمرير يبلغ 1 راديان/ثانية،  يجب تحديد تردد التوهين 3.01 ديسيبل، ويُستخدم هذا التردد لعكس مقياس كثيرات حدود تصميم المرشح الإهليلجي. ستكون النتيجة كثيرات حدود ذات توهين 3.01  ديسيبل عند تردد مُعَيَّر يبلغ 1 راديان/ثانية. يمكن استخدام طريقة نيوتن أو حل المعادلات مباشرةً باستخدام خوارزمية إيجاد الجذور لتحديد  تردد التوهين 3.01 ديسيبل.

قياس التردد باستخدام طريقة نيوتن

لوجي(s){\displaystyle G(s)}تُستخدم دالة نقل الساعة الرملية لإيجاد  التردد 3.01 ديسيبل، وωج{\displaystyle \omega _{c}}لإيجاد التردد المطلوب (3  ديسيبل)، يمكن استخدام الخطوات التالية.ωج{\displaystyle \omega _{c}}

  1. لوجي(s)جي(-s){\displaystyle G(s)G(-s)}غير متوفر بالفعل، اضربجي(s){\displaystyle G(s)}بواسطةجي(-s){\displaystyle G(-s)}للحصول علىجي(s)جي(-s){\displaystyle G(s)G(-s)}.
  2. ينفي جميع بنودsن{\displaystyle s^{n}}متى(ن+2){\displaystyle (n+2)}يقبل القسمة على4{\displaystyle 4}سيكون ذلكs2{\displaystyle s^{2}}،s6{\displaystyle s^{6}}،s10{\displaystyle s^{10}}وهكذا. سيتم استدعاء الدالة المعدلة جي2(s)جي2(-s){\displaystyle G_{2}(s)G_{2}(-s)}وهذا التعديل سيسمح باستخدام الأعداد الحقيقية بدلاً من الأعداد المركبة عند حساب قيمة متعددة الحدود ومشتقتها.ωأ{\displaystyle \omega _{a}}يمكن الآن استخدامها بدلاً من المركبجωأ{\displaystyle j\omega _{a}}
  3. حوّل قيمة التوهين المطلوبة إلى ديسيبل.أدب{\displaystyle A_{dB}}، إلى قيمة ربح حسابي تربيعي،بأرأناتح2{\displaystyle B_{arith}^{2}}باستخدامبأرأناتح2=10أدب/10{\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}على سبيل المثال،  يتم تحويل 3.010 ديسيبل إلى 0.5،  ويتم تحويل 1 ديسيبل إلى 0.79432823 وهكذا.
  4. احسب المعدل|جي2(s)جي2(-s)|{\displaystyle |G_{2}(s)G_{2}(-s)|}في طريقة نيوتن باستخدام القيمة الحقيقية،ωأ{\displaystyle \omega _{a}}. خذ القيمة المطلقة دائمًا.
  5. احسب المشتقة المعدلةجي2(ωأ)جي2(-ωأ){\displaystyle G_{2}(\omega _{a})G_{2}(-\omega _{a})}فيما يتعلق بالقيمة الحقيقية،ωأ{\displaystyle \omega _{a}}لا تأخذ القيمة المطلقة للمشتقة.

عند اكتمال الخطوات من 1 إلى 4، يمكن كتابة التعبير المتعلق بطريقة نيوتن على النحو التالي:

ωأ=ωأ-([جي2(ωأ)جي2(-ωأ)|-ب2)/(د[جي2(ωأ)جي2(-ωأ)]/دωأ){\displaystyle \omega _{a}=\omega _{a}-([G_{2}(\omega _{a})G_{2}(-\omega _{a})|-B^{2})/(d[G_{2}(\omega _{a})G_{2}(-\omega _{a})]/d\omega _{a})}

باستخدام قيمة حقيقية لـωأ{\displaystyle \omega _{a}}دون الحاجة إلى عمليات حسابية معقدة. حركةωأ{\displaystyle \omega _{a}}ينبغي تقييدها لمنعها من أن تصبح سالبة في المراحل المبكرة من التكرارات لزيادة الموثوقية. عند اكتمال التقارب،ωأ{\displaystyle \omega _{a}}يمكن استخدامه لـωج{\displaystyle \omega _{c}}يمكن استخدامها لتغيير حجم الأصلجي(s){\displaystyle G(s)}مقام دالة النقل. توهين المُعدَّلجي(s){\displaystyle G(s)}ستكون القيمة المطلوبة تقريبًا هي القيمة المطلوبة تمامًا عند 1 راديان/ثانية. إذا تم تنفيذ العملية بشكل صحيح، فلن يتطلب الأمر سوى بضع تكرارات لضبط التوهين عبر نطاق واسع من قيم التوهين المطلوبة لكل من المرشحات ذات الرتب الصغيرة والكبيرة جدًا.

قياس التردد مع إيجاد الجذر

منذ|جي(جωأ)|{\displaystyle |G(j\omega _{a})|}بما أن دالة التحويل لا تحتوي على أي معلومات عن الطور، فإن تحليلها المباشر لن ينتج عنه نتائج قابلة للاستخدام. ومع ذلك، يمكن تعديل دالة التحويل بضربها فيجي(-s){\displaystyle G(-s)}للقضاء على جميع القوى الشاذة لـجي(جωأ){\displaystyle G(j\omega a)}وهذا بدوره يجبرجي(جωأ){\displaystyle G(j\omega a)}أن يكون حقيقياً عند جميع الترددات، ثم إيجاد التردد الذي ينتج عنه مربع الانتباه المطلوب.

  1. لوجي(s)جي(-s){\displaystyle G(s)G(-s)}غير متوفر بالفعل، اضربجي(s){\displaystyle G(s)}بواسطةجي(-s){\displaystyle G(-s)}للحصول علىجي(s)جي(-s){\displaystyle G(s)G(-s)}.
  2. حوّل قيمة التوهين المطلوبة إلى ديسيبل.أدب{\displaystyle A_{dB}}، إلى قيمة ربح حسابي تربيعي،بأرأناتح2{\displaystyle B_{arith}^{2}}باستخدامبأرأناتح2=10أدب/10{\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}على سبيل المثال،  يتم تحويل 3.010 ديسيبل إلى 0.5،  ويتم تحويل 1 ديسيبل إلى 0.79432823 وهكذا.
  3. يجدP(S)=جينuم(S)جينuم(-S)-بأرأناتح2جيدهـن(S)جيدهـن(-S){\displaystyle P(S)=G_{num}(S)G_{num}(-S)-B_{arith}^{2}G_{den}(S)G_{den}(-S)}
  4. أوجد جذور P(S) باستخدام خوارزمية إيجاد الجذور.
  5. من بين مجموعة الجذور المذكورة أعلاه، اختر الجذر التخيلي الموجب لجميع مرشحات الرتبة، والجذر الحقيقي الموجب لمرشحات الرتبة الزوجية.ωج{\displaystyle \omega _{c}}.

تغيير مقياس دالة النقل

متىωج{\displaystyle \omega _{c}}بعد تحديد ذلك، يمكن قياس متعددة الحدود لدالة نقل الساعة الرملية على النحو التالي:

جي(s)الأصل=شمالننsنن++شمال2s2+شمال1s+شمال0دننsند++د2s2+د1s+د0 (متعددات حدود دالة النقل الأصلية غير المُقاسة)جي(s) مقياس=شمالنن(ωجs)نن++شمال2(ωجs)2+شمال1ωجs+شمال0دنن(ωجs)ند++د2(ωجs)2+د1ωجs+د0 (3.01 (ديسيبل عند 1 راديان/ثانية، متعددات حدود دالة النقل المُقاسة)ωج=|جي(s)| تردد التوهين 3.01 ديسيبل نن،ند= ترتيب البسط والمقام، على التواليشمال،د= معاملات البسط والمقام، على التوالي{\displaystyle {\begin{aligned}G(s)_{\text{orig}}&={\frac {N_{nn}s^{nn}+\dots +N_{2}s^{2}+N_{1}s+N_{0}}{D_{nn}s^{nd}+\dots +D_{2}s^{2}+D_{1}s+D_{0}}}{\text{ (original unscaled transfer function polynomials)}}\\G(s)_{\text{ scaled}}&={\frac {N_{nn}(\omega _{c}s)^{nn}+\dots +N_{2}(\omega _{c}s)^{2}+N_{1}\omega _{c}s+N_{0}}{D_{nn}(\omega _{c}s)^{nd}+\dots +D_{2}(\omega _{c}s)^{2}+D_{1}\omega _{c}s+D_{0}}}{\text{ (}}3.01{\text{ dB at 1 rad/sec scaled transfer function polynomials)}}\\\omega _{c}&=|G(s)|{\text{ 3.01 dB attenuation frequency }}\\nn,nd&={\text{ order of numerator and denominator, respectively}}\\N,D&={\text{ coefficients of numerator and denominator, respectively}}\\\end{aligned}}}

حتى تعديلات الطلبات

تُعاني مرشحات الساعة الرملية ذات الترتيب الزوجي من نفس القيود المتعلقة بالشبكات السلبية متساوية النهايات التي تُعاني منها مرشحات القطع الناقص الأخرى. وتُعالج التعديلات نفسها ذات الترتيب الزوجي التي تُحل مشكلة مرشحات القطع الناقص المشكلة نفسها التي تُعالجها مرشحات الساعة الرملية.

مراجع

  1. 1 2 3 رورابو، سي. بريتون (1 يناير 1993). دليل مصمم المرشحات الرقمية (  طبعة معاد طباعتها). بلو ريدج ساميت، بنسلفانيا، الولايات المتحدة الأمريكية: تاب بوكس، قسم ماكجرو هيل، الصفحات  93 إلى 95. ISBN 978-0830644315.
  2. بارمان، لاري د. (2001). تصميم وتحليل المرشحات التناظرية، منظور معالجة الإشارات . نورويل، ماساتشوستس، الولايات المتحدة الأمريكية: دار نشر كلوير الأكاديمية. الصفحات من 182 إلى 198. ISBN  0-7923-7373-1.
  3. ١ ٢ ٣ ٤ ٥ ملاحظات محاضرة تصميم المرشحات للدكتور بايرون بينيت، مؤرشفة بتاريخ ٢٣ أبريل ٢٠٢٤ في أرشيف الإنترنت ، ١٩٨٥، جامعة ولاية مونتانا . مؤرشفة بتاريخ ٢٨ مارس ٢٠٢٣ في أرشيف الإنترنت ، قسم الهندسة الكهربائية . مؤرشفة بتاريخ ١٥ أبريل ٢٠٢٤ في أرشيف الإنترنت ، بوزمان ، مونتانا، الولايات المتحدة الأمريكية.
  4. 1 2 سيدرا، عادل س.؛ براكيت، بيتر أو. (1978). نظرية وتصميم المرشحات: النشطة والسلبية . بيفرتون، أويغون، الولايات المتحدة الأمريكية: دار ماتريكس للنشر، ص 45-73 . ISBN  978-0916460143.
  5. 1 2 3 4 سال ، رودولف (يناير 1979). دليل تصميم المرشح (باللغتين الإنجليزية والألمانية) ( الطبعة الأولى). ميونيخ، ألمانيا: Allgemeine Elektricitais-Gesellschaft. الصفحات 25، 26، 56-61 ، 116، 117. ISBN   3-87087-070-2.
  6. بينيت، بايرون ج. (ديسمبر 1988). "تقنية جديدة لتوليف المرشحات - الساعة الرملية" . معاملات IEEE في الدوائر والأنظمة . 35 (12): 1469-1477 . doi : 10.1109/31.9910 عبر IEEE.
  7. ماثاي، جورج ل.؛ يونغ، ليو؛ جونز، إي إم تي (1984). مرشحات الميكروويف، وشبكات مطابقة المعاوقة، وهياكل الاقتران . 610 شارع واشنطن، ديدهام، ماساتشوستس، الولايات المتحدة الأمريكية: دار أرتيك هاوس للنشر (نُشر عام 1985). ص 44. ISBN  0-89006-099-1.{{cite book}}: CS1 maint: location ( link )