מסנן ספרתי

מסנן ספרתי הוא מונח בעיבוד אותות ספרתי, המתאר מערכת המבצעת פעולה ליניארית על דגימות של אות. במוצא המסנן מתקבלות דגימות חדשות כאשר כל דגימה היא סכום משוקלל של הדגימות שלפניה או אחריה. באמצעות המסנן ניתן להדגיש או להפחית תכונות של אות דגום. ישנו קשר הדוק בין המסננים הספרתיים לבין המסננים האנלוגיים, ובתנאים מסוימים ניתן לבצע סינון ספרתי על אות דגום במקום דגימה של אות אנלוגי שעבר סינון אנלוגי.

מימושים של מסנן ספרתי יכולים להיעשות בתוכנה על גבי מעבד כללי או על גבי מעבד אותות ספרתי (DSP), אשר מתוכנן בין השאר על מנת לתמוך בפעולות של סינון באמצעות מקבול של פעולות והוראות ייעודיות לביצוע של מכפלה וסכימה בפקודה אחת. מקובל גם לממש מסננים ספרתיים ב-FPGA או לבנות חומרה ייעודית על מנת לעמוד בדרישות של real time.

החשיבות של מסננים ספרתיים היא בתחומים רבים, למשל תקשורת ספרתית, עיבוד תמונה, עיבוד אותות קול, איכון, שערוך ועוד.

ייצוג

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

$y[n]=\sum _{m=-\infty }^{\infty }{b[m]x[n-m]}-\sum _{m=-\infty ,m\neq 0}^{\infty }{a[m]y[n-m]}$

כאשר $x[n]$ הוא אות הכניסה ו-$y[n]$ הוא אות המוצא, והמקדמים $a[n]$ ו-$b[n]$ הם מקדמי המסנן. במערכת פרקטית הסכימה אינה מתבצעת על גבי אינסוף מקדמים אלא על מספר סופי, והאינדקס הראשון של הסכימה הוא אפס, במקרה זה המסנן נקרא סיבתי. הסכום השני נקרא משוב (Feedback). שימו לב שצורה שקולה לרישום המשוואה היא באמצעות קונבולוציה:

$y[n]*a[n]=x[n]*b[n]$

כאשר $a[0]=1$.

כאשר $a[n]$ שווה לאפס לכל $n$ ו-$b[n]$ הוא בעל מספר סופי של מקדמים המסנן נקרא FIR - Finite Impulse Response כלומר, תגובה סופית להלם (ראו מסנן FIR), אחרת הוא נקרא IIR - Infinite Impulse Response (ראו מסנן IIR), כלומר תגובה אינסופית להלם. הכוונה להלם של קרונקר המוגדר כ-

$\delta [n]={\begin{cases}1,&{\text{if }}n{\text{=0}}\\0,&{\text{otherwise}}\end{cases}}$

נהוג לסמן את התגובה להלם כ-$h[n]$. אם מדובר בתגובה סופית להלם, כלומר $a[n]=0,\forall {n}$ אזי:

$h[n]=b[n],0\leq n\leq N-1$

במקרה ש-$a[n]$ אינו אפס, ניתן להציג את תגובת המערכת להלם כסדרה אינסופית של מקדמים הפועלים על אות הכניסה, ללא משוב.

$h[n]=b[n],n\geq 0$

ראו למטה דוגמה כיצד מוצאים תגובה אינסופית למערכת עם משוב.

את משוואת ההפרשים ניתן להמיר לייצוג באמצעות התמרת Z:

$Y(z)={\frac {B(z)}{A(z)}}X(z)$

כאשר $A(z)$נתון לפי:

$A(z)=\sum _{k=-\infty }^{\infty }{a[k]z^{-k}}$

ן-$B(z),X(z),Y(z)$ מוגדרים באופן דומה.

באופן דומה להתמרת ה-Z, ניתן להתייחס למקדמי המסנן כמקדמים של טור פורייה. בהקשר של מסננים ואותות דגומים הסיגנל בתדר ידוע כ-DTFT - Discrete Time Fourier Transform. למשל, עבור האות $x[n]$ נקבל:

$DTFT(x[n])=\sum _{n=-\infty }^{\infty }{x[n]e^{-j\omega n}}=X(e^{j\omega })$

