## Permutation Trellis Coded Multi-level FSK Signaling to Mitigate Primary User Interference in Cognitive Radio Networks

Raghed El-Bardan, Student Member, IEEE, Engin Masazade, Member, IEEE, Onur Ozdemir, Member, IEEE, Yunghsiang S. Han, Fellow, IEEE, and Pramod K. Varshney, Fellow, IEEE

## Abstract

We employ Permutation Trellis Code (PTC) based multi-level Frequency Shift Keying signaling to mitigate the impact of Primary Users (PUs) on the performance of Secondary Users (SUs) in Cognitive Radio Networks (CRNs). The PUs are assumed to be dynamic in that they appear intermittently and stay active for an unknown duration. Our approach is based on the use of PTC combined with multi-level FSK modulation to make the SU transmissions robust against PU interference and noise disturbances and help an SU improve its data rate by increasing its transmission bandwidth while operating at low power and not creating destructive interference for PUs. We evaluate system performance by obtaining an approximation for the actual Bit Error Rate (BER) using properties of the Viterbi decoder and carry out a thorough performance analysis in terms of BER and throughput. The results show that the proposed coded system achieves i) robustness by ensuring that SUs have stable throughput in the presence of heavy PU interference and ii) improved resiliency of SU links to interference in the presence of multiple dynamic PUs.

R. El-Bardan and P. K. Varshney are with the Department of Electrical Engineering and Computer Science, Syracuse University, Syracuse, NY 13244 USA.

- E. Masazade is with the Department of Electrical and Electronics Engineering, Yeditepe University, Turkey.
- O. Ozdemir is with Boston Fusion Corp., 1 Van de Graaff Drive, Suite 107, Burlington, MA 01803 USA.
- Y. S. Han is with the Department of Electrical Engineering at National Taiwan University of Science and Technology, Taipei, Taiwan. He visited Syracuse University during 2012-13.

Email: { raelbard, oozdemir, varshney } @syr.edu, engin.masazade@yeditepe.edu.tr, yshan@mail.ntpu.edu.tw

## Index Terms

Cognitive Radio Networks, Primary User Interference, Permutation Trellis Codes, FSK Systems

## I. INTRODUCTION

Today a large portion of the frequency spectrum assigned by the Federal Communications Commission (FCC) is used intermittently resulting in inefficient usage of the assigned frequency bands. In order to accommodate the increased demand for spectrum in wireless applications, Dynamic Spectrum Access (DSA) has been proposed as the new communication paradigm for wireless systems where the existing spectrum is utilized opportunistically [1], [2]. The radios using this technique are called Cognitive Radios (CRs). In a Cognitive Radio Network (CRN), Primary Users (PUs) are the licensed users who have the right to access their spectrum at any time. When Secondary Users (SUs) have data to transmit, they look for spectrum holes which are the frequency bands not being used by primary users and transmit their data over them [3], [4]. In order to ensure a non-degraded PU operation in the network, existing literature has considered several effective protection techniques, e.g., beamforming of SU signals and putting interference temperature limits to satisfy PU signal-to-interference plus noise constraint [5], [6]. In addition, there are scenarios in which SUs can co-exist with the PUs as long as they do not violate the quality of service requirements of the PUs [7], [8].

Orthogonal frequency division multiplexing (OFDM) has been suggested as a multi-carrier communication candidate for CR systems where the available spectrum is divided into subcarriers each of which carries a low rate data stream [9]-[13]. A typical approach for a CR using OFDM is to sense the PU activity over the sub-carriers and then adjust its communication parameters accordingly. The goal is to protect PUs as well as intended SU receivers from possible collisions resulting from the use of the same sub-carriers. Continuous spectrum sensing and reformation of wireless links may result in substantial performance degradation for SUs [14]. For instance, the throughput of the secondary system is affected by the time spent for channel sensing. When an SU spends more time on spectrum sensing, a smaller number of information bits will be transmitted over a smaller interval of time resulting in reduced system throughput. On the other hand, decreasing sensing duration may result in a larger probability of making incorrect decisions, thereby decreasing the throughput of both SU and PU. In this regard, the

authors in [15] numerically analyze the trade-offs between throughput and sensing duration in addition to coding blocklength and buffer constraints. In [16], the authors study the problem of designing a sensing slot duration to maximize the achievable throughput for the SUs under the constraint that the PUs are sufficiently protected.

Error correction coding (ECC) has been proposed for CRNs as one possible technique to better protect SU links and prevent degradation in their performance. In [10], it is assumed that an SU transmitter vacates the band once a PU is detected. Due to the sudden appearance of a PU, rateless codes have been considered to compensate for the packet loss in SU data which is transmitted through parallel subchannels. The authors in [11] consider the design of two efficient anti-jamming coding techniques for the recovery of lost transmitted packets via parallel channels, namely rateless and piecewise coding. Similar to the spectrum model defined in [10], a system's throughput and goodput are analyzed and a performance comparison is carried out between these two coding techniques. For an OFDM-based CRN presented in [12], SU transmitters and receivers continuously sense the spectrum, exchange information and decide on the available and unavailable portions of the frequency spectrum. Depending on frequency availability, an appropriate Reed-Solomon coding scheme is used to retrieve the bits transmitted over the unavailable portions of the frequency spectrum. The authors in [13] further explore LowDensity-Parity-Check (LDPC) codes in an OFDM scheme where a switching model is considered for dynamic and distributed spectrum allocation. They also analyze the effects of errors during PU detection on channel capacity and system performance. The switch is assumed to be open for each SU detecting a PU. When the switch is open, the channel is modeled as a binary erasure channel (BEC) and the cognitive transmitter continues to transmit its message allowing bits to be erroneous when received. Another major use of the ECC schemes is presented in [17] where the authors study the performance of cooperative relaying in cognitive radio networks using a rateless coding error-control mechanism. They assume that an SU transmitter participates in PU's transmission as a relay instead of vacating the band in order to reduce the channel access time by a PU. Since the use of rateless codes allows an SU receiver to decode data regardless of which packets it has received as long as enough encoded packets are received, these codes are very suitable for cooperative schemes. The authors in [18] propose an end-to-end hybrid ARQ scheme in CRNs consisting of unidirectional opportunistic links to reduce the number of retransmissions with a fixed throughput offset. Their error control approach is based on coded

cooperation among paths and amplify-and-forward relaying of packets within a path such that this hybrid ARQ works for CRNs even if some coded data are missing. The authors implement their scheme using convolutional codes combined with BPSK modulation.

In this paper, we investigate the error correction performance of Permutation Trellis Codes (PTCs) [19], [20] combined with multi-level Frequency Shift Keying modulation systems in CRNs. The proposed PTC based framework is quite general and is applicable to many systems beyond CRNs where interference is an issue. For example, the framework has been applied to power line communications [19], [20] where strong interferers are assumed to be always present as an extreme case. In this paper, we show a much wider applicability of the PTC based framework as an interference mitigation approach by applying it to CRNs where the special challenge is to guarantee reliable communication by SUs in the presence of intermittent PUs that stay active for an unknown duration without the need of detecting them. These dynamic PUs serve as the source of strong interference for SU communications. Another challenge lies in the susceptibility of SU transmissions in a CRN to intentional interference like jamming attacks. In such a case, it is extremely important to guarantee a minimum degree of operability of the network. In other words, additive white Gaussian noise (AWGN) is not the critical limitation here, but rather an unpredictable and/or intermittent interference, e.g., PU transmissions and jamming attacks. Therefore, the emphasis of the proposed PTC based framework is on robustness, rather than on data rate or bandwidth use. To the best of our knowledge, our work is the first to consider the use of PTCs in CRNs. The motivation for this proposed approach is to overcome the tremendous degrading effect of PU interference on an SU transmission and to provide a stable level of reliable information reception regardless of how severe or prolonged a dynamic PU activity is and without requiring an accurate detection of it. In our model, an SU transmits its own information using low power concurrently with PU transmissions without the need to relay PU's traffic [17] or to vacate the band [10]. Different from [12], we assume that no information exchange or spectrum negotiation takes place between the SU transmitter-receiver pair. Our proposed scheme is different from existing work where communication sessions are carried over parallel subchannels [10]-[13]. This is due to the fact that our scheme is based on multi-level FSK modulation with PTC using a single subchannel at a time. By using PTC, continuous channel sensing by the SU transmitter-receiver pair is no longer required, since an appropriate PTC can cope with high levels of PU interference on a given SU link. Thus, SUs

