{"id":112,"date":"2023-12-12T05:12:50","date_gmt":"2023-12-12T05:12:50","guid":{"rendered":"https:\/\/neilfoxman.com\/?page_id=112"},"modified":"2024-06-24T03:53:29","modified_gmt":"2024-06-24T03:53:29","slug":"discrete-fourier-series","status":"publish","type":"page","link":"https:\/\/neilfoxman.com\/?page_id=112","title":{"rendered":"Discrete Transforms"},"content":{"rendered":"\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"Intent\"><\/span>Intent<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p><a href=\"https:\/\/neilfoxman.com\/?page_id=273#Convolution_and_LTI_Systems\">Recall<\/a> we found that an input signal $x[n]$ fed into an LTI system with impulse response $h[n]$ generates an output $y[n]$ described by the convolution operation.<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>x[n] &amp;\\rightarrow y[n] = x[n] * h[n] \\\\<br>x[n] &amp;\\rightarrow y[n] = \\sum_{k=-\\infty}^{\\infty} x[k] h[n-k] \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>As it can be fairly hard to conceptualize how a system will shape an input to an output using convolution in the $n$ domain, it is commonly beneficial to transform the representation into a new domain that is easier to conceptualize.<\/p>\n\n\n\n<p>First, let&#8217;s imagine the simple case in discrete time where the input signal is merely some complex number raised to the $n$, that is<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>x[n] &amp;= z^n &amp;&amp;&amp;z \\in \\mathbb{C} \\\\<br>\\\\<br>y[n] &amp;= x[n] * h[n] \\\\<br>y[n] &amp;= h[n] * x[n] \\\\<br>y[n] &amp;= \\sum_{k=-\\infty}^{\\infty} h[k] x[n-k] \\\\<br>y[n] &amp;= \\sum_{k=-\\infty}^{\\infty} h[k] z^{n-k} \\\\<br>y[n] &amp;= z^n \\sum_{k=-\\infty}^{\\infty} h[k] z^{-k} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>The critical takeaway here is that the output is a scaled version of the input.  In this sense, $z^n$ is called an &#8220;eigenfunction&#8221; and the scaling factor $\\sum_{k=-\\infty}^{\\infty} h[k] z^{-k}$ is the &#8220;eigenvalue&#8221;.  This eigenvalue is dependent on both the original eigenfunction $z$ and the system impulse response, and is commonly expressed as<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>H(z) &amp;= \\sum_{k=-\\infty}^{\\infty} h[k] z^{-k} \\\\<br>\\\\<br>y[n] &amp;= H(z) z^n \\\\<br>y[n] &amp;= H(z) x[n] \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Also note that the variable of summation for $H(z)$ can be changed to $n$ to better envision the calculation for the impulse response i.e. the eigenvalue may also be calculated using<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>H(z) &amp;= \\sum_{n=-\\infty}^{\\infty} h[n] z^{-n} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Assuming an LTI system, if a new input signal is instead represented as a sum (series or integral) of $z$s, the output can be calculated as the sum of scaled $z$s.  That is<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>x[n] = z^n &amp;\\rightarrow y[n] = H(z) z^n \\\\<br>x[n] = \\sum_q a_q z_q^n &amp;\\rightarrow y[n] = \\sum_q a_q H(z_q) z_q^n \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Note that the eigenvalue is dependent on $z$, so if the input is expressed as a sum of $z$s, then $H(z)$ must be calculated for each $z$.<\/p>\n\n\n\n<p>The transformations here are all built on this principle and essentially just vary in what set of $z$s we can use to describe our input signal impulse response.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"Discrete_Time_Fourier_Series\"><\/span>Discrete Time Fourier Series<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>The Discrete Time Fourier Series is a method to represent any periodic signal as a series of periodic exponentials.<\/p>\n\n\n\n<p>Assume that a discrete signal may be represented as a series of periodic harmonics such that<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>x[n] &amp;= \\sum_k a_k e^{j k \\Omega_0 n} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p><a href=\"https:\/\/neilfoxman.com\/?page_id=251\">Recall<\/a> that there are only $N$ distinct harmonics due to discrete frequency periodicity. We can then represent this signal with<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>x[n] &amp;= \\sum_{k= \\langle N \\rangle} a_k e^{j k \\Omega_0 n} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>We can solve for Fourier Series coefficients, $a_k$ by scaling both sides with another complex exponential and summing both sides over one period with respect to $n$.<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br><br>x[n] e^{-j q \\Omega_0 n} &amp;= \\sum_{k= \\langle N \\rangle } a_k e^{j k \\Omega_0 n} e^{-j q \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k e^{j k \\Omega_0 n} e^{-j q \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k e^{j (k-q) \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;= \\sum_{k= \\langle N \\rangle } a_k \\sum_{n = \\langle N \\rangle} e^{j (k-q) \\Omega_0 n} \\\\<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>From how we defined the harmonics of $x[n]$ we know that $k \\in \\mathbb{Z}$.  Now if we further restrict our study to cases when $k &#8211; q \\in \\mathbb{Z} \\implies q \\in \\mathbb{Z}$, then we may note that the rightmost summation is 0 $ \\forall q, k: q  \\neq k$, but the exponential is 1 when $q=k$.<\/p>\n\n\n\n<p>If we know the LHS must be nonzero, then<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br><br>q &amp;= k \\\\<br>\\\\<br>\\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;= \\sum_{k= q} a_k \\sum_{n = \\langle N \\rangle} 1 &amp;&amp;(k = q) \\land (k, q \\in \\mathbb{Z}) \\\\<br><br>\\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;= a_q N &amp;&amp;(k = q) \\land (k, q \\in \\mathbb{Z}) \\\\<br><br>a_q &amp;= \\frac{1}{N} \\sum_{n = \\langle N \\rangle} x[n] e^{-j q \\Omega_0 n} &amp;&amp;(k = q) \\land (k, q \\in \\mathbb{Z})<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>and then replacing $q$ with $k$ for notation consistency, we have the synthesis and analysis equations:<\/p>\n\n\n\n<p>$$<br>\\boxed{<br>\\begin{align}<br>x[n] &amp;= \\sum_{k = \\langle N \\rangle} a_k e^{j k \\Omega_0 n} \\\\<br>a_k &amp;= \\frac{1}{N} \\sum_{n = \\langle N \\rangle} x[n] e^{-j k \\Omega_0 n}<br>\\end{align}<br>}<br>$$<\/p>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"Discrete_Time_Fourier_Transform\"><\/span>Discrete Time Fourier Transform<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>The Fourier Transform is an extension of the Fourier series that may represent finite aperiodic sequences in addition to periodic sequences.<\/p>\n\n\n\n<p>Consider the finite duration signal $x[n]$ where $x[n] = 0$ when $-N_1 \\leq n \\leq N_2$. Then imagine periodic version of that signal $\\tilde{x}[n]$ that repeats every $N \\geq N_1 + N_2$.<\/p>\n\n\n\n<figure data-wp-context=\"{&quot;imageId&quot;:&quot;69d31fed24bca&quot;}\" data-wp-interactive=\"core\/image\" data-wp-key=\"69d31fed24bca\" class=\"wp-block-image size-full wp-lightbox-container\"><img loading=\"lazy\" decoding=\"async\" width=\"810\" height=\"574\" data-wp-class--hide=\"state.isContentHidden\" data-wp-class--show=\"state.isContentVisible\" data-wp-init=\"callbacks.setButtonStyles\" data-wp-on--click=\"actions.showLightbox\" data-wp-on--load=\"callbacks.setButtonStyles\" data-wp-on-window--resize=\"callbacks.setButtonStyles\" src=\"http:\/\/neilfoxman.com\/wp-content\/uploads\/2023\/12\/image-20220708221432703.png\" alt=\"\" class=\"wp-image-207\" srcset=\"https:\/\/neilfoxman.com\/wp-content\/uploads\/2023\/12\/image-20220708221432703.png 810w, https:\/\/neilfoxman.com\/wp-content\/uploads\/2023\/12\/image-20220708221432703-300x213.png 300w, https:\/\/neilfoxman.com\/wp-content\/uploads\/2023\/12\/image-20220708221432703-768x544.png 768w\" sizes=\"auto, (max-width: 810px) 100vw, 810px\" \/><button\n\t\t\tclass=\"lightbox-trigger\"\n\t\t\ttype=\"button\"\n\t\t\taria-haspopup=\"dialog\"\n\t\t\taria-label=\"Enlarge\"\n\t\t\tdata-wp-init=\"callbacks.initTriggerButton\"\n\t\t\tdata-wp-on--click=\"actions.showLightbox\"\n\t\t\tdata-wp-style--right=\"state.imageButtonRight\"\n\t\t\tdata-wp-style--top=\"state.imageButtonTop\"\n\t\t>\n\t\t\t<svg xmlns=\"http:\/\/www.w3.org\/2000\/svg\" width=\"12\" height=\"12\" fill=\"none\" viewBox=\"0 0 12 12\">\n\t\t\t\t<path fill=\"#fff\" d=\"M2 0a2 2 0 0 0-2 2v2h1.5V2a.5.5 0 0 1 .5-.5h2V0H2Zm2 10.5H2a.5.5 0 0 1-.5-.5V8H0v2a2 2 0 0 0 2 2h2v-1.5ZM8 12v-1.5h2a.5.5 0 0 0 .5-.5V8H12v2a2 2 0 0 1-2 2H8Zm2-12a2 2 0 0 1 2 2v2h-1.5V2a.5.5 0 0 0-.5-.5H8V0h2Z\" \/>\n\t\t\t<\/svg>\n\t\t<\/button><\/figure>\n\n\n\n<p>Recall from the Fourier series that this periodic signal $\\tilde{x}[n]$ may be represented as<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>\\tilde{x}[n] &amp;= \\sum_{k=\\langle N \\rangle} a_k e^{j k \\Omega_0 n} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>and so we then have Fourier Series coefficients that can be described by<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br><br>a_k &amp;= \\frac{1}{N} \\sum_{n=\\langle N \\rangle} \\tilde{x}[n] e^{- j k \\Omega_0 n} \\\\<br>a_k &amp;= \\frac{1}{N} \\sum_{n=-N_2}^{N_1} \\tilde{x}[n] e^{- j k \\Omega_0 n} \\\\<br>a_k &amp;= \\frac{1}{N} \\sum_{n=-N_2}^{N_1} x[n] e^{- j k \\Omega_0 n} \\\\<br>a_k &amp;= \\frac{1}{N} \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j k \\Omega_0 n} \\\\<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Define the &#8220;envelope function&#8221; $X$ as<\/p>\n\n\n\n<p>$$<br>X(e^{j \\Omega}) = \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j \\Omega n} \\\\<br>$$<\/p>\n\n\n\n<p>Using this definition of $X$ we may rewrite the Fourier Series coefficients as<\/p>\n\n\n\n<p>$$<br>a_k = \\frac{1}{N} X(e^{j k \\Omega_0})<br>$$<\/p>\n\n\n\n<p>and then the synthesis equation can be rewritten as<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br><br>\\tilde{x}[n] &amp;= \\sum_{k=\\langle N \\rangle} \\left[ \\frac{1}{N} X(e^{j k \\Omega_0}) \\right] e^{j k \\Omega_0 n} \\\\<br><br>\\tilde{x}[n] &amp;= \\sum_{k=\\langle N \\rangle} \\left[ \\frac{\\Omega_0}{2 \\pi} X(e^{j k \\Omega_0}) \\right] e^{j k \\Omega_0 n} \\\\<br><br>\\tilde{x}[n] &amp;= \\frac{1}{2 \\pi} \\sum_{k=\\langle N \\rangle} X(e^{j k \\Omega_0}) e^{j k \\Omega_0 n} \\Omega_0 \\\\<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Now if we extend the &#8220;period&#8221; of our signal until $N \\rightarrow \\infty$, we essentially have a finite duration signal, and we may note the phenomena<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>\\Omega_0 \\rightarrow d \\Omega \\\\<br>k \\Omega_0 \\rightarrow \\Omega \\\\<br>\\tilde{x} \\rightarrow x \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>So that the Analysis and Synthesis equation become, for any signal,<\/p>\n\n\n\n<p>$$<br>\\boxed{<br>\\begin{align}<br>x[n] &amp;= \\frac{1}{2 \\pi}\\int_{2 \\pi} X(e^{j \\Omega}) e^{j \\Omega n} d \\Omega \\\\<br>X(e^{j \\Omega}) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j \\Omega n} \\\\<br>\\end{align}<br>}<br>$$<\/p>\n\n\n\n<p>Note that $X$ is essentially a function of $\\Omega$ and is periodic every $2 \\pi$ due to the frequency periodicity of discrete signals as shown below.<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>X(e^{j (\\Omega+2 \\pi)}) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j (\\Omega+2 \\pi) n} \\\\<br><br>X(e^{j (\\Omega+2 \\pi)}) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j \\Omega n} e^{- j \\cdot 2 \\pi \\cdot n} \\\\<br><br>X(e^{j (\\Omega+2 \\pi)}) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] e^{- j \\Omega n} \\\\ \\end{align}<br>$$<\/p>\n\n\n\n<p>In summary, these equations suggest that<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>Any finite discrete signal (periodic or aperiodic) may be represented by a summation of complex exponentials in the frequency domain $\\Omega$<\/li>\n\n\n\n<li>Unique complex exponentials used to describe a signal reside in a frequency range of width $2 \\pi$<\/li>\n\n\n\n<li>The amplitude of a signal&#8217;s constituent complex exponentials are described by the function $\\frac{1}{2 \\pi} X(e^{j \\Omega})$ which may, itself, be complex.<\/li>\n<\/ol>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"Z_Transform\"><\/span>Z Transform<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>Note that for the Discrete Fourier Transform, we have been using the analysis function defined by<\/p>\n\n\n\n<p>$$<br>X(e^{j \\Omega})<br>$$<\/p>\n\n\n\n<p>Which for any $\\Omega \\in \\mathbb{R}$ is defined on a unit circle in the imaginary plane and repeats every $2 \\pi$. If we further generalize the independent variable of this function to accept any complex input, this would be described by<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>X(r e^{j \\Omega})<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>This quantity is often represented by $z = r e^{j \\Omega}$, and in the special case where $r=1$ we have the Fourier Transform.  Re-doing the same analysis as above, but for a general complex variable $z$, we have (in &#8220;fast motion&#8221;)<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>\\tilde{x}[n] &amp;= \\sum_{k= \\langle N \\rangle} a_k (r e^{j k \\Omega_0})^n \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0 })^{-n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k (r e^{j k \\Omega_0 })^n (r e^{j q \\Omega_0 })^{-n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0 })^{-n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k r^n e^{j k \\Omega_0 n} r^{-n} e^{-j q \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0 })^{-n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k r^{n-n} e^{j (k-q) \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0 })^{-n} &amp;= \\sum_{n = \\langle N \\rangle} \\sum_{k= \\langle N \\rangle } a_k e^{j (k-q) \\Omega_0 n} \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n](r e^{j q \\Omega_0 })^{-n} &amp;= \\sum_{k= \\langle N \\rangle } a_k \\sum_{n = \\langle N \\rangle} e^{j (k-q) \\Omega_0 n} \\\\<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>We effectively get the same implication as before: for the kth harmonic, $k, q \\in \\mathbb{Z} \\implies$<\/p>\n\n\n\n<p>$$<br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0})^{-n} = \\cases{<br>\\begin{align}<br>\\sum_{k= \\langle N \\rangle } a_k \\sum_{n = \\langle N \\rangle} 1 &amp;&amp; k = q \\\\<br>0 &amp;&amp; k \\neq q \\\\<br>\\end{align}<br>}<br>$$<\/p>\n\n\n\n<p>Assuming the LHS is nonzero, we can say<\/p>\n\n\n\n<p>\\begin{align}<br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0})^{-n} &amp;= \\sum_{k=q} a_k \\sum_{n = \\langle N \\rangle} 1 \\\\<br><br>\\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0})^{-n} &amp;= a_q N \\\\<br><br>a_q &amp;= \\frac{1}{N} \\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j q \\Omega_0})^{-n} \\\\<br><br>a_k &amp;= \\frac{1}{N} \\sum_{n = \\langle N \\rangle} \\tilde{x}[n] (r e^{j k \\Omega_0})^{-n} \\\\<br>\\end{align}<\/p>\n\n\n\n<p>Assuming $\\forall n \\in \\langle N \\rangle, \\tilde{x}[n] = x[n]$, we can expand the region of integration such that<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>a_k &amp;= \\frac{1}{N} \\sum_{n=\\langle N \\rangle} \\tilde{x}[n] (r e^{j k \\Omega_0})^{-n} \\\\<br>a_k &amp;= \\frac{1}{N} \\sum_{n=-\\infty}^{\\infty} x[n] (r e^{j k \\Omega_0})^{-n} \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Now we define<\/p>\n\n\n\n<p>$$<br>X(re^{j \\Omega}) = \\sum_{n=-\\infty}^{\\infty} x[n] (r e^{j \\Omega})^{-n} \\\\<br>$$<\/p>\n\n\n\n<p>And expanding upon our earlier findings<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>a_k &amp;= \\frac{1}{N} X(r e^{j k \\Omega_0}) \\\\<br>\\\\<br>\\tilde{x}[n] &amp;= \\sum_{k= \\langle N \\rangle} a_k (r e^{j k \\Omega_0})^n \\\\<br><br>\\tilde{x}[n] &amp;= \\sum_{k=\\langle N \\rangle} \\left[ \\frac{1}{N} X(r e^{j k \\Omega_0}) \\right] (r e^{j k \\Omega_0})^n \\\\<br><br>\\tilde{x}[n] &amp;= \\sum_{k=\\langle N \\rangle} \\left[ \\frac{\\Omega_0}{2 \\pi} X(r e^{j k \\Omega_0}) \\right] (r e^{j k \\Omega_0})^n \\\\<br><br>\\tilde{x}[n] &amp;= \\frac{1}{2 \\pi} \\sum_{k=\\langle N \\rangle} X(r e^{j k \\Omega_0}) \\cdot (r e^{j k \\Omega_0})^n \\Omega_0 \\\\<br>\\\\<br>N &amp;\\rightarrow \\infty \\\\<br>\\Omega_0 &amp;\\rightarrow d \\Omega \\\\<br>k \\Omega_0 &amp;\\rightarrow \\Omega \\\\<br>\\tilde{x} &amp;\\rightarrow x \\\\<br>\\\\<br>x[n] &amp;= \\frac{1}{2 \\pi}\\int_{2 \\pi} X(r e^{j \\Omega}) \\cdot (r e^{j \\Omega})^n d \\Omega \\\\<br><br>\\end{align}<br>$$<\/p>\n\n\n\n<p>To put in terms of $z$, the above equation can be conceptualized as an integral along a contour with fixed $r$ and varying $\\omega$ along one full counterclockwise revolution.<\/p>\n\n\n\n<p>$$<br>\\begin{align}<br>z &amp;= r e^{j \\Omega} \\\\<br>dz &amp;= j r e^{j \\Omega} d \\Omega \\\\<br>dz &amp;= j z d \\Omega \\\\<br>d \\Omega &amp;= \\frac{1}{j z} dz \\\\<br>\\\\<br>x[n] &amp;= \\frac{1}{2 \\pi} \\oint_{\\Omega = \\langle 2\\pi \\rangle} X(z) z^n \\frac{1}{j z} dz \\\\<br>x[n] &amp;= \\frac{1}{2 \\pi j} \\oint_{\\Omega = \\langle 2\\pi \\rangle} X(z) z^{n-1} dz \\\\<br>\\end{align}<br>$$<\/p>\n\n\n\n<p>Integration along counterclockwise closed circular contour with fixed $r$ chosen from within ROC.  As this approach is unusual, partial fraction expansion and transform tables are generally preferred.  In closing, the Analysis and Synthesis equations for the z-Transform are then<\/p>\n\n\n\n<p>$$<br>\\boxed{<br>\\begin{align}<br>X(z) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] z^{-n} \\\\<br>x[n] &amp;= \\frac{1}{2 \\pi j} \\oint_{\\Omega = \\langle 2\\pi \\rangle} X(z) z^{n-1} dz \\\\<br>\\end{align}<br>}<br>$$<\/p>\n\n\n\n<p>Note that the real part of $z$ in the analysis equation suggests that the summation does not always converge, therefore a region of Convergence or ROC must be defined for each transform.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"Region_of_Convergence\"><\/span>Region of Convergence<span class=\"ez-toc-section-end\"><\/span><\/h3>\n\n\n\n<p>\\begin{align}<br>X(z) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] z^{-n} \\\\<br>X(r e^{j \\Omega}) &amp;= \\sum_{n=-\\infty}^{\\infty} x[n] r^{-n} e^{-j \\Omega n} \\\\<br>\\end{align}<\/p>\n\n\n\n<p>Notes:<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>When $|z| = r = 1$, then we effectively get the Fourier Transform. Convergence in that case is determined solely by $x[n]$. Otherwise, convergence depends on both $x[n]$ and $r$.<\/li>\n\n\n\n<li>ROC must never contain poles (doesn&#8217;t converge at those locations)<\/li>\n\n\n\n<li>If $x[n]$ is finite duration and absolutely integrable, then the ROC is the entire $z$ plane except possibly $z=0$ and\/or $z \\rightarrow \\infty$<\/li>\n\n\n\n<li>If $X(z)$ is rational, then ROC is bounded by poles, or extends to infinity.<\/li>\n\n\n\n<li>If the signal is right sided and circle $|z| = r_0$ is in the ROC, then $|r_1| &gt; |r_0|$ is also in the ROC.<br><br>Note that $r_1^{-n}$ will converge faster than $r_0^{-n}$<br><br>Given a right sided signal that starts at $n=N_1$, the ROC of $z$ will include infinity when $N_1 \\geq 0$ but will not include infinity when $N_1 &lt; 0$ because then there will be some terms with positive exponent in the summation.  The case when $N_1 &lt; 0$ is further analyzed with<br><br>$$<br>\\begin{align}<br>X(z) &amp;= \\sum_{n=N_1}^{\\infty} x[n] z^{-n} &amp;&amp;N_1 &lt; 0 \\\\<br>X(z) &amp;= \\sum_{n=N_1}^{-1} x[n] z^{-n} + \\sum_{n=0}^{\\infty} x[n] z^{-n} &amp;&amp;N_1 &lt; 0 \\\\<br>X(z) &amp;= \\sum_{n=1}^{|N_1|} x[n] z^{n} + \\sum_{n=0}^{\\infty} x[n] z^{-n} &amp;&amp;N_1 &lt; 0 \\\\<br>\\end{align}<br>$$<br><br>Note that the leftmost term would then include values where $z=\\infty$ would be raised to a positive number which does not converge.<br><\/li>\n\n\n\n<li>If the signal is left sided and circle $|z| = r_0$ is in the ROC, then all values $0 &lt; |z| &lt; r_0$ are also in the ROC.  ROC of $z$ will include 0 when $N_2 \\leq 0$ but will not include 0 when $N_2 &gt; 0$.<br><br>Consider the left-sided signal where the rightmost nonzero value is at $n=N_2$.  The z-Transform is<br><br>$$<br>\\begin{align}<br>X(z) &amp;= \\sum_{n=-\\infty}^{N_2} x[n] z^{-n} \\\\<br>X(z) &amp;= \\sum_{n=-\\infty}^{0} x[n] z^{-n} + \\sum_{n=1}^{N_2} x[n] z^{-n} \\\\<br>\\end{align}<br>$$<br><br>ROC does not include $r=0$ because the rightmost summation would put 0 in the denominator.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Intent Recall we found that an input signal $x[n]$ fed into an LTI system with impulse response $h[n]$ generates an output $y[n]$ described by the convolution operation. $$\\begin{align}x[n] &amp;\\rightarrow y[n] = x[n] * h[n] \\\\x[n] &amp;\\rightarrow y[n] = \\sum_{k=-\\infty}^{\\infty} x[k] h[n-k] \\\\\\end{align}$$ As it can be fairly hard to conceptualize how a system will shape [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"parent":251,"menu_order":1,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-112","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/pages\/112","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/neilfoxman.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=112"}],"version-history":[{"count":170,"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/pages\/112\/revisions"}],"predecessor-version":[{"id":1191,"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/pages\/112\/revisions\/1191"}],"up":[{"embeddable":true,"href":"https:\/\/neilfoxman.com\/index.php?rest_route=\/wp\/v2\/pages\/251"}],"wp:attachment":[{"href":"https:\/\/neilfoxman.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=112"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}