עיבוד אותות דיגיטלי

זהו הקורס השנים-עשר בהנדסת חשמל באוניברסיטה: עיבוד אותות דיגיטלי — האחרון במסלול החשמל. עיבוד אותות** — הסיום של מסלול ההנדסת חשמל. EE06 נבנה המרתים החיות במחשב, המסננים שתוכלו לשלוח למעשה, ו המנגנון הרב-קצבי שמעביר בין קצבים. המנגנון רב-קצבי שמעביר בין קצבים.

• יחידות 12.1–12.3 מעצבות מסננים: FIR באמצעות חלונות, IIR באמצעות טרנספורמציה בינארית
• יחידות 12.4–12.5 מעצבות פילטרים: FIR באמצעות חלונות, IIR באמצעות טרנספורמציה דו-ליניארית transform.
• Units 12.6–12.7 change rates and build a small verified pipeline.

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

ספקטרום, כל שיפור מהירות – בסופו של דבר הוא תוכנית הפועלת על רשימת מספרים. כל התנהגות התוכנית זו נקבעת על ידי רשימה אחת: תגובת הדחיקה. מספרים. כל התנהגותו של התוכנית נקבעת על ידי רשימה אחת: ה- תגובת הדחיפה.

הרעיון, עם תמונה. מערכת דיסקרטית היא LTI כאשר היא מלאה התכנסות $h[n]$ — החלק את הקדמים התמזגות $y[n] = \sum_{k} h[k]\cdot x[n-k]$ — הזז את הנקודות על פני הקלט וסכום אלכסון אחד בכל פעם.

סיבתיות פירושה $h[n] = 0$ עבור $n < 0$ (אין תגובה לפני כניסה קלט); יציבות BIBO פירושה כי $\sum |h[n]|$ מתכנס. כל הטענה בקורס הזה מצטמצמת לאחת משתי המשפטים הללו.

$x[n] = e^{j\omega n}$ — המעריך המרוכב הוא דרכו דרך סכום ההתמזגות: כל איבר לוכד את הגורם $h[k]e^{j\omega(n-k)} = e^{j\omega n}h[k]e^{-j\omega k}$, ולכן התוצאה היא $e^{j\omega n}\sum_{k}h[k]e^{-j\omega k}$ — האקספוננציאל היא פונקציית מובן, וערך המובן $H(e^{j\omega}) = \sum h[n]e^{-j\omega n}$ הוא תגובת התדר. תגובת התדר היא פשוט DTFT של הקוצבים: יחיד 12.2 נשען בדיוק על כך.