no longer suffer from the huge overhead created by continuous sensing of the spectrum.

We consider a similar system model as devised in [19], [20] and carry out a thorough performance analysis in terms of BERs and throughput for dynamically varying interference. In [19], [20], the authors were interested in the design of PTCs that mitigate the impact of permanent interferences that exist in power line communications. In their work, they conducted a simulation study and presented results on the performance of different PTCs with respect to the probability of an element of the code matrix being in error in the presence of different noises on one or more channels. They also investigated the choice of distance increasing/conservative/reducing mappings for a given PTC with respect to narrowband interference that always exists on one band in the presence of background noise in terms of BER. The authors in [19], [20] provide a BER performance analysis via simulations and not analytically.

The analytical BER analysis of a PTC coded M-FSK system is imperative to determine its link quality. In other words, the BER estimate is a useful tool for cross-layer design. Given that the PU stays active once it starts transmitting, the exact bit error rate (BER) is derived using an exhaustive search [21]. Hence, the analytical evaluation of BER becomes computationally prohibitive for large codes. Rather than using an exhaustive search, in this paper, we develop an approximation of BER using the properties of the Viterbi decoder. Also, in our previous work [21], a very special case is considered where a PU stays active once it starts transmission. In this paper, we consider a more practical scenario where a 2-state Markov chain is employed to model the dynamic PU activity in terms of alternating On-Off periods. It should be noted that the work on [19] and [20] also considered an extreme case where the authors assumed that the interference was always present and did not consider dynamically varying interference. In the more general framework presented in this paper, the performance results provided, e.g., BER, are intuitive and more useful in comparison to the worst-case guaranteed performance presented in [21].

Under the same resources such as energy per bit, transmission bandwidth, and transmission time, we compare the throughput performance of the proposed communication system with a simple uncoded opportunistic M-FSK system and a convolutionally coded OFDM system using BPSK modulation. The simulation results show that the proposed system outperforms the aforementioned systems in the presence of heavy PU interference. We also compare the BER performance of the proposed communication system with an LDPC code coupled with an M -

FSK modulation system for static and dynamic activities of a PU. It is shown that the proposed scheme outperforms the other one. Moreover, we present results on the degree of resiliency to interference that a PTC can guarantee for the SU link in the presence of interference by multiple dynamic PUs.

To summarize, the main contributions of the paper are:

- · We consider a more practical scenario in which we represent the cognitive radio network using a dynamic PU channel occupancy model including the case of multiple PUs.
- · To mitigate interference created by these dynamic PUs, we propose the use of PTC-based multi-level FSK signaling that incurs no overhead costs as no channel sensing is required before SU transmissions.
- · Using the proposed PTC based framework, the SU transmissions achieves robustness against the dynamic PUs' interference and noise disturbances and, at the same time, the resiliency of SU links improves.
- · A thorough analytical performance analysis of the proposed scheme is presented in terms of BER and throughput for static and dynamic activities of multiple PUs.

The rest of the paper is organized as follows. In Section II, we provide the system model and introduce the concept of PTC. In Section III, we introduce hard-decision Viterbi decoding of the system. We also approximate the actual BER through the use of a specified finite number of paths forming the Viterbi trellis. Section V presents numerical results demonstrating the improvement in the performance of the PTC-based multi-level FSK system, validating the suggested BER approximation, supporting the effectiveness of the proposed communication scheme, as well as demonstrating the efficiency of PTC in the presence of multiple dynamic PU interference. Finally, we conclude the paper and address some future research directions in Section VI.

## II. SYSTEM MODEL

We consider a simple cognitive radio network as shown in Fig. 1, which consists of a PU transmitter-receiver pair and an SU transmitter-receiver pair. The SU pair is assumed to be located within the transmission range of the PU. The transmission range shown in Fig. 1 is the maximum distance covered by a PU transmission such that the signal to interference plus noise ratio (SINR) at the PU receiver equals a minimum threshold value, SINR ∗ PU . In this model, we assume a free space path loss model and a Line-Of-Sight (LOS) AWGN channel. The power in

Fig. 1: A simple cognitive radio network.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-92946945-f681-401b-a4ba-bfeb750b02db.jpg)

the transmitted signal is P t , so the received power is given in [22] by

$$\Phi _ { \Phi } = \ P _ { T } \cdot \left ( \frac { \sqrt { G _ { l } } \lambda } { 4 \pi d } \right ) ^ { 2 },$$

where √ G l is the product of the transmit and receive antenna field gains in the LOS direction and λ is the signal wavelength defined as the ratio of the speed of light to the center frequency of the band in operation, C f j . In this model, we assume omnidirectional antennas so √ G l = 1 . Both licensed and unlicensed users coexist and they can simultaneously operate over the same band. One further assumption made here is that interference, created by an SU, that deteriorates the QoS of the PU is negligible compared to the received PU power. This is due to the fact that the transmission power of PU is much larger than the transmission power of SU, i.e., P SU T glyph[lessmuch] P PU T . In this approach, the SINR at the PU receiver [23] given in (2) satisfies SINR ≥ SINR ∗ PU , that is,

$$S I N R \ = \ \frac { P _ { R } ^ { P U } } { N _ { 0 } + P _ { I } ^ { S U } } \geq S I N R _ { P U } ^ { * },$$

where P SU I is the SU interference at the PU receiver and SINR ∗ PU is the minimum threshold above which a PU transmission is received successfully. It is assumed that the SU transmitter knows the location of its corresponding receiver possibly by means of extra signaling. An overview of the signal processing model that combines the PTC scheme with the multi-level FSK communication system is provided in Fig. 2. The information bits, m , are loaded into a rate R = m n convolutional encoder in parallel, e.g., R = 1 2 with a two-stage shift register and generator 7, 5 (octal) [19]. It is well-known that due to the incorporation of a channel coding scheme in a

Fig. 2: Block diagram of the Coded multi-level FSK System

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-aa2c86a6-2e11-46a1-81fc-3b4e5746a11c.jpg)

communication system, the overall transmission energy per bit is less than that without a coding scheme. However, due to its error correction capabilities, the purpose of convolutional encoding is to improve the performance of the system.

For each of the convolutional coded output symbols, the PTC encoder assigns a certain code matrix which is transmitted over both time and frequency. We present an example for illustration, in which the coded symbol '01' is mapped onto '213'. This permutation code matrix is to be transmitted in both time and frequency domains resulting in a 3 × 3 ( H × H ) binary code matrix as given below.

$$\mathbf T _ { i } = \left [ \begin{array} { c c c } 0 & 1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{array} \right ],$$

where 1 ≤ i ≤ M M . = 2 m denotes the number of symbols and H defines the number of frequency bands as well as the number of time steps used in transmitting the outputs of the encoder. In this paper, we define the PTC-based multi-level FSK system as the H -FSK system since the PTC-based multi-level FSK modulated system leads to an H x H binary code matrix. According to (3), transmission takes place on f 2 , f 1 , and f 3 (first, second and third columns of T i ) corresponding to the time steps T s , 2 T s , and 3 T s respectively. Tables I and II present the symbol (with 1 bit and 2 bits) mappings onto the corresponding unique permutation code matrices respectively [19]. For mapping tables for larger values of H , we refer the reader to [19]. For M different symbols, T i denotes the set of transmitted code matrices such that T i ∈ { T 1 , . . . , T M }

TABLE I: Mapping of symbols into permutation code matrices ( M = H = 2)

| Label   |   Symbol |   Permutation code matrix |
|---------|----------|---------------------------|
| T 1     |        0 |                        12 |
| T 2     |        1 |                        21 |

TABLE II: Mapping of symbols into permutation code matrices ( M = 4, H = 3)

