KonkurStudy
همه کارت‌ها

سیستم‌های LTI، پاسخ ضربه و پاسخ پله، و کانولوشن

یک سیستم خطی تغییرناپذیر با زمان به‌طور کامل با پاسخ ضربه‌اش hh توصیف می‌شود (یا به‌طور معادل با پاسخ پله‌اش، که مشتق یا تفاضل مرتبهٔ اول آن همان hh است)، خروجی آن کانولوشن ورودی با hh است، و نمایی همیشگی este^{st} تابع ویژهٔ آن است.

کارت مفهوم کنکور ارشد مهندسی برق

اهمیت
۰/۵

۱. ایدهٔ اصلی: یک تابع کل سیستم را توصیف می‌کند

یک سیستم پیوسته‌زمان y=T{x}y = \mathcal{T}\{x\} خطی است اگر از اصل برهم‌نهی پیروی کند — T{a x1+b x2}=a T{x1}+b T{x2}\mathcal{T}\{a\,x_1 + b\,x_2\} = a\,\mathcal{T}\{x_1\} + b\,\mathcal{T}\{x_2\} — و تغییرناپذیر با زمان است اگر تأخیر در ورودی صرفاً خروجی را به همان اندازه به تأخیر بیندازد: x(t−t0)→y(t−t0)x(t-t_0) \to y(t-t_0). سیستمی که هر دو ویژگی را دارد LTI نامیده می‌شود.

دلیل اهمیت محوری سیستم‌های LTI خاصیت غربالی است: هر ورودی را می‌توان به‌صورت برهم‌نهی ضربه‌های انتقال‌یافته نوشت،

x(t)=∫−∞∞x(τ) δ(t−τ) dτ.x(t) = \int_{-\infty}^{\infty} x(\tau)\,\delta(t-\tau)\,d\tau.

عملگر را اعمال کنید. خطی‌بودن، T\mathcal{T} را به درون انتگرال می‌برد (وزن‌های x(τ)x(\tau) نسبت به tt ثابت‌اند)؛ تغییرناپذیری با زمان می‌گوید پاسخ به δ(t−τ)\delta(t-\tau) همان پاسخ به δ(t)\delta(t) است که به اندازهٔ τ\tau انتقال یافته است. پاسخ ضربه را چنین تعریف کنید

h(t)  =  T{δ(t)}.h(t) \;=\; \mathcal{T}\{\delta(t)\}.

آنگاه پاسخ به δ(t−τ)\delta(t-\tau) برابر h(t−τ)h(t-\tau) است، و

y(t)=∫−∞∞x(τ) h(t−τ) dτ=(x∗h)(t).y(t) = \int_{-\infty}^{\infty} x(\tau)\,h(t-\tau)\,d\tau = (x * h)(t).

پس h(t)h(t) به‌تنهایی پاسخ سیستم به هر ورودی را تعیین می‌کند. دانستن hh یعنی شناختن سیستم. همین مطلب دربارهٔ هر یک زوج ورودی–خروجی معلوم که بتوان hh را از آن بازیابی کرد نیز درست است — رایج‌ترین نمونه پاسخ پله است (بخش ۲.۳).

۲. روابط کلیدی

۲.۱ کانولوشن (قانون ورودی–خروجی)

y(t)=(x∗h)(t)=∫−∞∞x(τ) h(t−τ) dτ,y[n]=∑k=−∞∞x[k] h[n−k].y(t) = (x * h)(t) = \int_{-\infty}^{\infty} x(\tau)\,h(t-\tau)\,d\tau,\qquad y[n] = \sum_{k=-\infty}^{\infty} x[k]\,h[n-k].

کانولوشن جابه‌جایی‌پذیر (x∗h=h∗xx*h = h*x)، شرکت‌پذیر ((x∗h1)∗h2=x∗(h1∗h2)(x*h_1)*h_2 = x*(h_1*h_2) — اتصال سری) و توزیع‌پذیر (x∗(h1+h2)x*(h_1+h_2) — اتصال موازی) است. ضربه عضو همانی آن است: x∗δ=xx * \delta = x، و x(t)∗δ(t−t0)=x(t−t0)x(t) * \delta(t-t_0) = x(t-t_0) (کانولوشن با ضربهٔ انتقال‌یافته فقط سیگنال را انتقال می‌دهد).