דוגמה מפורשת 1 — שולחן הזז. $h = [1, 2, 1]$ (מסנן חלק עם 3 קוצבים על צעד $x = [1, 1, 1, 1, \ldots]$**: הזז והוסף — $y = 1, 3, 4, 4, \ldots$ הערך היציב הוא $\sum h = 4$ ✓ (צעד דרכו בכל FIR מתייצב בסכום הקוצבים — זהו הרווח DC $H(1)$). פולס סופי $[1,1,1]$ נותן $1, 3, 4, 3, 1$, המדמה את הנקודות: הפלט קורא את הנקודות קדימה ואז אחורה. יציאה קוראת את הקוצבים קדימה ואחורה.

דוגמה מפורטת 2 — האח הגנרטיבי. הממוצע הנע $y[n] = (y[n-1] + x[n])/2$ הוא בעל תגובת פולס אינסופית $h[n] = 0.5^{n+1}$: תגובת השלב שלו מחושבת ל-$0.5, 0.75, 0.875, 0.9375$ — מחצית המרחק ל-1 בכל צעד, האקספוננציאל $1 - 2^{-(n+1)}$ ✓, כי הקוטב היחיד שלו יושב ב-$z = 0.5$. השווו את שני מסנני החלקה: הזיכרון של ה-FIR הוא בדיוק 3 דגימות; זה הרקורסיבי מעולם לא שוכח לחלוטין. $|pole| < 1$ היא יציבות — יחידה 12.5 עושה קריירה על כך.

טעויות נפוצות. הפוך את האות הלא נכון בטבלה (ההתכתבות היא קומוטטיבית, אך טבלת יד חייבת להישאר עקבית אחרת אלכסונות ישקר). ציטוט $\sum h$ כ"הגברת" מבלי לומר שהוא הגברה בתדר DC רק. וציפייה שמסנן רקורסיבי יתייצב — הוא מתקרב, לעולם מגיע.

סיכום. LTI פירונו קונבולוציה עם $h[n]$; זמניות היא $h[n]=0$ עבור סיכום. LTI פירונו התכתבות עם $n<0$; סיבה היא $\sum|h| < \infty$ עבור אקספוננטים הם פונקציות מובן עם ערך מובן $H(e^{j\omega})$.

אקספוננציאלים הם פונקציות אופייניות עם ערך אופייני $N$. הספקטרום חייב להפוך לרשימה סופית של מספרים. איך רשימה זו מופיעה לך חייב להיות רשימה סופית של מספרים. איך רשימה זו מייצגת אותך

הרעיון, עם תמונה. ה-DTFT $X(e^{j\omega}) = \sum_{n} x[n]e^{-j\omega n}$ הוא מחזורי ב-$\omega$ עם מחזור $2\pi$. דגימתו ב-$N$ נקודות מספקת את ה- DFT — והדגימה בתדר משמעות that the time signal מתועדת כמחזורית: עובדה זו בלבד מסבירה את הדליפה.

רזולוציה היא החלטה בזמן: התאים הם $\Delta f = f_{s}/N$ רוחבים, לכן רזולוציה של 1 הרץ ב-$f_{s} = 1$ קילוהרץ משמעות $N = 1000$ — שנייה של נתונים, ללא יוצאי מן הכלל.

הגזירה. דגימת ה-DTFT ב-$\omega_{k} = 2\pi k/N$ היא כפילה שלה במסרק תדרים; לפי הדואליות האות בזמן הופך למחזורי עם מחזור $N$. אם הגוון שלך יושב בדיוק בתא, העתקים המחזוריים מתחברים בצורה חלקה — שני שיאים נקיים. אם הוא יושב בין תאים, העתקים אינם מסכימים והחוסר הסכמה מפזר אנרגיה לכל תא. עיטור הנתונים עם $w[n]$ מבצע קונבולוציה של הספקטרום האמיתי עם $W(e^{j\omega})$: צלעות הלוח הן רק $-13$ dB למטה, ולכן הריכוך גלוי; חלון המהמינג מדכא את צלעות הלוח ל $-43$ דציבל, ומשיג ספקטרום נקי על מחיר לבה ראשית רחבה יותר.

דוגמה פתורה 1 — הפיזור, בכמות. DFT בן 8 נקודות של קוסין בדיוק בתא 1 ($f = f_{s}/8$): $|X| = 4$ בתאים 1 ו-7, אפס בשאר המקומות ✓ (כל 8 הדגימות מתארות מחזור שלם. הזז את הגוון חצי תא ($f = 1.5\cdot f_{s}/8$): הגדלים מחושבים כ-$2.40, 2.80, 1.19, 1.0, \ldots$ בכל תא ✓ — גוון בין תאים נראה לעין בכל מקום, והאות הקטן הסמוך לאות גדול נעלם בתוך ה- החלקים הצדדיים שלו. זהו העלות האמיתית של דליפה.

דוגמה מופתרת 2 — תקציב הרזולוציה. ב-$f_{s} = 1$ kHz: $N = 100$ נותן תאים ב-10 הרץ — שני גוונים במרחק 6 הרץ הם תלישה אחת, אולם ככל שהגרף נקי יותר. $N = 1000$ נותן תאים ב-1 הרץ לשנייה אחת של נתונים. אפס-פיזור ל-4000 נקודות מצייר את אותה תלישה עם 40× דיו more: האפס-פיזור ממערבל את אותו ספקטרום, אך אינו פותר כלום חדש. הרזולוציה נקבעה כשבחרת כמה זמן למדוד.

טעויות נפוצות. קריאת מספרי תאים כהרץ — המפה היא $f = k\cdot f_{s}/N$. האמינה שאפס-ריפוי רוכש רזולוציה. וכן שכחת שה-DFT מניחה שה-$N$ דגימות שלך חוזרות לנצח: השבר הגבול בין העתקה לעתקה הוא בעצמו אות רחב-פס.

סיכום. $X[k]$ = DTFT מדוגם ב-$2\pi k/N$ ⇔ תקופתי בזמן עם תקופה $N$; גוונים מחוץ לתאים דולפים דרך החלקים הצדדיים של החלון (מלבן $-13$ dB, המינג $-43$ dB); $\Delta f = f_{s}/N$ רוכש בעלות זמן מדידה, לעולם לא בעלות ריפוי.

למה זה חשוב. DFT עולה ב-$N^{2}$ פעולות מורכבות — המרת נקודה 1024 היא מיליון כפלות. ה-FFT הוא האלגוריתם שהפך את עיבוד אות דיגיטלי לנושא קיים.

הרעיון, עם תמונה. פצו $N$ נקודות DFT לשתי DFTs בעלות $N/2$ נקודות ועם $N/2$ כפילות שילוב — והמרוכזות. ה פרפר היא יחידת השילוב: $\log_{2}N$ שלבי DFTs בעלות $N/2$ נקודות, כל אחת מכפילה-וסכמה.

המסקע. חלקו $x[n]$ לדגימות זוגיות ואי-זוגיות: $X[k] = \sum_{even} x[n]W_{N}^{kn} + W_{N}^{k}\sum_{odd}x[n]W_{N}^{kn}$ עם גורם הפנייה $W_{N} = e^{-j2\pi/N}$. כל סכום הוא DFT בעלת $N/2$ נקודות (Danielson–Lanczos); הזרוע האי-זוגית זקוקה רק ל גורם $W_{N}^{k}$ לשילוב מחדש. המרוכזות בעומק $\log_{2}N$. ה פירוד זוגי/אי-זוגי מתבצע בחלוקה הפוכה דרך ביטי הכתובת — ולכן סדר כניסה בהפכת ביטים עבור הגרסת החלוקה בזמן.

דוגמה מפורטת 1 — הספירה. ספירת פעולות כפילה-וסכמה מרוכבת פעולות: ישירות $N^{2} = 1024^{2} = 1\,048\,576$; FFT $N\log_{2}N = 1024\times 10 = 10\,240$ — 102× פחות ✓. ב- $N = 4096$: $4096^{2}$ מול $4096\times 12$ = 341× ✓. המרווח מתרחב עם $N$: ה-FIT הוא $O(N\log N)$, והכלל כולו של החישוב משתנה.

דוגמה מפורטת 2 — המנתח הספקטרלי. ⍴ bins בתדר 1 Hz $f_{s} = 1$ ק"הצ זקוק ל-$N \ge 1000$: בחר $N = 1024$ (חזקה של שתיים), ומספק $\Delta f = 1000/1024 = 0.977$ הרץ ✓ בדיוק מספיק. בעשר ננושניות לכל פעולה, המרתק נגמר בתוך כ-$100\ \mu s$ — בזמן אמת מול רקורד 1.024. ערוץ 100 ב-DFT הוא תדר $100\times 0.977 = 97.66$ הרץ ✓ — ההתאמה לתא-הרץ היא $k\cdot f_{s}/N$, ושכחון ממנה הוא השגיאה הקלאסית במעבדה.

טעויות נפוצות. הנחה ש-$N$ חייבת להיות חזקה של שתיים (מערכי שורש מעורבים ועיבוד Bluestein מטפלים בכל $N$, מקרה החזקה של שתיים הוא פשוט ה- של שתיים הוא פשוט ה fastest). ציפייה שה-FFT ישנה את התשובה — היא מחשבת את אותו DFT, רק מהר יותר. וקריאת תאגרים כתדרים, שוב.

סיכום. חלוקה זוגי/אי-זוגי + גורמי עיוות, $\log_{2}N$ שלבים, $N\log_{2}N$ פעולות: מהירות פי 102 ב-$N=1024$ הגדלה עם $N$; המתמטיקה היא DFT, והמהנדסה היא הפרפר.

למה זה חשוב. המסנן FIR הוא הכלי העיקרי שניתן להוכיח עליו דברים: יציבות מוחלטת, ובנוסף — עם טפס — בדיוק פאז ליניארי. תדרים בפס עבר שטוח שומרים על הפאזים היחסיים שלהם; סינון משרעת עדיין יכול לשנות גל.

כלל עיוות לגל. הרעיון, עם איור. התחל מתגובה ההילוך תגובה (סינק), חתוך אותה ל-$N$ טפרות עם חלון $w[n]$, וה- צורת הספקטרלית של החלון היא גורל הפילטר: רוחב האנבה הראשי הלוב האופייני שלו הופך לחגורת המעבר שלך, והלובות הצדדיות הופכות

הגדרות עם מספרים. צלעות לצידי מקסימליות של החלונים עצמם: רבועית $-13$ dB, Hann $-31$ ד"ב, Hamming $-43$ ד"ב, Blackman $-58$ ד"ב. רוחב מעבר $\approx 3.3\cdot f_{s}/N$ עבור חלון המינג. טפסים סימטריים אמיתיים $h[n] = h[N-1-n]$ נותנים פאז ליניארי עם עיכוב קבוצתי $(N-1)/2$ דגימות בכל מקום שבו התגובה אינה אפס. הפאז עשוי לקפוץ ב- $\pi$ כאשר גורם המשרעת הממשי משנה סימן; פאז ועיכוב קבוצתי אינם מוגדרים באפסי התגובה.

הגזירה. חיתוך $h[n]$ הוא כפל ב$w[n]$; ב תדר התגובה האמיתית מתחלקת (convolves) עם $W(e^{j\omega})$. לכן התגובה המופעלת היא הקיר הבלוק האידיאלי שמורח על ידי החלון: גלימת מעבר, מעבר בעובי של לוב main אחד, וסרן עכבה שהעיכוב שלהם תלוי בקונבולוציה, בתדר הגזירה ובמספר הטפסים. ערך צלע צדדי בחלון אינו מבטיח פס חסום למסנן הנובע ממנו. יש להשתמש בטבלת החלונים וההערכה של רוחב המעבר לבחור עיצוב התחלתי, ואז לחשב את תגובתו. (משפט האלטרנטציה שבבסיס Parks–McClellan הופך את העיצוב באקו-גלימה לאופטימלי כאשר חלונים הם פשוט נוחים.)

דוגמה מפורשת 1 — הגדרת גודל הפילטר. FIR Hamming בתדר $f_{s} = 1$ kHz חייב לעבור עד 200 Hz ולהתעכב החל מ-250 Hz. לבחור 225 Hz כהגזרה אידיאלית בתוך רוחב המעבר 50 Hz: $N = 3.3\cdot 1000/50 = 66$ ✓ — לקחת 67, אי-זוגי, עבור Type-I סימטרי עם עיכוב $(67-1)/2 = 33$ דגימות (33 ms). בתדר $N = 31$ אותו חלון נותן $3.3\cdot 1000/31 = 106$ Hz מעבר: שני העיצובים סוגרים את המפרט; ההשהייה היא עלות להקלטה, ולא הערה שולית.

דוגמה מפורטת 2 — לוודא את העיצוב. יצירת 67 טפסים מתוך $h[n] = 0.45\operatorname{sinc}(0.45(n-33))w[n]$, שבה $\operatorname{sinc}(x)=\sin(\pi x)/(\pi x)$ וערךו באפס הוא 1. נרמל את סכום הקצוות ל-1. הערך $H(f)=\sum_{n=0}^{66}h[n]e^{-j2\pi fn/1000}$. על רשת 0.025 הרץ, פס העברת החישוב (0–200 הרץ) נע בין $-0.034$ ל-$+0.0045$ דציבלים, והרמה המקומית הגדולה ביותר בפס חוסם (250–500 הרץ) היא כ-$-52.3$ דציבלים. בדיקות הדגימה הללו עומדות ביעדים של גליגול של ±0.05 דציבלים ורמת פס חוסם של עד $-43$ דציבל; רשת דקה יותר או בדיקות חיפוש קיצוניים בין הדגימות. אם עיצוב אחר מניב $-37$ דציבלים, הוא לא עומד ביעד זה: בדוק את הקצוות, תדר החיתוך, החלון, מספר הקצוות ואת פס המדידה במקום לנחש סיבה יחידה.

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

סיכום. שכן את ה-sinc האידיאלי: טבלת האזניות המדריכה לבחירת החלון, סרטון החוסם, $3.3f_{s}/N$ קובע את המעבר, סימטריה קובעת פאזה ליניארית בהשהייה $(N-1)/2$; ואז מדוד את התגובה מול כל המספר המובטח.

למה זה חשוב. אותו סף 100 הרץ שמחייב FIR ב-67 טפרות נקודות הקצפה, עולה לביוקוואד IIR 5 מקדם — כי הוא יורש ארבעים שנות ידע מעצבי אנלוגיים דרך החלפה אחת.

הרעיון, עם תמונה. המרת בייליאר $s = \frac{2}{T}\cdot\frac{1-z^{-1}}{1+z^{-1}}$ ממיפה את כל האנלוגי חצי המישור השמאלי לתוך עיגול היחידה: כל מסנן אנלוגי יציב ופך למסנן דיגיטלי יציב, אוטומטית.

המחיר הוא עיוות תדר עיוות: ציר האנלוג $-\infty \ldots \infty$ מתעגל פעם אחת סביב עיגול היחידה, בצורה דחוסה, כך שכל תדר קריטי חייב להיות מותאם מראש $\omega_{a} = \frac{2}{T}\tan(\omega_{d}T/2)$ לפני השימוש.

הנגזרת. הצבת $s = j\omega_{a}$ בממיר: הפתרון הוא $z = \frac{1+sT/2}{1-sT/2}$, שהמגנווד שלו הוא 1 בדיוק כאשר $\mathrm{Re}\,s = 0$ — ציר הדמיון האנלוגי ממופה לעיגול היחידה, חצי המישור השמאלי פנימה באופן מוחלט. הקשר הטנגנסיאלי הוא ה עיוות: בקרבת תדר אפס $\tan\theta \approx \theta$ (עיוות זניח), אך תדר דיגיטלי המתקרב לניקוויסט מתכווץ לכיוון אינסוף אנלוגי. טבלאות בטרוורת', צ'בישב ואלפטיות מעוברות בשלמותן כי רק ציר התדר כופף — הגיאומטריה הרב-מודעית נותרה intact.

דוגמה פתורה 1 — הבייקוואד. מסנן נמוך בטרוורת', דרגת 2 100 הרץ, $f_{s} = 1$ קילו הרץ: מתאים מראש $\omega_{a} = 2000\tan(\pi\cdot 100/1000) = 649.8$ רדיאן/שנייה ✓. עם $s$ מנורמל על ידי $\omega_{a}$, קבוע ההצבה הוא $K = 1/\tan(0.1\pi) = 3.078$ והפרוטוטיפוס $H(s) = 1/(s^{2}+1.4142s+1)$ מתמזג ל-

דוגמה פתורה 2 — ודא, ואז תשלם. רווח DC: המכנה סוכם ל-$14.82 - 16.94 + 6.12 = 4.00$, ולכן $H(1) = 4/4 = 1.000$ ✓ (Butterworth שטוח ב-DC). התגובה בדיוק ב-100 Hz נחשבת ל- $-3.0$ dB ✓ — ה-prewarp הנחיש את נקודת ה-$-3$ dB בדיוק ב- התדר המובטח, שהוא העניין כולו של prewarp. הקוטבים יושבים ב-$|z| = 0.643$, בתוך מעגל היחידה כי המפה מבטיח את זה ✓. אל תעבור על העיוות המוקדם — עצב את הפילטר האנלוגי ב- $2\pi \times 100 = 628.3$ רדיאן/שנייה — ונקודת ה-$-3$ ד"ב תיפול ב-96.9 הרץ; set $T = 1$ as well ($K = 1$) and it cuts at 250 Hz, reading a complacent $-0.05$ dB at 100 Hz while missing the specification entirely — the המספרים תמיד מראים את ההבדל. ה 5 מקדמים מול ה-67 taps של ה-FIR, אך הפזה היא לא ליניארית בסמוך לגבול החיתוך ובנקודה קבועה, המקדמים של המכנה הם הרגישים ביותר — עיגול של $6.12$ מזיז את גבול החיתוך במידה מורגשת. יעילות, פזה ורגישות נמצאים במסחר באותה שורה.

Common mistakes. Skipping the prewarp (the cutoff lands low — 96.9 Hz for 100 Hz at 1 kHz is small, 2.86 kHz for 4 kHz at 10 kHz is fatal). הגדרת $T = 1$ בעוד שציטוט התדרים האנלוגיים עדיין ב-Hz. בדיקת יציבות על מזל במקום לסמוך על המפה. והתקוויות לפזה ליניארית — IIR השלב הקרוב לקוטבים הוא הפוך למישורי.

סיכום. דו-ליניאר: LHP → עיגול יחיד (ללא יציבות), תדרים מעוות על ידי $\frac{2}{T}\tan(\cdot)$ (עיוות מוקדם לכל תדר קריטי); הביקוואד של בטרבורת מתמקם כ- $H(z) = \frac{1+2z^{-1}+z^{-2}}{14.82 - 16.94z^{-1} + 6.12z^{-2}}$ (dc 1.000, $-3.0$ dB ב-100 Hz, מקטבים ב-0.643).

למה זה חשוב. קצב הדגימה הוא תקציב חישוב: מסנן הפועל ב- 48 kHz עולה שש פעמים יותר מאותו מסנן ב-8 kHz. רב-קצבי הוא האמנות של ביצוע העבודה שם שהיא זולה.

הרעיון, עם תמונה. צמצום ב-$M$ (שמירה על כל $M$-ת דגימה) חייב קודם לכן פס-נמוך מתחת ל-$f_{s}/(2M)$ — אחרת הספקטרום המתקפל יוצר אליאסים, לא ניתן לשחזור. השלמה ב-$L$ (הכנסת $L-1$ אפסים בין דגימות) חייבת לאחר מכן פס-נמוך מתחת לתדר נייקווסט הישן כדי לנטר את התמונות.

הגזירה. הכנסת אפסים ב-$L$ בזמן גורמת ל- $X_{up}(e^{j\omega}) = X(e^{j\omega L})$ — אותו ספקטרום מוחזר $L$ פעמים סביב המעגל: אלו הן התמונות, והן בדיוק ממשי כמו האות. שמירת דגימה אחת ב-1 ב-$M$ מעטפת $M$ העתקות ספקטרליות על עצמן: כל מה שקיים מעל תדר נייקווסט החדש מתקפל בחזרה אל הלהקה. שני מסנני מעבר נמוך אינם ביורוקרטיה; כל אחד מוסיף ודאות מתמטית. פוליפאזה היא אז אלגברה טהורה: ה- מסנן האינטרפולציה $L$ מחלק את הקבוצות הפאזיות של כל כניסה ל-$L$, כך שמסנן אחד ב-$f_{s}L$ הופך ל-$L$ תת-מסננים ב-$f_{s}$ — זהים בערכי האריתמטיקה, ו-$1/L$ מכפילים.

דוגמה פתורה 1 — קול דיבור ב-8 kHz ל-48 kHz. אינטרפולציה ב-6, ולאחר מכן מסנן low-pass at 4 kHz (the old Nyquist) removes the images at 8k, 16k, פוליפייז מארגן מחדש אותו ל-40 פילטרים משניים כל אחד פועל ב-48 קילוהרץ: ה- Polyphase reorganises it into 6 subfilters each running at 8 kHz: the same arithmetic, $1/6$ the multiplies ✓ — and this is not an הערכה, זו אותה חישוב שנעשתה בצורה מתוזמנת שונה.

Worked example 2 — one converter, not two. Resampling 48 kHz audio up to 72 kHz needs $L = 3$, $M = 2$ ($48\times 3/2 = 72$ ✓). Two passes would דורשים ממיר הרחבה, מסנן, ממיר צמצום ומסנן נוסף. ממיר rational converter does it: upsample by 3 to 144 kHz, apply a single מעבר נמוך בתדר 24 קילוהרץ (תדר ניקווי המכניס עצמו — הוא חוסם את $\uparrow 3$ התמונות ואת ה$\downarrow 2$ האליאסים בו-זמנית), שומר על כל זוגי טעויות נפוצות. דצימה לפני סינון (הalias then is not

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

סיכום. $\downarrow M$ דורש מסנן מעבר נמוך ב$f_{s}/2M$ תחילה; $\uparrow L$ דורש מסנן מעבר נמוך בתדר נייקוויסט הישן לאחר מכן; פולא-פייז מספק את אותו חישוב עבור $1/L$ הכופלים; ממירים רציונליים $L/M$ מתמזגים למסנן אחד בקצב הביניים.

למה זה חשוב. הקורס סוגר את העיגול שבו הפעולה פותחת: עם מפרט, מימוש ומדידה — והפער האמיתי בין השלושה.

הרעיון, עם איור. בנה גלאי טון כזרם עיבוד, כל שלב מאומת לבד לפני שקיים הבא: מסנן אנטי-אלייאזিং → חלון → FFT → סף. תקלה שנמצאת בשלב יחיד עולה שעה אחת; נמצאת במערכת המאוחדת, שבוע של חצייה.

השיטה. כתוב את המפרט תחילה — גלאי טון: 1 זריזות בתוך 1 שניות, זיהוי ב$-10$ dB SNR — ביצע בזרמי שלב מאומתים, ודווח ביצועים נמדדים מול כל שורה, כולל כישלונות.

דוגמה מפורשת 1 — המפרט מחזיר אש. 1 זריזות בתדר $f_{s} = 1$ ק"ה נדרשים $N \ge 1000$ דגימות — כלומר זהו כבר $1.0$ שניות של נתונים, וה- המפרט קובע תחת 1 שניות השהייה. הסתירה מופיעה בעת עיצוב, שזה כאשר אתה רוצה אותה: $N = 1024$ מספק 0.977 הרץ ✓ אך רשומה של 1.024 שניות ✏; $N = 1000$ מספק 1.000 הרץ — בדיוק על הקו. טבלת השוליים רשומה את ההתחלפות; אין שום דבר שקט.

דוגמה פותרת 2 — האם ניתן לזהות אפילו $-10$ dB? אנרגיה קוהרנטית של סינוסיד נופלת בתיבה אחת בעוד שרעש לבן מתפשט בכל $N/2$: הרווח העיבוד בתיבה זו הוא $10\log_{10}(N/2)$ — ב $N = 1024$ זהו $27$ dB ✓. יחס אות לרעש כניסה של $-10$ dB הופך לכ- $+17$ dB בתיבת החלטה: ניתן לזיהוי עם שוליים — חלון עולה במעט dB מהרווח הזה ועדיין עובר — והמספר הגיע מהמרת הפוריה, לא מהתקווה.

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

סיכום. הגדר מפרט תחילה, בדוק שלבים בנפרד, הסתמך על השוליים; רווח עיבוד $10\log_{10}(N/2)$ קובע את ההיתכנות.