הסיבה לכך שכותבים $X$ בתלות ב-$e^{j\omega }$ היא כדי להדגיש את מחזוריות ה-$2\pi$ שלו.

דוגמה

נניח כי מקדמי המסנן הם $b[0]=1,b[1]=-0.5,a[1]=-0.25$ אז משוואת ההפרשים היא:

$y[n]=x[n]-0.5x[n-1]+0.25y[n-1]$

אם נסתכל על האיבר הקודם במוצא נקבל:

$y[n-1]=x[n-1]-0.5x[n-2]+0.25y[n-2]$

אם נציב את המשוואה השנייה בראשונה נקבל:

$y[n]=x[n]-0.5x[n-1]+0.25(x[n-1]-0.5x[n-2]+0.25y[n-2])$

$y[n]=x[n]-0.25x[n-1]-0.125x[n-2]+0.0675y[n-2]$

וכך ניתן להמשיך הלאה, כאשר נשים לב כי המקדם של $y[n-k]$ יתקרב לאפס ככל שנמשיך. על מנת לקבל נוסחה סגורה, נסתכל על התמרת Z של המערכת:

$Y(z)={\frac {1-0.5z^{-1}}{1-0.25z^{-1}}}X(z)=X(z)(1-0.5z^{-1})\sum _{k=0}^{\infty }0.25^{k}z^{-k}$

$Y(z)=X(z)(\sum _{k=0}^{\infty }0.25^{k}z^{-k}-2\sum _{k=1}^{\infty }0.25^{k}z^{-k})=X(z)(1-\sum _{k=1}^{\infty }0.25^{k}z^{-k})$

ואם נחזור לתחום הזמן נקבל:

$y[n]=x[n]-\sum _{k=1}^{\infty }x[n-k]0.25^{k}$

מימוש

עבור מסנני FIR מדובר בסכימה פשוטה למדי, למשל על מימוש של המסנן $h[n]={\frac {1}{3}}\delta [n]+{\frac {1}{3}}\delta [n-1]+{\frac {1}{3}}\delta [n-2]$.

ישנם שני סוגים נפוצים למימוש של מסנן IIR:

Direct Form 1

מימוש DIRECT FORM 1 של מסנן IIR מסדר 2

זהו מימוש המבוסס על חישוב מפורש של משוואת ההפרשים של המסנן. עבור סכימה זו נדרשים 2N השהיות ו-2N סכימות, כאשר N הוא סדר הפולינום של התמרת ה-Z של המערכת (בהנחה שסדר המונה וסדר המכנה זהה).

Direct Form 2

מימוש DIRECT FORM 2 של מסנן IIR מסדר 2

מימוש זה מבוסס על הפעלת המכנה של התמרת ה-Z של המסנן ולאחר מכן הפעלת המונה על התוצאה. עבור סכימה זו נדרשים N השהיות ו-2N סכימות.

הקשר לסינון אנלוגי

נניח שאות רציף עובר מסנן אנלוגי ולאחר מכן דגימה:

$y(t)=\textstyle \int \limits _{-\infty }^{\infty }\displaystyle h(\tau )x(t-\tau )d\tau$

$y_{s}[n]=y(nT)$

לפי משפט הדגימה מתקיים הקשר הבא:

$Y_{s}(e^{j\omega })={\frac {1}{T}}\sum _{k=-\infty }^{\infty }Y({\frac {\omega -2\pi k}{T}})={\frac {1}{T}}\sum _{k=-\infty }^{\infty }X({\frac {\omega -2\pi k}{T}})H({\frac {\omega -2\pi k}{T}})$

אם נניח ש-$H$ ו-$X$ חסומי סרט ובנוסף הדגימה היא בקצב נייקוויסט אזי אין חפיפות בתדר (Aliasing) בין האיברים שבסכימה. במקרה כזה אם נגדיר $x_{s}[n]=x(nT),h_{s}[n]=h(nT)$,

ֿאזי:

$Y_{s}(e^{j\omega })={\frac {1}{T}}\sum _{k=-\infty }^{\infty }X({\frac {\omega -2\pi k}{T}})H({\frac {\omega -2\pi k}{T}})={\frac {1}{T}}\sum _{k=-\infty }^{\infty }X({\frac {\omega -2\pi k}{T}})\sum _{k=-\infty }^{\infty }H({\frac {\omega -2\pi k}{T}})=TH_{s}(e^{j\omega })X_{s}(e^{j\omega })$

כאשר השוויון השני נכון מכיוון שתנאי נייקוויסט מתקיים, כלומר:

$X({\frac {\omega -2\pi k}{T}})H({\frac {\omega -2\pi m}{T}})=0,\forall \omega ,k\neq m$

מכיוון שמכפלה בתדר שקולה לקונבולוציה בזמן אז בתחום הזמן מקבלים:

$y_{s}[n]=T\sum _{k=-\infty }^{\infty }x_{s}[n-k]h_{s}[k]$

המשמעות היא שניתן להחליף סינון אנלוגי ולאחריו דגימה של אות המוצא בדגימה ולאחריה סינון ספרתי כאשר המסנן הספרתי הוא דגימות של המסנן האנלוגי.

מסננים נפוצים

מסנן מעביר נמוכים (Lowpass Filter)

הכוונה במסנן מעביר נמוכים היא למסנן המעביר תדרים נמוכים מ-$\omega _{c}$ ומאפס את כל התדרים שמעליו, כאשר $\omega _{c}$ ידוע כתדר הקיטעון. מסנן מעביר נמוכים אידיאלי נתון לפי:

$\mathrm {Lowpass} _{ideal}={\begin{cases}1,&{\text{if }}{|\omega |<\omega _{c}}\\0,&{\text{otherwise}}\end{cases}}$

התגובה בזמן של המסנן האידיאלי היא אינסופית, ולכן במימוש מעשי מגדירים 'קטע מעבר' (transition band) בו המסנן יורד מ-1 ל-0 בצורה מתונה יותר.

מסנן מעביר גבוהים (Highpass Filter) ומסנן מעביר פס (Bandpass Filter)

באופן דומה למסנן מעביר הנמוכים מסנן המעביר גבוהים מעביר רק תדרים גבוהים מ-$\omega _{c}$, בעוד מסנן מעביר פס מעביר תדרים בין $\omega _{1}$ ל-$\omega _{2}$.

גרסאות ספרתיות למסנן אנלוגיים

כפי שהראינו קודם ניתן להשתמש בדגימות של מסנן אנלוגי לתכנון מסנן ספרתי, ולכן למסננים אנלוגיים רבים יש גרסה ספרתית. למשל מסנן צ'בישב או מסנן באטרווארת'.

מסננים מסתגלים

מסנן מסתגל הוא מסנן אשר מקדמיו משתנים בזמן בהתאם לקריטריון אופטימיזציה כלשהו. לדוגמה, נניח שאות ידוע, $x(t)$, עובר סינון אנלוגי ואנו מעוניינים לשערך את המסנן באופן ספרתי מתוך דגימות של אות המוצא $y(nT)$, למשל בשיטת Least Mean Squares.

ניתן להתחיל מניחוש התחלתי $h_{estimated}^{(0)}[n]=\delta [n]$, ולחשב את השגיאה הרגעית:

$e^{(0)}=y(nT)-h_{estimated}^{(0)}[n]*x(nT)$

כאשר $*$ מסמן קונבולציה בזמן בדיד.

כעת ניתן לעדכן את מקדמי המסנן, למשל באופן הבא:

$h_{estimated}^{(1)}[n]=h_{estimated}^{(0)}[n]+\mu e^{(0)}x(nT)$

כך ניתן להמשיך ולעדכן את השערוך באופן איטרטיבי.

ראו גם

• מסנן מתואם
• מסנן קלמן
• אקולייזר
• תורת הדגימה
• מסנן Root-raised-cosine
• ADC
• DAC

קריאה נוספת

• Oppenheim, Alan V., and Ronald W. Schafer. Digital Signal Processing. Prentice Hall, 1975. ISBN: 9780132146357.
• Simon Haykin: Adaptive Filter Theory, Prentice Hall, 2002

הערך באדיבות ויקיפדיה העברית, קרדיט,
רשימת התורמים
רישיון cc-by-sa 3.0

29894111מסנן ספרתי