۲.۲ توابع ویژه و تابع سیستم

برای یک سیستم LTI، نمایی مختلط همیشگی este^{st} (که برای همهٔ ttها تعریف شده و با u(t)u(t) روشن نمی‌شود) یک تابع ویژه است:

y(t)=∫−∞∞h(τ) es(t−τ) dτ=est∫−∞∞h(τ) e−sτ dτ⏟H(s)=H(s) est.y(t) = \int_{-\infty}^{\infty} h(\tau)\,e^{s(t-\tau)}\,d\tau = e^{st}\underbrace{\int_{-\infty}^{\infty} h(\tau)\,e^{-s\tau}\,d\tau}_{H(s)} = H(s)\,e^{st}.

H(s)H(s) تابع سیستم (تابع تبدیل) است، یعنی تبدیل لاپلاس hh. برای سیستمی که با یک معادلهٔ دیفرانسیل با ضرایب ثابت ∑kak y(k)=∑kbk x(k)\sum_k a_k\,y^{(k)} = \sum_k b_k\,x^{(k)} داده شده است، جایگزینی ddt→s\tfrac{d}{dt}\to s مستقیماً H(s)=∑kbksk∑kakskH(s) = \dfrac{\sum_k b_k s^{k}}{\sum_k a_k s^{k}} را می‌دهد. با قرار دادن s=jωs = j\omega پاسخ فرکانسی H(jω)H(j\omega) به دست می‌آید، و ejωt→H(jω) ejωte^{j\omega t} \to H(j\omega)\,e^{j\omega t}.

حالت خاص s=0s=0 همان بهرهٔ DC است: ورودی ثابت x(t)=cx(t)=c که در تمام زمان‌ها اعمال شده باشد همان نمایی همیشگی c e0tc\,e^{0t} است، پس خروجی مقدار ثابت زیر است

y(t)=c H(0),H(0)=∫−∞∞h(t) dt(discrete time: H(z)∣z=1=∑nh[n]),y(t) = c\,H(0), \qquad H(0) = \int_{-\infty}^{\infty} h(t)\,dt \quad\Big(\text{discrete time: } H(z)\big\vert_{z=1} = \sum_{n} h[n]\Big),

که هرگاه آن انتگرال (مجموع) همگرا باشد معتبر است، به‌ویژه برای هر سیستم پایدار BIBO.

۲.۳ پاسخ پله و چگونگی بازیابی hh از آن

پاسخ پله ss خروجی سیستم به ازای ورودی پلهٔ واحد است. کانولوشن hh با پله یک انتگرال انباشته (مجموع انباشته) است:

s(t)=(u∗h)(t)=∫−∞th(τ) dτ,s[n]=(u∗h)[n]=∑k=−∞nh[k].s(t) = (u*h)(t) = \int_{-\infty}^{t} h(\tau)\,d\tau, \qquad s[n] = (u*h)[n] = \sum_{k=-\infty}^{n} h[k].

با وارون‌کردن این روابط، پاسخ ضربه از پاسخ پلهٔ اندازه‌گیری‌شده یا داده‌شده به دست می‌آید:

h(t)=ds(t)dt,h[n]=s[n]−s[n−1].h(t) = \frac{ds(t)}{dt}, \qquad h[n] = s[n] - s[n-1].

رابطهٔ گسسته چیزی نیست جز اعمال خطی‌بودن و تغییرناپذیری با زمان بر δ[n]=u[n]−u[n−1]\delta[n] = u[n] - u[n-1]. هر دو رابطه نمونه‌هایی از یک قاعدهٔ کلی‌اند: سیستم LTI با مشتق‌گیری، تفاضل‌گیری، انتگرال‌گیری و انتقال جابه‌جا می‌شود. اگر x→yx \to y، آنگاه