| Label   |   Symbol |   Permutation code matrix |
|---------|----------|---------------------------|
| T 1     |       00 |                       231 |
| T 2     |       01 |                       213 |
| T 3     |       10 |                       132 |
| T 4     |       11 |                       123 |

has the following general form:

$$\mathbf T _ { i } = \left [ \begin{array} { c c c c } q _ { 1, 1 } & q _ { 1, 2 } & \dots & q _ { 1, H } \\ \colon & \ddots & \dots & \colon \\ q _ { H, 1 } & \dots & \dots & q _ { H, H } \end{array} \right ],$$

where q j,k ∈ { 0 1 , } denotes the ( j, k ) binary element in the i -th transmitted code matrix, j indicates the output of the detector for frequency f j at time step k in the code matrix. At each time step, the element where 1 is present in the code matrix is modulated and transmitted. At a given time, only one frequency band is used for transmission. Using the H -FSK scheme, the transmitted signal, assumed to be sufficiently narrowband, over j -th frequency and k -th time step where j, k ∈ { 1 2 , , . . . , H } , is given by:

$$\dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

where E s = P SU T · T s is the transmitted signal energy per information symbol and T s is the symbol duration. In this model, we assume a strong LOS path between a stationary SU transmitter and a stationary receiver. Therefore, channel noise is modeled as Additive White Gaussian

Noise (AWGN) with zero-mean and variance N 0 2 . Since path loss is a function of the operating frequency, we need to ensure that the received signal power over the different frequency channels is the same, P SU R ( f 1 ) = . . . = P SU R ( f H ) = P SU R , in order to satisfy a fixed signal-to-noise ratio, E /N s 0 at the receiver. Therefore, we adjust each SU's transmitting power over each frequency, P SU T ( f j ) , according to the free space path loss model given in (1). Therefore, the SU transmitter needs to know the location of its corresponding receiver. At the input of the SU demodulator, the received signal is given by:

$$x _ { j } ( t ) = \begin{cases} \ s _ { j } ^ { r } ( t ) + i _ { P U } ^ { r } ( t ) + w ( t ), & \text{if $\exists\,PU$ on $f_{j}$,} \\ \ s _ { j } ^ { r } ( t ) + w ( t ), & \text{otherwise,} \end{cases}$$

where s r j ( ) = t √ 2 E /H r s T s cos(2 πf t j + ) θ , i r PU ( ) = t √ 2 I PU /H T s cos(2 πf t j + φ ) , and E r s glyph[defines] P SU R T s is defined as the symbol energy at the receiver. I PU glyph[defines] P PU I T s and P PU I are defined as the interference energy per coded symbol and the interference power due to the PU transmitter at the SU receiver, respectively, and w t ( ) represents the channel noise at the receiver. It should be mentioned that a PU transmitter can modulate the transmitted signal using any scheme and is not limited to the use of FSK-modulated signals. The fact that the PU power is very high compared to that of the SU implies that the received signal is dominated by the PU signal which, in turn, shapes the behavior of the system. Without loss of generality, we assume that the PU uses an FSK-modulated signal to transmit its information similar to the SU.

At the receiver side, non-coherent detection is employed using a bank of H quadrature receivers so that each consists of two correlation receivers corresponding to the in-phase and quadrature components of the signal. The in-phase component of the signal received, x I j,k , is given by:

$$x _ { I _ { j, k } } & = \\ & \begin{cases} \frac { E _ { s } ^ { r } } { H } \cos \theta + \frac { I _ { P W } } { H } \cos \phi + w, & \text{if PU exists on $f_{j}$ at time $k$,} \\ \frac { E _ { s } ^ { r } } { H } \cos \theta + w, & \text{otherwise,} \end{cases}$$

where θ and φ are uniformly distributed over [0, 2 π ], i.e., θ , φ ∼ U (0 , 2 π ) , and denote the random phase components of the SU and the PU signals, respectively. The noise term w in (6) is modeled as AWGN, i.e., w ∼ N (0 , N 0 2 ) . E r s is the received symbol energy over a given signaling interval.

Similarly, x Q j,k is the received signal's quadrature component defined as follows:

$$& x _ { Q _ { j, k } } = \\ & \left \{ \begin{array} { c c c c c } \frac { E ^ { r } } { H } \sin \theta + \frac { I _ { P W } } { H } \sin \phi + w, & \text{if} \, PU \, \text{exists on} \, f _ { j } \, \text{at} \, \text{time} \, k, \\ & \frac { E ^ { r } } { H } \sin \theta + w, & \text{otherwise}. \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. \end{array} \right \}$$

The envelope of each quadrature receiver over frequency j and time step k , l j,k , is defined as the square root of the sum of the squared in-phase and quadrature components of the correlator output as

$$l _ { j, k } = \sqrt { x _ { I _ { j, k } } ^ { 2 } + x _ { Q _ { j, k } } ^ { 2 } }.$$

At the receiver, a hard decision decoding scheme is used where the envelope value of each of the H quadrature receivers is compared to a threshold value, l th . The threshold value we use in this paper is the same as that used by the authors in [19], namely l th = 0 6 . √ E r s .

The received code matrix R i is of the form,

$$\cdots \dots \dots \dots \dots \dots, \\ \mathbf R _ { i } = \left [ \begin{array} { c c c c } b _ { 1, 1 } & b _ { 1, 2 } & \dots & b _ { 1, H } \\ \colon & \ddots & \dots & \colon \\ b _ { H, 1 } & \dots & \dots & b _ { H, H } \end{array} \right ],$$

where b j,k ∈ { 0 1 , } and can be determined from,

$$b _ { j, k } \ = \ \begin{cases} \ 1, & l _ { j, k } \geq l _ { t h }, \\ 0, & \text{otherwise}. \end{cases}$$

## III. BIT ERROR RATE ANALYSIS

The link quality between two SUs is determined by how much PU interference it can tolerate. We employ BER as the QoS metric to characterize successful communication for an SU. We consider the system model presented in Fig. 2 and obtain an approximation for BER. This system model is the same as the one considered in [19], [20] where the PTC encoder includes a convolutional encoder and the corresponding decoder uses Viterbi algorithm for decoding.

Here, we assume that the PTC encoder uses a rate 1 2 convolutional code whose output is converted to a symbol which, in turn, is mapped onto a permutation matrix to be transmitted. It is important to point out that the exact BER analysis of the proposed system is computationally prohibitive. Hence, we approximate BER using some properties of the Viterbi decoder.

Fig. 3: Finite State Machine of a ( R = 1 2 ) convolutional coded PTC system

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-3e0724ea-8a14-4807-867e-21c1e98e60de.jpg)

Fig. 4: A ( R = 1 2 ) convolutional code trellis for our PTC-based system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-2b78e5dd-d04e-4460-8417-8346f02a9e45.jpg)

The encoder of the binary convolutional code with a given rate is viewed as a finite state machine as shown in Fig. 3. This state diagram, corresponding to the PTC shown in Table II, results in the trellis presented in Fig. 4. In order to optimally decode the received information, we employ the Viterbi algorithm which reconstructs the maximum-likelihood path given the input information sequence. The Viterbi decoder gives as an output the path in the trellis that has the minimum overall Hamming distance with respect to the received sequence. In this regard, an error event may occur if for a transmitted path in the trellis the overall number of differences with the demodulator outputs is larger than or equal to that of a competing path.

In order to compute an upper-bound on the probability of bit error, all error events have to be taken into account. It is important to note that the number of error events is proportional to the number of input bits to be transmitted. In other words, the number of competing paths and the number of error events increase as more input bits are serially loaded into the convolutional encoder. Then, the probability of bit error of a convolutional code is upper-bounded according to [24], [25] as

$$P _ { e } \leq \sum _ { d = d _ { f r e e } } ^ { \infty } a _ { d } P _ { 2 } ( d ),$$

where d free is the free distance of the convolutional code, a d is the number of paths that differ by d bits from the transmitter codeword, and P 2 ( d ) is the probability that the decoded path differs by d bits from the transmitted codeword. In particular, a d and P 2 ( d ) are independent of the transmitted codeword. In fact, they can be calculated assuming that an all-zero codeword is transmitted. Knowing that all the code matrices to be sent include H non-zero elements, we next prove that there is a similar upper-bound on the system given in Fig. 2.

