מערכת ליניארית בלתי תלויה בזמן

מערכת ליניארית בלתי תלויה בזמן (באנגלית: Linear time-invariant system, הידועה גם בתור מערכת LTI) היא מושג במתמטיקה שימושית, שיש לו יישומים בתחומים כגון מעגלים חשמליים, עיבוד אותות, תורת הבקרה, ספקטרוסקופיית NMR , סייסמולוגיה ועוד.

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

== סקירה ==
הגדרת המאפיינים של כל מערכת LTI הן ליניאריות וקביעות (אי-תלות) בזמן.
* ליניאריות – הקשר בין הכניסה למוצא המערכת היא העתקה ליניארית: אם עבור הכניסה <math>x_1(t)\,</math> המוצא הוא <math>y_1(t)\,</math>, ועבור הכניסה <math>x_2(t)\,</math> המוצא הוא <math>y_2(t)\,</math>, אז עבור הצירוף הליניארי של הכניסות <math>a_1 x_1(t)+ a_2 x_2(t)\,</math> המוצא יהיה <math>a_1 y_1(t) + a_2y_2(t)\,</math>, עם אותם קבועים <math>a_1,a_2</math>. אפשר להרחיב תכונה זאת למספר רב יותר של כניסות, מוכפלות במספרים ממשיים <math>c_1, c_2, \ldots, c_k</math>,
*: כך שהכניסה <math>\sum_k c_k\,x_k(t)</math> מייצרת יציאה <math>\sum_k c_k\,y_k(t)\,</math>. (נוסחה 1)
*: בפרט:
*: הכניסה <math>\int_{-\infty}^{\infty} c_{\omega}\,x_{\omega}(t) \, \operatorname{d}\omega</math> מייצרת יציאה <math>\int_{-\infty}^{\infty} c_{\omega}\,y_{\omega}(t) \, \operatorname{d}\omega\,</math>
*: כאשר <math>c_{\omega}</math> ו־<math>x_{\omega}</math> הם סקלרים וכניסות, בהתאמה, אשר משתנים לאורך רצף אינדקסים <math>\omega</math>. לכן אם פונקציית כניסה יכולה להיות מיוצגת על ידי צירוף ליניארי של פונקציות כניסה, כפי שמוצג בדוגמה, אז פונקציית המוצא המתאימה יכולה להיות מיוצגת על ידי צירוף לילנארי של פונקציות יציאה מתאימות, באותה העתקה ליניארית שבכניסה.
* קביעות בזמן – עבור כניסה בהזזה של T שניות, המוצא יהיה זהה למעט עיכוב זמן של T שניות. כלומר, אם עבור הכניסה <math>x(t)</math> המוצא הוא <math>y(t)</math>, אז עבור הכניסה <math>x(t-T)</math> המוצא הוא <math>y(t-T)</math>. לפיכך, המערכת קבועה בזמן כי המוצא אינו תלוי בזמן שהכניסה נכנסת.

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

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

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

== מערכות רציפות ==
=== תגובה להלם וקונבולוציה ===

לכל מערכת LTI, ניתן להגדיר פונקציית תגובה להלם, שמבטאת את המוצא של המערכת עבור כניסה של פונקציית הלם <math>x(t) = \delta(t)</math> (הדלתא של דיראק), ומסומנת באות h: <math>y(t)=h(t)</math>.

ההתנהגות של מערכת ליניארית וקבועה בזמן עבור כניסה כללית <math>x(t) </math> ואות מוצא <math>y(t)</math> מתוארת על ידי אינטגרל הקונבולוציה עם התגובה להלם:
:<math>y(t) = (x * h)(t)\quad \stackrel{\mathrm{def}}{=}\ \int\limits_{-\infty}^{\infty} x(t-\tau)\cdot h(\tau) \, \mathrm{d}\tau\quad = \int\limits_{-\infty}^\infty x(\tau)\cdot h(t-\tau) \,\mathrm{d}\tau</math>,
כאשר המעבר נעשה על ידי שימוש בתכונת הקומוטטיביות.

