Signal Analysis


Quadrature Phase-shift Keying (QPSK) is a widely used method of transferring digital data by changing or modulating the phase of a carrier signal. In QPSK digital data is represented by 4 points around a circle which correspond to 4 phases of the carrier signal. These points are called symbols. Figure 1 shows this mapping. Gray coding is used in this mapping so that no two adjacent symbols differ by more than 1 bit. This helps in reducing the Bit Error Rate (BER).

Grey Coded Constelation
QPSK Constellation

To model QPSK, we need:
-Input Data to transmit/modulate with carrier
-Noisy Medium to transfer the data
-Demodulation of transmitted data.
-Comparison of original signal, and the demodulated signal to calculate the Bit Error Rate (BER)

For the input signal we can use the following command:
x_signal = randint(1000,1,4);

To modulate the data against the grey coded constellation, you can use the genqammod matlab function.
y = genqammod(x_signal,constell_gray);

I then made a for loop, each iteration of the loop represents the Signal to Noise Ratio (SNR).
for SNR=0:2:10
y_noisy = awgn(y,SNR,'measured');

From the below graphs, you can see that the higher the SNR, the better the received signal.
SNR = 6
QPSK SNR = 6

SNR = 10
QPSK with SNR of 10

QPSK is simple and works fine for most applications. To increase bandwidth, you can use more than 4 symbols to represent your data. For instance, digital cable companies such as charter and Comcast, use 256-QAM. 256 QAM looks similar to the constellation above, except there are 64 symbols in each quadrant arranged in a rectangular fashion. The symbols are determined by modulating the Phase, and Amplitude of the carrier signal.

To demonstrate how effective QAM is, I modulated the original data input using 4-QAM instead of QPSK. The graph is shown below:
4-QAM
4_QAM

Spent some time on the problem for homework in my Communication Theory Course. The problem is a Fourier Transform of a damped sinusoidal wave. I plan on adding Matlab Graphs in a couple days…

Here, g(t) is the signal in the time domain. G(f) is the transform to the frequency domain. Enjoy!

g(t)=e^{-t}sin(2\pi{f_{c}}t)u(t)

G(f)=\int_{\small-\infty}^{\small\infty}e^{-t}sin(2\pi{f_{c}}t)u(t)e^{-j2\pi{f}t}dt

G(f)=\int_{\small0}^{\small\infty}e^{\small-t}\frac{e^{j2\pi{f_{c}}t}-e^{-j2\pi{f_{c}}t}}{2j}e^{-j2\pi{f}t}dt

G(f)=\frac{1}{2j}\int_{\small0}^{\small\infty}e^{j2\pi{f_{c}}t-t-j2\pi{f}t}-e^{-j2\pi{f_{c}}t-t-j2\pi{f}t}dt

G(f)=\frac{1}{2j}\biggl[\frac{e^{-t(-j2\pi{f_{c}}+1+j2\pi{f})}}{-j2\pi{f_{c}}+1+j2\pi{f}}\biggr]_0^{\infty}-\frac{1}{2j}\biggl[\frac{e^{-t(j2\pi{f_{c}}+1+j2\pi{f})}}{j2\pi{f_{c}}+1+j2\pi{f}}\biggr]_0^{\infty}

G(f)=\frac{1}{2j}\biggl[\frac{1}{j2\pi{f_{c}}-1-j2\pi{f}}\biggr]-\frac{1}{2j}\biggl[\frac{1}{-j2\pi{f_{c}}-1-j2\pi{f}}\biggr]

G(f)=\frac{1}{2j}\biggl[\frac{-j2\pi{f_{c}}-1-j2\pi{f}-j2\pi{f_{c}}+1+j2\pi{f}}{(j2\pi{f_{c}}-1-j2\pi{f})(-j2\pi{f_{c}}-1-j2\pi{f})}\biggr]

G(f)=\frac{1}{2j}\biggl[\frac{j4\pi{f_{c}}}{4\pi^2f_{c}^2+4\pi^2f_{c}f+1+j2\pi{f}-4\pi^2f_{c}f+j2\pi{f}-4\pi^2{f}^2}\biggr]

G(f)=\frac{2\pi{f_{c}}}{1+4\pi^2({f_{c}^2}-f^2)+j4\pi{f}}

Graph Damped sin with Frequency of 30Hz
Damped Sinusoidal Wave

Graph of Frequency Response
Damped Sinusoidal Wave Frequency Response