x′(t)→y′(t),x[n]−x[n−1]→y[n]−y[n−1],∫−∞tx→∫−∞ty.x'(t) \to y'(t), \qquad x[n]-x[n-1] \to y[n]-y[n-1], \qquad \int_{-\infty}^{t} x \to \int_{-\infty}^{t} y .

بنابراین پاسخ معلوم به هر ورودی را می‌توان بدون آنکه اصلاً hh را پیدا کنیم، به پاسخ به مشتق آن، تفاضل آن، یا هر ترکیبی از نسخه‌های انتقال‌یافته‌اش تبدیل کرد.

مقدار نهایی پاسخ پله همان بهرهٔ DC است: s(∞)=∫h=H(0)s(\infty) = \int h = H(0) (در حالت گسسته: ∑nh[n]\sum_n h[n]).

۲.۴ علّی‌بودن و پایداری

هر دو مستقیماً از روی hh خوانده می‌شوند:

  • علّی   ⟺  h(t)=0\iff h(t)=0 برای t<0t<0 (خروجی نمی‌تواند پیش از ورودی ظاهر شود).
  • پایدار BIBO   ⟺  ∫−∞∞∣h(t)∣ dt<∞\iff \displaystyle\int_{-\infty}^{\infty}\lvert h(t)\rvert\,dt < \infty (مطلقاً انتگرال‌پذیر).

۲.۵ انرژی خروجی

انرژی یک سیگنال عبارت است از

Ey=∫−∞∞∣y(t)∣2 dt,Ey=∑n=−∞∞∣y[n]∣2.E_y = \int_{-\infty}^{\infty} \lvert y(t)\rvert^{2}\,dt, \qquad E_y = \sum_{n=-\infty}^{\infty} \lvert y[n]\rvert^{2}.

انرژی خطی نیست: انرژی y=x∗hy = x*h برابر حاصل‌ضرب یا مجموع انرژی‌های xx و hh نیست. ابتدا yy را بیابید، سپس مجذور کنید و جمع بزنید. دو نتیجه که ارزش به‌خاطرسپردن دارند، برای a>0a>0 و ∣r∣<1\lvert r\rvert<1:

∫0∞(e−at)2 dt=12a,∑n=0∞∣rn∣2=11−∣r∣2.\int_{0}^{\infty} \big(e^{-at}\big)^{2}\,dt = \frac{1}{2a}, \qquad \sum_{n=0}^{\infty} \big\lvert r^{n}\big\rvert^{2} = \frac{1}{1-\lvert r\rvert^{2}} .

۳. روش — محاسبهٔ یک کانولوشن

برای محاسبهٔ دستی y(t)=(x∗h)(t)y(t) = (x*h)(t):

  1. متغیر ظاهری τ\tau را انتخاب کنید و x(τ)x(\tau) و h(t−τ)h(t-\tau) را بنویسید. معمولاً وارون‌کردن و انتقال‌دادن سیگنال ساده‌تر / کوتاه‌تر آسان‌تر است (جابه‌جایی‌پذیری این انتخاب را به شما می‌دهد).
  2. یکی از سیگنال‌ها را قرینه کنید: h(τ)→h(−τ)h(\tau) \to h(-\tau). آن را به اندازهٔ tt انتقال دهید: h(−τ)→h(t−τ)h(-\tau) \to h(t-\tau).
  3. tt را از −∞-\infty تا +∞+\infty بلغزانید؛ در هر tt از حاصل‌ضرب x(τ)h(t−τ)x(\tau)h(t-\tau) فقط روی ناحیهٔ هم‌پوشانی انتگرال بگیرید.
  4. نقاط شکست را دنبال کنید — مقادیری از tt که در آن‌ها لبه‌های جلویی/عقبی دو تکیه‌گاه به هم می‌رسند. این نقاط پاسخ را به بازه‌هایی تقسیم می‌کنند؛ انتگرال هر بازه را جداگانه حل کنید.
  5. آزمون‌های درستی: تکیه‌گاه x∗hx*h از (شروع xx + شروع hh) تا (پایان xx + پایان hh) امتداد دارد؛ مدت آن مجموع دو مدت است. مساحت‌ها در هم ضرب می‌شوند: ∫(x∗h) dt=(∫x)(∫h)\int (x*h)\,dt = \big(\int x\big)\big(\int h\big).