עבור מערכת ליניארית, O צריכה לקיים את נוסחה 1:

:<math>O_t\left\{\int\limits_{-\infty}^\infty c_{\tau}\ x_{\tau}(u) \, \mathrm{d}\tau ;\ u\right\} = \int\limits_{-\infty}^\infty c_\tau\ \underbrace{y_\tau(t)}_{O_t\{x_\tau\}} \, \mathrm{d}\tau
</math>, (נוסחה 2)

מדרישת האי-תלות בזמן, היא מקיימת גם:

:<math>
O_t\{x(u-\tau);\ u\}\ \stackrel{\quad}{=}\ y(t-\tau)
\ \stackrel{\text{def}}{=}\ O_{t-\tau}\{x\}
</math>. (נוסחה 3)

כעת נוכל לרשום את תגובת ההלם כ- <math>\textstyle h(t) \ \stackrel{\text{def}}{=}\ O_t\{\delta(u);\ u\}</math>.

באופן דומה:
:<math>
h(t-\tau)\ \stackrel{\text{def}}{=}\ O_{t-\tau}\{\delta(u);\ u\}
= O_t\{\delta(u-\tau);\ u\}
</math>

כאשר הצורה של צד ימין של נוסחה 2 נשמרת עבור המקרה <math>\textstyle c_\tau = x(\tau)</math> ו- <math>\textstyle x_\tau(u) = \delta(u-\tau)</math>. נוסחה 2 מאפשרת את השלב הבא:

:<math>
(x * h)(t) = O_t\left\{\int\limits_{-\infty}^\infty x(\tau)\cdot \delta(u- \tau) \, \mathrm{d}\tau;\ u \right\}
= O_t\left\{x(u);\ u \right\}
\ \stackrel{\text{def}}{=}\ y(t).
</math>

לסיכום, קלט הפונקציה, <math>\textstyle \{x\}</math>, יכול להיות מיוצגת על ידי רצף של פונקציות הלם מוזזות בזמן בצירוף ליניארי. תכונת הליניאריות של המערכת מאפשרת לייצג את תגובת המערכת על ידי רצף של תגובות הלם מתאימות, מסוכמות באותה הדרך, ותכונת הקביעות בזמן מאפשרת לייצג את הצירוף על ידי אינטגרל הקונבולוציה.

=== אקספוננטים כפונקציות עצמיות ===
פונקציה עצמית היא פונקציה שעבורה פעולת האופרטור היא מכפלה של אותה פונקציה בקבוע. זאת אומרת, f היא פונקציה עצמית של <math>\hat\mathcal{H}</math>, אם
:<math>\hat\mathcal{H}f = \lambda f,</math>
כאשר <math>\lambda</math> הוא קבוע מספרי, שנקרא ערך עצמי.

פונקציות האקספוננט <math>A e^{s t}</math>, כאשר <math>A, s \in \mathbb{C}</math>, הן פונקציות עצמיות של אופרטור ליניארי וקבוע בזמן.

ניתן להוכיח זאת באופן הבא: נניח שהכניסה היא <math>x(t) = A e^{s t}</math>. היציאה של המערכת היא הקונבולוציה עם התגובה להלם <math>\int_{-\infty}^\infty h(t - \tau) A e^{s \tau}\, \operatorname{d} \tau</math>

אשר לפי תכונת הקומוטטיביות של הקונבולוציה, שקולה ל:

:<math>\begin{align}
\overbrace{\int_{-\infty}^\infty h(\tau) \, A e^{s (t - \tau)} \, \operatorname{d} \tau}^{\mathcal{H} f}
&= \int_{-\infty}^\infty h(\tau) \, A e^{s t} e^{-s \tau} \, \operatorname{d} \tau \\
&= A e^{s t} \int_{-\infty}^{\infty} h(\tau) \, e^{-s \tau} \, \operatorname{d} \tau \\
&= \overbrace{\underbrace{A e^{s t}}_{\text{Input}}}^{f} \overbrace{\underbrace{H(s)}_{\text{Scalar}}}^{\lambda},
\end{align}</math>

