KonkurStudy
همه کارت‌ها

فوریهٔ زمان‌گسسته: DTFS، DTFT و DFT

سه ابزار فوریهٔ زمان‌گسسته — DTFS برای دنباله‌های متناوب، DTFT برای دنباله‌های نامتناوب و DFT برای بلوک‌های با طول محدود — همراه با زوج‌های تبدیل (از جمله هستهٔ دیریکله)، خواص (انتقال، وارون‌سازی زمانی، مزدوج‌گیری، مقیاس‌دهی، مشتق‌گیری در حوزهٔ فرکانس، تقارن، پارسوال)، تقارن‌های ضرایب DTFS، ارتباط قطار ضربه با سری فوریهٔ زمان‌پیوسته، و روش‌هایی برای وارون‌کردن طیف‌های مرکب، محاسبهٔ انتگرال‌های طیفی به‌صورت نمونه و دسته‌بندی سیستم‌ها از روی یک رابطهٔ حوزهٔ فرکانس.

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

اهمیت
۰/۵

۱. ایدهٔ اصلی

یک سیگنال زمان‌گسسته x[n]x[n] را می‌توان به‌صورت مجموعی از نمایی‌های مختلط ejΩne^{j\Omega n} نوشت. اینکه از کدام ابزار فوریه استفاده کنید فقط به دامنهٔ x[n]x[n] بستگی دارد:

سیگنال متغیر فرکانس ابزار
متناوب، با دورهٔ تناوب NN گسسته، NN هارمونیک DTFS — ضرایب aka_k
نامتناوب، نامحدود Ω\Omega پیوسته، متناوب با دورهٔ 2π2\pi DTFT — X(ejΩ)X(e^{j\Omega})
طول محدود NN NN نمونه از DTFT DFT — X[k]X[k]

تنها واقعیتی که همه‌چیز را سامان می‌دهد: چون nn عددی صحیح است، ejΩn=ej(Ω+2π)ne^{j\Omega n}=e^{j(\Omega+2\pi)n}، پس هر طیف زمان‌گسسته نسبت به Ω\Omega متناوب با دورهٔ 2π2\pi است. DTFS روی NN نقطه از یک دوره قرار دارد؛ DTFT یک منحنی پیوسته روی یک دوره است؛ DFT همان منحنی را در NN نقطهٔ هم‌فاصله نمونه‌برداری می‌کند. این‌ها یک نظریهٔ واحدند که با سه دقت تفکیک مختلف دیده می‌شود.

واقعیت سامان‌دهندهٔ دوم به همان اندازهٔ اولی ارزش دارد: مجموع تحلیل یک سری توانی برحسب e−jΩe^{-j\Omega} است که ضرایب آن همان نمونه‌ها هستند. هرگاه بتوان طیفی را به‌صورت ∑ncne−jΩn\sum_n c_n e^{-j\Omega n} نوشت، دنباله به‌سادگی x[n]=cnx[n]=c_n است — هیچ انتگرالی لازم نیست. بیشتر سؤال‌های «سخت» DTFT به‌محض بازنویسی طیف در این قالب فرو می‌ریزند.

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

۲.۱ سه زوج تبدیل

DTFT (نامتناوب). تحلیل و سنتز:

X(ejΩ)=∑n=−∞∞x[n] e−jΩn,x[n]=12π∫2πX(ejΩ) ejΩn dΩ.X(e^{j\Omega}) = \sum_{n=-\infty}^{\infty} x[n]\,e^{-j\Omega n}, \qquad x[n] = \frac{1}{2\pi}\int_{2\pi} X(e^{j\Omega})\,e^{j\Omega n}\,d\Omega.

انتگرال سنتز روی هر بازه‌ای به طول 2π2\pi گرفته می‌شود. با قراردادن n=0n=0 اتحاد مفید x[0]=12π∫2πX(ejΩ) dΩx[0]=\tfrac{1}{2\pi}\int_{2\pi}X(e^{j\Omega})\,d\Omega به دست می‌آید، و نیز X(ej0)=∑nx[n]X(e^{j0})=\sum_n x[n].

DTFS (متناوب، با دورهٔ تناوب NN). با Ω0=2π/N\Omega_0 = 2\pi/N:

x[n]=∑k=⟨N⟩ak e jkΩ0n,ak=1N∑n=⟨N⟩x[n] e−jkΩ0n.x[n] = \sum_{k=\langle N\rangle} a_k\, e^{\,jk\Omega_0 n}, \qquad a_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n]\,e^{-jk\Omega_0 n}.

ضرایب خودشان نسبت به kk متناوب با دورهٔ NN هستند: ak+N=aka_{k+N}=a_k. فقط NN تا از آن‌ها متمایزند. دو مقدار رایگان به دست می‌آید: a0=1N∑n=⟨N⟩x[n]a_0=\tfrac1N\sum_{n=\langle N\rangle}x[n] (میانگین) و ∑k=⟨N⟩ak=x[0]\sum_{k=\langle N\rangle}a_k=x[0].

DFT (طول محدود NN). X[k]X[k] همان DTFT نمونه‌برداری‌شده در Ωk=2πk/N\Omega_k = 2\pi k/N است (معادلاً N⋅akN\cdot a_k مربوط به گسترش متناوب):

X[k]=∑n=0N−1x[n] e−j2πkn/N,x[n]=1N∑k=0N−1X[k] e j2πkn/N.X[k] = \sum_{n=0}^{N-1} x[n]\,e^{-j 2\pi kn/N}, \qquad x[n] = \frac{1}{N}\sum_{k=0}^{N-1} X[k]\,e^{\,j 2\pi kn/N}.

۲.۲ زوج‌های DTFT که باید در نگاه اول شناخت