Now, let c be a codeword of the convolutional code and P ( c ) be the corresponding expanded codeword comprising the permutation code matrices mapped from c . The length of P ( c ) is ( L + m H ) 2 , where L denotes the size of packet to be transmitted and L + m the number of branches in the trellis. In order to apply (11), we need to prove that the upper-bound given in (11) is independent of the transmitted expanded codeword. Let P e ( P ( c )) denote the probability of bit error when P ( c ) is the transmitted codeword. Knowing that P ( c ) is a one-to-one mapping from the codeword c , by the argument of (11) given in [24], [25], we have

$$P _ { e } ( \mathcal { P } ( \mathbf c ) ) \quad \leq \sum _ { d = d _ { f r e e } ^ { * } } ^ { \infty } a _ { d } ^ { \mathcal { P } ( \mathbf c ) } P _ { 2 } ^ { \mathcal { P } ( \mathbf c ) } ( d ),$$

where d ∗ free is the free distance of the expanded code, and a P ( c ) d and P P ( c ) 2 ( d ) are the number of paths and the probability that the paths differ by d bits from the transmitted codeword P ( c ) , respectively.

Proposition 1: The bound in (12) is independent of the transmitted codeword if the noise caused by the PU in a subchannel will affect all bits in this channel that are generated by one or more trellis branches.

Proof: See proof in Appendix.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-6cddd62b-4014-4611-85f5-7c34c96ff8c1.jpg)

The fact that we have no information regarding the length of the transmitted bits sequence means that the calculation of the upper bound, in (11), involves an infinite sum over all possible error events. In other words, an infinite number of trellis' paths have to be taken into consideration which will render the process of evaluating the error bound computationally prohibitive for applications such as link adaptation. Therefore, it is important to approximate BER using a suitably selected finite number of trellis' paths.

For every convolutional code, there exists a transfer function that uniquely determines the number of paths by which a received sequence is d bits different from an assumed all-zero transmitted bit sequence. For instance, the transfer function of the mapped permutation code in the proposed H -FSK scheme based on Fig. 3, i.e. H = 3, is given by T D ( ) where

$$T ( D ) = D ^ { 1 6 } + 2 \cdot D ^ { 2 0 } + 4 \cdot D ^ { 2 4 } + 8 \cdot D ^ { 2 8 } + \cdots$$

Equivalently, we are saying that there is only one path, { a → b, b → c, c → a } , with an input sequence { 1 0 0 , , } and permutation code matrix sequence { 123 132 123 , , } that has 16 bits different from the assumed all-zero transmitted bit sequence ( d H (231,123) + d H (231,132) + d H (231,123) = 6 + 4 + 6), 2 paths that have 20 bits difference, and so on. It is most likely that the first few error events (paths) in the trellis dominate all the remaining ones. In this manner, the approximate BER can be written as

$$\hat { P } _ { e } \quad \approx \sum _ { d = d _ { f r e e } ^ { * } } ^ { d _ { f r e e } ^ { * } + z } a _ { d } ^ { \mathcal { P } ( \mathbf c ) } P _ { 2 } ^ { \mathcal { P } ( \mathbf c ) } ( d ),$$

where z ∈ N and z + 1 denotes the number of paths involved in the BER approximation for which numerical results are presented in Section V.

In order to compute P P ( c ) 2 ( d ) we denote by T u and R u the transmitted and received code matrices at stage u in a trellis composed of V stages. Then, P P ( c ) 2 ( d ) can be expressed as,