כאשר הסקלר
:<math>H(s)\ \stackrel{\text{def}}{=}\ \int_{-\infty}^\infty h(t) e^{-s t} \, \operatorname{d} t</math>
תלוי רק בפרמטר s, ולא במשתנה t.

תגובת המערכת היא העתקה ליניארית של הכניסה. בפרט, עבור כל <math>A,s\in\Complex</math> היציאה של המערכת היא המכפלה של הכניסה <math>Ae^{st}</math> והקבוע <math>H(s)</math>. לפיכך, <math>Ae^{st}</math> היא פונקציה עצמית של מערכת ליניארית בלתי תלויה בזמן (LTI), והערך העצמי המתקבל הוא <math>H(s)</math>.

====דוגמה====
נאמר כי <math>v(t) = e^{i \omega t}</math> אקספוננט מדומה (<math>i^2\ \stackrel{\text{def}}{=} -1</math>) והגרסה המוזזת בזמן שלו הוא: <math>v_a(t) = e^{i \omega (t+a)}</math> .

<math>H[v_a](t) = e^{i\omega a} H[v](t)</math> כהעתקה ליניארית של האקספוננט לפי הקבוע <math>e^{i \omega a}</math> .

<math>H[v_a](t) = H[v](t+a)</math> לפי קביעות בזמן לפי <math>H</math> .

אז, <math>H[v](t+a) = e^{i \omega a} H[v](t)</math> . עבור <math>t = 0</math> ונקבל:<math display="block">H[v](\tau) = e^{i\omega \tau} H[v](0)</math>

כלומר, עבור אקספוננט מדומה <math>e^{i \omega \tau}</math> ככניסה, נקבל אקספוננט מדומה עם אותו תדר במוצא.

===התמרת פורייה והתמרת לפלס===

תכונת הפונקציה העצמית של אקספוננטים שימושית מאוד כדי לנתח את המערכת בכלים אנילטיים. במסגרת ניתוח זה פותחו שתי התמרות המאפשרות להתמיר את פונקציית הכניסה והמוצא למרחב אחר שבו נוח יותר לעבוד. התמרת לפלס:
:<math> \mathcal{L}\{x(t)\}\ \stackrel{\text{def}}{=}\ \int_{-\infty}^\infty x(t) e^{-s t} \, \operatorname{d} t</math>
כלומר התמרת לפלס של פונקציה היא הפריסה שלה על גבי הפונקציות העצמיות.

בפרט עבור פונקציית התגובה להלם:
:<math>H(s)\ \stackrel{\text{def}}{=}\ \mathcal{L}\{h(t)\}\ \stackrel{\text{def}}{=}\ \int_{-\infty}^\infty h(t) e^{-s t} \, \operatorname{d} t</math>
היא בדיוק הדרך למצוא את הערכים העצמיים מתוך התגובה להלם.

התמרה נוספת נוצרת כאשר מתמקדים בפונקציות סינוסואידיות טהורות (למשל פונקציות אקספוננטים מדומים מהצורה <math>e^{j \omega t}</math> כאשר <math>\omega \in \mathbb{R}</math> ו- <math>j^2\ \stackrel{\text{def}}{=} -1</math>). התמרת פורייה <math>H(j \omega) = \mathcal{F}\{h(t)\}</math> מביאה את הערכים העצמיים עבור פונקציות סינוסואידיות מדומות טהורות. שתי הפונקציות <math>H(s)</math> ו- <math>H(j\omega)</math> נקראות "תגובת המערכת", "פונקציית המערכת" או פונקציית המעבר.

התמרת לפלס בדרך כלל משמשת אותות חד-צדדיים, למשל, אותות ששווים לאפס עבור כל <math>t<t_0</math>, כאשר <math>t_0</math> ערך קבוע כלשהו.

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

:<math>y(t) = (h*x)(t)\ \stackrel{\text{def}}{=}\ \int_{-\infty}^\infty h(t - \tau) x(\tau) \, \operatorname{d} \tau\ =\ \mathcal{L}^{-1}\{H(s)X(s)\}</math>

