תגובת הלם סופית

בעיבוד אותות, מסנן תגובת הלם סופית (FIR) הוא מסנן שתגובת ההלם שלו (או התגובה לכל קלט באורך סופי) היא בעלת אורך סופי, מכיוון שהיא מתייצבת לאפס בזמן סופי. זאת בניגוד למסנני תגובת הלם אינסופית (IIR), אשר עשויים להיות בעלי משוב פנימי ועשויים להמשיך להגיב ללא הגבלה (בדרך כלל מתכלה).

תגובת ההלם (כלומר, הפלט בתגובה לכניסת דלתא של קרונקר) של מסנן FIR בזמן בדיד מסדר N נמשכת בדיוק $N+1$ דגימות (מהאלמנט הראשון שאיננו אפס ועד האלמנט האחרון שאינו אפס) לפני שהוא מתייצב לאפס.

מסנני FIR יכולים להיות בזמן בדיד או בזמן רציף, דיגיטליים או אנלוגיים.

הגדרה

מסנן FIR בצורה ישירה בזמן דיסקרטי בסדר N . החלק העליון הוא קו השהייה של N עם N + 1 יציאות. כל השהיה של יחידה הוא אופרטור z ^-1 בסימון Z-transform .

מסנן FIR בזמן דיסקרטי בצורת סריג בסדר N. כל השהיה של יחידה הוא אופרטור z ⁻¹ בסימון Z-transform.

עבור מסנן FIR בזמן דיסקרטי סיבתי בסדר N, כל ערך בפלט הוא סכום משוקלל של ערכי הקלט העדכניים ביותר :

${\begin{aligned}y[n]&=b_{0}x[n]+b_{1}x[n-1]+\cdots +b_{N}x[n-N]\\&=\sum _{i=0}^{N}b_{i}\cdot x[n-i],\end{aligned}}$

כאשר:

• ${\textstyle x[n]}$ הוא אות הכניסה,
• ${\textstyle y[n]}$ הוא אות הפלט,
• ${\textstyle N}$ הוא סדר הסינון; במסנן מסדר ${\textstyle N}$ יש ${\textstyle N+1}$ ביטויים באגף ימין
• ${\textstyle b_{i}}$ הוא הערך של תגובת ההלם ברגע ה-i עבור ${\textstyle 0\leq i\leq N}$ של מסנן FIR מסדר N. אם המסנן הוא מסנן FIR אז ${\textstyle b_{i}}$ הוא גם מקדם של המסנן.

חישוב זה ידוע גם בתור קונבולוציה בדידה.

האיברים ${\textstyle x[n-i]}$ מכונים בדרך כלל s, המבוסס על מבנה של קו השהייה מוצמד שבמימושים רבים או בדיאגרמות בלוקים מספק את הכניסות המושהות לפעולות הכפל.

תגובת ההלם של המסנן כפי שהוגדר אינה אפס לאורך זמן סופי. כולל אפסים, תגובת ההלם היא הרצף האינסופי :

$h[n]=\sum _{i=0}^{N}b_{i}\cdot \delta [n-i]={\begin{cases}b_{n}&0\leq n\leq N\\0&{\text{otherwise}}.\end{cases}}$

אם מסנן FIR אינו סיבתי, טווח הערכים שאינם אפס בתגובת ההלם שלו יכול להתחיל לפני $n=0$, עם הנוסחה מוגדרת ומוכללת.

תכונות

למסנן FIR יש מספר מאפיינים שימושיים שלעיתים הופכים אותו לעדיף על מסנן של תגובה אינסופית (IIR). מסנני FIR:

• אינו דורש משוב. משמעות הדבר היא ששגיאות עיגול אינן מורכבות על ידי איטרציות מסוכמות. אותה שגיאה יחסית מתרחשת בכל חישוב. זה גם הופך את היישום לפשוט יותר.
• הם יציבים, מכיוון שהפלט הוא סכום של מספר סופי של כפולות סופיות של ערכי הקלט, ולכן אינו יכול להיות גדול מ ${\textstyle \sum |b_{i}|}$ כפול הערך במקסימלי בקלט.
• ניתן לעצב בקלות לשלב ליניארי על ידי הפיכת רצף המקדם לסימטרי. מאפיין זה רצוי לפעמים עבור יישומים Phase-sensitive, למשל תקשורת נתונים, סייסמולוגיה, מסנני מוצלבים ומאסטרינג.