$$\text{ matrices at stage } u \text{ in a trellis composed of } \mathcal { V } \text{ stages. Then, } P _ { 2 } ^ { \rho } \stackrel { \nu } { \langle } ( d ) \text{ can be expressed as, } \\ P _ { 2 } ^ { \mathcal { P } ( e ) } ( d ) \text{ } & = \prod _ { u = 1 } ^ { \mathcal { V } } P \left ( \mathbf R ^ { u } | \mathbf T ^ { u } \right ) \\ & = \prod _ { u = 1 } ^ { \mathcal { V } } j = 1 \, k = 1 } ^ { \mathbf V } P \left ( b _ { j, k } ^ { u } | q _ { j,$$

where q u j,k and b u j,k are the ( j, k ) binary elements given by (4) and (9) respectively in the u -th stage of the trellis and PU j and ¯ PU j denote the presence and absence of PU activity over channel f j respectively. Here we note that P PU ( j ) and P PU ( ¯ j ) depend on channel f j 's occupancy model which, in turn, is discussed in Section IV-1.

In each stage of the trellis, the absence and presence of PU need to be taken into consideration when computing the likelihoods in (15) that also depend on the code matrix to be transmitted.

## A. Computation of likelihoods in the absence of PU

Suppose there is no PU transmission over frequency band j at time step k . If the SU is active, i.e. transmitting a bit '1', then the demodulator in-phase component x I j,k ∼ N ( √ E r s H cos θ, N 0 2 ) and the demodulator quadrature component x Q j,k ∼ N ( √ E r s H sin θ, N 0 2 ) . The fact that x I j,k and x Q j,k are statistically independent random variables with non-zero means, then implies l j,k ∼ Rice ( √ E r s H , √ N 0 2 ) . Accordingly, P ( b j,k = 0 | q j,k = 1 , P U ¯ j ) can be computed as

$$\cdot \quad \cdot \quad & \cdot \\ P \left ( b _ { j, k } = 0 | q _ { j, k } = 1, \bar { P } U _ { j } \right ) \ & = \ \ P ( l _ { j, k } < l _ { t h } | q _ { j, k } = 1, \bar { P } U _ { j } ) \\ & = \ \ F _ { L _ { j, k } } ( l _ { t h } ) \\ & = \ 1 - Q _ { 1 } \left ( \sqrt { 2 \frac { E ^ { r } s / H } { N _ { 0 } } }, 0. 6 \sqrt { 2 \frac { E ^ { r } } { N _ { 0 } } } \right ), \quad \quad \quad$$

where F L j,k ( l th ) is the cumulative distribution function of l j,k evaluated at l th , Q 1 ( v, w ) is the Marcum's Q-function defined as [26],

$$Q _ { 1 } ( v, w ) = \int _ { w } ^ { \infty } x \exp \left \{ - \frac { x ^ { 2 } + v ^ { 2 } } { 2 } \right \} I _ { 0 } ( v x ) \, d x$$

and I 0 ( vx ) is the zeroth order modified Bessel function. P ( b j,k = 1 | q j,k = 1 , P U ¯ j ) , in this case, is nothing but the complement of (16) and is given by

$$P \left ( b _ { j, k } = 1 | q _ { j, k } = 1, \, \bar { P } \bar { U } _ { j } \right ) \ & = \ \, P ( l _ { j, k } \geq l _ { t h } | q _ { j, k } = 1, \, \bar { P } \bar { U } _ { j } ) \\ & = \ 1 - F _ { L _ { j, k } } ( l _ { t h } ) \\ & = \ Q _ { 1 } \left ( \sqrt { 2 \frac { E _ { s } ^ { r } / H } { N _ { 0 } } }, 0. 6 \sqrt { 2 \frac { E _ { s } ^ { r } } { N _ { 0 } } } \right ). \\.$$

If the SU is not transmitting any information on subchannel j at time step k , then, x I j,k ∼ N (0 , N 0 2 ) and x Q j,k ∼ N (0 , N 0 2 ) , and we have l j,k ∼ Rayleigh ( √ N 0 2 ) . Then, P b ( j,k = 1 | q j,k =

0 , P U ¯ j ) and P b ( j,k = 0 | q j,k = 0 , P U ¯ j ) can be computed as,

$$P ( b _ { j, k } = 1 | q _ { j, k } = 0, \bar { P U } _ { j } ) \ & = \ P ( l _ { j, k } \geq l _ { t h } | q _ { j, k } = 0, \bar { P U } _ { j } ) \\ & = \ 1 - F _ { L _ { j, k } } ( l _ { t h } ) \\ & = \ \exp \left ( \ - 0. 3 6 \frac { E ^ { r } } { N _ { 0 } } \right )$$

and

$$P ( b _ { j, k } = 0 | q _ { j, k } = 0, \bar { P U } _ { j } ) \ & = \ P ( l _ { j, k } < l _ { t h } | q _ { j, k } = 0, \bar { P U } _ { j } ) \\ & = \ F _ { L _ { j, k } } ( l _ { t h } ) \\ & = \ 1 - \exp \left ( - \ 0. 3 6 \frac { E _ { r } ^ { r } } { N _ { 0 } } \right ).$$

## B. Computation of likelihoods in the presence of PU

In the presence of PU activity and an active SU, i.e., the SU is transmitting a '1', the signal received at the input of the demodulator has two significant terms besides noise. One corresponds to the actual signal energy of the transmitted symbol and the other is created by the PU activity. It is assumed that the SU transmitting power is sufficiently small so that QoS of the PU session is maintained. Therefore, the SU signal is negligible as compared to that of the PU. In this case, x I j,k ∼ N ( √ I PU H cos φ, N 0 2 ) and x Q j,k ∼ N ( √ I PU H sin φ, N 0 2 ) . Thus, l j,k ∼ Rice ( √ I PU H , √ N 0 2 ) . P b ( j,k = 0 | q j,k = 1 , P U j ) and P b ( j,k = 1 | q j,k = 1 , P U j ) , in this scenario, are calculated according to,

$$P ( b _ { j, k } = 0 | q _ { j, k } = 1, P U _ { j } ) & \ = \ \ P ( l _ { j, k } < l _ { t h } | q _ { j, k } = 1, P U _ { j } ) \\ & \ = \ \ F _ { L _ { j, k } } ( l _ { t h } ) \\ & \ = \ \ 1 - Q _ { 1 } \Big ( \sqrt { 2 \frac { I _ { P U } / H } { N _ { 0 } } }, 0. 6 \sqrt { 2 \frac { E _ { s } ^ { r } } { N _ { 0 } } } \Big ) \quad \quad ( 2 1 ) \\. \, \text{a}$$

and

$$P ( b _ { j, k } = 1 | q _ { j, k } = 1, P U _ { j } ) \ & = \ \ P ( l _ { j, k } \geq l _ { t h } | q _ { j, k } = 1, P U _ { j } ) \\ & = \ 1 - F _ { L _ { j, k } } ( l _ { t h } ) \\ & = \ Q _ { 1 } \left ( \sqrt { 2 \frac { I _ { P U } / H } { N _ { 0 } } }, 0. 6 \sqrt { 2 \frac { E _ { s } ^ { r } } { N _ { 0 } } } \right ).$$

If the SU does not transmit on subchannel j at time step k , the output of the demodulator has the noisy signal received from the PU. In that case, and similar to the analysis approach followed earlier, x I j,k ∼ N ( √ I PU H cos φ, N 0 2 ) and x Q j,k ∼ N ( √ I PU H sin φ, N 0 2 ) . Consequently, l j,k ∼ Rice ( √ I PU H , √ N 0 2 ) . So, the following probabilities P b ( j,k = 0 | q j,k = 0 , P U j ) and P b ( j,k = 1 | q j,k = 0 , P U j ) are the same as P b ( j,k = 0 | q j,k = 1 , P U j ) and P b ( j,k = 1 | q j,k = 1 , P U j ) , respectively.

Once the probabilities presented above are evaluated, we can calculate P ( R T i ′ | i ) for each stage of the V stages in the trellis and, eventually, approximate the BER of the proposed scheme as given in (14).

## IV. THROUGHPUT ANALYSIS

So far, we have computed the BER approximately and evaluated the performance of the H -FSK system in the presence of PU interference. We have also proposed the adoption of the PTC-based multi-level FSK in the transmission strategy of SU in order to provide resilience to spontaneous narrowband PU interference. In other words, heavy redundancy introduced by the PTC ensures successful communication sessions of the SUs in the presence of PU interference. In order to further explore the performance of the proposed H -FSK communication system and validate its effectiveness, we perform a throughput analysis based on the approximate BER provided in Section III. Accordingly, an SU, adopting the proposed H -FSK communication scheme, is expected to provide a level of reliable information reception which is better than what the SU could achieve otherwise under heavy PU interference.

1) Licensed Channel Dynamic Occupancy Model: Different from [21] where the licensed channel was assumed to be always occupied by the PU, we model the channel available for PU active transmissions and its occupancy as a 2-state Markov chain as shown in Fig. 5. In this model, the channel has two alternative states denoted by On (Busy) and Off (Idle). This assumption is practical in the sense that the licensed spectrum occupancy will experience alternating On-Off periods rather than the assumption of the channel being always occupied (PU being always present once it starts transmission during an SU transmission). The periods during which the licensed spectrum's state is idle or busy, also known as holding times, are geometrically distributed with known independent parameters r and p as shown in Fig. 5. An On state represents the state in which the licensed band is occupied by a PU resulting in degraded

Fig. 5: Two-state Markov chain model for PU activity over licensed spectrum

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-435db259-48cd-4667-91ef-5a627849b615.jpg)

BER performance for SU transmissions, while an Off (idle) state represents the state in which the primary band is idle. Accordingly, the steady state probabilities of the primary channel's state being Off or On are given respectively as

$$P _ { c h } ^ { O f f } = \frac { r } { r + p }$$

and

$$P _ { c h } ^ { Q n } = \frac { p } { r + p }.$$

2) Throughput Analysis: In order to examine the effectiveness of the proposed H -FSK approach, we conduct a throughput analysis that monitors the average throughput of the SU communication session in the presence of PU activity modeled as a 2-state Markov chain. Knowing that the SU transmitter has no information regarding the presence or absence of the PU, we determine the unlicensed throughput for SU communications as a function of the probability of a licensed channel being occupied by a PU. Such analysis illustrates the efficiency of the proposed scheme in ensuring a stable throughput regardless of any PU activity over the licensed spectrum. Accordingly, we present the approximate expected throughput ( T e ) of the H -FSK communication scheme based on the approximate BER given as

$$T _ { e } \approx \left ( 1 - P \hat { E } R \right ) \times R _ { p } \times L,$$

where ˆ PER = 1 -(1 -ˆ ) P e L is the approximate packet error rate (PER), ˆ P e is the approximate BER computed with z +1 paths knowing the presence of dynamic PU interference as described in (14) and (15), and R p is the rate at which information packets are sent. Throughput, in our

model, is defined as the number of correctly received packets per unit time. For simplicity, we do not consider header and overhead bits when computing the link throughput.

## V. NUMERICAL RESULTS AND DISCUSSION

## A. Approximate BER with Single Permanent PU Interference

For the coded H -FSK communication system shown in Fig. 2, the following values of the physical parameters are assumed. The transmit power of the SU is varied between [25 µW, 4 mW ] and the noise power density is selected as N 0 = 2 5 . × 10 -14 . These values are chosen in such a way so as to not create destructive interference to the PU. We further assume that the lowest frequency band in the available frequency spectrum, f 1 , is selected as 56 MHz, and the bandwidth spacing between any two subchannels is selected as 6 MHz. The PU is assumed to be operating at f 2 and has a total transmitting power of 1 MW. The parameters' values used in the simulation are chosen in accordance with the IEEE 802.22 standard [27]. The distances between the PU transmitter and SU receiver and that between the SU transmitter-receiver pair are selected as 10 m, in order to analyze the quality of the SU link under a high interference scenario.

In Fig. 6, we fix H = 3 and compare the simulated BER with its approximations using different number of paths used in the approximation, i.e., z +1 . It is clear from the figure that z = 3 approximates BER well for H = 3 . In Fig. 7, we fix SNR = 7 dB and compare absolute error of the BER approximation for different values of H . It is again clear that z = 3 , i.e., using only z +1 = 4 paths of the convolutional code's trellis, provides good BER approximation for different sizes of PTC. This is due to the fact that the first few error events with a few errors in the trellis dominate all the remaining transitions. This result confirms what we have mentioned earlier regarding the domination of a small number of paths in the trellis over all the remaining paths that may evolve due to incoming bits. The approximate BER we provide in the section could be used as a QoS metric in an adaptive scenario where certain parameters of the system such as H could be adapted in an online manner using this approximation. In Table III, we present results on the performance of exhaustive search 1 compared to the proposed