===דוגמאות===
*דוגמה פשוטה לכך היא אופרטור הנגזרת:
** <math> \frac{\operatorname{d}}{\operatorname{d}t} \left( c_1 x_1(t) + c_2 x_2(t) \right) = c_1 x'_1(t) + c_2 x'_2(t) </math> (כלומר האופרטור ליניארי)
** <math> \frac{\operatorname{d}}{\operatorname{d}t} x(t-\tau) = x'(t-\tau) </math> (כלומר האופרטור קבוע בזמן)
:כאשר מתמירים את הנגזרת בהמרת לפלס, הנגזרת מותמרת להתמרת הפונקציה המקורית מוכפלת במשתנה הלפלס, s.
::<math> \mathcal{L}\left\{\frac{\operatorname{d}}{\operatorname{d}t}x(t)\right\} = s X(s) </math>
:הפשטות של התמרת הנגזרת מעידה על למה התמרת לפלס יעילה על פונקציות מעבר.
*עוד דוגמה פשוטה היא אופרטור הממוצע
::<math> \mathcal{A}\left\{x(t)\right\}\ \stackrel{\text{def}}{=}\ \int_{t-a}^{t+a} x(\lambda) \, \operatorname{d} \lambda. </math>
:מליניאריות האינטגרל,
::<math>\begin{align}
\mathcal{A} \{c_1 x_1(t) + c_2 x_2(t)\}
&= \int_{t-a}^{t+a} ( c_1 x_1(\lambda) + c_2 x_2(\lambda)) \, \operatorname{d} \lambda\\
&= c_1 \int_{t-a}^{t+a} x_1(\lambda) \, \operatorname{d} \lambda + c_2 \int_{t-a}^{t+a} x_2(\lambda) \, \operatorname{d} \lambda\\
&= c_1 \mathcal{A}\{x_1(t)\} + c_2 \mathcal{A} \{x_2(t) \},
\end{align}</math>
::ליניארי. בנוסף, מפני ש-
::<math>\begin{align}
\mathcal{A}\left\{x(t-\tau)\right\}
&= \int_{t-a}^{t+a} x(\lambda-\tau) \, \operatorname{d} \lambda\\
&= \int_{(t-\tau)-a}^{(t-\tau)+a} x(\xi) \, \operatorname{d} \xi\\
&= \mathcal{A}\{x\}(t-\tau),
\end{align}</math>
::קבוע בזמן. למען האמת, <math>\mathcal{A}</math> יכולה להיות מיוצגת כקונבולוציה עם פונקציית המלבן <math>\Pi(t)</math>. כלומר,
::<math> \mathcal{A}\left\{x(t)\right\} = \int_{-\infty}^\infty \Pi\left(\frac{\lambda-t}{2a}\right) x(\lambda) \, \operatorname{d} \lambda, </math>
:כאשר פונקציית המלבן היא:
::<math>\Pi(t)\ \stackrel{\text{def}}{=}\ \begin{cases} 1 &\text{if } |t| \leq \frac{1}{2},\\ 0 &\text{if } |t| > \frac{1}{2}.\end{cases}</math>

===תכונות ===
חלק מהתכונות החשובות ביותר של מערכת הן סיבתיות ויציבות.

====סיבתיות====
מערכת נקראת "סיבתית" כאשר מוצא המערכת נקבע אך ורק על פי ערכי הכניסה בעבר ובהווה, ולא על פי ערכי הכניסה בעתיד.

עבור התגובה להלם, מערכת כזאת מתקבלת רק כאשר לכל t שלילי, התגובה מתאפסת, משום שפונקציית ההלם שהיא הכניסה, גם היא מתאפסת בתחום זה. בייצוג מתמטי:
:<math>\forall t < 0 : h(t) = 0 \Rightarrow \text{causal system}</math>

כאשר <math>h(t)</math> הוא התגובה להלם.

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

====יציבות====