میان‌بر — ورودی‌ها/هسته‌های ضربه‌ای. کانولوشن با ∑ici δ(t−ti)\sum_i c_i\,\delta(t-t_i) صرفاً ∑ici x(t−ti)\sum_i c_i\,x(t-t_i) است — هر ضربه را با یک نسخهٔ مقیاس‌شده و انتقال‌یافته از سیگنال دیگر جایگزین کنید. ورودی گسسته‌ای که به‌صورت نمودار میله‌ای رسم شده همیشه به همین شکل است: آن را به‌صورت x[n]=∑kx[k] δ[n−k]x[n]=\sum_k x[k]\,\delta[n-k] بخوانید و بنویسید y[n]=∑kx[k] h[n−k]y[n]=\sum_k x[k]\,h[n-k]. سپس نمونه به نمونه محاسبه کنید؛ نسخه‌های انتقال‌یافته اغلب از اندیسی به بعد یکدیگر را حذف می‌کنند و خروجی‌ای با طول متناهی باقی می‌ماند.

میان‌بر — ورودی‌های پله‌مانند (نیمه‌نامتناهی). وقتی یک سیگنال پلهٔ واحد است، پله‌های واحد فقط حدود انتگرال‌گیری را تعیین می‌کنند:

ورودی x(t)x(t) خروجی y(t)y(t) برحسب پاسخ پله
u(t)u(t) ∫−∞th(τ) dτ\int_{-\infty}^{t} h(\tau)\,d\tau s(t)s(t)
u(t−t0)u(t-t_0) ∫−∞t−t0h(τ) dτ\int_{-\infty}^{t-t_0} h(\tau)\,d\tau s(t−t0)s(t-t_0)
u(−t)u(-t) ∫t∞h(τ) dτ\int_{t}^{\infty} h(\tau)\,d\tau H(0)−s(t)H(0)-s(t)
u(t−a)−u(t−b)u(t-a)-u(t-b) ∫t−bt−ah(τ) dτ\int_{t-b}^{t-a} h(\tau)\,d\tau s(t−a)−s(t−b)s(t-a)-s(t-b)

سطر u(−t)u(-t) از u(−t)=1−u(t)u(-t) = 1 - u(t) نتیجه می‌شود: ثابت 11 مقدار H(0)H(0) را تولید می‌کند و پله s(t)s(t) را (این کار نیازمند همگرایی ∫h\int h است). برای یک سیستم علّی داریم s(t)=0s(t)=0 برای t<0t<0، پس خروجی به ازای هر t<0t<0 برابر ثابت H(0)H(0) است و تنها از t=0t=0 شروع به تغییر می‌کند: ورودی از ازل یک سطح DC بوده، سیستم در حالت ماندگار DC خود قرار دارد، و پس از برداشته‌شدن ورودی، خروجی آزادانه میرا می‌شود.

دو زوج کانولوشن با پله که پیوسته ظاهر می‌شوند (برای a>0a>0، r≠1r\neq 1):

e−atu(t)∗u(t)=1−e−ata u(t),rnu[n]∗u[n]=1−rn+11−r u[n].e^{-at}u(t) * u(t) = \frac{1-e^{-at}}{a}\,u(t), \qquad r^{n}u[n] * u[n] = \frac{1-r^{n+1}}{1-r}\,u[n].

۴. مثال‌های کوچک حل‌شده

۴.۱ دو پالس متناهی