دنبالهٔ x[n]x[n] DTFT یعنی X(ejΩ)X(e^{j\Omega})
δ[n−n0]\delta[n-n_0] e−jΩn0e^{-j\Omega n_0}
anu[n]a^{n}u[n], ∣a∣<1\lvert a\rvert<1 11−a e−jΩ\dfrac{1}{1-a\,e^{-j\Omega}}
11 برای ∣n∣≤M\lvert n\rvert\le M، وگرنه 00 (2M+12M+1 عدد یک، متقارن حول مبدأ) sin⁡((2M+1)Ω/2)sin⁡(Ω/2)=1+2∑m=1Mcos⁡(mΩ)\dfrac{\sin\left((2M+1)\Omega/2\right)}{\sin(\Omega/2)} = 1+2\sum_{m=1}^{M}\cos(m\Omega)
11 برای 0≤n≤L−10\le n\le L-1، وگرنه 00 (LL عدد یک، علّی) e−jΩ(L−1)/2 sin⁡(LΩ/2)sin⁡(Ω/2)e^{-j\Omega(L-1)/2}\,\dfrac{\sin(L\Omega/2)}{\sin(\Omega/2)}
L−∣n∣L-\lvert n\rvert برای ∣n∣≤L−1\lvert n\rvert\le L-1، وگرنه 00 (مثلثی) sin⁡2(LΩ/2)sin⁡2(Ω/2)\dfrac{\sin^{2}(L\Omega/2)}{\sin^{2}(\Omega/2)}
sin⁡(Wn)πn\dfrac{\sin(Wn)}{\pi n}, 0<W<π0<W<\pi 11 برای ∣Ω∣≤W\lvert\Omega\rvert\le W، و 00 برای W<∣Ω∣≤πW<\lvert\Omega\rvert\le\pi (پایین‌گذر ایده‌آل)
11 (به‌ازای همهٔ nnها) 2π∑lδ(Ω−2πl)2\pi\sum_{l}\delta(\Omega-2\pi l)
ejΩ0ne^{j\Omega_0 n} 2π∑lδ(Ω−Ω0−2πl)2\pi\sum_{l}\delta(\Omega-\Omega_0-2\pi l)

نسبت sin⁡(LΩ/2)/sin⁡(Ω/2)\sin(L\Omega/2)/\sin(\Omega/2) همان هستهٔ دیریکله (سینک متناوب) است — همتای زمان‌گسستهٔ تابع سینک. این نسبت، تبدیلِ رشته‌ای از LL عدد یک است؛ مقدار آن در Ω=0\Omega=0 برابر LL است. سطر مثلثی از همین نتیجه می‌شود: مثلث L−∣n∣L-\lvert n\rvert کانولوشن یک پالس به طول LL با وارونِ زمانی خودش است، پس تبدیل آن مجذور اندازهٔ هسته است.

۲.۳ خواص DTFT

هر سطر یک عمل زمانی را به یک عمل طیفی نگاشت می‌کند.

حوزهٔ زمان حوزهٔ فرکانس
x[n−n0]x[n-n_0] e−jΩn0 X(ejΩ)e^{-j\Omega n_0}\,X(e^{j\Omega}) (تأخیر خالص)
ejΩ0nx[n]e^{j\Omega_0 n}x[n] X(ej(Ω−Ω0))X(e^{j(\Omega-\Omega_0)}) (مدولاسیون)
x[−n]x[-n] X(e−jΩ)X(e^{-j\Omega}) (وارون‌سازی زمانی)
x∗[n]x^{*}[n] X∗(e−jΩ)X^{*}(e^{-j\Omega})
x∗[−n]x^{*}[-n] X∗(ejΩ)X^{*}(e^{j\Omega}) (مزدوج‌گیری از طیف = مزدوج‌گیری و وارون‌سازی در زمان)
n x[n]n\,x[n] j dX(ejΩ)dΩ\displaystyle j\,\frac{dX(e^{j\Omega})}{d\Omega} (مشتق‌گیری در حوزهٔ فرکانس)
x[n/m]x[n/m] (درج صفر، بیش‌نمونه‌برداری با ضریب mm) X(ejmΩ)X(e^{jm\Omega}) (مقیاس‌دهی/فشرده‌سازی فرکانسی)
x[n]∗h[n]x[n]*h[n] X(ejΩ) H(ejΩ)X(e^{j\Omega})\,H(e^{j\Omega})
x[n] g[n]x[n]\,g[n] 12π\tfrac{1}{2\pi} برابرِ کانولوشن متناوب X,GX,G
x[n]∗x∗[−n]x[n]*x^{*}[-n] (خودهمبستگی) ∣X(ejΩ)∣2\lvert X(e^{j\Omega})\rvert^{2}

تقارن سیگنال حقیقی. اگر x[n]x[n] حقیقی باشد، X(ejΩ)=X∗(e−jΩ)X(e^{j\Omega}) = X^{*}(e^{-j\Omega}) (تقارن مزدوج): ∣X∣|X| زوج و ∠X\angle X فرد است. برای یک سیگنال متناوب حقیقی، معادل آن در DTFS عبارت است از ak∗=a−k=aN−ka_k^{*}=a_{-k}=a_{N-k}.

قسمت‌های حقیقی و موهومی یک طیف. چون X∗(ejΩ)↔x∗[−n]X^{*}(e^{j\Omega})\leftrightarrow x^{*}[-n]، گرفتن قسمت حقیقی یا موهومی یک طیف، دنباله را به قسمت‌های مزدوج‌متقارن و مزدوج‌پادمتقارنش تجزیه می‌کند:

Re[X(ejΩ)]  ⟷  x[n]+x∗[−n]2,j Im[X(ejΩ)]  ⟷  x[n]−x∗[−n]2.\mathrm{Re}\left[X(e^{j\Omega})\right]\;\longleftrightarrow\;\frac{x[n]+x^{*}[-n]}{2}, \qquad j\,\mathrm{Im}\left[X(e^{j\Omega})\right]\;\longleftrightarrow\;\frac{x[n]-x^{*}[-n]}{2}.