1 Using this technique, decisions on the transmitted code matrices are made by the permutation trellis decoder which decides in favor of the transmitted code matrix, i.e., among 2 m possible code matrices, which has the minimum Hamming distance with respect to the received one. In order to determine all the decisions for every possible code matrix, a brute force method is used to compare each possible received code matrix to every code matrix. This method needs roughly 2 H 2 + m comparisons [21].

Fig. 6: BER approximations for PTC with H = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-206bb307-8b38-43ee-a31e-7a492f4b460b.jpg)

Fig. 7: | BER -BER apx | for different PTC ( H = 2, 3, and 4) and SNR = 7 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-862acf66-d360-406f-8114-5acee9aed0b1.jpg)

approximation approach in this paper in terms of the time needed to compute BER for different PTC. As shown in Table III, it takes much shorter time to compute the BER approximately as compared to its exact calculation.

TABLE III: A comparison of the time needed to compute BER between exhaustive search [21] and approximation approach for different PTC ( H = 2,3, and 4)

| PTC   | Exhaustive Search Approach   | Approximate Approach   |
|-------|------------------------------|------------------------|
| H = 2 | 0 . 33 s                     | 0 . 31 s               |
| H = 3 | 42 . 21 s                    | 7 . 24 s               |
| H = 4 | 26274 . 99 s                 | 245 . 32 s             |

## B. Throughput Analysis of H -FSK

For comparison purposes, we first consider a competing scenario where an SU does not employ any coding technique in its spectrum interweaved transmission strategy which, in turn, is assumed to be based on M -FSK modulated signals. In this scenario, an SU makes use of the CR channel sensing feature in order to determine if there is an active PU in the network or not. Depending on the channel sensing measurements of the licensed spectrum, a decision is made on whether the SU should vacate the band or not. In the case where the SU detects a white space, it adapts the transmission parameters to utilize the full spectrum available and achieve the maximum throughput. On the other hand, when a PU is detected an SU vacates the band and adjusts its transmission parameters accordingly by employing a modified FSK modulation scheme supporting the transmission of M ′ symbols where M &lt; M ′ . In order to gain more insight regarding the throughput performance comparison of the PTC-based multi-level FSK with the uncoded-opportunistic (adaptive) M -FSK system, we assume that packet transmissions take place over an observation window of length T seconds divided into x time slots. For each time slot, the steady state probability of a channel being occupied by an active PU, given in (24), is evaluated. Knowing each interval's time and the time required for transmitting a single packet, the incorrectly received packets can be determined. As a result, the unlicensed network's throughput is computed and analyzed. We further assume that the available spectrum consists of four frequency bands. With the use of the PTC-based multi-level FSK approach, an SU does not need to sense the channel or adapt its wireless links. In the uncoded opportunistic system, an SU senses the licensed frequency band at the beginning of a time slot. Assuming zerocost perfect detection of a licensed user's activity, an SU decides on which M -FSK modulated

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-9e937ddf-537c-4056-9291-6eefc3799203.jpg)

P Ch

Fig. 8: Throughput Comparison of H -FSK and uncoded opportunistic M -FSK as well as convolutionally coded BPSK schemes for SNR = 4 dB, packet size = 256 bits, and R p = 100 packets/sec.

signals are to be transmitted. In the case where an SU detects a spectrum hole, it utilizes the full available spectrum (4 bands) using a 4-FSK modulation scheme. When a PU is detected, an SU vacates the licensed band and utilizes the remaining spectrum (2 frequency bands out of the 4 bands assumed to comprise the spectrum) using BFSK modulated signals. Different from [28], the licensed channel can change its state at any instant in a time slot. In Fig. 8, we present the average throughput of H -FSK and that of an uncoded opportunistic M -FSK system as a function of P On ch for H = M = 4 where the packet size L = 256 bits, R p = 100 packets/sec, and SNR is 4 dB. Fig. 8 illustrates that the average throughput in H -FSK is constant regardless of P On ch values while the average throughput of the uncoded opportunistic M -FSK decreases as P On ch increases and it matches the theoretical approximation given in (25). As shown, there exists a probability p ∗ 1 ≈ 0 7 . such that for P On ch ≥ p ∗ 1 the H -FSK approach outperforms the uncoded opportunistic scheme in terms of throughput performance besides BER. Thus, the use of the proposed H -FSK communication scheme for values of P On ch less than p ∗ 1 is not beneficial because higher throughput values are obtained using the uncoded opportunistic scheme. In other

words, the H -FSK system is effective for situations where the PU activity is relatively high. In practice, the throughput evaluated for the uncoded M -FSK system is expected to be lower than what is shown in Fig. 8 due to channel sensing intervals. It should be noted that the proposed H -FSK communication scheme does not need the channel sensing mechanism that might not be available for some applications.

Along these lines, we compare the average throughput of the proposed H -FSK scheme to a rate1 2 convolutionally coded BPSK modulated OFDM system. We assume that the latter transmits a BPSK modulated signal in parallel over the available frequency bands in a given time window similar to what is considered in the previous comparison setting. For fair comparison, neither scheme carries out sensing, i.e., they transmit at all times. It is worth noting that the throughput performance of this BPSK OFDM system is investigated assuming the same transmission bandwidth, power, and time frame as with the H -FSK system. As shown in Fig. 8, the average throughput of the convolutionally coded BPSK system also decreases as P On ch increases. Similar to the previous throughput comparison, there exists a probability p ∗ 2 ≈ 0 55 . such that for P On ch ≥ p ∗ 2 the H -FSK system outperforms the convolutionally coded BPSK scheme in terms of average SU's throughput. Therefore, the use of the proposed H -FSK communication scheme for values of P On ch less than p ∗ 2 is not beneficial because higher throughput values would be obtained using the other coded scheme.

## C. SU Link's Resiliency in the Presence of Multiple PUs

In this section, we further carry out our performance analysis and examine the resiliency of the network (SU's link) to interference provided by the proposed PTC based framework in the presence of multiple PUs when i) PUs always transmit and ii) PUs are dynamic and their activities vary dynamically. It is worth noting here that the PUs are located 10 m away from the SU receiver in order to analyze the quality of the SU link under a high interference scenario.

1) Static Channel Occupancy: First, we consider multiple PUs that are always active in the network during SU communication. In Fig. 9, we plot the approximate BER performance of the SU link for different PTCs in the presence of multiple PUs. The results obtained in the presence of multiple PUs that are always ON are similar to those presented in [19], [20]. These results show that the use of larger PTC (higher value of H ) adds more robustness to SU communications, i.e., BER decreases as the value of H increases. It is also shown in these results that the approximate

Fig. 9: Approximate BER performance of different PTCs with respect to PUs that are always ON and transmitting.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-d37e6f65-bc4a-4b93-9db5-d697c71be264.jpg)

Fig. 10: Approximate BER versus Simulated BER for a PTC with H = 2 , 3 , and 4 and P ON ch = 0 35 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-6f2cc69e-68d5-4dc3-a104-21c364e419b1.jpg)

BER performance of an SU link degrades as the number of PUs joining the network increases. Of course, this is achieved at the expense of larger overhead when larger values of H are used.

2) Dynamic Channel Occupancy: The case of static PUs that always transmit is more of a pessimistic scenario as presented above. In this paper, we extend the work in [19], [20] and consider a more practical scenario to model the intermittent dynamic activities of PUs in the network. Using a 2-state Markov chain to model channels' occupancy by PUs, we provide a more

Fig. 11: Approximate BER performance of a PTC ( H = 4) with respect to dynamically varying PUs activities and PUs that are always ON.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-ff1ae0ae-dbf7-4c2d-b9fc-9949993cc92d.jpg)

