خط تأخير رقمي

تمثيل مخطط الكتلة القياسي لخط التأخير ذي العدد الصحيح M. [ 1 ]

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

تُستخدم خطوط التأخير الرقمية على نطاق واسع في طرق محاكاة الصوتيات في الغرف ، والآلات الموسيقية ، ووحدات المؤثرات الصوتية . ويُظهر توليف الموجات الرقمية كيف يمكن استخدام خطوط التأخير الرقمية كطرق لتوليف الصوت لمختلف الآلات الموسيقية، مثل الآلات الوترية وآلات النفخ .

إذا احتوى خط التأخير على قيمة غير صحيحة أصغر من واحد، فإنه ينتج عنه خط تأخير كسري (يُسمى أيضًا خط تأخير مُستكمل أو مرشح تأخير كسري). تُستخدم سلسلة من خط تأخير صحيح ومرشح تأخير كسري بشكل شائع لنمذجة مرشحات التأخير العشوائية في معالجة الإشارات الرقمية . [ 2 ] يُعد مخطط داتورو تطبيقًا قياسيًا في الصناعة للمرشحات الرقمية باستخدام خطوط التأخير الكسرية. [ 3 ]

نظرية

يُشتق خط التأخير القياسي ذو التأخير الصحيح من تحويل Z لإشارة زمنية منفصلةx{\displaystyle x}تأخر بسببم{\displaystyle M}العينات [ 4 ] :

y[ن]=x[ن-م]{\displaystyle y[n]=x[nM]}Z{\displaystyle {\xrightarrow[{}]{\mathcal {Z}}}}Y(z)=z-محم(z)X(z).{\displaystyle Y(z)=\overbrace {z^{-M}} ^{H_{M}(z)}X(z).}

في هذه الحالة،z-م=حم(z){\displaystyle z^{-M}=H_{M}(z)}مرشح تأخير الأعداد الصحيحة مع:

