This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
סיכום–עיבוד אותות ספרתי התמרת פורייה בזמן רציף.1
שיווין פרסבל.2
הדלתא של דיראק.3
קונבולוציה עם הלם.4
סכום פואסון.5הקשר שבין התמרת פורייה של שרשרת דגימות עם טור פורייה.5.1
)3/6, תרגול 26/5הורדת קצב דגימה (הרצאה .17M ונרצה להוריד את כמות הדגימות פי נתון אות דגום .17.1
כדי למנוע קיפול בתדר ותדר קטעון 1 עם הגבר LPFנעביר דרך מסנן .17.2
נדגום את האות הבדיד בתדר שזמן המחזור שלו הוא .17.3
)3/6, תרגול 2/6העלאת קצב דגימה (הרצאה .18L ונרצה להעלות את כמות הדגימות פי נתון אות דגום .18.1
נדגום את האות הבדיד בתדר שזמן המחזור שלו הוא .18.2
אפסים ביניהןL-1לוקחים את הדגימות המקוריות ושמים .18.3
לא נעביר דרך מסנן שהתגובה להלם שלו היא .18.4
זהו אינטרפולטור לא טוב, כיוון שכדי לקבל דיוק גבוה צריך לחשב הרבה איברים.18.4.1
לכן נעביר דרך מסנן שהתגובה להלם שלו היא.18.5
3
)3/6, תרגול 2/6 של קצב דגימה (הרצאה ישינוי רציונאל.19Lמעלים את קצב הדגימה ב-.19.1מעבירים דרך מסנן לביצוע אינטרפולציה.19.2מעבירים דרך מסנן להורדת תדרים גבוהים (ניתן להפוך את שני המסננים לאחד).19.3Mהורדת קצב הדגימה ב-.19.4
4
20.FT
21.DTFTעבור סדרה אינסופית.21.1מחזור הפונקציה הוא .21.2
22.DFT 17/6, תרגול 2/6 (הרצאה(Nעבור סדרה סופית באורך .22.1 מבצעיםDFTבכל חישוב כל סדרת איברי ה- .22.2
הכפלות מרוכבות.22.2.1 חיבורים מרוכבים.22.2.2
הערה.22.3
)10/6שוויון פרסבל לסדרה סופית (הרצאה .22.4
)26/5 (התמרת סדרה בדידה אין-סופית) (הרצאה DTFT (התמרת האות הרציף) ל-FTהקשר בין .23
)2/6 (התמרת סדרה בדידה אין-סופית) (הרצאה DTFT (התמרת סדרה בדידה סופית) ל-DFTהקשר בין .24)N איברים (סדרה סופית באורך Nאם נתונה סדרה אינסופית, אשר מורכבת מאפסים בלבד פרט לרצף של .24.1נדגום את הסדרה האינסופית ע"י הצבה.24.2
ומקבלים סדרה סופית.24.3 דגימותNניתן לראות את הסדרה כחלוקה של מעגל היחידה ל-.24.3.1
כיוון שעתהDFT אך תשפיע על ה-DTFT לא תשפיע על ה-M לסדרה סופית באורך Nשינוי אורך הסדרה הסופית באורך .24.4N ולא ל-Mמעגל היחידה יחולק ל-
)2/6הסטה מעגלית של סדרה (הרצאה .25Nנתונה סדרה באורך .25.1Nהסטה מעגלית של הסדרה שקולה להסטת האיברים וקיפול הסדרה על מעגל שהיקפו .25.2
)N+1 ימוקם עם האיבר ה 1 לסדרה (האיבר ה-Nחישוב ההסטה נעשה על בסיס מודולו .25.2.1 שמאלהN-1 ימינה שקול להסטה ב 1כלומר הסטה ב-.25.3DFTהסטה מעגלית של סדרה תחת .25.4
זהויות תחת הסטה מעגלית (עבור סדרה ממשית).25.5
5
)2/6קונבולוציה מעגלית (הרצאה .26 (חייבות להיות זהות באורכן)N סדרות באורך 2נתונות .26.1N=M+L סדרות באורך 2) , ניתן לרפד באפסים לקבלת M והשנייה L סדרות שונות באורכן (אחת באורך 2כדי להשוות בין .26.2
זהויות שכדאי לדעת.26.3
)10/6 באמצעות קונבולוציה מעגלית (הרצאה תחישוב קונבולוציה ליניארי.27 סדרות שונות באורכן2נתונות .27.1
27.1.1.x באורך M27.1.2.y באורך LN=M+Lנרפד באפסים כל סדרה כך ששתיהן יהיו באורך .27.1.3 לכל סדרה בנפרדDFTנחשב .27.1.4 זה בזה (זוגות של אינדקסים זה עם זה)DFTנכפיל את תוצאות ה-.27.1.5 נקודות הראשונות של הקונבולוציה הליניאריתN ב-DTFT על התוצאה ונקבל את הIDFTנבצע .27.1.6
)10/6 מתוך סדרה סופית (הרצאה DTFTחישוב .28Nנתונה סדרה סופית באורך .28.1
עד0 רגילה עם אינדקס מ-DTFT ע"י נוסחת כל האיברים מתאפסים, ניתן לחשב את Nמכיוון שהחל מהמקום ה-.28.2
N-1x[n] של IDFT נציב במשוואה את חישוב הx[n]במקום .28.3 :DTFT לבין ה-DFT נמצא יחס בין סדרת ה-–לאחר פישוט המשוואה .28.4
OVERLAP-ADDשיטת .32 בין סדרה קצרה לסדרה ארוכה מאוד (או אינסופית)תהרצון לחשב קונבולוציה ליניארי.32.1Nניתן לחלק סדרה אינסופית לאינסוף סדרות סופיות באורך .32.2
תהיהM באורך h לבין סידרה סופית x ת ציקלית בין סדרה אינסופיהקונבולוצי.32.3
(אשר יכולות להיות שונות באורכן) סדרות סופיות2קונבולוציה ציקלית בין חישובניתן לחשב כל חלק קונבולוציה בנפרד לפי .32.4
מבצעים חיבור של מול קונבולוציות 2 איברים בין כל M-1מכיוון שקיימת חפיפה בגודל .32.5
איברים חופפים.
)1/7 (תרגול OVERLAP-SAVEשיטת .33 . כאשרM באורך h ציקלית עם סדרה סופית ה וביצוע קונבולוציNניתן לחלק סדרה אינסופית לאינסוף סדרות סופיות באורך .33.1
N>M הסדרות השוות2, ונבצע קונבולוציה ציקלית בין N כך שתהיה באורך hנוסיף אפסים ל-.33.2 האיברים הראשונים של הקונבולוציה הציקלית מכילה קיפול של איברים, ולכן נזרוק אותםM-1בכל איטרציה נקבל ש.33.3 נשאיר ונחבר את הסדרות המתקבלותN-(M-1) = N-M+1את שאר האיברים, .33.4
)1/7אלגוריתם גרצל (תרגול .34 לדעת את ערכו), ללא חישוב כל הסדרה– מסוימת (לדעת אם יש משהו בתדר, ואם יש DFTחישוב נקודת .34.1
35.ZoomFFT 1/7 (תרגול()k0-k+1 עד k0 (כלומר מ-k רק עבור קטע מסויים באורך DTFTחישוב נקודת .35.135.2.N ,יהיה אורך הסדרה k יהיה מספר התדרים שמעוניינים בהם, לכן L=N/k
8
36.FFT 16/6 בזמן (הרצאה(שיטה זו חוסכת פעולות חישוב כיוון שהיא משתמשת בחישובים שכבר נעשו.36.1, כך שניתן לחצות את הסדרות עד לקבלת זוגות איברים2מניחים כי אורך הסדרה הוא חזקה של .36.2 סדרות: סדרת האיברים במקום הזוגי וסדרת האיברים במקום האי-זוגי2בדרך כלל חוצים ל-.36.3, ניתן לרפד באפסים את אורך הסדרה.2אם האורך אינו חזקה של .36.4 היא FFT מקדמים/תדרים כאשר סיבוכיות ה-N עבור חישוב היא DFTסיבוכיות ה-.36.5 שלהDFTכניסת התרשים היא הסדרה בזמן, ומוצאה הוא ה-.36.6הצבת כניסה מול יציאה הוא ע"י התאמת הביטוי הבינארי כך שסיביות הכניסה הן היפוך סדר סיביות היציאה.36.7
011 = 3 מול 110 = 6למשל, .36.7.1 הכפלות2 ובכך לבצע הכפלה יחידה במקום 1-ה"פרפר" הוא ביצוע הכפלה מרוכבת לפני הפיצול, ולאחר מכן ניתן להכפיל ב .36.8
מרוכבותבצומת מתבצע חיבור הקווים ועל הקו מתבצע כפל בקבוע.36.8.1
37.FFT 16/6 בתדר (הרצאה( שלהIDFTכניסת התרשים היא הסדרה בתדר, ומוצאה הוא ה-.37.1
9
)7/7 ניתוח התדר (הרצאה – חלונותחלון ריבועי.38
SINCהאות בתדר יהיה .38.1
Bartletחלון .39חלון משולש.39.1
Hannחלון .40N-2 המקדמים הקיצוניים מתאפסים בחלון זה, לכן סדר המסנן הוא 2 –חשוב .40.1
Hammingחלון .41אונות הצד קבועות.41.1
Blackmanחלון .42
Kaiser-Besselחלון .43 (היחס בין רוחב החלון לגובה האונה) ישפיע בעיקר על גובה אונות הצדβהמשתנה .43.1 ישפיע בעיקר על רוחב האונה המרכזיתMהמשתנה .43.2
10
Nאורך החלון = Mסדר המסנן =
רוחב אונההגדרה בזמן של החלוןמרכזית
חלון
-13.5dB0.090.7521מלבני
-41dB0.00220.01953Hamming
-31dB0.00630.05544Hann(raised
cosine)
-57dB0.00020.001774blackman
זוגי או א"ז( Nמשולש )הגדרה בזמן תלויה אם בתדר קל להיווכח כי זוהי ההתמרה של חלון מלבני בריבוע
-25dB0.050.4525bartlett
תלוי בפרמטר
0פונקצית בסל מסדר
Kaiser
11
)7/7 (הרצאה מדידת תדר של סינוסים/קוסינוסים.44נתון אות.44.1
Tsנדגום אותו בתדר .44.2
נכפיל אותו בחלון.44.3
לסדרת הדגימותDFTנבצע .44.4
ונקבל כאשר w = wmנביט בתדר מסויים .44.5
תהיה זניחהwm גדול מחצי רוחב האונה הראשית, אז ההפרעה לאות בתדר אם .44.6
)15/7, תרגול 21/7 (הרצאה FIR מסנן .61 אין קטביםFIRלמסנן .61.1
), למעט יוצאי דופןהאפסים מגיעים ברביעיות (.61.2
GLPטבלת תכונות .61.3IIIIIIIVטיפוסאי-זוגיזוגיאי-זוגיזוגיסדר
h[n]סימטריה של סביב מרכז המסנן
אנטי סימטריאנטי סימטריסימטריסימטרי
אנטי סימטריאנטי סימטריסימטריסימטריA(w)סימטריה של A(w)מחזוריות של
00
00שרירותישרירותי
שרירותי00שרירותי
הפרמטרים מחושבים לפי:Kaiserעבור חלון .61.3.1 מקבלים חלון מלבניכאשר .61.3.1.1
IRT בשיטת FIR תכנון מסנן .61.4 הרצויהתגובת המסנןמגדירים את .61.4.1
את טיפוס המסנן הרצוי מתוך אוסף הדרישות הנתונות GLP מתוך טבלת הבוחרים.61.4.2
לפי הבחירה שנעשה נוכל להגדיר את.61.4.2.1
תגובת התדר של המסנן האידיאלי.61.4.3
את מחשבים.61.4.4
As ואת Apנחשב את הנחתת האות ע"י חישוב .61.4.5
גבוה יותר מערכים שמצאנוAs נמוך יותר וערך Apמתוך טבלת החלונות נבחר חלון שנותן ערך .61.4.5.1נחשב את סדר המסנן.61.4.6
w = wp-wsנחשב את רוחב התדר הנחוץ .61.4.6.1
מהביטויMלפי המסנן שבחרנו, נדרוש שרוחב האונה המרכזית שלו יהיה קטן מרוחב התדר הנחוץ, ונחלץ את .61.4.6.2, עלינו להתאים את סדר המסנן הרצוי לטיפוס הנבחר (זוגי\אי-זוגי)GLPלפי טיפוס המסנן שבחרנו בטבלת ה.61.4.6.3
הגדרת איברי המסנן.61.4.7n = 0,1,…,Mקוצצים את סדרת איברי המסנן האידיאלי כך ש: .61.4.7.1
מציאת פונקצית המסנן.61.4.8 בפונקצית החלון שבחרנוMמציבים את .61.4.8.1כופלים את פונקצית החלון במסנן האידיאלי שבחרנו.61.4.8.2