מערכת נקראת יציבה במובן BIBO אם לכל כניסה חסומה יתקבל מוצא סופי (כלומר, חסום גם כן). בייצוג מתמטי, אם כל כניסה שמקיימת:
:<math>\exists k(\ \|x(t)\|_{\infty}<k)</math>

מתקבל מוצא שמקיים:
:<math>\exists m(\ \|y(t)\|_{\infty}<m)</math>
:(כלומר, המקסימום של <math>|x(t)|</math> סופית גוררת שהמקסימום של <math>|y(t)|</math> סופית גם כן), לכן המערכת יציבה.
תנאי הכרחי ומספיק הוא שהתגובה להלם, <math>h(t)</math>, נמצאת במרחב L<sup>1</sup> (בעלת נורמת <math>L^1</math> סופית):
:<math>\exists k(\ \|h(t)\|_1 = \int_{-\infty}^\infty |h(t)| \, \operatorname{d}t <k)</math>

במישור התדר, תחום ההתכנסות חייבת לכלול את הציר המדומה <math>s=j\omega</math>.

לדוגמה, מסנן מעביר תדרים נמוכים (LPF) אידיאלי עם תגובה להלם של פונקציית sinc היא לא יציבה במובן BIBO, מפני שנורמת L<sup>1</sup> שלה לא סופית. לכן, עבור כניסה חסומה לא הכרחי שנקבל מוצא חסום. בפרט, עבור הכניסה אפס ל-<math>t < 0\,</math> ושווה לסינוסואיד בתדר cutoff ל־<math>t > 0\,</math>, אז המוצא לא יהיה חסום לכל זמן מלבד אלה שחוצים את האפס.

==מערכות בדידות==
כמעט לכל מערכת בזמן רציף יש חלק משלים במערכת בזמן בדיד.

===מערכות בדידות ממערכות רציפות===

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

אם <math>x(t)</math> הוא אות אנלוגי רציף, ממיר אנלוגי לדיגיטלי יהפוך אותו לאות דיגיטלי DT:

:<math>x[n] \ \stackrel{\text{def}}{=}\ x(nT), \qquad\, n \in \mathbb{Z}</math>

כאשר T זה זמן הדגימה. חשוב להגביל את טווח התדרים של אות הכניסה לייצוג נאמן באות בדיד, כך שמשפט הדגימה יבטיח שאף מידע של האות הרציף לא יאבד. אות בדיד יכול להכיל טווח תדרים רק עד <math>\frac{1}{2T}</math>. תדרים אחרים "יתקפלו" לאותו טווח.

===תגובה להלם וקונבולוציה===
יהי <math>\{x[m-k];\ m\}</math> ייצוג של הרצף <math>\{x[m-k]|\forall m \in \Z \}</math>.
יהי הייצוג הקצר יותר <math>\{x\}\,</math> תייצג <math>\{x[m];\ m\}</math>.

מערכת בדידה הופכת רצף קלט <math>\{x\}</math>, לרצף פלט <math>\{y\}</math>. באופן כללי, כל אלמנט של הפלט תלוי בכל רכיב הקלט. נייצג את אופרטור ההעתקה על ידי <math>O</math>, ועתה אנו יכולים לכתוב:

:<math>y[n] \ \stackrel{\text{def}}{=}\ O_n\{x\}</math>