accurate characterization of system performance and verify the approximate BER in the presence of a dynamically transmitting PU with simulated BER for a given PTC and known parameters of the channel occupancy model. In Fig. 10, we present the approximate BER performance of the proposed framework as a function of E s N 0 for different PTC schemes, i.e., H = 2 3 , , and 4 , and P ON ch = 0 35 . . Both analytical and simulation results are presented which match each other quite well.

In addition, in Fig. 11, we further examine the approximate BER performance of SU communications for a given PTC, i.e., H = 4 , in the presence of one, two, and three PUs that i) are static in nature and always transmit and ii) dynamic in nature. The results obtained in Fig. 11 clearly show that there exists a gap between the performance of the systems in both cases. Our analysis based on a more practical occupancy model is able to predict SU communication performance more accurately and which is better than the pessimistic case of PUs being always ON. For example, the approximate BER value in the presence of one dynamic PU is lower than when a PU is always transmitting. This is intuitive because only a fraction of the transmitted information packets are affected by interference from the dynamic PU when

Fig. 12: A BER comparison between the proposed PTC scheme with H = 4 and a 1 / 2 rate ( 64800 32400 , ) LDPC encoding with a 4 -FSK modulation system with respect to static ( P On ch = 1) and dynamic ( P On ch = 0.75) activities of a PU.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0259v2-figure-b0616ff2-0082-4228-bada-e4aa34c9aae6.jpg)

it is transmitting compared to the static PU case where all the transmitted information packets experience permanent interference. In fact, this is also observed for the case of 2 and 3 PUs being active in the network.

Moreover, we compare the BER performance of the proposed communication system when H = 4 with another system that is based on a 1 / 2 rate ( 64800 32400 , ) LDPC encoding given in Matlab toolbox with a 4 -FSK modulation system. Fig. 12 illustrates how the proposed communication system outperforms the aforementioned system for both static and dynamic activities of a PU.

## VI. CONCLUSION

In this paper, we employed a PTC based framework to mitigate the impact of PUs modeled using a practical dynamic channel occupancy model in CRNs. We, further, computed the SU link's BER approximately which was shown to be quite accurate. This approximation allows one to use BER as a QoS metric to determine the link quality of an SU link for applications such as link adaptation. Furthermore, in order to assess the effectiveness of the proposed PTC based

multi-level FSK communication scheme, we compared the performance to that of an uncoded opportunistic M -FSK system, a coded M -FSK system, and a coded BPSK modulated system deploying parallel transmissions. Based on the comparisons described in Section V, we showed that the proposed scheme outperforms the latter two under relatively heavy PU interference. We also presented results that exhibit the resiliency of an SU link to interference for PTCs in the presence of multiple dynamic PUs activities.

In this work, we evaluated the performance of PTC codes for a single SU link. Optimal code assignments in a multiple SU scenario will be considered as a future research direction. We also plan to consider 1) more realistic channel models that account for shadowing and fading effects; and 2) more effective decoding techniques, e.g., soft-decision decoding. Based on the promising results we have obtained in a CRN, we intend to generalize the proposed scheme to tackle wireless networks subject to jamming attacks and malicious nodes' activities.

## APPENDIX

## PROOF OF PROPOSITION 1

Let P ( c 1 ) be the transmitted codeword and R ( ) i be an arbitrary received vector. Assume that P ( c 2 ) is another possible codeword. If there exists a received vector R ( j ′ ) such that d H ( R ( j ′ ) , P ( c 2 )) = d H ( R ( ) i , P ( c 1 )) and P ( R ( j ′ ) |P ( c 2 )) = P ( R ( ) i |P ( c 1 )) , then the upper-bound given in (12) is the same no matter what the transmitted codeword is ( P ( c 2 ) or P ( c 1 ) ). Let s 1 = ( s 1 ,glyph[lscript]H 2 , s 1 ,glyph[lscript]H 2 +1 , . . . , s 1 ( , glyph[lscript] +1) H 2 -1 ) be a portion of P ( c 1 ) generated from the corresponding ( glyph[lscript] +1) th output branch of the convolutional code and s 2 be that of P ( c 2 ) for the same bit positions. Since the number of non-zero elements in s 1 and s 2 is H , we can assume that they occupy the following bit indices o 1 0 , , o 1 1 , , . . . , o 1 ,H -1 in s 1 and o 2 0 , , o 2 1 , , . . . , o 2 ,H -1 in s 2 , respectively. Let r 1 = ( r 1 ,glyph[lscript]H 2 , r 1 ,glyph[lscript]H 2 +1 , . . . , r 1 ( , glyph[lscript] +1) H 2 -1 ) be the corresponding received vector of s 1 . Note that r 1 here is in the form of an arbitrary received code matrix for each output branch of the convolutional code and is mapped to a permutation trellis code matrix. Then, R ( j ′ ) can be constructed from R ( ) i as provided below.

Define a one-to-one mapping π glyph[lscript] from indices { glyph[lscript]H 2 , glyph[lscript]H 2 +1 , . . . , ( glyph[lscript] +1) H 2 - } 1 to itself as

$$\mu _ { \ell } ( q ) = \begin{cases} \lambda _ { \ell } \\ c \end{cases}$$

$$\pi _ { \ell } ( q ) = \begin{cases} \ o _ { 2, i }, & \text{if $q=0_{1,i}$ for $0\leq i\leq H-1$,} \\ \ o _ { 1, i }, & \text{if $q=0_{2,i}$ for $0\leq i\leq H-1$,} \\ \ q, & \text{otherwise.} \end{cases}$$

It is easy to see that the above mapping rule switches s 1 to s 2 and vice versa. We can write π ( P ( c 1 )) = P ( c 2 ) and π ( P ( c 2 )) = P ( c 1 ) . At the same time, the mapping rule also switches r 1 ,o 1 ,i with r 1 ,o 2 ,i for 0 ≤ i ≤ H -1 when applied to r 1 . Then, the combination of glyph[lscript] mappings in π with respect to R i , results in

$$\mathbf R ^ { ( j ^ { \prime } ) } = \pi ( \mathbf R ^ { ( i ) } ) = ( r _ { 1, \pi _ { 0 } ( 0 ) }, r _ { 1, \pi _ { 0 } ( 1 ) }, \dots, r _ { 1, \pi _ { 0 } ( H ^ { 2 } - 1 ) }, r _ { 1, \pi _ { 1 } ( H ^ { 2 } ) }, \dots, r _ { 1, \pi _ { L + m - 1 } ( ( L + m - 1 ) H ^ { 2 } - 1 ) } ).$$

## Accordingly,

$$d _ { H } ( \mathbf R ^ { ( i ) }, \mathcal { P } ( \mathbf c _ { 1 } ) ) = d _ { H } ( \pi ( \mathbf R ^ { ( i ) } ) ), \pi ( \mathcal { P } ( \mathbf c _ { 1 } ) ) ) = d _ { H } ( \mathbf R ^ { ( j ^ { \prime } ) } ), \mathcal { P } ( \mathbf c _ { 2 } ) )$$

for the one-to-one mapping conserves the Hamming distance of the convolutional code. In the case where noise for each received bit is independent, we have

$$P ( { \mathbf R } ^ { ( i ) } | { \mathcal { P } } ( { \mathbf c } _ { 1 } ) ) = P ( \pi ( { \mathbf R } ^ { ( i ) } ) | \pi ( { \mathcal { P } } ( { \mathbf c } _ { 1 } ) ) ) = P ( { \mathbf R } ^ { ( j ^ { \prime } ) } | { \mathcal { P } } ( { \mathbf c } _ { 2 } ) )$$

and since noise of the received bits in a subchannel occupied by the PU is not independent, P ( R ( ) i |P ( c 1 )) is not always equal to P ( R ( j ′ ) |P ( c 2 )) . However, if we assume that the noise caused by the PU in a subchannel will affect all bits in this channel that are generated by one or more trellis branches, then

$$P ( { \mathbf R } ^ { ( i ) } | { \mathcal { P } } ( { \mathbf c } _ { 1 } ) ) = P ( \pi ( { \mathbf R } ^ { ( i ) } ) | \pi ( { \mathcal { P } } ( { \mathbf c } _ { 1 } ) ) ) = P ( { \mathbf R } ^ { ( j ^ { \prime } ) } | { \mathcal { P } } ( { \mathbf c } _ { 2 } ) ).$$