{||=1=0دب،كسب صفر ديسيبل=-ωم،طور خطي مع ω=2πوتيs أين تيs فترة أخذ العينات بالثواني [s].{\displaystyle {\begin{cases}|\centerdot |=1=0dB,&{\text{كسب صفر ديسيبل}}\\\measuredangle =-\omega M,&{\text{طور خطي مع }}\omega =2\pi fT_{s}{\text{ حيث }}T_{s}{\text{ هي فترة أخذ العينات بالثواني }}[s].\end{cases}}}

مرشح المجال الزمني المنفصل لتأخير الأعداد الصحيحةم{\displaystyle M}باعتباره التحويل العكسي لزيتا لـحم(z){\displaystyle H_{M}(z)}وهو أمر تافه، لأنه عبارة عن دافع تم إزاحته بواسطةم{\displaystyle M}[ 5 ] :

حم[ن]={1،ل ن=م0،ل نم.{\displaystyle h_{m}[n]={\begin{cases}{\text{1}},&{\text{for }}n=M\\0,&{\text{for }}n\neq M.\end{cases}}}

يُعدّ العمل في مجال الزمن المتقطع مع التأخيرات الكسرية أقل بساطة. في أبسط صوره النظرية، يُعرَّف خط التأخير ذو التأخير الكسري العشوائي بأنه خط تأخير قياسي ذو تأخيردR{\displaystyle D\in \mathbb {R} }والتي يمكن نمذجتها على أنها مجموع مكون صحيحمZ{\displaystyle M\in \mathbb {Z} }ومكون كسريدR{\displaystyle d\in \mathbb {R} }وهو أصغر من عينة واحدة:

خط التأخير (الجزئي) -Z{\displaystyle {\mathcal {Z}}}اِختِصاص

هذا هوZ{\displaystyle {\mathcal {Z}}}تمثيل المجال لمشكلة تصميم مرشح رقمي غير تافهة : الحل هو أي مرشح في المجال الزمني يمثل أو يقارب التحويل العكسي Z لـحد(z){\displaystyle H_{D}(z)}[ 2 ]

حلول تصميم المرشحات

حل ساذج

يُمكن الحصول على الحل الأسهل من الناحية المفاهيمية عن طريق أخذ عينات من حل مجال الزمن المستمر، وهو أمر بديهي لأي قيمة تأخير. بالنظر إلى إشارة زمنية مستمرةx{\displaystyle x}تأخر بسببدR{\displaystyle D\in \mathbb {R} }عينات، أو τ=دتيs{\displaystyle \tau =DT_{s}}ثوانٍ [ 6 ] :

y(ت)=x(ت-د){\displaystyle y(t)=x(tD)}F{\displaystyle {\xrightarrow[{}]{\mathcal {F}}}}Y(ω)=هـ-جωدحأنادهـأل(ω)X(ω).{\displaystyle Y(\omega )=\overbrace {e^{-j\omega D}} ^{H_{ideal}(\omega )}X(\omega ).}

في هذه الحالة،هـ-جωد=حأنادهـأل(ω){\displaystyle e^{-j\omega D}=H_{مثالي}(\أوميغا )}هو مرشح التأخير الجزئي في المجال الزمني المستمر مع:

{||=1=0دب،كسب صفر ديسيبل=-ωد،الطور الخطيτزر=-ددω=د،تأخير المجموعة المستمرτصح=-ω=-د،تأخير طور ثابت.{\displaystyle {\begin{cases}|\centerdot |=1=0dB,&{\text{كسب صفر ديسيبل}}\\\measuredangle =-\omega D,&{\text{طور خطي}}\\\tau _{gr}=-{d\measuredangle \over {d\omega }}=D,&{\text{تأخير مجموعة ثابت}}\\\tau _{ph}=-{\measuredangle \over {\omega }}=-D,&{\text{تأخير طور ثابت.}}\end{cases}}}

الحل البسيط للمرشح المأخوذ عينات منهحأنادهـأل[ن]{\displaystyle h_{ideal}[n]}هو التحويل العكسي لفورييه المأخوذ عينات منهحأنادهـأل(ω){\displaystyle H_{ideal}(\omega )}، مما ينتج عنه مرشح IIR غير سببي على شكل دالة جيبية أساسيةsأنانج(){\displaystyle sinc()}تم تغييرها بواسطةد{\displaystyle D}[ 6 ] :

حأنادهـأل[ن]=F-1[حأنادهـأل(ω)]=12π-π+πهـجωدهـجωندω=sأنانج(ن-د)=sأنان(π(ن-د))π(ن-د){\displaystyle h_{ideal}[n]={\mathcal {F}}^{-1}[H_{ideal}(\omega )]={1 \over {2\pi }}\int \limits _{-\pi }^{+\pi }e^{j\omega D}e^{j\omega n}d\omega =sinc(n-D)={sin(\pi (n-D)) \over {\pi (n-D)}}}

مجال الزمن المستمرsأنانج{\displaystyle sinc}يتم إزاحتها بمقدار التأخير الجزئي بينما تتم محاذاة أخذ العينات دائمًا مع المستوى الديكارتي، لذلك:

  • عندما يكون التأخير عددًا صحيحًا من العيناتدشمال{\displaystyle D\in \mathbb {N} }، العينة التي تم تحويلهاsأنانج{\displaystyle sinc}يتحول إلى نبضة منزاحة تمامًا كما هو الحال في الحل النظري.
  • عندما يكون التأخير عددًا كسريًا من العيناتدR{\displaystyle D\in \mathbb {R} }، العينة التي تم تحويلهاsأنانج{\displaystyle sinc}ينتج مرشح IIR غير سببي، وهو أمر غير قابل للتنفيذ عمليًا.
رسوم متحركة لحركة التزامن
يُحصل على خط التأخير الكسري المثالي عن طريق أخذ عينات من تحويل فورييه العكسي لمرشح التأخير الكسري في مجال الزمن المستمر. لاحظ كيف يتحول هذا الأمر، عند قيمة تأخير صحيحة، إلى نبضات مُزاحة بسيطة. إن تأخير إشارة مُعينة باستخدام هذا المرشح يُعادل من الناحية النظرية إعادة أخذ عينات من مصدرها التناظري بنفس فترة أخذ العينات ولكن مع إزاحة محاذاة العينات بمقداردR{\displaystyle D\in \mathbb {R} }لاحظ أيضًا أن الصورة لا تُظهر سوى عدد قليل من العينات حول الصفر، ولكن يتم تعريف IIR غير السببي لعدد لا نهائي من العينات في كلا اتجاهي المحور السيني.

حل FIR السببي المقتطع

إن الحل الأسهل من الناحية النظرية والقابل للتطبيق هو الحذف السببي للحل البسيط المذكور أعلاه. [ 7 ]

حτ[ن]={sأنانج(ن-د)ل 0نشمال0خلاف ذلكأينشمال-12<د<شمال+12وشمالهذا هو ترتيب المرشح.{\displaystyle h_{\tau }[n]={\begin{cases}sinc(n-D)&{\text{for }}0\leq n\leq N\\0&{\text{otherwise}}\end{cases}}\;\;\;\;\;{\text{where}}\;\;\;\;\;{N-1 \over {2}}<D<{N+1 \over {2}}\;\;\;\;\;{\text{and}}\;\;\;\;\;N\;{\text{is the order of the filter.}}}

لكن تقليص استجابة النبضة قد يتسبب في عدم الاستقرار، والذي يمكن التخفيف منه بعدة طرق:

  • يتم تطبيق نافذة على استجابة النبضة المقتطعة، وبالتالي تنعيمها. لاحظ أنه في هذه الحالة، يجب إضافة إزاحة أخرى.ل{\displaystyle L}من أجل محاذاة النافذة وsأنانج(){\displaystyle sinc()}وتوفير ترشيح متناظر [ 7 ] [ 8 ] .

    حτ[ن]={w(ن-د)sأنانج(ن-د)ل لنل+شمال0خلاف ذلكأينل={رouند(د)-شمال2حتى شمالد-شمال-12للفردي شمال{\displaystyle h_{\tau }[n]={\begin{cases}w(n-D)sinc(n-D)&{\text{for }}L\leq n\leq L+N\\0&{\text{otherwise}}\end{cases}}\;\;\;\;\;{\text{where}}\;\;\;\;\;L={\begin{cases}round(D)-{N \over {2}}&{\text{for even }}N\\\lfloor D\rfloor -{N-1 \over {2}}&{\text{for odd }}N\end{cases}}}

  • طريقة المربعات الصغرى العامة (GLS): [ 2 ] تقوم بتعديل استجابة التردد بشكل متكرر عن طريق تحديد نافذة لتصميم خطأ التكامل للمربعات الصغرى، مما يقلل من خطأ التكامل التربيعي بين استجابات التردد المثالية والمقتطعة للمرشح، والمحددة على النحو التالي:

هـلS=12π-απαπw(ω)|حدترuنجأتهـد(هـجω)-حدأناد(هـجω)|2دωأين 0<α1 هو معلمة عرض نطاق التمرير{\displaystyle E_{LS}={1 \over {2\pi }}\int \limits _{-\alpha \pi }^{\alpha \pi }w(\omega )|H_{D}^{truncated}(e^{j\omega })-H_{D}^{id}(e^{j\omega })|^{2}d\omega \;\;\;\;\;{\text{where }}0<\alpha \leq 1{\text{ is the passband width parameter}}}

  • مُستكمِل لاغرانج (مُرشِّح التأخير الكسري ذو التسطيح الأقصى): [ 9 ] يُضيف قيود "التسطيح" إلى المشتقات N الأولى لخطأ التكامل التربيعي الأدنى. وتكتسب هذه الطريقة أهمية خاصة لأنها تمتلك حلاً مغلقاً.
تمثيل تخطيطي لصيغة لاغرانج للاستيفاء. [ 10 ]

حد[ن]=ك=0،كنشمالد-كن-كأين0نشمال{\displaystyle h_{D}[n]=\prod _{k=0,\;k\neq n}^{N}{D-k \over {n-k}}\;\;\;\;\;{\text{where}}\;\;\;\;\;0\leq n\leq N}

فيما يلي شرح موسع للصيغة المذكورة أعلاه، يعرض نتائج التصفية حسب الترتيب حتىشمال=3{\displaystyle N=3}:

توسيع صيغة لاغرانج Interpolator [ 7 ]
حτ[0]{\displaystyle h_{\tau }[0]}حτ[1]{\displaystyle h_{\tau }[1]}حτ[2]{\displaystyle h_{\tau }[2]}حτ[3]{\displaystyle h_{\tau }[3]}
ن = 11-د{\displaystyle 1-D}د{\displaystyle D}--
ن = 2(د-1)(د-2)2{\displaystyle {(D-1)(D-2) \over {2}}}-د(د-2){\displaystyle -D(D-2)}د(د-1)2{\displaystyle {D(D-1) \over {2}}}-
ن = 3-(د-1)(د-2)(د-3)6{\displaystyle -{(D-1)(D-2)(D-3) \over {6}}}د(د-2)(د-3)2{\displaystyle {D(D-2)(D-3) \over {2}}}-د(د-1)(د-3)2{\displaystyle -{D(D-1)(D-3) \over {2}}}د(د-1)(د-2)6{\displaystyle {D(D-1)(D-2) \over {6}}}

حل تقريبي للطور باستخدام خوارزمية IIR ذات المرور الكامل

يتمثل نهج آخر في تصميم مرشح IIR من الرتبةشمال{\displaystyle N}مع بنية تحويل Z التي تجبرها على أن تكون تمريرة كاملة مع الاستمرار في تقريبهاد{\displaystyle D}التأخير [ 7 ] :

حد(z)=z-شمالأ(z)أ(z-1)=أشمال+أشمال-1z-1+...+أ1z-(شمال-1)+z-شمال1+أ1z-1+...+أشمال-1z-(شمال-1)+أشمالz-شمالوالذي لديه{||=1=0دب0دب يكسبحد(z)=-شمالω+2أ(z)=-دωالقيمة المطلوبة للتأخير د{\displaystyle H_{D}(z)={z^{-N}A(z) \over {A(z^{-1})}}={a_{N}+a_{N-1}z^{-1}+...+a_{1}z^{-(N-1)}+z^{-N} \over {1+a_{1}z^{-1}+...+a_{N-1}z^{-(N-1)}+a_{N}z^{-N}}}\;\;\;\;\;{\text{which has}}\;\;\;\;\;{\begin{cases}|\centerdot |=1=0dB&0dB{\text{ gain}}\\\measuredangle _{H_{D}(z)}=-N\omega +2\measuredangle _{A(z)}=-D\omega &{\text{desired value for delay }}D\end{cases}}}

الأصفار والأقطاب الموضوعة بشكل متبادل لـأ(z) و أ(z-1){\displaystyle A(z){\text{ and }}A(z^{-1})}على التوالي، تسطيح التردد||{\displaystyle |\centerdot |}الاستجابة ، بينما الطور هو دالة لطورأ(z){\displaystyle A(z)}وبالتالي، تصبح المشكلة هي تصميم مرشح FIRأ(z){\displaystyle A(z)}أي إيجاد معاملاتهأك{\displaystyle a_{k}}كدالة لـ D (لاحظ أنأ0=1{\displaystyle a_{0}=1}دائماً)، بحيث تقترب المرحلة من القيمة المطلوبة على أفضل وجهحد(z)=-دω{\displaystyle \measuredangle _{H_{D}(z)}=-D\omega }[ 7 ]

الحلول الرئيسية هي:

  • التقليل التكراري لخطأ الطور المربع الأدنى، [ 2 ] والذي يُعرَّف على النحو التالي:

هـلS=12π-ππw(ω)|-دωأناد-(-شمالω+2أ(z))حΔحد|2دω{\displaystyle E_{LS}={1 \over {2\pi }}\int \limits _{-\pi }^{\pi }w(\omega )|\underbrace {\underbrace {-D\omega } _{\measuredangle _{ID}}-\underbrace {(-N\omega +2\measuredangle _{A(z)})} _{\measuredangle _{H}}} _{\Delta \measuredangle _{H_{D}}}|^{2}d\omega }

  • التقليل التكراري لخطأ تأخير الطور في المربعات الصغرى ، [ 2 ] والذي يُعرَّف على النحو التالي:

هـلS=12π-ππw(ω)|Δحدω|2{\displaystyle E_{LS}={1 \over {2\pi }}\int \limits _{-\pi }^{\pi }w(\omega )|{{\Delta \measuredangle _{H_{D}}} \over {\omega }}|^{2}}

  • مرشح تمرير منخفض ذو أقطاب كاملة من نوع ثيران مع تأخير مجموعة مسطح إلى أقصى حد . [ 11 ] ينتج عن هذا حل مغلق لإيجاد المعاملات.أك{\displaystyle a_{k}}للتأخير الإيجابيد>0{\displaystyle D>0}:

أك=(-1)ك(شمالك)ل=0شمالد+لد+ك+لأين(نك)=شمال!ك!(شمال-ك)!{\displaystyle a_{k}=(-1)^{k}{\binom {N}{k}}\prod _{l=0}^{N}{D+l \over {D+k+l}}\;\;\;\;\;{\text{where}}\;\;\;\;\;{\binom {n}{k}}={N! \over {k!(N-k)!}}}

فيما يلي شرح موسع للصيغة أعلاه يعرض المعاملات الناتجة من الرتبة حتىشمال=3{\displaystyle N=3}:

توسيع صيغة معاملات مرشح تمرير الترددات المنخفضة ذي جميع الأقطاب من نوع ثيران [ 7 ]
أ0{\displaystyle a_{0}}أ1{\displaystyle a_{1}}أ2{\displaystyle a_{2}}أ3{\displaystyle a_{3}}
ن = 11-د-1د+1{\displaystyle -{D-1 \over {D+1}}}--
ن = 21-2د-2د+1{\displaystyle -2{D-2 \over {D+1}}}(د-1)(د-2)(د+1)(د+2){\displaystyle {(D-1)(D-2) \over {(D+1)(D+2)}}}-
ن = 31-3د-3د+1{\displaystyle -3{D-3 \over {D+1}}}3(د-2)(د-3)(د+1)(د+2){\displaystyle 3{(D-2)(D-3) \over {(D+1)(D+2)}}}-(د-1)(د-2)(د-3)(د+1)(د+2)(د+3){\displaystyle -{(D-1)(D-2)(D-3) \over {(D+1)(D+2)(D+3)}}}

التاريخ التجاري

خط التأخير الرقمي Eventide DDL 1745

استُخدمت خطوط التأخير الرقمية لأول مرة عام 1973 لتعويض سرعة الصوت في الهواء، وذلك لتوفير أوقات تأخير مناسبة لأبراج مكبرات الصوت البعيدة في مهرجان "سمر جام" الموسيقي في واتكينز غلين بنيويورك، والذي حضره 600 ألف شخص. وقد وفرت شركة "إيفنتايد كلوك ووركس" التي تتخذ من مدينة نيويورك مقرًا لها ، أجهزة تأخير رقمية، كل منها قادر على تأخير الصوت لمدة 200 مللي ثانية. وُضعت أربعة أبراج مكبرات صوت على بُعد 60 مترًا من المسرح، مع تأخير إشارتها بمقدار 175 مللي ثانية لتعويض سرعة الصوت بين مكبرات الصوت الرئيسية وأبراج التأخير. كما وُضعت ستة أبراج أخرى على بُعد 120 مترًا من المسرح، مما استلزم تأخيرًا قدره 350 مللي ثانية، بينما وُضعت ستة أبراج أخرى على بُعد 180 مترًا من المسرح، مع تأخير قدره 525 مللي ثانية. احتوت كل وحدة من وحدات Eventide DDL 1745 على مئة شريحة مسجل إزاحة سعة 1000 بت ومحول رقمي تناظري مصمم خصيصًا ، وبلغت تكلفتها 3800 دولار ( ما يعادل 29248 دولارًا في عام 2025 ). [ 12 ] [ 13 ] 

انظر أيضاً

مراجع

  1. "خط تأخير العينة M" . ccrma.stanford.edu . تم الاطلاع عليه بتاريخ 2023-07-06 .
  2. 1 2 3 4 5 لاكسو، تيمو آي؛ فاليماكي، فيسا؛ كارجالاينن، ماتي أ؛ Laine, Unto K. (يناير 1996), "تقسيم تأخير الوحدة [تصميم FIR/جميع مرشحات المرور]"، مجلة IEEE لمعالجة الإشارات ، المجلد. 13، لا. 1، الصفحات من 30 إلى 60، بيب كود : 1996ISPM...13...30L ، دوى : 10.1109/79.482137   
  3. سميث، جوليوس أو.؛ ​​لي، نيلسون (5 يونيو 2008)، "النمذجة الصوتية الحاسوبية مع التأخير الرقمي" ، مركز أبحاث الحاسوب في الموسيقى والصوتيات ، تم الاطلاع عليه بتاريخ 21 أغسطس 2007
  4. "خطوط التأخير" . ccrma.stanford.edu . تم الاطلاع عليه بتاريخ 2023-07-06 .
  5. "مقدمة في المرشحات الرقمية مع تطبيقات صوتية" . ccrma.stanford.edu . تم الاطلاع عليه بتاريخ 2023-07-06 .
  6. 1 2 "الاستيفاء المثالي المحدود النطاق (Sinc)" . ccrma.stanford.edu . تم الاطلاع عليه بتاريخ 2023-07-06 .
  7. 1 2 3 4 5 6 فاليماكي، فيسا (1998). "نمذجة الوقت المنفصل للأنابيب الصوتية باستخدام مرشحات التأخير الجزئي" .
  8. Harris, F.J. (1978). "On the use of windows for harmonic analysis with the discrete Fourier transform". Proceedings of the IEEE. 66 (1): 51–83. doi:10.1109/proc.1978.10837. ISSN 0018-9219. S2CID 426548.
  9. Hermanowicz, E. (1992). "Explicity [sic] formulas for weighting coefficients of maximally flat tunable FIR delays". Electronics Letters. 28 (20): 1936. doi:10.1049/el:19921239.
  10. Smith, Julius (5 September 2022). "Explicit Formula for Lagrange Interpolation Coefficients". ccrma.
  11. Thiran, J.-P. (1971). "Recursive digital filters with maximally flat group delay". IEEE Transactions on Circuit Theory. 18 (6): 659–664. doi:10.1109/TCT.1971.1083363. ISSN 0018-9324.
  12. Nalia Sanchez (July 29, 2016), "Remembering the Watkins Glen Festival", Eventide Audio, retrieved February 20, 2020
  13. "DDL 1745 Digital Delay". Eventide Audio. Retrieved 2023-07-22.

Further reading