دو پالس مستطیلی با ارتفاع واحد را کانوالو کنید. فرض کنید x(t)=h(t)=1x(t)=h(t)=1 روی [0,1][0,1] و 00 در جاهای دیگر. آنگاه y=x∗hy=x*h دارای تکیه‌گاه [0,2][0,2] است و

y(t)=∫−∞∞x(τ) h(t−τ) dτ=(length of the overlap of [0,1] with [t−1,t]).y(t) = \int_{-\infty}^{\infty} x(\tau)\,h(t-\tau)\,d\tau = \text{(length of the overlap of } [0,1] \text{ with } [t-1,t]).
  • برای 0≤t≤10 \le t \le 1: پالس‌ها روی [0,t][0,t] هم‌پوشانی دارند، پس y(t)=ty(t)=t.
  • برای 1≤t≤21 \le t \le 2: روی [t−1,1][t-1,1] هم‌پوشانی دارند، پس y(t)=1−(t−1)=2−ty(t) = 1-(t-1) = 2-t.
  • در غیر این صورت y(t)=0y(t)=0.

نتیجه پالس مثلثی y(t)=max⁡(0, 1−∣t−1∣)y(t)=\max(0,\,1-\lvert t-1\rvert) است — با پهنای 22 (مجموع دو پهنا)، قلهٔ 11 در t=1t=1، و مساحت کل 1=(1)(1)1 = (1)(1) که با آزمون حاصل‌ضرب مساحت‌ها سازگار است.

۴.۲ پلهٔ چپ‌سو به یک سیستم علّی

فرض کنید h(t)=3e−2tu(t)h(t) = 3e^{-2t}u(t) و x(t)=u(−t)x(t) = u(-t). ورودی را وارون کنید: x(t−τ)=u(τ−t)x(t-\tau) = u(\tau - t)، پس

