Given a vector $\mathbf{x} =\begin{bmatrix}x_0 & x_1 & \ldots & x_{N-1}\end{bmatrix} \in\mathbb{C}^N$ and its DFT $\mathbf{X} =\begin{bmatrix}X_0 & X_1 & \ldots & X_{N-1}\end{bmatrix}$ show that
As explained before, unfortunately the denominator of the transfer function does not admit a simple factorization. We can find the roots with a numerical package and we have
\begin{align*}
z_1 &= -0.4534\\
z_{2,3}&= 0.2267 \pm 1.4677j
\end{align*}
Since the magnitude of the complex roots is greater than one, the system is unstable.
\vspace{2em}
{\em This was the intended exercise:} If the feedback branch had had a gain of $+1$ as intended, the transfer function would have been
with $f_0=300$Hz, is raw-sampled into the discrete time sequence $x[n]= x_c(nT_s)$ with $T_s =5\cdot10^{-4}$s. Plot the DTFT of $x[n]$, $X(e^{j\omega})$.
\begin{sol}
\Solution
The sampling frequency is $F_s =1/T_s =(1/5)\cdot10^4$Hz $=2000$Hz. Since the continuous-time signal is just a linear combination of pure sinusoids, we can easily determine its highest frequency which is $F_{\max} =4\cdot300=1200$Hz. Since $F_{\max} > F_s/2$ we will have aliasing. With this in mind, the easiest way to solve the problem is by working in the time domain, although of course we can operate also in the frequency domain.
\vspace{1em}
\textbf{Working in the time domain:} by replacing the value for $T_s$ and we have
which yields the same plot as before. Note that we had to use the scaling property of the Dirac delta, $\delta(t/\alpha)\equiv\alpha\delta(t)$, which is perhaps not widely known. Hence the preference for the time-domain approach.
The channel has three usable bands with the following characteristics:
\begin{itemize}\setlength\itemsep{0pt}
\item band $W_1$, from 500Hz to 900Hz
\item band $W_2$, from 1000Hz to 1600Hz
\item band $W_3$, from 1900Hz to 2000Hz
\end{itemize}
To transmit the data you can use one or more configurable passband data transmitters. For each transmitter you can set the following parameters:
\begin{enumerate}\setlength\itemsep{0pt}
\item[1-] the center frequency for the trasmission band (in Hz)
\item[2-] the number of symbols per second (i.e. the symbol rate or Baud rate)
\item[3-] the number of bits per symbol (4, 8, 16, or 32)
\end{enumerate}
Each transmitter can also adjust its gain so that its transmitted signal reaches the maximum power allowed by the channel's power constraint. Because of different noise levels across the spectrum the resulting SNR levels for each subband are
\begin{itemize}\setlength\itemsep{0pt}
\item band $W_1$: SNR of 16dBs
\item band $W_2$: SNR of 10dBs
\item band $W_3$: SNR of 20dBs
\end{itemize}
The operational characteristic for the data transmitters is shown in this chart, where each curve shows the attainable probability of error as a function of the SNR; each curve corresponds to a different bit-per-symbol operation mode:
\dspLegendStyle(22,-.4){{linecolor=black}{4 bits per symbol}{linecolor=black,linestyle=dashed,dash=15pt\space3pt}{8 bits per symbol}{linecolor=black,linestyle=dashed,dash=5pt\space4pt}{16 bits per symbol}{linecolor=black,linestyle=dashed,linestyle=dotted}{32 bits per symbol}}
\end{dspPlot}
\end{center}
We want the transmission scheme to operate at an overall probability of error of $10^{-5}$ or less. Determine the maximum achievable transmission rate (in bits per second).
\begin{sol}
\Solution
Since we have three independent subbands, we will use a separate transmitter for each one. The number of symbols per second for each transmitter will be equal to the (positive) width of each subband while the trasmission center frequency will be equal to each subband's center frequency
\begin{itemize}
\item band $W_1$: 400 symbols per second, center frequency $f_1=700$Hz
\item band $W_2$: 600 symbols per second, center frequency $f_1=1300$Hz
\item band $W_3$: 100 symbols per second, center frequency $f_1=1950$Hz
\end{itemize}
To compute the total throughput we need to determine how many bits per symbol can be sent over each subchannel at the given reliability figure. Since we need to operate at most at $P_e=10^{-5}$ we can find the required minimum SNR for each operation curve by looking at the intersection with the line $P_e=10^{-5}$:
\item to transmit at 4 bits per symbol we need at least $\approx12$dB SNR
\item to transmit at 8 bits per symbol we need at least $\approx15$dB SNR
\item to transmit at 16 bits per symbol we need at least $\approx18$dB SNR
\item to transmit at 32 bits per symbol we need at least $\approx21$dB SNR
\end{itemize}
therefore:
\begin{itemize}
\item band $W_1$: SNR of 16dBs $\Rightarrow$ we can transmit at 8 bits per symbol
\item band $W_2$: SNR of 10dBs $\Rightarrow$ we cannot transmit since all possible rates will have $P_e>10^{-5}$ at this SNR
\item band $W_3$: SNR of 20dBs $\Rightarrow$ we can transmit at 16 bits per symbol
\end{itemize}
The final capacity of the transmission scheme is therefore $R = W_1\cdot8+ W_3\cdot16=4800$bps.
\end{sol}
\end{exercise}
\begin{exercise}{Quantization}
We have seen that one of the fundamental ingredients of JPEG compression is the {\em deadzone} quantizer, i.e. a quantizer which has a quantization interval centered around zero. To see the effects of the deadzone quantizer on SNR consider the following problem.
\vspace{1ex}
Assume $x[n]$ is an i.i.d. discrete-time signal whose values are over the $[-1, 1]$ interval. Consider the following uniform 2-bit quantizers for the interval:
\multido{\n=0+1}{3}{\pscircle*[linecolor=gray](!0 \n\space 1 sub 2 mul 3 div){2pt}}
\psplot[plotpoints=800,linewidth=2pt,linecolor=gray]{-0.99}{0.99}{x 3 mul 2 div round 2 mul 3 div}
\end{pspicture}
\end{tabular}
\end{center}
Both quantizers operate at two bits per sample but the deadzone quantizer "wastes" a fraction of a bit since it has only 3 quantization intervals instead of 4; for a uniformly distributed input over $[-1, 1]$, therefore, the Mean Square Error of the deadzone quantizer will be larger than the MSE of the standard quantizer.
Assume now that the probability distribution for each input sample is the following:
\[
P[x[n]=\alpha]=\begin{cases}
0 & \mbox{if $|\alpha| > 1$} \\
p & \mbox{if $|\alpha| =0$} \\
(1-p)/2 & \mbox{otherwise}
\end{cases}
\]
In other words, each sample is either zero with probability $p$ or drawn from a uniform distribution over the $[-1, 1]$ interval; we can express this distribution as a pdf like so:
\[
f(x)=\frac{1-p}{2} + p\delta(x)
\]
where the Dirac delta encodes the higher probability of the zero value.
Determine the minimum value of $p$ for which it is better to use the deadzone quantizer, i.e. the value of $p$ for which the MSE of the deadzone quantizer is smaller that the MSE of the uniform quantizer.
\vspace{1ex}
{\footnotesize\em Hint: remember that the formula for the MSE of a scalar quantizer over the $[-1, 1]$ interval (under the hypotheses of iid samples with distribution $f(x)$) is
For a uniform quantizer with $M$ quantization levels (as the ones in this exercise) if the input distribution is uniform (i.e. $f(x)=1/2$) the above simplifies to:
The number of quantization levels in the two quantizers are $M=M_n=4$ for the normal 2-bit quantizer and $M=M_d=3$ for the deadzone quantizer. Let's compute the MSE for the normal quantizer using the composite pdf for the input