We can reason this from the fact that bit-switching by π on R ( ) i does not move any bit out of the row it is located in.

## ACKNOWLEDGMENT

This material is based on research sponsored in part by CASE: The Center for Advanced Systems and Engineering, an NYSTAR center for advanced technology at Syracuse University. Part of Han's work is supported by National Science Council, Taiwan, R.O.C., under grants NSC 101-2221-E-011 -069 -MY3 and NSC 99-2221-E-011 -158 -MY3.

## REFERENCES

- [1] I. F. Akyildiz, W.-Y. Lee, M. C. Vuran, and S. Mohanty, 'Next generation/dynamic spectrum access/cognitive radio wireless networks: A survey,' Computer Networks , vol. 50, no. 13, pp. 2127-2159, 2006.
- [2] Q. Zhao and A. Swami, 'A decision-theoretic framework for opportunistic spectrum access,' IEEE Trans. Wireless Comm. , vol. 14, no. 4, pp. 14-20, Aug. 2007.

- [3] S. Haykin, 'Cognitive radio: brain-empowered wireless communications,' IEEE J. Sel. Areas in Comm. , vol. 23, no. 2, pp. 201-220, Feb. 2005.
- [4] J. Mitola and G. Q. Maguire, 'Cognitive radio: making software radios more personal,' IEEE Personal Comm. , vol. 6, no. 4, pp. 13-18, Aug. 1999.
- [5] Y.-S. Kwon, H.-S. Kim, J.-H. Yoo, and J.-H. Chung, 'A non-interfering cognitive radio system for spectrum sharing,' in Proc. Int. Symp. on Comm. and Information Technology , Sep. 2009, pp. 664-665.
- [6] Y. Wu and Z. Yang, 'Spectrum sharing under interference temperature and SINR constraints,' in Proc. Int. Conf. on Wireless Comm., Networking and Mobile Computing , Oct. 2008, pp. 1-4.
- [7] Federal Communications Commission, 'Establishment of interference temperature metric to quantify and manage interference and to expand available unlicensed operation in certain fixed mobile and satellite frequency bands,' Et Docket 03-289, Notice of Inquiry and Proposed Rulemaking , 2003.
- [8] T. C. Clancy, 'Formalizing the interference temperature model,' Wiley J. on Wireless Comm. and Mobile Computing , vol. 7, no. 9, pp. 1077-1086, Nov. 2007.
- [9] H. Mahmoud, T. Yucek, and H. Arslan, 'OFDM for cognitive radio: merits and challenges,' IEEE Trans. Wireless Comm. , vol. 16, no. 2, pp. 6-15, Apr. 2009.
- [10] F. Shayegh and M. R. Soleymani, 'Rateless codes for cognitive radio in a virtual unlicensed spectrum,' in Proc. IEEE Sarnoff Symp. , May 2011, pp. 1-5.
- [11] G. Yue and X. Wang, 'Anti-jamming coding techniques with application to cognitive radio,' IEEE Trans. Wireless Comm. , vol. 8, no. 12, pp. 5996-6007, Dec. 2009.
- [12] J. Ko and C. Kim, 'Communication method between spectrum heterogeneous secondary users in OFDM-based cognitive radio,' Electronics Letters , vol. 47, no. 14, pp. 827-829, Jul. 2011.
- [13] C. Nicola, H. Mercier, and V. K. Bhargava, 'Error-correcting codes for dynamic spectrum allocation in cognitive radio systems,' in Proc. Int. Symp. on Signals, Syst. and Electron. , Aug. 2007, pp. 247-250.
- [14] S. Wang, Y. Wang, J. P. Coon, and A. Doufexi, 'Energy-efficient spectrum sensing and access for cognitive radio networks,' IEEE Trans. Veh. Technol. , vol. 61, no. 2, pp. 906-912, Feb. 2012.
- [15] G. Ozcan and M. C. Gursoy, 'Throughput of cognitive radio systems with finite blocklength codes,' in Proc. IEEE Annual Conf. on Inf. Sciences and Systems (CISS) , Mar. 2012, pp. 1-6.
- [16] Y.-C. Liang, Y. Zeng, E. Peh, and A. T. Hoang, 'Sensing-throughput tradeoff for cognitive radio networks,' IEEE Trans. Wireless Comm. , vol. 7, no. 4, pp. 1326-1337, Apr. 2008.
- [17] B. Shahrasbi and N. Rahnavard, 'Rateless-coding-based cooperative cognitive radio networks: Design and analysis,' in Proc. Annual IEEE on Comm. Society Conf. on Sensor, Mesh and Ad Hoc Comm. and Networking (SECON) , June 2011, pp. 224-232.
- [18] W.-C. Ao and K.-C. Chen, 'End-to-end HARQ in cognitive radio networks,' in Proc. IEEE Wireless Comm. and Networking Conf. (WCNC) , Apr. 2010, pp. 1-6.
- [19] H. C. Ferreira, A. J. H. Vinck, T. G. Swart, and I. De Beer, 'Permutation trellis codes,' IEEE Trans. Communications , vol. 53, no. 11, pp. 1782-1789, Nov. 2005.
- [20] T. G. Swart, I. De Beer, H. C. Ferreira, and A. J. H. Vinck, 'Simulation results for permutation trellis codes using M-ary FSK,' in Proc. Int. Symp. on Power Line Comm. and its Applications , Apr. 2005, pp. 317-321.
- [21] R. El Bardan, E. Masazade, O. Ozdemir, and P. K. Varshney, 'Performance of permutation trellis codes in cognitive radio networks,' in Proc. IEEE Sarnoff Symp. , May 2012, pp. 1-6.

- [22] A. Goldsmith, Wireless Communications . New York, NY, USA: Cambridge University Press, 2005.
- [23] M. Elalem, L. Zhao, and Z. Liao, 'Interference mitigation using power control in cognitive radio networks,' in Proc. IEEE Vehicular Tech. Conf. , May 2010.
- [24] J. G. Proakis, Digital Communications . McGraw -Hill, 2000.
- [25] S. Lin and D. J. Costello, Jr., Error Control Coding: Fundamentals and Applications . Prentice-Hall, 1983.
- [26] M. K. Simon, Probability Distributions Involving Gaussian Random Variables: A Handbook for Engineers, Scientists and Mathematicians . Springer-Verlag New York, 2006.
- [27] C. Stevenson, G. Chouinard, Z. Lei, W. Hu, S. J. Shellhammer, and W. Caldwell, 'IEEE 802.22: The first cognitive radio wireless regional area network standard,' IEEE Comm. Mag. , vol. 47, no. 1, pp. 130-138, Jan. 2009.
- [28] Q. Zhao, B. Krishnamachari, and K. Liu, 'On myopic sensing for multi-channel opportunistic access: structure, optimality, and performance,' IEEE Trans. Wireless Comm. , vol. 7, no. 12, pp. 5431-5440, Dec. 2008.
Raghed El-Bardan Engin Masazade Onur Ozdemir Yunghsiang S. Han Pramod K. Varshney
Permutation Trellis Coded Multi-level FSK Signaling to Mitigate Primary User Interference in Cognitive Radio Networks
We employ Permutation Trellis Code (PTC) based multi-level Frequency Shift
Keying signaling to mitigate the impact of Primary Users (PUs) on the
performance of Secondary Users (SUs) in Cognitive Radio Networks (CRNs). The
PUs are assumed to be dynamic in that they appear intermittently and stay
active for an unknown duration. Our approach is based on the use of PTC
combined with multi-level FSK modulation so that an SU can improve its data
rate by increasing its transmission bandwidth while operating at low power and
not creating destructive interference for PUs. We evaluate system performance
by obtaining an approximation for the actual Bit Error Rate (BER) using
properties of the Viterbi decoder and carry out a thorough performance analysis
in terms of BER and throughput. The results show that the proposed coded system
achieves i) robustness by ensuring that SUs have stable throughput in the
presence of heavy PU interference and ii) improved resiliency of SU links to
interference in the presence of multiple dynamic PUs.