החיסרון העיקרי של מסנני FIR הוא שנדרש הרבה יותר כוח חישוב במעבד בהשוואה למסנן IIR בעל חדות או סלקטיביות דומים, במיוחד כאשר יש צורך בחיתוכים בתדר נמוך (ביחס לקצב הדגימה). עם זאת, מעבדי אותות דיגיטליים רבים מספקים תכונות חומרה מיוחדות כדי להפוך מסנני FIR יעילים בערך כמו IIR עבור יישומים רבים.

תגובת תדר

השפעת המסנן על הרצף $x[n]$ מתואר בתחום התדר על ידי משפט הקונבולציה :

$\underbrace {{\mathcal {F}}\{x*h\}} _{Y(\omega )}=\underbrace {{\mathcal {F}}\{x\}} _{X(\omega )}\cdot \underbrace {{\mathcal {F}}\{h\}} _{H(\omega )}$  ו  $y[n]=x[n]*h[n]={\mathcal {F}}^{-1}{\big \{}X(\omega )\cdot H(\omega ){\big \}},$

כאשר האופרטורים ${\mathcal {F}}$ ו ${\mathcal {F}}^{-1}$ בהתאמה מייצגים התמרת פורייה בזמן הבדיד (DTFT) ואת ההופכי שלו. לכן, הפונקציה המורכבת, המכפלת $H(\omega )$ היא תגובת התדר של המסנן. זה מוגדר על ידי סדרת פורייה :

$H_{2\pi }(\omega )\ \triangleq \sum _{n=-\infty }^{\infty }h[n]\cdot \left({e^{i\omega }}\right)^{-n}=\sum _{n=0}^{N}b_{n}\cdot \left({e^{i\omega }}\right)^{-n},$

כאשר כתב המשנה שנוסף מציין מחזוריות של 2π. $\omega$ מייצג תדר ביחידות מנורמלות (רדיאנים/דגימה). ההחלפה $\omega =2\pi f,$, משנה את יחידות התדר $(f)$ למחזורים/דגימה והמחזוריות ל-1. כאשר לרצף x[n] יש קצב דגימה ידוע, $f_{s}$ דגימות/שנייה, ההחלפה $\omega =2\pi f/f_{s}$ משנה את יחידות התדירות $(f)$ למחזורים/שנייה (הרץ) והמחזוריות ל $f_{s}.$ הערך $\omega =\pi$ מתאים לתדירות של $f={\tfrac {f_{s}}{2}}$ הרץ $={\tfrac {1}{2}}$ מחזורים/מדגם, שהוא תדר Nyquist.

$H_{2\pi }(\omega )$ ניתן לבטא גם במונחים של טרנספורמציה Z של תגובת דחף המסנן:

${\widehat {H}}(z)\ \triangleq \sum _{n=-\infty }^{\infty }h[n]\cdot z^{-n}.$

$H_{2\pi }(\omega )=\left.{\widehat {H}}(z)\,\right|_{z=e^{j\omega }}={\widehat {H}}(e^{j\omega }).$

תכנון מסנן

מסנן FIR מתוכנן על ידי מציאת המקדמים וסדר הסינון העומדים במפרטים מסוימים, שיכולים להיות בתחום הזמן (לדוגמה מסנן מתואם) ו/או בתחום התדר (הנפוץ ביותר). מסננים מתואמים מבצעים קורלציה צולבת בין אות הכניסה לצורת פולס ידועה. קונבולוציית ה-FIR היא מתאם צולב בין אות הכניסה לבין עותק הפוך בזמן של תגובת ההלם. לכן, תגובת ההלם של המסנן המתואם "מעוצבת" על ידי דגימת צורת ההלם הידועה ושימוש בדגימות הללו בסדר הפוך כמקדמי המסנן.

כאשר רצויה תגובת תדר מסוימת, נפוצות מספר שיטות עיצוב שונות:

1. שיטת תכנון חלונות
2. שיטת דגימת תדר
3. שיטת מזעור MSE (ממוצע שגיאה בריבוע).
4. שיטת Parks–McClellan (מוכרת גם כשיטת equiripple, אופטימלית או מינימקס). אלגוריתם החלפת רמז משמש בדרך כלל למציאת קבוצת מקדמים אופטימלית equiripple. כאן המשתמש מציין תגובת תדר רצויה, פונקציית שקלול לשגיאות מתגובה זו וסדר סינון N. לאחר מכן האלגוריתם מוצא את סט ${\textstyle N+1}$ מקדמים שממזערים את הסטייה המקסימלית מהאידיאל. באופן אינטואיטיבי, זה מוצא את המסנן הקרוב ככל האפשר לתגובה הרצויה בהינתן זאת בלבד ${\textstyle N+1}$ ניתן להשתמש במקדמים. שיטה זו קלה במיוחד בפועל שכן לפחות טקסט אחד כולל תוכנית שלוקחת את המסנן הרצוי ואת N, ומחזירה את המקדמים האופטימליים.
5. ניתן לעצב מסנני Equiripple FIR באמצעות אלגוריתמי DFT גם כן. האלגוריתם הוא איטרטיבי באופיו. ה-DFT של עיצוב מסנן ראשוני מחושב באמצעות אלגוריתם FFT (אם אומדן ראשוני אינו זמין, ניתן להשתמש ב-h[n]=delta[n]). בתחום הפורייה, או תחום ה-DFT, תגובת התדר מתוקנת בהתאם למפרט הרצוי, ולאחר מכן מחשבים את ה-DFT ההפוך. בתחום הזמן, רק ה-N המקדמים הראשונים נשמרים (המקדמים האחרים מוגדרים לאפס). לאחר מכן, התהליך חוזר על עצמו באופן איטרטיבי: ה-DFT מחושב פעם נוספת, תיקון מוחל בתחום התדר וכן הלאה.

חבילות תוכנה כגון MATLAB, GNU Octave, Scilab ו־SciPy מספקות דרכים נוחות ליישם את השיטות השונות הללו.

שיטת עיצוב חלונות

בשיטת עיצוב החלונות, מתכננים תחילה מסנן IIR אידיאלי ולאחר מכן חותכים את תגובת ההלם האינסופי על ידי הכפלתה עם פונקציית חלון באורך סופי. התוצאה היא מסנן תגובת הלם סופית שתגובת התדר שלו שונה מזו של מסנן IIR. הכפלת ההלם האינסופי בפונקציית החלון בתחום הזמן גורמת לכך שתגובת התדר של ה-IIR מתפתלת עם טרנספורמציה פורייה (או DTFT) של פונקציית החלון. אם האונה הראשית של החלון צרה, תגובת התדר המרוכבת נשארת קרובה לזו של מסנן IIR האידיאלי.

התגובה האידיאלית היא לרוב מלבנית, וה-IIR המקביל הוא פונקציית sinc. התוצאה של קונבולוציית תחום התדר היא שהקצוות של המלבן מתחדדים, ואדוות מופיעות בפס המעבר ובפס העצירה. בעבודה לאחור, אפשר לציין את השיפוע (או הרוחב) של האזור המחודד (רצועת מעבר) ואת גובה האדוות, ובכך לגזור את הפרמטרים של תחום התדר של פונקציית חלון מתאימה. המשך אחורה לתגובת הלם יכול להיעשות על ידי איטרציה של תוכנית עיצוב מסנן כדי למצוא את סדר הסינון המינימלי. שיטה נוספת היא להגביל את ערכת הפתרון למשפחה הפרמטרית של חלון קייזר, המספקת יחסים בצורה סגורה בין הפרמטרים של תחום הזמן והתחום התדר. באופן כללי, שיטה זו לא תשיג את סדר הסינון המינימלי האפשרי, אך היא נוחה במיוחד עבור יישומים אוטומטיים הדורשים עיצוב סינון דינמי, תוך כדי תנועה.

שיטת עיצוב החלונות מועילה גם ליצירת מסנני חצי פס יעילים, מכיוון שפונקציית ה-sinc המתאימה היא אפס בכל נקודת דגימה אחרת (למעט המרכזית). המוצר עם פונקציית החלון אינו משנה את האפסים, כך שכמעט מחצית מהמקדמים של תגובת הדחף הסופית הם אפס. יישום מתאים של חישובי FIR יכול לנצל את המאפיין הזה כדי להכפיל את יעילות המסנן.

ראו גם

• מסנן (אלקטרוניקה)
• התמרת Z
• תומך (מתמטיקה)
• תגובת הלם אינסופית

קישורים חיצוניים

מדיה וקבצים בנושא תגובת הלם סופית בוויקישיתוף