נשים לב שאם ההתמרה עצמה משתנה עם n, רצף הפלט הוא קבוע בלבד, והמערכת אינה מעניינת. במערכת טיפוסית, [y[n "תלויה בדרך כלל מהאלמנטים של x שהאינדקסים שלהם קרובים ל- n.

עבור מקרה מיוחד של הדלתא של קרונקר, <math>x[m]=\delta[m]</math>, רצף המוצא יהיה התגובה להלם:

:<math>h[n] \ \stackrel{\text{def}}{=}\ O_n\{\delta[m];\ m\}.\,</math>

עבור מערכת ליניארית, <math>O</math> צריכה לקיים את:
:<math>
O_n\left\{\sum_{k=-\infty}^{\infty} c_k\cdot x_k[m];\ m\right\} =
\sum_{k=-\infty}^{\infty} c_k\cdot O_n\{x_k\}.
\,</math> (נוסחה 4)

והתנאי לקביעות בזמן:
:<math>
\begin{align}
O_n\{x[m-k];\ m\}\ &\stackrel{\quad}{=}\ y[n-k]\\
&\stackrel{\text{def}}{=}\ O_{n-k}\{x\}.\,
\end{align}
</math> (נוסחה 5)

במערכת כזאת, התגובה להלם, <math>\{h\},\,</math>, מתארת את המערכת לגמרי. כלומר, עבור כל רצף קלט, רצף הפלט יכול להיות מחושב במונחים של קלט התגובה להלם. כדי לראות איך זה נעשה, יש לשקול את הזהות:

:<math>
x[m] \equiv \sum_{k=-\infty}^{\infty} x[k]\cdot \delta[m-k],
\,</math>

אשר מייצגת את <math>\{x\}\,</math> כסכום של דלתאות במושקלים.
לכן,

:<math>
\begin{align}
y[n] = O_n\{x\}
&= O_n\left\{\sum_{k=-\infty}^\infty x[k]\cdot \delta[m-k];\ m \right\}\\
&= \sum_{k=-\infty}^\infty x[k]\cdot O_n\{\delta[m-k];\ m\},\,
\end{align}
</math>
כאשר עוררנו את נוסחה 4 עבור המקרה <math>c_k=x[k]</math> ו-<math>x_k[m]=\delta[m-k]</math>.

ובגלל נוסחה 5, נוכל לכתוב:

:<math>
\begin{align}
O_n\{\delta[m-k];\ m\}\ &\stackrel{\quad}{=}\ O_{n-k}\{\delta[m];\ m\}\\
&\stackrel{\text{def}}{=}\ h[n-k].\,
\end{align}
</math>

לכן:

:{|
|-
|<math>= \sum_{k=-\infty}^{\infty} x[k]\cdot h[n-k]\,</math>
|<math>y[n]\,</math>
|-
|
|<math>= \sum_{k=-\infty}^{\infty} x[n-k]\cdot h[k],</math> (קומוטטיביות)
|} אשר מניבה את הנוסחה המוכרת עבור קונבולוציה בדידה.

=== אקספוננטים כפונקציות עצמיות ===
פונקציה עצמית היא פונקציה שעבורה הפלט של האופרטור היא מכפלה של אותה פונקציה בקבוע. זאת אומרת, f היא פונקציה עצמית של <math>\mathcal{H}</math>, אם
:<math>\mathcal{H}f = \lambda f,</math>
כאשר <math>\lambda</math> הוא קבוע, שנקרא בהקשר זה ערך עצמי.

את פונקציות האקספוננט <math>z^n = e^{sT n}</math>, כאשר <math>n \in \mathbb{Z}</math>, הם פונקציות עצמיות של אופרטור ליניארי, וקבוע בזמן. <math>T \in \mathbb{R}</math> הוא קצב הדגימה, וכן <math>z=e^{sT}</math>, <math>z</math>, <math>s\in\Complex</math>. הוכחה פשוטה ממחישה את המושג הזה.

נניח שהכניסה היא <math>x[n] = \,\!z^n</math> . היציאה של המערכת עבור תגובה להלם <math>h[n]</math> יהיה:
:<math>\sum_{m=-\infty}^{\infty} h[n-m] \, z^m</math>

כאשר לפי תכונת הקומוטטיביות של הקונבולוציה, שקולה ל:

:<math>\sum_{m=-\infty}^{\infty} h[m] \, z^{(n - m)} = z^n \sum_{m=-\infty}^{\infty} h[m] \, z^{-m} = z^n H(z)</math>

כאשר הסקלר
:<math>H(z)\ \stackrel{\text{def}}{=}\ \sum_{m=-\infty}^\infty h[m] z^{-m}</math>
תלוי רק בפרמטר z.

אז <math>z^n</math> היא פונקציה עצמית של מערכת LTI מפני שתגובת המערכת שקולה לקלט כפול הקבוע <math>H(z)</math>.

===התמרת Z והתמרת פורייה הבדידה===
תכונת הפונקציה העצמית של אקספוננטים שימושית מאוד כדי להבין יותר מערכות LTI. התמרת Z:
:<math>H(z) = \mathcal{Z}\{h[n]\} = \sum_{n=-\infty}^\infty h[n] z^{-n}</math>
היא בדיוק הדרך למצוא את הערכים העצמיים מתוך התגובה להלם. בפרט, נתעניין בפונקציות סינוסואידיות טהורות (למשל אקספוננטים מרוכבים מהצורה <math>e^{j \omega n}</math> כאשר <math>\omega \in \mathbb{R}</math> ו- <math>j^2\ \stackrel{\text{def}}{=} -1</math>). התמרת פורייה הבדידה (DTFT) <math>H(j \omega) = \mathcal{F}\{h[n]\}</math> מביאה את הערכים העצמיים עבור פונקציות סינוסואידיות מדומות טהורות. שתי הפונקציות <math>H[z]</math> ו- <math>H(e^{j\omega})</math> נקראות "תגובת המערכת", "פונקציית המערכת" או פונקציית המעבר.

התמרת Z בדרך כלל מופעלת על אותות חד-צדדיים, למשל, אות ששווה לאפס עבור כל <math>t<t_0</math>, כאשר <math>t_0</math> ערך קבוע כלשהו.

בגלל תכונת הקונבולוציה של שתי ההתמרות, הקונבולוציה שמביאה את מוצא המערכת יכולה להיות מיוצגת על ידי מכפלה במישור התדר, עבור מערכות אשר קיימות במישור התדר:
:<math>y[n] = (h*x)[n] = \sum_{m=-\infty}^\infty h[n-m] x[m] = \mathcal{Z}^{-1}\{H(z)X(z)\}</math>

===דוגמאות===
*דוגמה פשוטה לכך היא אופרטור העיכוב (delay):
*:<math>D\{x[n]\} \ \stackrel{\text{def}}{=}\ x[n-1]</math>
**<math> D \left( c_1 \cdot x_1[n] + c_2 \cdot x_2[n] \right) = c_1 \cdot x_1[n-1] + c_2\cdot x_2[n-1] = c_1\cdot Dx_1[n] + c_2\cdot Dx_2[n]</math> (כלומר האופרטור ליניארי)
**<math> D\{x[n-m]\} = x[n-m-1] = x[(n-1)-m] = D\{x\}[n-m]\,</math> (כלומר האופרטור קבוע בזמן)
:כאשר מתמירים את הנגזרת בהמרת Z, הנגזרת מותמרת להתמרת הפונקציה המקורית מוכפלת במשתנה, z<sup>−1</sup>. כלומר,
::<math> \mathcal{Z}\left\{Dx[n]\right\} = z^{-1} X(z). </math>
*עוד דוגמה פשוטה היא אופרטור הממוצע
::<math> \mathcal{A}\left\{x[n]\right\}\ \stackrel{\text{def}}{=}\ \sum_{k=n-a}^{n+a} x[k]</math>
:מליניאריות הסכום,
::<math>\begin{align}
\mathcal{A}\left\{c_1 x_1[n] + c_2 x_2[n] \right\}
&= \sum_{k=n-a}^{n+a} \left( c_1 x_1[k] + c_2 x_2[k] \right)\\
&= c_1 \sum_{k=n-a}^{n+a} x_1[k] + c_2 \sum_{k=n-a}^{n+a} x_2[k]\\
&= c_1 \mathcal{A}\left\{x_1[n] \right\} + c_2 \mathcal{A}\left\{x_2[n] \right\},
\end{align}</math>
::ליניארי. בנוסף, מפני ש-
::<math>\begin{align}
\mathcal{A}\left\{x[n-m]\right\}
&= \sum_{k=n-a}^{n+a} x[k-m]\\
&= \sum_{k'=(n-m)-a}^{(n-m)+a} x[k']\\
&= \mathcal{A}\left\{x\right\}[n-m],
\end{align}</math>
::קבוע בזמן.

===תכונות===
ניתן לתאר באופן מלא את הקשר שבין הקלט לפלט במערכת LTI בדידה על ידי תגובת ההלם של המערכת <math>h[n]</math>. חלק מהתכונות החשובות ביותר של מערכת הן סיבתיות ויציבות. בניגוד למערכות בזמן רציף, מערכות בדידות לא סיבתיות יכולות להתממש. זה טריוויאלי כדי להפוך מערכת עם תגובה להלם סופית סיבתי על ידי הוספת עיכובים (Delays). זה אפשרי אפילו לעשות מערכות שלא סיבתיות עם תגובה להלם אינסופית.

====סיבתיות====
מערכת בדידה נקראת "סיבתית" כאשר מוצא המערכת נקבעת אך ורק על פי ערכי הכניסה בעבר ובהווה.

תנאי הכרחי ומספיק למערכת סיבתית:

:<math>\forall n(n < 0 \Rightarrow h[n] = 0)</math>

כאשר <math>h[n]</math> הוא התגובה להלם. באופן כללי, לא ניתן לדעת האם המערכת היא סיבתית מתוך התמרת Z, מפני שההתמרה ההפוכה אינה יחידה. רק כאשר התחום התכנסות ידוע, ניתן לדעת האם המערכת סיבתית.

====יציבות====
מערכת נקראת יציבה במובן BIBO אם לכל כניסה חסומה מתקבל מוצא סופי (כלומר, חסום גם כן). בייצוג מתמטי, אם כל כניסה שמקיימת:
:<math>\exists k (\ \|x[n]\|_{\infty} < k)</math>
מתקבל מוצא שמקיים:
:<math>\exists k(\ \|y[n]\|_{\infty} < k)</math>
:(כלומר, המקסימום של <math>|x[n]|</math> סופית גוררת שהמקסימום של <math>|y[n]|</math> סופית גם כן), לכן המערכת יציבה. תנאי הכרחי ומספיק הוא שהתגובה להלם, <math>h[n]</math>, נמצאת בL<sup>1</sup> (בעלת נורמת<math>L^1</math> סופית):

:<math>\exists k(\|h[n]\|_1\ \stackrel{\text{def}}{=}\ \sum_{n = -\infty}^\infty |h[n]| < k)</math>

במישור התדר, תחום ההתכנסות חייבת לכלול את מעגל היחידה (כלומר, המקום הגאומטרי המקיים <math>|z|=1</math> עבור <math>z</math> מרוכב)

== ראו גם ==
* תגובת הלם
* פונקציית גרין

== לקריאה נוספת ==
<div class="mw-content-ltr">
* 
* 
* 
* 
</div>

==קישורים חיצוניים==
* [http://www.tedpavlic.com/teaching/osu/ece209/support/circuits_sys_review.pdf ECE 209: Review of Circuits as LTI Systems] – קדימון קצר על ניתוח מתמטי של מערכות LTI חשמליות
* [http://www.tedpavlic.com/teaching/osu/ece209/lab3_opamp_FO/lab3_opamp_FO_phase_shift.pdf ECE 209: Sources of Phase Shift] – הסבר אינטואיטיבי של המקור של הזזות בפאזה בשתי מערכות LTI חשמליות משותפות
*[http://www.ece.jhu.edu/~cooper/courses/214/signalsandsystemsnotes.pdf JHU 520.214 Signals and Systems course notes] – קורס מקוצר על תורת מערכת ה־LTI, הולם להוראה עצמית
*[http://www.etti.unibw.de/labalive/tutorial/lti/ LTI system example: RC low-pass filter] – תגובת ההגבר ותגובת הפאזה

==הערות שוליים==

קטגוריה:הנדסת חשמל
קטגוריה:עיבוד אותות
קטגוריה:עיבוד אותות ספרתי
קטגוריה:מתמטיקה שימושית