y(t)=∫0∞3e−2τ u(τ−t) dτ=∫max⁡(0, t)∞3e−2τ dτ=32 e−2max⁡(0, t).y(t) = \int_{0}^{\infty} 3e^{-2\tau}\,u(\tau-t)\,d\tau = \int_{\max(0,\,t)}^{\infty} 3e^{-2\tau}\,d\tau = \tfrac{3}{2}\,e^{-2\max(0,\,t)} . y(t)={32,t≤032 e−2t,t>0y(t) = \begin{cases} \tfrac{3}{2}, & t \le 0 \\ \tfrac{3}{2}\,e^{-2t}, & t > 0 \end{cases}

با جدول بررسی کنید: H(0)=∫0∞3e−2tdt=32H(0) = \int_0^\infty 3e^{-2t}dt = \tfrac32 و s(t)=32(1−e−2t)u(t)s(t) = \tfrac32\big(1-e^{-2t}\big)u(t)، پس H(0)−s(t)H(0)-s(t) دقیقاً همان دو شاخهٔ بالا را می‌دهد. برای زمان‌های منفی خروجی تخت است، پس برای نمونه y(−1)−y(−4)=0y(-1)-y(-4)=0، در حالی که y(1)−y(−1)=32(e−2−1)y(1)-y(-1) = \tfrac32\big(e^{-2}-1\big).

۴.۳ از پاسخ پله تا پاسخ به یک ورودی نمودار میله‌ای

یک سیستم LTI گسسته دارای پاسخ پلهٔ s[n]=(3−2(13)n)u[n]s[n] = \big(3 - 2(\tfrac13)^{n}\big)u[n] است. پاسخ آن به x[n]=3δ[n]−δ[n−1]x[n] = 3\delta[n] - \delta[n-1] و انرژی این پاسخ را بیابید.

پاسخ ضربه. h[n]=s[n]−s[n−1]h[n] = s[n]-s[n-1]، که نمونهٔ اول را باید جداگانه بررسی کرد زیرا s[−1]=0s[-1]=0:

  • n=0n=0: h[0]=s[0]−0=3−2=1h[0] = s[0] - 0 = 3-2 = 1.
  • n≥1n\ge 1: h[n]=2(13)n−1−2(13)n=4(13)nh[n] = 2(\tfrac13)^{n-1} - 2(\tfrac13)^{n} = 4(\tfrac13)^{n}.

فرم بستهٔ 4(13)n4(\tfrac13)^n در n=0n=0 مقدار 44 می‌دهد که نادرست است — پس h[n]=4(13)nu[n]−3δ[n]h[n] = 4(\tfrac13)^{n}u[n] - 3\delta[n]. بررسی: ∑nh[n]=1+4⋅1/31−1/3=3=s(∞)\sum_n h[n] = 1 + 4\cdot\tfrac{1/3}{1-1/3} = 3 = s(\infty).

خروجی. y[n]=3h[n]−h[n−1]y[n] = 3h[n] - h[n-1]:

  • y[0]=3⋅1=3y[0] = 3\cdot 1 = 3.
  • y[1]=3⋅43−h[0]=4−1=3y[1] = 3\cdot\tfrac43 - h[0] = 4 - 1 = 3.
  • n≥2n\ge 2: 3⋅4(13)n−4(13)n−1=03\cdot 4(\tfrac13)^{n} - 4(\tfrac13)^{n-1} = 0.

پس y[n]=3δ[n]+3δ[n−1]y[n] = 3\delta[n] + 3\delta[n-1] و Ey=32+32=18E_y = 3^{2}+3^{2} = 18.

بررسی با تبدیل. H(z)=1+z−11−13z−1H(z) = \dfrac{1+z^{-1}}{1-\tfrac13 z^{-1}} و X(z)=3−z−1=3(1−13z−1)X(z) = 3 - z^{-1} = 3\big(1-\tfrac13 z^{-1}\big)، پس Y(z)=3(1+z−1)Y(z) = 3\big(1+z^{-1}\big): صفر ورودی در z=13z=\tfrac13 قطب سیستم را حذف می‌کند و به همین دلیل دنبالهٔ هندسی نامتناهی ناپدید می‌شود.

۵. گونه‌ها و حالت‌های خاص

  • تلهٔ نمایی همیشگی. اگر ورودی este^{st} بدون u(t)u(t) باشد، سراغ معادلهٔ دیفرانسیل نروید: خروجی صرفاً H(s) estH(s)\,e^{st} است — نه پلهٔ روشن‌شدنی در کار است و نه مدهای طبیعی (ورودی‌صفر). ضریب u(t)u(t) نشانه‌ای است که باید حالت گذرا و مدهای طبیعی را نیز در نظر گرفت.
  • ورودی‌ای که به‌جای وصل‌شدن قطع می‌شود. سطح ثابت cc که از t=−∞t=-\infty اعمال شده و در t=t0t=t_0 برداشته می‌شود، یعنی x(t)=c u(t0−t)x(t)=c\,u(t_0-t)، خروجی y(t)=c [H(0)−s(t−t0)]y(t) = c\,\big[H(0) - s(t-t_0)\big] را می‌دهد. برای یک سیستم علّی و پایدار، این خروجی پیش از t0t_0 برابر c H(0)c\,H(0) و پس از آن یک میرایی آزاد است. اگر سیستم علّی نباشد، خروجی پیش از t0t_0 شروع به تغییر می‌کند، و اگر ∫h\int h واگرا باشد (مثلاً انتگرال‌گیر، h=uh=u) خروجی وجود ندارد.
  • حذف قطب و صفر. در H(s)H(s) ممکن است یک صفر یک قطب را حذف کند (مثلاً s+2s(s+2)=1s\tfrac{s+2}{s(s+2)}=\tfrac{1}{s}) و مرتبهٔ مؤثر را کاهش دهد — همیشه پیش از محاسبهٔ H(s)H(s) آن را ساده کنید. همین اتفاق میان یک ورودی و سیستم نیز رخ می‌دهد: اگر X(z)X(z) (یا X(s)X(s)) در محل یکی از قطب‌های HH صفر داشته باشد، مد متناظر در خروجی حضور ندارد و یک پاسخ با طول نامتناهی می‌تواند به چند نمونه فروبکاهد.
  • سایر زوج‌های ورودی–خروجی معلوم. اگر پاسخ y1y_1 به ورودی‌ای مانند x1x_1 داده شده باشد و ورودی جدید ترکیبی از انتقال‌ها، مشتق‌ها یا تفاضل‌های x1x_1 باشد، همان ترکیب را بر y1y_1 اعمال کنید. گذر از hh تنها زمانی لازم است که ورودی جدید ارتباطی با x1x_1 نداشته باشد. در زمان پیوسته، مشتق‌گیری از پاسخ پله‌ای که در t=0t=0 پرش دارد یک ضربه تولید می‌کند: ddt[g(t)u(t)]=g′(t)u(t)+g(0) δ(t)\tfrac{d}{dt}\big[g(t)u(t)\big] = g'(t)u(t) + g(0)\,\delta(t).
  • همبستگی در برابر کانولوشن. همبستگی متقابل همان کانولوشن است با یکی از سیگنال‌ها به‌صورت قرینه‌شده: Rxy(τ)=∫x(t) y(t−τ) dt=(x(τ)∗y(−τ))R_{xy}(\tau)=\int x(t)\,y(t-\tau)\,dt = \big(x(\tau)*y(-\tau)\big). برای یک سیگنال انرژی حقیقی، خودهمبستگی Rv(τ)=v(τ)∗v(−τ)R_v(\tau)=v(\tau)*v(-\tau) زوج است و در τ=0\tau=0 بیشینه می‌شود، جایی که برابر انرژی است: Rv(0)=EvR_v(0)=E_v. برای مجموعی از ضربه‌ها v=∑kδ(t−tk)v=\sum_k \delta(t-t_k)، به ازای هر تفاضل تأخیر دوبه‌دو ti−tjt_i-t_j یک ضربهٔ خروجی قرار می‌گیرد، و جملهٔ τ=0\tau=0 وزنی برابر تعداد ضربه‌ها انباشته می‌کند.
  • ضرب‌کننده ≠ LTI. ضرب x(t)x(t) در یک سیگنال ثابت و معین g(t)g(t) (حامل، پنجره یا قطار ضربه) سیستم را خطی نگه می‌دارد، اما تنها زمانی تغییرناپذیر با زمان است که gg ثابت باشد — ضرب‌کنندهٔ وابسته به زمان، سیستم خطی را متغیر با زمان می‌کند، و چنین سیستمی توصیف یکتایی با پاسخ ضربه / کانولوشن ندارد.
  • زمان گسسته همه‌چیز را با مجموع‌ها بازتاب می‌دهد؛ zn→H(z) znz^n \to H(z)\,z^n و شرط پایداری ∑n∣h[n]∣<∞\sum_n \lvert h[n]\rvert < \infty می‌شود.

تلهٔ رایج: هنگام تفاضل‌گیری از پاسخ پله، نمونهٔ اول را دستی حساب کنید. فرم بسته‌ای که برای n≥1n\ge 1 به دست می‌آید لزوماً در n=0n=0 برقرار نیست، زیرا s[−1]=0s[-1]=0 از فرمول درون u[n]u[n] به دست نمی‌آید. به همین ترتیب، عبارتی برای y(t)y(t) که به ازای t>0t>0 به دست آمده (مانند یک نمایی میرا) نباید در زمان‌های منفی محاسبه شود — هر بازه از یک کانولوشن فرمول خودش را دارد.

نکتهٔ کنکوری: ابزار را متناسب با ورودی انتخاب کنید. برای نمایی همیشگی este^{st} کانولوشن نگیرید — تابع تبدیل را یک‌بار محاسبه کنید: y=H(s)esty=H(s)e^{st}؛ ورودی ثابت همان حالت s=0s=0 است. برای پله یا پله‌ای که قطع می‌شود، از hh بین حدودی که پله‌ها تحمیل می‌کنند انتگرال بگیرید، یا از s(t)s(t) و H(0)−s(t)H(0)-s(t) استفاده کنید. برای پالس‌های با تکیه‌گاه متناهی، وارون کنید و بلغزانید و نقاط شکست را دنبال کنید. برای مجموعی از ضربه‌ها (هسته، سیگنال یا هر نمودار میله‌ای)، کانولوشن/همبستگی به نسخه‌های انتقال‌یافته و مقیاس‌شده فرومی‌کاهد — نیازی به انتگرال نیست، و پیش از جمع‌زدن هر سری نامتناهی بررسی کنید که آیا نسخه‌ها یکدیگر را حذف می‌کنند یا نه. بودن یا نبودن u(t)u(t) به شما می‌گوید که اصلاً مدهای طبیعی ظاهر می‌شوند یا خیر.

فرمول‌های کلیدی

Convolution (continuous time)
y(t)=(x∗h)(t)=∫−∞∞x(τ) h(t−τ) dτy(t) = (x * h)(t) = \int_{-\infty}^{\infty} x(\tau)\,h(t-\tau)\,d\tau
Convolution (discrete time)
y[n]=(x∗h)[n]=∑k=−∞∞x[k] h[n−k]y[n] = (x * h)[n] = \sum_{k=-\infty}^{\infty} x[k]\,h[n-k]
Sifting property of the impulse
∫−∞∞x(τ) δ(t−τ) dτ=x(t)\int_{-\infty}^{\infty} x(\tau)\,\delta(t-\tau)\,d\tau = x(t)
Eigenfunction property
est  ⟶  H(s) est,H(s)=∫−∞∞h(t) e−st dte^{st} \;\longrightarrow\; H(s)\,e^{st}, \qquad H(s) = \int_{-\infty}^{\infty} h(t)\,e^{-st}\,dt
DC gain (response to a constant input)
x(t)=c  ⟶  y(t)=c H(0),H(0)=∫−∞∞h(t) dtx(t)=c \;\longrightarrow\; y(t)=c\,H(0), \qquad H(0)=\int_{-\infty}^{\infty} h(t)\,dt
Step response (continuous time)
s(t)=∫−∞th(τ) dτ,h(t)=ds(t)dts(t) = \int_{-\infty}^{t} h(\tau)\,d\tau, \qquad h(t) = \frac{ds(t)}{dt}
Step response (discrete time)
s[n]=∑k=−∞nh[k],h[n]=s[n]−s[n−1]s[n] = \sum_{k=-\infty}^{n} h[k], \qquad h[n] = s[n]-s[n-1]
Response to a left-sided step
u(−t)  ⟶  ∫t∞h(τ) dτ=H(0)−s(t)u(-t) \;\longrightarrow\; \int_{t}^{\infty} h(\tau)\,d\tau = H(0) - s(t)
Exponential convolved with a step
e−atu(t)∗u(t)=1−e−ata u(t),rnu[n]∗u[n]=1−rn+11−r u[n]e^{-at}u(t) * u(t) = \frac{1-e^{-at}}{a}\,u(t), \qquad r^{n}u[n] * u[n] = \frac{1-r^{n+1}}{1-r}\,u[n]
BIBO stability
∫−∞∞∣h(t)∣ dt<∞\int_{-\infty}^{\infty} |h(t)|\,dt < \infty
Signal energy
Ey=∫−∞∞∣y(t)∣2 dt,Ey=∑n=−∞∞∣y[n]∣2E_y = \int_{-\infty}^{\infty} |y(t)|^{2}\,dt, \qquad E_y = \sum_{n=-\infty}^{\infty} |y[n]|^{2}
Deterministic autocorrelation (energy signal)
Rv(τ)=∫−∞∞v(t) v(t−τ) dt=(v(τ)∗v(−τ))R_v(\tau) = \int_{-\infty}^{\infty} v(t)\,v(t-\tau)\,dt = \big(v(\tau) * v(-\tau)\big)
#lti-systems#impulse-response#step-response#convolution#eigenfunction#system-function#correlation#causality-stability

سوالات تمرینی · 8/16

همهٔ 16 سؤال →