به‌طور دوگان، Re{x[n]}↔12[X(ejΩ)+X∗(e−jΩ)]\mathrm{Re}\{x[n]\}\leftrightarrow\tfrac12\left[X(e^{j\Omega})+X^{*}(e^{-j\Omega})\right]. برای x[n]x[n] حقیقی این روابط به «قسمت زوج ↔\leftrightarrow قسمت حقیقی» و «قسمت فرد ↔\leftrightarrow j ×j\,\times قسمت موهومی» ساده می‌شوند.

تغییر آفینِ آرگومان فرکانس. برای هر عدد صحیح ناصفر cc (منفی هم مجاز است) و هر عدد صحیح n0n_0:

G(ejΩ)=e−jΩn0 X(ejcΩ)⟺g[n]={x[n−n0c],(n−n0) a multiple of c0,otherwiseG(e^{j\Omega}) = e^{-j\Omega n_0}\,X(e^{jc\Omega}) \quad\Longleftrightarrow\quad g[n] = \begin{cases} x\left[\dfrac{n-n_0}{c}\right], & (n-n_0) \text{ a multiple of } c \\ 0, & \text{otherwise} \end{cases}

همین یک خط، تأخیر (c=1c=1)، وارون‌سازی (c=−1c=-1)، انبساط (c=m>0c=m>0) و هر ترکیبی از آن‌ها را در بر دارد. اثبات آن با جایگذاری مجموع تحلیل انجام می‌شود: e−jΩn0∑kx[k]e−jΩcke^{-j\Omega n_0}\sum_k x[k]e^{-j\Omega ck} نمونهٔ x[k]x[k] را در مکان n=ck+n0n=ck+n_0 قرار می‌دهد.

پارسوال (انرژی/توان).

∑n∣x[n]∣2=12π∫2π∣X(ejΩ)∣2 dΩ(DTFT),1N∑n=⟨N⟩∣x[n]∣2=∑k=⟨N⟩∣ak∣2(DTFS).\sum_{n} |x[n]|^{2} = \frac{1}{2\pi}\int_{2\pi} |X(e^{j\Omega})|^{2}\,d\Omega \quad(\text{DTFT}),\qquad \frac{1}{N}\sum_{n=\langle N\rangle}|x[n]|^{2} = \sum_{k=\langle N\rangle}|a_k|^{2}\quad(\text{DTFS}).

صورت دوسیگنالی (تعمیم‌یافته) عبارت است از ∑nx[n] y∗[n]=12π∫2πX(ejΩ) Y∗(ejΩ) dΩ\sum_n x[n]\,y^{*}[n]=\tfrac{1}{2\pi}\int_{2\pi}X(e^{j\Omega})\,Y^{*}(e^{j\Omega})\,d\Omega.

انتگرال یک طیف در برابر یک هارمونیک. معادلهٔ سنتز را وارونه بخوانید — انتگرال روی یک دوره، یک نمونه است:

∫2πX(ejΩ) ejΩk dΩ=2π x[k],∫2πX(ejΩ)cos⁡(kΩ) dΩ=π(x[k]+x[−k]),\int_{2\pi}X(e^{j\Omega})\,e^{j\Omega k}\,d\Omega = 2\pi\,x[k], \qquad \int_{2\pi}X(e^{j\Omega})\cos(k\Omega)\,d\Omega = \pi\left(x[k]+x[-k]\right), ∫2πX(ejΩ)sin⁡(kΩ) dΩ=−jπ(x[k]−x[−k]).\int_{2\pi}X(e^{j\Omega})\sin(k\Omega)\,d\Omega = -j\pi\left(x[k]-x[-k]\right).

۲.۴ خواص DTFS

فرض کنید x[n]↔akx[n]\leftrightarrow a_k و y[n]↔bky[n]\leftrightarrow b_k، هر دو با دورهٔ تناوب NN و Ω0=2π/N\Omega_0=2\pi/N.

حوزهٔ زمان ضرایب
x[n−n0]x[n-n_0] e−jkΩ0n0 ake^{-jk\Omega_0 n_0}\,a_k
ejMΩ0n x[n]e^{jM\Omega_0 n}\,x[n] ak−Ma_{k-M} (انتقال در kk)
x[−n]x[-n] a−ka_{-k}
x∗[n]x^{*}[n] a−k∗a^{*}_{-k}
x[n]−x[n−1]x[n]-x[n-1] (1−e−jkΩ0)ak\left(1-e^{-jk\Omega_0}\right)a_k
∑r=⟨N⟩x[r] y[n−r]\sum_{r=\langle N\rangle}x[r]\,y[n-r] (کانولوشن متناوب) N ak bkN\,a_k\,b_k
x[n] y[n]x[n]\,y[n] ∑l=⟨N⟩al bk−l\sum_{l=\langle N\rangle}a_l\,b_{k-l}

انتقال به اندازهٔ نیم‌دوره (NN زوج). با قراردادن n0=−N/2n_0=-N/2 در سطر انتقال داریم x[n+N/2]↔ejkπak=(−1)kakx[n+N/2]\leftrightarrow e^{jk\pi}a_k=(-1)^{k}a_k. بنابراین

x[n+N/2]=−x[n]  ⟺  ak=0 for all even k,x[n+N/2]=+x[n]  ⟺  ak=0 for all odd k.x[n+N/2]=-x[n]\;\Longleftrightarrow\;a_k=0\ \text{for all even } k, \qquad x[n+N/2]=+x[n]\;\Longleftrightarrow\;a_k=0\ \text{for all odd } k.

حالت اول پادتقارن نیم‌موج است (فقط هارمونیک‌های فرد)؛ حالت دوم صرفاً می‌گوید دورهٔ تناوب واقعی N/2N/2 است.

دوگانی تُنُکی–تناوب. فرض کنید MM مقسوم‌علیه NN باشد. طبق سطر مدولاسیون، ak−M=aka_{k-M}=a_k برای همهٔ kkها یعنی (ejMΩ0n−1)x[n]=0\left(e^{jM\Omega_0 n}-1\right)x[n]=0 برای همهٔ nnها، پس x[n]x[n] فقط در جاهایی می‌تواند ناصفر باشد که Mn/NMn/N عددی صحیح است:

ak+M=ak  ∀k  ⟺  x[n]=0 unless n is a multiple of N/M.a_{k+M}=a_k\ \ \forall k\;\Longleftrightarrow\;x[n]=0 \text{ unless } n \text{ is a multiple of } N/M .

گزارهٔ قرینه نقش‌ها را عوض می‌کند: x[n+L]=x[n]x[n+L]=x[n] با LL مقسوم‌علیه NN (دورهٔ تناوب واقعیِ کوتاه‌تر) ⟺\Longleftrightarrow ak=0a_k=0 مگر آنکه kk مضربی از N/LN/L باشد. متناوب‌بودن در یک حوزه به معنای تُنُک‌بودن در حوزهٔ دیگر است.

۲.۵ از دنبالهٔ متناوب به قطار ضربهٔ زمان‌پیوسته

نمونه‌های یک دنباله با دورهٔ تناوب NN را روی ضربه‌هایی به فاصلهٔ TsT_s از هم قرار دهید:

y(t)=∑n=−∞∞x[n] δ(t−nTs).y(t)=\sum_{n=-\infty}^{\infty}x[n]\,\delta(t-nT_s).

آنگاه y(t)y(t) متناوب است با T=NTsT=NT_s و ω0=2π/(NTs)\omega_0=2\pi/(NT_s)، و هر دوره دقیقاً NN ضربه دارد. طبق خاصیت غربالی،

ck=1T∫Ty(t) e−jkω0t dt=1NTs∑n=⟨N⟩x[n] e−jk2πNn=akTs.c_k=\frac{1}{T}\int_{T}y(t)\,e^{-jk\omega_0 t}\,dt =\frac{1}{NT_s}\sum_{n=\langle N\rangle}x[n]\,e^{-jk\frac{2\pi}{N}n} =\frac{a_k}{T_s}.

پس ضرایب سری فوریهٔ زمان‌پیوسته همان دنباله‌اند که در 1/Ts1/T_s ضرب شده است — مستقل از NN — و برخلاف ضرایب یک سیگنال زمان‌پیوستهٔ معمولی، میرا نمی‌شوند: ck+N=ckc_{k+N}=c_k. (آزمون درستی: x[n]≡1x[n]\equiv1 در مضارب NN دارای ak=1a_k=1 است که همان ضرایب آشنای قطار ضربه یعنی 1/Ts1/T_s را می‌دهد.) همتای نامتناوب آن Y(jω)=X(ejωTs)Y(j\omega)=X(e^{j\omega T_s}) است: تبدیل زمان‌پیوستهٔ قطار ضربهٔ وزن‌دار همان DTFT با Ω=ωTs\Omega=\omega T_s است.

۳. روش‌ها

۳.۱ استخراج ساختار از یک رابطهٔ طیفی

بسیاری از سؤال‌های آزمون یک رابطهٔ ورودی–خروجی در حوزهٔ فرکانس می‌دهند و می‌پرسند سیستم از چه رده‌ای است، یا ادعایی دربارهٔ تقارن/ضرایب مطرح می‌کنند. روال مطمئن:

  1. جمله‌به‌جمله به حوزهٔ زمان ترجمه کنید. هر عمل طیفی را با همتای آن در جدول خواص جایگزین کنید: e−jkΩ ⁣→e^{-j k\Omega}\!\to تأخیر به اندازهٔ kk؛ ddΩ ⁣→\tfrac{d}{d\Omega}\!\to ضرب در nn؛ X(ejmΩ) ⁣→X(e^{jm\Omega})\!\to بیش‌نمونه‌برداری x[n/m]x[n/m]؛ حاصل‌ضرب تبدیل‌ها →\to کانولوشن.
  2. رابطهٔ حوزهٔ زمان y[n]=T{x}[n]y[n]=T\{x\}[n] را صریحاً بنویسید.
  3. از روی همان رابطه دسته‌بندی کنید.
    • خطی‌بودن: آیا yy تابعکی خطی (پیرو اصل جمع آثار) از xx است؟ بهره‌های معلومِ وابسته به اندیس مانند (n−3)(n-3) خطی‌بودن را حفظ می‌کنند.
    • تغییرناپذیری با زمان: آیا انتقال ورودی همان انتقال را در خروجی ایجاد می‌کند؟ بهره‌های وابسته به اندیس و بیش/کم‌نمونه‌برداری آن را از بین می‌برند؛ تأخیر ثابت e−jkΩe^{-jk\Omega} آن را از بین نمی‌برد.
    • علّی‌بودن: آیا y[n]y[n] هرگز به x[m]x[m] با m>nm>n نیاز دارد؟ بیش‌نمونه‌برداری x[n/2]x[n/2] برای ساختن y[−2]y[-2] به x[−1]x[-1] نیاز دارد — غیرعلّی.
  4. برای ادعاهای تقارن/حقیقی‌بودن، با آزمونه‌های ساده امتحان کنید. یک انتگرال یا قید اسکالرِ تنها روی XX، یک معادله برای بی‌نهایت نمونه است — بسیار ضعیف‌تر از آنکه یک ردهٔ تقارن کامل را تحمیل کند. δ[n]\delta[n]، jδ[n]j\delta[n] و یک ضربهٔ انتقال‌یافته δ[n−k]\delta[n-k] را امتحان کنید تا «لزوماً حقیقی / زوج / فرد» نقض شود.
  5. برای مقیاس‌دهی دورهٔ تناوب در DTFS، همان سیگنال را با دورهٔ تناوب مفروضِ طولانی‌تر دوباره محاسبه کنید: ضرب دورهٔ تناوب در عدد صحیح mm، هارمونیک kk را به هارمونیک mkmk می‌نگارد و بین آن‌ها صفر می‌گذارد.

۳.۲ وارون‌کردن یک طیف مرکب با جای‌گذاری نمونه‌ها

وقتی طیفی چند عمل را روی هم انباشته است — یک ضریب فاز، آرگومان مقیاس‌شده یا منفی‌شده، مزدوج، قسمت حقیقی — زنجیره‌کردن خواص از حفظ همان جایی است که خطاهای علامت و ترتیب رخ می‌دهد. به‌جای آن:

  1. ابتدا Re/Im\mathrm{Re}/\mathrm{Im} را حذف کنید: 2 Re[G]=G+G∗2\,\mathrm{Re}[G]=G+G^{*} و 2j Im[G]=G−G∗2j\,\mathrm{Im}[G]=G-G^{*}.
  2. مجموع تحلیل را در جملهٔ درونی جایگذاری کنید و یک نمایی واحد e−jΩ(⋅)e^{-j\Omega(\cdot)} را جمع‌وجور کنید. نما به شما می‌گوید هر نمونهٔ x[k]x[k] در کدام اندیس nn می‌نشیند — این همان قاعدهٔ آفینِ §2.3 است.
  3. جملهٔ مزدوج را با G∗(ejΩ)↔g∗[−n]G^{*}(e^{j\Omega})\leftrightarrow g^{*}[-n] مدیریت کنید: در gg به‌جای nn مقدار −n-n بگذارید و مزدوج بگیرید.
  4. در اندیس‌های خواسته‌شده مقدار بگیرید، و به یاد داشته باشید هر جا اندیس عدد صحیح نباشد g[n]=0g[n]=0 است.

۳.۳ محاسبهٔ یک انتگرال مثلثاتی به‌صورت یک نمونه

انتگرالِ روی یک دوره از (نسبت سینوس‌ها) ×\times (یک کسینوس یا سینوس) تقریباً هیچ‌گاه برای انتگرال‌گیری مستقیم طرح نشده است.

  1. تبدیل را بشناسید. sin⁡(LΩ/2)/sin⁡(Ω/2)\sin(L\Omega/2)/\sin(\Omega/2) رشته‌ای از LL عدد یک است؛ مجذور آن مثلث L−∣n∣L-\lvert n\rvert است؛ حاصل‌ضرب طیف‌ها کانولوشن دنباله‌هاست.
  2. دنبالهٔ r[n]r[n] را صریحاً بنویسید — فهرست کوتاهی از اعداد.
  3. پاسخ را از §2.3 بخوانید: ∫2πRcos⁡(kΩ) dΩ=π(r[k]+r[−k])\int_{2\pi}R\cos(k\Omega)\,d\Omega=\pi\left(r[k]+r[-k]\right)، ∫2πR dΩ=2πr[0]\int_{2\pi}R\,d\Omega=2\pi r[0]، و ∫2π∣X∣2dΩ=2π∑∣x[n]∣2\int_{2\pi}\lvert X\rvert^{2}d\Omega=2\pi\sum\lvert x[n]\rvert^{2}.

۳.۴ بازسازی یک دنبالهٔ متناوب از روی قیدهای ضرایب

  1. هر شرط را به گزاره‌ای دربارهٔ نمونه‌ها یا دربارهٔ اینکه کدام aka_kها صفر می‌شوند تبدیل کنید، با استفاده از §2.4: تغییر علامت در نیم‌دوره هارمونیک‌های زوج را حذف می‌کند؛ متناوب‌بودن aka_k نسبت به kk، دنبالهٔ x[n]x[n] را تُنُک می‌کند؛ حقیقی‌بودن a−k=ak∗a_{-k}=a_k^{*} را می‌دهد.
  2. مجهول‌های باقی‌مانده را بشمارید — پس از شرط‌های تُنُکی و تقارن فقط چند نمونه باقی می‌ماند.
  3. از شرط(های) عددی استفاده کنید؛ مجموع تحلیل را فقط روی همان نمونه‌ها بنویسید و حل کنید.
  4. فقط پس از آن خواستهٔ مسئله را محاسبه کنید (یک مجموع، یک انرژی، یک نمونهٔ خاص). نمونه‌هایی که صفرند همچنان در هر مجموعی از یک تابع از x[n]x[n] شرکت دارند.

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

(الف) دسته‌بندی یک سیستم. با فرض Y(ejΩ)=e−j2ΩX(ejΩ)+dX(ejΩ)dΩ\displaystyle Y(e^{j\Omega}) = e^{-j2\Omega}X(e^{j\Omega}) + \frac{dX(e^{j\Omega})}{d\Omega}، سیستم را دسته‌بندی کنید.

ترجمه. جملهٔ اول: تأخیر خالص به اندازهٔ 22، پس x[n−2]x[n-2]. جملهٔ دوم: dXdΩ↔−j n x[n]\tfrac{dX}{d\Omega}\leftrightarrow -j\,n\,x[n] (زیرا n x[n]↔jdXdΩn\,x[n]\leftrightarrow j\tfrac{dX}{d\Omega}). در نتیجه

y[n]=x[n−2]−j n x[n].y[n] = x[n-2] - j\,n\,x[n].

دسته‌بندی. هر دو جمله نسبت به xx خطی‌اند ⇒\Rightarrow خطی. بهرهٔ −jn-jn به اندیس nn وابسته است ⇒\Rightarrow تغییرپذیر با زمان (ورودی را انتقال دهید؛ ضریب −jn-jn دیگر منطبق نمی‌شود). هر نمونهٔ به‌کاررفته، یعنی x[n−2]x[n-2] و x[n]x[n]، اندیسی ≤n\le n دارد ⇒\Rightarrow علّی. پس: خطی، تغییرپذیر با زمان، علّی. (توجه کنید که جملهٔ تأخیر به‌تنهایی LTI می‌بود؛ آنچه تغییرناپذیری را از بین برد مشتق فرکانسی است.)

(ب) طیف مرکب. با فرض Y(ejΩ)=2 Re[ej4ΩX(e−j2Ω)]Y(e^{j\Omega})=2\,\mathrm{Re}\left[e^{j4\Omega}X(e^{-j2\Omega})\right]، مقادیر y[0]y[0]، y[1]y[1] و y[2]y[2] را بیابید.

قرار دهید G=ej4ΩX(e−j2Ω)G=e^{j4\Omega}X(e^{-j2\Omega})، پس Y=G+G∗Y=G+G^{*}. با جایگذاری مجموع تحلیل،

G(ejΩ)=ej4Ω∑kx[k] ej2Ωk=∑kx[k] e−jΩ(−4−2k),G(e^{j\Omega})=e^{j4\Omega}\sum_{k}x[k]\,e^{j2\Omega k}=\sum_{k}x[k]\,e^{-j\Omega(-4-2k)},

پس x[k]x[k] در n=−4−2kn=-4-2k می‌نشیند، یعنی g[n]=x[−2−n2]g[n]=x\left[-2-\tfrac n2\right] برای nn زوج و 00 برای nn فرد (قاعدهٔ آفین با c=−2c=-2, n0=−4n_0=-4). جملهٔ مزدوج برابر است با g∗[−n]=x∗[−2+n2]g^{*}[-n]=x^{*}\left[-2+\tfrac n2\right]. بنابراین برای nn زوج،

y[n]=x[−2−n2]+x∗[−2+n2],y[n]=x\left[-2-\tfrac n2\right]+x^{*}\left[-2+\tfrac n2\right],

که نتیجه می‌دهد y[0]=x[−2]+x∗[−2]=2 Re{x[−2]}y[0]=x[-2]+x^{*}[-2]=2\,\mathrm{Re}\{x[-2]\}، y[2]=x[−3]+x∗[−1]y[2]=x[-3]+x^{*}[-1] و y[1]=0y[1]=0. توجه کنید که دو اندیس همواره نسبت به یک مرکز مشترک (در اینجا −2-2) قرینهٔ یکدیگرند — آزمونی سریع برای هر پاسخ پیشنهادی.

(ج) انتگرالی که به‌صورت نمونه خوانده می‌شود. محاسبه کنید: I=∫−ππsin⁡2(5Ω/2)sin⁡2(Ω/2)cos⁡(2Ω) dΩ\displaystyle I=\int_{-\pi}^{\pi}\frac{\sin^{2}(5\Omega/2)}{\sin^{2}(\Omega/2)}\cos(2\Omega)\,d\Omega.

عامل اول انتگرال‌ده، تبدیلِ مثلث r[n]=5−∣n∣r[n]=5-\lvert n\rvert, ∣n∣≤4\lvert n\rvert\le4 است، یعنی {1,2,3,4,5,4,3,2,1}\{1,2,3,4,5,4,3,2,1\}. بنابراین

I=π(r[2]+r[−2])=π(3+3)=6π.I=\pi\left(r[2]+r[-2]\right)=\pi(3+3)=6\pi .

بدون کسینوس، همین انتگرال برابر 2πr[0]=10π2\pi r[0]=10\pi می‌شد که همان پارسوال برای پنج عدد یک است: 2π⋅52\pi\cdot5.

(د) دنباله از روی قیدهای ضرایب. یک دنبالهٔ حقیقی دورهٔ تناوب N=12N=12 دارد، با ak+4=aka_{k+4}=a_k، x[n+6]=−x[n]x[n+6]=-x[n] و a1=1−j3a_1=\tfrac{1-j}{3}. دنبالهٔ x[n]x[n] و انرژی آن در یک دوره را بیابید.

ak+4=aka_{k+4}=a_k با M=4M=4 ایجاب می‌کند x[n]=0x[n]=0 مگر آنکه nn مضربی از N/M=3N/M=3 باشد، و x[0],x[3],x[6],x[9]x[0],x[3],x[6],x[9] باقی می‌مانند. شرط نیم‌دوره می‌دهد x[6]=−x[0]x[6]=-x[0] و x[9]=−x[3]x[9]=-x[3] — دو مجهول. با Ω0=π/6\Omega_0=\pi/6:

a1=112[x[0]+x[3]e−jπ/2−x[0]e−jπ−x[3]e−j3π/2]=x[0]−j x[3]6.a_1=\frac{1}{12}\left[x[0]+x[3]e^{-j\pi/2}-x[0]e^{-j\pi}-x[3]e^{-j3\pi/2}\right]=\frac{x[0]-j\,x[3]}{6}.

برابرقراردادن با 1−j3\tfrac{1-j}{3} با نمونه‌های حقیقی می‌دهد x[0]=2x[0]=2, x[3]=2x[3]=2، پس یک دوره برابر {2,0,0,2,0,0,−2,0,0,−2,0,0}\{2,0,0,2,0,0,-2,0,0,-2,0,0\} و انرژی 4⋅4=164\cdot4=16 است. وارسی با پارسوال: هارمونیک‌های زوج صفر می‌شوند و هر یک از شش هارمونیک فرد دارای ∣ak∣2=29\lvert a_k\rvert^{2}=\tfrac29 است، پس ∑∣ak∣2=43=1612\sum\lvert a_k\rvert^{2}=\tfrac43=\tfrac{16}{12}. ✓

(ه) قطار ضربه. دنباله‌ای با دورهٔ تناوب N=5N=5 ضرایب aka_k دارد. ضرایب سری فوریهٔ y(t)=∑nx[n] δ(t−0.5n)y(t)=\sum_n x[n]\,\delta(t-0.5n) را بیابید.

T=NTs=2.5T=NT_s=2.5 و ck=ak/Ts=2akc_k=a_k/T_s=2a_k. دورهٔ تناوب NN هیچ نقشی در ضریب مقیاس ندارد؛ فقط TT و دورهٔ تناوب ckc_k نسبت به kk را تعیین می‌کند.

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

  • مقیاس‌دهی فرکانسی X(ejmΩ)X(e^{jm\Omega}). متناظر با بیش‌نمونه‌برداری (درج صفر) است: x[n/m]x[n/m] برای nn مضرب mm، وگرنه 00. طیف را فشرده/تکرار می‌کند و ذاتاً تغییرپذیر با زمان است — تلهٔ کلاسیکِ «به نظر LTI می‌آید ولی نیست».
  • آرگومان منفی X(e−jmΩ)X(e^{-jm\Omega}). وارون‌سازی و انبساط: x[−n/m]x[-n/m] روی مضارب mm. در ترکیب با یک ضریب فاز، به‌جای اعمال انتقال و وارون‌سازی با ترتیبی حدسی، از قاعدهٔ آفین استفاده کنید — e−jΩn0X(e−jΩ)↔x[n0−n]e^{-j\Omega n_0}X(e^{-j\Omega})\leftrightarrow x[n_0-n]، نه x[−n−n0]x[-n-n_0].
  • مزدوج در فرکانس در برابر مزدوج در زمان. X∗(ejΩ)↔x∗[−n]X^{*}(e^{j\Omega})\leftrightarrow x^{*}[-n] اما X∗(e−jΩ)↔x∗[n]X^{*}(e^{-j\Omega})\leftrightarrow x^{*}[n]. برای x[n]x[n] حقیقی، اولی فقط وارون‌سازی زمانی است.
  • طولانی‌کردن دورهٔ تناوب در DTFS. برای یک سیگنال حقیقی با دورهٔ تناوب اصلی NN، محاسبهٔ ضرایب با دورهٔ mNmN می‌دهد bmk=akb_{mk}=a_k و در بقیهٔ جاها صفر؛ هارمونیک واقعی را با حل mk=(new index)mk = (\text{new index}) بخوانید، سپس برای درایه‌های مزدوج‌شده ak∗=aN−ka_k^{*}=a_{N-k} را اعمال کنید.
  • طول پالس زوج در برابر فرد. یک پالس متقارن حول مبدأ به طول فرد 2M+12M+1 نیاز دارد؛ طول زوج LL را نمی‌توان حول یک عدد صحیح متقارن کرد و فاز خطی e−jΩ(L−1)/2e^{-j\Omega(L-1)/2} را به همراه دارد. زوجِ مجذور هسته/مثلث برای هر LL برقرار است، زیرا فاز در ∣⋅∣2\lvert\cdot\rvert^{2} حذف می‌شود.
  • پادتقارن نیم‌موج به NN زوج نیاز دارد. برای NN فرد نه انتقال نیم‌دوره‌ای وجود دارد و نه گزارهٔ «فقط هارمونیک‌های فرد».
  • DFT در برابر DTFS. X[k]=NakX[k]=N a_k مربوط به گسترش متناوب — همان اعداد، با نرمال‌سازی متفاوت. FFT صرفاً یک DFT سریع با O(Nlog⁡N)O(N\log N) است.
  • کانولوشن دایره‌ای در برابر خطی. ضرب دو DFT کانولوشن دایره‌ای به پیمانهٔ NN می‌دهد؛ برای بازیابی کانولوشن خطی معمولی، تا طول ≥N1+N2−1\ge N_1+N_2-1 صفر اضافه کنید.
  • بهرهٔ ثابت در برابر بهرهٔ وابسته به اندیس. c⋅x[n]c\cdot x[n] (با cc ثابت) LTI است؛ g[n] x[n]g[n]\,x[n] (با g[n]g[n] وابسته به اندیس) خطی ولی تغییرپذیر با زمان است.
  • ضریب مقیاس قطار ضربه. ck=ak/Tsc_k=a_k/T_s بر فاصلهٔ ضربه‌ها یعنی TsT_s تقسیم می‌شود، نه بر NN: ضریب 1/T=1/(NTs)1/T=1/(NT_s) سری زمان‌پیوسته با مجموع خام ∑x[n]e−jkΩ0n=Nak\sum x[n]e^{-jk\Omega_0 n}=N a_k روبه‌رو می‌شود و NN حذف می‌گردد.

نکتهٔ کنکوری: یک ضریب فاز ثابت e−jkΩe^{-jk\Omega} که در یک DTFT ضرب شود فقط یک تأخیر است — سیستم را LTI نگه می‌دارد. اما دو عملِ ddΩ\tfrac{d}{d\Omega} (که به وزنِ اندیسیِ n x[n]n\,x[n] تبدیل می‌شود) و X(ejmΩ)X(e^{jm\Omega}) (که به بیش/کم‌نمونه‌برداری تبدیل می‌شود) هر دو تغییرناپذیری با زمان را از بین می‌برند و خطی‌بودن را حفظ می‌کنند، و عملِ نمونه‌برداری می‌تواند علّی‌بودن را هم نقض کند. و هرگز اجازه ندهید یک قید اسکالر/انتگرالیِ تنها روی X(ejΩ)X(e^{j\Omega}) شما را قانع کند که سیگنالی «لزوماً» حقیقی، زوج یا فرد است — یک معادله نمی‌تواند یک ردهٔ تقارن کامل را مشخص کند؛ یک مثال نقض سریع (δ[n]\delta[n], jδ[n]j\delta[n], δ[n−2]\delta[n-2]) معمولاً مسئله را فیصله می‌دهد.

نکتهٔ کنکوری: پیش از انتگرال‌گیری یا زنجیره‌کردن خواص، طیف را به‌صورت ∑ncne−jΩn\sum_n c_n e^{-j\Omega n} بازنویسی کنید. نسبت سینوس‌ها رشته‌ای از یک‌هاست، مجذور آن یک مثلث است، انتگرال در برابر cos⁡(kΩ)\cos(k\Omega) برابر π(r[k]+r[−k])\pi\left(r[k]+r[-k]\right) است، و آرگومان مقیاس‌شده/منفی‌شده صرفاً هر نمونه را به n=ck+n0n=ck+n_0 جابه‌جا می‌کند.

تلهٔ رایج: «متناوب نسبت به kk» و «صفر به‌ازای برخی kkها» دو شرط متفاوت با معانی متفاوت در حوزهٔ زمان هستند — اولی x[n]x[n] را تُنُک می‌کند، دومی تقارن یا دورهٔ تناوب آن را مقید می‌کند. و وقتی معلوم شد دنباله‌ای عمدتاً صفر است، آن صفرها همچنان در مجموع‌هایی مانند ∑nf(x[n])\sum_n f(x[n]) به حساب می‌آیند، هرگاه f(0)≠0f(0)\ne0.

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

DTFT analysis
X(ejΩ)=∑n=−∞∞x[n] e−jΩnX(e^{j\Omega}) = \sum_{n=-\infty}^{\infty} x[n]\,e^{-j\Omega n}
DTFT synthesis
x[n]=12π∫2πX(ejΩ) ejΩn dΩx[n] = \frac{1}{2\pi}\int_{2\pi} X(e^{j\Omega})\,e^{j\Omega n}\,d\Omega
DTFS pair (period N)
x[n]=∑k=⟨N⟩akejk(2π/N)n,ak=1N∑n=⟨N⟩x[n] e−jk(2π/N)nx[n] = \sum_{k=\langle N\rangle} a_k e^{jk(2\pi/N)n},\quad a_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n]\,e^{-jk(2\pi/N)n}
DFT
X[k]=∑n=0N−1x[n] e−j2πkn/NX[k] = \sum_{n=0}^{N-1} x[n]\,e^{-j2\pi kn/N}
Frequency differentiation
n x[n]  ⟷  j dX(ejΩ)dΩn\,x[n] \;\longleftrightarrow\; j\,\frac{dX(e^{j\Omega})}{d\Omega}
Frequency scaling (upsampling)
x[n/m]  ⟷  X(ejmΩ)x[n/m] \;\longleftrightarrow\; X(e^{jm\Omega})
Time shift
x[n−n0]  ⟷  e−jΩn0X(ejΩ)x[n-n_0] \;\longleftrightarrow\; e^{-j\Omega n_0}X(e^{j\Omega})
Real-signal symmetry
X(ejΩ)=X∗(e−jΩ),ak∗=aN−kX(e^{j\Omega}) = X^{*}(e^{-j\Omega}),\qquad a_k^{*}=a_{N-k}
Parseval (DTFT)
∑n∣x[n]∣2=12π∫2π∣X(ejΩ)∣2 dΩ\sum_{n} |x[n]|^{2} = \frac{1}{2\pi}\int_{2\pi} |X(e^{j\Omega})|^{2}\,d\Omega
Time reversal / spectrum conjugation
x[−n]  ⟷  X(e−jΩ),x∗[−n]  ⟷  X∗(ejΩ)x[-n] \;\longleftrightarrow\; X(e^{-j\Omega}),\qquad x^{*}[-n] \;\longleftrightarrow\; X^{*}(e^{j\Omega})
Real part of a spectrum
Re[X(ejΩ)]  ⟷  12(x[n]+x∗[−n])\mathrm{Re}\left[X(e^{j\Omega})\right] \;\longleftrightarrow\; \tfrac12\left(x[n]+x^{*}[-n]\right)
Affine frequency argument
e−jΩn0X(ejcΩ)  ⟷  x[n−n0c] on multiples of c, 0 otherwisee^{-j\Omega n_0}X(e^{jc\Omega}) \;\longleftrightarrow\; x\left[\tfrac{n-n_0}{c}\right]\ \text{on multiples of } c,\ 0 \text{ otherwise}
Rectangular pulse (Dirichlet kernel)
∑n=−MMe−jΩn=sin⁡((2M+1)Ω/2)sin⁡(Ω/2)\sum_{n=-M}^{M} e^{-j\Omega n} = \frac{\sin\left((2M+1)\Omega/2\right)}{\sin(\Omega/2)}
Triangle / squared kernel
L−∣n∣ (∣n∣≤L−1)  ⟷  sin⁡2(LΩ/2)sin⁡2(Ω/2)L-\lvert n\rvert\ (\lvert n\rvert\le L-1) \;\longleftrightarrow\; \frac{\sin^{2}(L\Omega/2)}{\sin^{2}(\Omega/2)}
Integral against a cosine
∫2πX(ejΩ)cos⁡(kΩ) dΩ=π(x[k]+x[−k])\int_{2\pi} X(e^{j\Omega})\cos(k\Omega)\,d\Omega = \pi\left(x[k]+x[-k]\right)
DTFS time shift
x[n−n0]  ⟷  e−jk(2π/N)n0 akx[n-n_0] \;\longleftrightarrow\; e^{-jk(2\pi/N)n_0}\,a_k
Half-wave antisymmetry (N even)
x[n+N/2]=−x[n]  ⟺  ak=0 for even kx[n+N/2]=-x[n] \;\Longleftrightarrow\; a_k=0 \text{ for even } k
Sparsity–periodicity duality
ak+M=ak  ⟺  x[n]=0 unless n is a multiple of N/Ma_{k+M}=a_k \;\Longleftrightarrow\; x[n]=0 \text{ unless } n \text{ is a multiple of } N/M
Impulse-train Fourier series
y(t)=∑nx[n] δ(t−nTs):T=NTs,ck=akTsy(t)=\sum_{n} x[n]\,\delta(t-nT_s):\quad T=NT_s,\qquad c_k=\frac{a_k}{T_s}
#dtft#dtfs#dft#dirichlet-kernel#time-reversal-conjugation#half-wave-symmetry#impulse-train#system-classification

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

همهٔ 17 سؤال →