<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.0 20120330//EN" "http://jats.nlm.nih.gov/publishing/1.0/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.0" specific-use="sps-1.6" xml:lang="en" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
	<front>
		<journal-meta>
			<journal-id journal-id-type="publisher-id">rium</journal-id>
			<journal-title-group>
				<journal-title>Revista Ingenierías Universidad de Medellín</journal-title>
				<abbrev-journal-title abbrev-type="publisher">Rev. ing. univ. Medellín</abbrev-journal-title>
			</journal-title-group>
			<issn pub-type="ppub">1692-3324</issn>
			<publisher>
				<publisher-name>Universidad de Medellín</publisher-name>
			</publisher>
		</journal-meta>
		<article-meta>
			<article-id pub-id-type="doi">10.22395/rium.v15n29a8</article-id>
			<article-categories>
				<subj-group subj-group-type="heading">
					<subject>Articles</subject>
				</subj-group>
			</article-categories>
			<title-group>
				<article-title>Spectrum Sensing Framework based on Blind Source Separation for Cognitive Radio Environments</article-title>
				<trans-title-group xml:lang="es">
					<trans-title>Método para sensado del espectro basado en separación ciega de fuentes para ambientes de radio cognitiva</trans-title>
				</trans-title-group>
			</title-group>
			<contrib-group>
				<contrib contrib-type="author">
					<name>
						<surname>Sepúlveda-Cano</surname>
						<given-names>Lina María</given-names>
					</name>
					<xref ref-type="aff" rid="aff1">**</xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Quiza-Montealegre</surname>
						<given-names>Jhon Jair</given-names>
					</name>
					<xref ref-type="aff" rid="aff2">***</xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Gil-Taborda</surname>
						<given-names>Camilo</given-names>
					</name>
					<xref ref-type="aff" rid="aff3">****</xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Gómez</surname>
						<given-names>Jorge Andrés</given-names>
					</name>
					<xref ref-type="aff" rid="aff4">*****</xref>
				</contrib>
			</contrib-group>
			<aff id="aff1">
				<institution content-type="original"> ** Electronic Engineer. MSc.In Telecommunications Engineering, PhD in Engineering - Automatic. Researcher of the ArkadiusGroup , University of Medellín. Carrera 87 # 30-65, Medellín, Colombia. E-mail: lmsepulveda@ udem.edu.co </institution>
				<institution content-type="normalized">Universidad de Medellín</institution>
				<institution content-type="orgname">University of Medellín</institution>
				<country country="CO">Colombia</country>
				<email>lmsepulveda@ udem.edu.co</email>
			</aff>
			<aff id="aff2">
				<institution content-type="original">*** Electronic Engineer. Specialist in Teleinformatics.MSc.In Telecommunications Engineering. Researcher of the Arkadius Group, University of Medellín .Carrera 87 # 30-65, Medellín, Colombia. E-mail: jhquiza@udem. edu.co </institution>
				<institution content-type="normalized">Universidad de Medellín</institution>
				<institution content-type="orgname">University of Medellín</institution>
				<country country="CO">Colombia</country>
				<email>jhquiza@udem. edu.co</email>
			</aff>
			<aff id="aff3">
				<institution content-type="original">**** Electronic Engineer .MSc. en Multimedia and Communications, UC3M. PhD en Multimedia and Communica tions, UC3M., Medellín, Colombia. E-mail: cgil@udem.edu.co </institution>
				<institution content-type="orgname">Multimedia and Communications, UC3M</institution>
				<country country="CO">Colombia</country>
				<email>cgil@udem.edu.co</email>
			</aff>
			<aff id="aff4">
				<institution content-type="original">***** Electronic Engineer. MSc.In Engineering - Automation.Researcher of the group Bioengineering and Opto electronics (ByO). Polytechnic University of Madrid . Madrid, España. E-mail: jorge.gomez.garcia@upm.es</institution>
				<institution content-type="orgname">Polytechnic University of Madrid</institution>
				<country country="ES">España</country>
				<email>jorge.gomez.garcia@upm.es</email>
			</aff>
			<pub-date pub-type="epub-ppub">
				<season>Jun-Dec</season>
				<year>2016</year>
			</pub-date>
			<volume>15</volume>
			<issue>29</issue>
			<fpage>129</fpage>
			<lpage>140</lpage>
			<history>
				<date date-type="received">
					<day>14</day>
					<month>12</month>
					<year>2015</year>
				</date>
				<date date-type="accepted">
					<day>02</day>
					<month>05</month>
					<year>2016</year>
				</date>
			</history>
			<permissions>
				<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by-nc/4.0/" xml:lang="en">
					<license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License</license-p>
				</license>
			</permissions>
			<abstract>
				<title><italic>Abstract</italic> </title>
				<p>The efficient use of spectrum has become an active research area due to its scarcity and underutilization. In a spectrum sharing scenario as cognitive radio (CR), the vacancy of licensed frequency bands could be detected by a secondary user through spectrum sensing techniques. Usually, this sensing approaches are performed with prior knowledge of the channel features. In the present work, a blind spectrum sensing approach based on independent component analysis (ICA) and singular spectrum analy sis (SSA) is proposed. The approach is tested and compared with other outcomes. Results show that the proposed scheme is capable of detecting most of the sources with low time consumption, which is a remarkable aspect for online applications with demanding time issues. </p>
			</abstract>
			<trans-abstract xml:lang="es">
				<title><italic>Resumen</italic> </title>
				<p>El uso eficiente del espectro se ha convertido en un área de investigación activa, debido a la escasez de este recurso y a su subutilización. En un escenario en el que el espectro es un recurso compartido como en la radio cognitiva (CR), los espacios sin uso dentro de las bandas de frecuencias con licencia podrían ser detectados y posteriormente utilizados por un usuario secundario a través de técnicas de detección y sensado del espe ctro. Generalmente, estas técnicas de detección se utilizan a partir de un conocimiento previo de las características de canal. En el presente trabajo se propone un enfoque de detección ciega del espectro basado en análisis de componentes independientes (ICA) y análisis de espectro singular (SSA). La técnica de detección se valida a través de simulación, y su des empeño se compara con metodologías propuestas por otros autores en la literatura. Los resultados muestran que el sistema propuesto es capaz de detectar la mayoría de las fuentes con bajo consumo de tiempo, un aspecto que cabe resaltar para aplicaciones en línea con exigencias de tiempo. </p>
			</trans-abstract>
			<kwd-group xml:lang="en">
				<title>Key words:</title>
				<kwd>spectrum sensing</kwd>
				<kwd>blind source separation</kwd>
				<kwd>cognitive radio</kwd>
				<kwd>ica</kwd>
				<kwd>ssa</kwd>
			</kwd-group>
			<kwd-group xml:lang="es">
				<title>Palabras clave:</title>
				<kwd>sensado del espectro</kwd>
				<kwd>separación ciega de fuentes</kwd>
				<kwd>radio cognitiva</kwd>
				<kwd>ica</kwd>
				<kwd>ssa</kwd>
			</kwd-group>
			<counts>
				<fig-count count="0"/>
				<table-count count="1"/>
				<equation-count count="13"/>
				<ref-count count="18"/>
				<page-count count="12"/>
			</counts>
		</article-meta>
	</front>
	<body>
		<sec sec-type="intro">
			<title>INTRODUCTION</title>
			<p>Spectrum sensing is one of the key elements in the cognitive radio (cr) communication, and it is the first step that must be followed before allowing unlicensed users to access vacant licensed channels <xref ref-type="bibr" rid="B1"><sup>1</sup></xref>. In a cr environment, frequency bands that are legally assigned to primary users are exploited by secondary users when primary users are idle <xref ref-type="bibr" rid="B2"><sup>2</sup></xref>; nevertheless, if the primary user needs to establish the communication again, the secondary user has to vacate the channel and seek another one. Detecting the presence of primary signal (spectrum sensing) is a challenging task in cr <xref ref-type="bibr" rid="B3"><sup>3</sup></xref>. Several primary users will perform different modulation schemes, data rates and transmission powers in the presence of both variable propagation environment and interference generated by other secondary users <xref ref-type="bibr" rid="B4"><sup>4</sup></xref>. The spectrum sensing techniques most commonly used are: Energy Detection (ed) <xref ref-type="bibr" rid="B5"><sup>5</sup></xref>, Cyclostationary Detection <xref ref-type="bibr" rid="B6"><sup>6</sup></xref> and Matched Filtering <xref ref-type="bibr" rid="B7"><sup>7</sup></xref>. Other techniques proposed to improve the poor performance of ed in noisy environ ments are the Maximum to Minimum Eigenvalue method (mme) and the Covariance method (cov) <xref ref-type="bibr" rid="B3"><sup>3</sup></xref>. All of these techniques require some previous knowledge that is not always available, regarding the received signals, e.g. the bandwidth, frequency range, modulation scheme, power, among others. Additionally, any sensing spectrum technique should be fast and effective, attempting to reduce processing time and misdetection. </p>
			<p>This paper presents a technique for spectrum sensing based on the independent component analysis (ica) for blind source separation. The proposed scheme includes a single-to-multichannel conversion through singular spectrum analysis (ssa) as an alternative for multi-channelization. </p>
			<p>The rest of this paper is organized as follows: first, the spectrum sensing tech niques are introduced; then, the proposed scheme for spectrum sensing is described in detail. Lastly, the effectiveness of the proposal is depicted and compared to the other outcomes by using artificial signals with different signal to noise relation (snr) levels, followed by a discussion of the results and conclusion. </p>
		</sec>
		<sec>
			<title>1. BACKGROUND</title>
			<sec>
				<title>1.1 Spectrum Sensing</title>
				<p>The spectrum sensing goal in cr is centered in hypothesis (Eq. (1,2)):</p>
				<p>
					<disp-formula id="e1">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e1.png"/> (<label>1</label>)</disp-formula>
				</p>
				<p>
					<disp-formula id="e2">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e2.png"/> (<label>2</label>)</disp-formula>
				</p>
				<p>where <italic>x</italic> (<italic>n</italic>) is the signal received by the secondary user; <italic>s</italic> (<italic>n</italic>) is the signal trans mitted by the primary user; <italic>ƞ</italic> (<italic>n</italic>) is an additive white Gaussian noise (AWGN); and <italic>H</italic>
 <sub>0</sub> is the null hypothesis that indicates absence of licensed user signal. </p>
			</sec>
			<sec>
				<title>1.2 Energy Detection (ed)</title>
				<p>In the ED technique, the signal power is measured in a particular frequency band along the time axis. The decision parameter for edis shown in Eq. (3):</p>
				<p>
					<disp-formula id="e3">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e3.png"/>(<label>3</label>)</disp-formula>
				</p>
				<p>where <italic>N</italic> is the number of samples in the analysis window. In this case, given a threshold <italic>ϒ</italic>, if <italic>E</italic>&gt;ϒ, the hypothesis <italic>H</italic>
 <sub>1</sub> is accepted, and if <italic>E</italic> &lt; ϒ the hypothesis <italic>H</italic>
 <sub>0</sub> is accepted. </p>
			</sec>
			<sec>
				<title>1.3 Cyclostationary Detection (cd)</title>
				<p>The cyclostationary features are produced by the signal periodicity (or by the periodic ity of its statistic measures). The cyclic spectral analysis implies second order trans formations of a function and its respective spectral representation <xref ref-type="bibr" rid="B8"><sup>8</sup></xref>. The distinction between primary and secondary users is possible only if different cyclic frequencies are involved <xref ref-type="bibr" rid="B9"><sup>9</sup></xref>. </p>
				<p>The Spectral Correlation Density (scd) of a signal could be calculated by Eq. (4) (10):</p>
				<p>
					<disp-formula id="e4">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e4.png"/> (<label>4</label>)</disp-formula>
				</p>
				<p>where<inline-graphic xlink:href="1692-3324-rium-15-29-00129-i005.png"/>is the Fourier Transform; <italic>α</italic> is the cyclic frequency and<inline-graphic xlink:href="1692-3324-rium-15-29-00129-i006.png"/>is the Cyclic Autocorrelation Function (CAF) of <italic>x</italic>, defined in Eq. (5):</p>
				<p>
					<disp-formula id="e5">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e5.png"/>(<label>5</label>)</disp-formula>
				</p>
				<p>where <italic>x</italic>
 <sup>*</sup> is the conjugate version of <italic>x</italic>. The caf exhibits peaks when the cyclic frequency <italic>_</italic> has the same value of the fundamental f requencies of the t ransmitted signal <xref ref-type="bibr" rid="B6"><sup>6</sup></xref>. </p>
				<p>The frequency-smoothing method is used to estimate the scd as recommended in <xref ref-type="bibr" rid="B11"><sup>11</sup></xref>
				</p>
				<p>
					<disp-formula id="e6">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e6.png"/>(<label>6</label>)</disp-formula>
				</p>
				<p>Where:</p>
				<p>
					<disp-formula id="e7">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e7.png"/>(<label>7</label>)</disp-formula>
				</p>
				<p>where <italic>N</italic> is the number of samples in the analysis window; <italic>s T</italic> is the time-sampling increment; <italic>s F</italic> is the frequency-sampling increment; and () <italic>Wv</italic>is the frequency smooth ing function centered at 0 =<italic>v</italic>of width <italic>L</italic>.</p>
			</sec>
		</sec>
		<sec>
			<title>2. PROPOSED SCHEME</title>
			<sec>
				<title>2.1 Blind Source Separation</title>
				<p>Blind source separation techniques allow separating different signals without previous knowledge about the signal features. The most common technique for source separa tion is ica.</p>
				<p>In ica, the received signal is assumed as a linear mixture of statistically indepen dent sources that could be separated maximizing the non-Gaussianity. This technique requires a mimo system with <italic>p</italic> sources and <italic>m</italic> receptors (<italic>p ≤ m</italic>), their combined out puts build the matrix of observable data <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i010.png"/>, where xi (j) is the j-th sample (or time instant) of the signal in the i-th receptor, and N is the number of samples in the analysis window. The principal idea in ica, is to get the original sources S from the matrix of observable data X as in Eq. (8):</p>
				<p>
					<disp-formula id="e8">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e8.png"/> (<label>8</label>)</disp-formula>
				</p>
				<p>where <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i012.png"/>is the mixture matrix and <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i013.png"/> , where <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i014.png"/> is the <italic>j</italic>-th sample (or time instant) of the signal transmitted by the <italic>l</italic>-th source, and <italic>N</italic> is the number of samples in the analysis window.</p>
				<p>If the matrix <bold>
 <italic>A</italic> 
</bold> is known, the source data could be reconstructed by its pseudo-inverse matrix. If <bold>
 <italic>A</italic> 
</bold> is unknown, it is necessary to calculate <bold>
 <italic>A</italic> 
</bold> and 𝑺 from the obser vation matrix <bold>
 <italic>X</italic> 
</bold> as in Eq: (9), in such a way that the estimated components <bold>
 <italic>ŝ</italic>
</bold> must be the linear combination of the data with the maximum independence <xref ref-type="bibr" rid="B12"><sup>12</sup></xref>.</p>
				<p>
					<disp-formula id="e9">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e9.png"/>(<label>9</label>)</disp-formula>
				</p>
				<p>where <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i016.png"/>is a matrix that maximizes the statistical independence.</p>
			</sec>
			<sec>
				<title>2.2 ms - ica</title>
				<p>In this work, the algorithm proposed by Molgedey (<xref ref-type="bibr" rid="B13">13</xref>), ms - ica, is chosen (Algorithm 1 (14)) due to its performance in online applications (<xref ref-type="bibr" rid="B15">15</xref>). </p>
				<p>Algorithm 1. ms - ica</p>
				<p>Input: Data matrix X </p>
				<p>Output: Estimated source matrix Ŝ</p>
				<p>1. Center the observation matrix variance <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i017.png"/>is the expectation operator. </p>
				<p>2. Whitening: transform X with non correlated components and unit variance </p>
				<p>a) Calculate the covariance matrix <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i018.png"/>
				</p>
				<p>b) Find the eigenvalues <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i019.png"/>where V is the eigenvector matrix and Ó is a diagonal matrix of eigenvalues</p>
				<p>c) Find the transformation matrix <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i020.png"/>
				</p>
				<p>d) Calculate the data matrix projection <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i021.png"/>
				</p>
				<p>3. Fix the lag τ</p>
				<p>4. Solve the eigenvalue problem for the ratio matrix <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i022.png"/> is the autocorrelation matrix of S at lag τ</p>
				<p>5. Compute the matrix <bold>
 <italic>A</italic>
</bold> as the eigenvalues matrix of <bold>
 <italic>Q</italic>
</bold> (16). </p>
				<p>6. Find the mixture matrix <bold>
 <italic>W</italic> 
</bold> =<bold>
 <italic>A</italic>
</bold> 
 <sup>-1</sup></p> <p>7. Calculate the source matrix <bold>
 <italic>ŝ</italic>
</bold> by Eq. (9) </p>
				<p>The lag parameter τ could be calculated as proposed in <xref ref-type="bibr" rid="B16"><sup>16</sup></xref>: a) compute <bold>
 <italic>ŝ</italic> 
</bold> as in algorithm 1 for τ =1; b) find the minimum value of <italic>δ</italic> (τ ) in Eq. (10) by iterating the parameter τ, where <sub>
 <sup>
 <italic>cs</italic> 
</sup> 
</sub> are the normalized autocorrelation functions; and c) recalculate <bold>
 <italic>ŝ</italic> 
</bold> after <italic>q</italic> iterations.</p>
				<p>
					<disp-formula id="e10">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e10.png"/>(<label>10</label>)</disp-formula>
				</p>
			</sec>
			<sec>
				<title>2.3 Single-to-Multichannel Conversion</title>
				<p>Aiming to find a matrix representation for a single-channel data, the procedure pro posed by (17) through ssa was applied as follows: </p>
				<p>1. Calculate the autocorrelation function of <italic>x</italic> (<italic>n</italic>) as in Eq. (11), where <italic>l</italic> = 0,…,1 =…<italic>L</italic>− and <italic>L</italic> is the analysis window length. A suitable value for <italic>L</italic> is an important issue in ssa. In this work, the value is fixed as one more than the number of peaks in the periodogram of the received signal <italic>x</italic> (<italic>n</italic>), aiming to find an approximate number of expected sources.</p>
				<p>
					<disp-formula id="e11">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e11.png"/>(<label>11</label>)</disp-formula>
				</p>
				<p>2. Find the Toeplitz matrix <sub>
 <sup>
 <bold>
 <italic>To</italic>
</bold> 
</sup> 
</sub> of the autocorrelation function. </p>
				<p>3. Calculate de eigendecomposition of <sub>
 <sup>
 <bold>
 <italic>To</italic> 
</bold> 
</sup> 
</sub> as <sub>
 <sup>
 <bold>
 <italic>To</italic>
</bold> 
</sup> 
</sub> = <bold>
 <italic>V</italic>
</bold> 
 <bold>Ó</bold>
 <sub>
 <sup>
 <bold>
 <italic>VT</italic>
</bold> 
</sup> 
</sub> 
 <bold>.</bold> Sort eigenvalues (<bold>Ó</bold>) in descending order and preserve the most informative components. </p>
				<p>4. Create the Hankel matrix <bold>
 <italic>TH</italic>
</bold> of <italic>x</italic> (<italic>n</italic>), with <italic>K</italic> = <italic>N</italic> −1 + <italic>L</italic> (Eq. (12)).</p>
				<p>
					<disp-formula id="e12">
						<graphic xlink:href="1692-3324-rium-15-29-00129-e12.png"/>(<label>12</label>)</disp-formula>
				</p>
				<p>5. Project <bold>
 <italic>H T</italic> 
</bold> into the subspace of eigenvectors of <inline-graphic xlink:href="1692-3324-rium-15-29-00129-i026.png"/> The main idea is to define the temporal behaviors. </p>
				<p>6. Reconstruct <bold>
 <italic>X</italic> 
</bold> as a sum of the projections obtained by the inverse Hankel procedure of <bold>
 <italic>Y j</italic>
</bold> .</p>
				<p>7. Calculate algorithm 1 with the following modifications: </p>
				<p>a. Choose two subsets, i.e., <italic>j</italic> = 2.</p>
				<p>b. Use the first subset as the input of the algorithm. </p>
				<p>c. Divide the remaining subset into two subsets. </p>
				<p>d. Go to step b. and repeat the process.</p>
			</sec>
		</sec>
		<sec>
			<title>3. EXPERIMENTAL SET-UP</title>
			<p>The different schemes for spectrum sensing are tested with an artificial mixture signal defined by Eq. (13).</p>
			<p>
				<disp-formula id="e13">
					<graphic xlink:href="1692-3324-rium-15-29-00129-e13.png"/>(<label>13</label>)</disp-formula>
			</p>
			<p>The length <italic>N</italic> of <italic>x</italic> (<italic>n</italic>) with <italic>n</italic> =1,…,<italic>N</italic> is fixed in 4,000 bins, the sampling frequency <sub>
 <sup>
 <italic>fs</italic> 
</sup> 
</sub> = 10 <italic>kHz</italic>, and the parameter <italic>L</italic>= 512 for the cyclostationary analysis (smoothed window length). The white Gaussian noise <italic>ƞ</italic>(<italic>n</italic>) is a variable parameter, with SNR between 10 −and 20 <italic>dB</italic>.</p>
			<p>In order to check the technique’s efficiency, some variations in the number of original sources and noise power (<sub>
 <sup>
 <italic>Pƞ</italic>
</sup> 
</sub> ) are made. The amplitude of the sources is <italic>ai</italic> = 1 and the phase <italic>Φi</italic> =0 . For ED, a system with six fir filters is implemented, each one with order 100, Hamming window, bandwidth 520<italic>Hz</italic>, a center frequency of <sub>
 <sup>
 <italic>fi</italic> 
</sup> 
</sub> and 10<italic>Hz</italic> of overlapping; then, the procedure recommended by (18) is followed for each channel with 1,000 Monte Carlo simulations, aiming to determine whether a channel is available or not. For CD, the parameter <italic>μ</italic> is fixed by heuristic tuning and it determines when a peak of the caf corresponds to a fundamental frequency or not. </p>
			<p>The elapsed time that the algorithms took to run is measured in an Intel<sup>
 <bold>®</bold>
</sup> Core™ i5-2,400 with 4 GB of ram memory and Python<sup>
 <bold>®</bold>
</sup> of 64 bits.</p>
		</sec>
		<sec sec-type="results|discussion">
			<title>4. RESULTS AND DISCUSSION</title>
			<p>
				<xref ref-type="table" rid="t1">Table 1</xref> summarizes the results for the three tested schemes. The technique most sensi tive to noise is ed, with the disadvantage that if the frequency bands are unknown, the algorithm is unable of detect different frequency components. The mean time elapsed for this technique is 97.85 ± 1.02 s, and tends to increase with the number of channels. The cd and ms-ica’s are noise-robust algorithms. The time elapsed for cd is around 106.94 ± 1.77s that could be compensated by its precision. The cd sensing technique generates some false positives, but the rate of false negatives is zero. If the intention of the secondary user is to search for vacancy without disturbing the primary user, and the time is not a critical issue in the application, this technique could be a convenient solution. The only inconvenient detected is the parameter μ, because the algorithm performance depends on the choice of this parameter.</p>
			<p>
				<table-wrap id="t1">
					<label>Table 1</label>
					<caption>
						<title>Spectrumsensing test results</title>
					</caption>
					<graphic xlink:href="1692-3324-rium-15-29-00129-gt1.jpg"/>
					<graphic xlink:href="1692-3324-rium-15-29-00129-gt1.jpg"/>
					<table-wrap-foot>
						<fn id="TFN1">
							<p>Source: authors</p>
						</fn>
					</table-wrap-foot>
				</table-wrap>
			</p>
			<p>The proposed ms-ica technique exhibits a suitable compromise between time and performance. The algorithm is fast enough to sense without causing long interrup tions in transmissions. The major disadvantage is the misdetection of some sources. To overcome this inconvenient, it is recommended to run the algorithm a couple of times and increase the amount of information about the frequency occupancy. The remarkable aspect is the elapsed time (0.15 ± 0.1 s), especially for online applications with demanding time issues. </p>
		</sec>
		<sec sec-type="conclusions">
			<title>5. CONCLUSION</title>
			<p>A technique for spectrum sensing based on blind source separation is explored. The methodology includes a single-to-multichannel procedure through singular spectrum analysis (ssa), aiming to simulate a mimo communication system. In this work, three approaches were tested: Energy Detection, Cyclostationary Analysis and the proposed scheme with ms-ica. The methodology is proved on a collection of artificial mixture signals with different frequencies and noise power. The proposed scheme is able to detect most of the sources, even with high noise presence, with low time consumption. The most accurate technique is Cyclostationary Analysis but it is more demanding in terms of time than the others. As future work, it would be interesting to deal with other kinds of modulations, codifications and channel models.</p>
		</sec>
		<sec>
			<title>6. ACKNOWLEDGMENTS</title>
			<p>Thisworkissupportedby “Validación de algoritmos de aprendizaje automático usados en sistemas de telecomunicaciones de radio cognitiva (Código No. 737)” financedby Universidad de Medellín and by “Ayudas para la realización del doctorado” (RR01/2011) of the Universidad Politécnica de Madrid y TEC2012-38630-C04-01 del Ministerio de Educación de España” </p>
		</sec>
	</body>
	<back>
		<ref-list>
			<title>REFERENCES</title>
			<ref id="B1">
				<label>1</label>
				<mixed-citation>1. M. Sarijari, A. Marwanto, N. Fisal, S. K. S. Yusof, R. Rashid &amp; M. Satria, “Energy detection sensing based on gnu radio and usrp: An analysis study”, in <italic>Proceedings of the 2009 IEEE 9th Malaysia International Conference on Communications</italic>, Kuala Lumpur, Malaysia, Dec. 15-17, 2009, pp. 338-342. Available: <ext-link ext-link-type="uri" xlink:href="https://pdfs.semanticscholar.org/2319/11ca01bab559aed1 7bd745bb6ae4c97a8d70.pdf">https://pdfs.semanticscholar.org/2319/11ca01bab559aed1 7bd745bb6ae4c97a8d70.pdf</ext-link>.</mixed-citation>
				<element-citation publication-type="confproc">
					<person-group person-group-type="author">
						<name>
							<surname>Sarijari</surname>
							<given-names>M</given-names>
						</name>
						<name>
							<surname>Marwanto</surname>
							<given-names>A</given-names>
						</name>
						<name>
							<surname>Fisal</surname>
							<given-names>N</given-names>
						</name>
						<name>
							<surname>Yusof</surname>
							<given-names>S. K. S</given-names>
						</name>
						<name>
							<surname>Rashid</surname>
							<given-names>R</given-names>
						</name>
						<name>
							<surname>Satria</surname>
							<given-names>M</given-names>
						</name>
					</person-group>
					<source>Energy detection sensing based on gnu radio and usrp: An analysis study</source>
					<conf-name>Proceedings of the 2009 IEEE 9th Malaysia International Conference on Communications</conf-name>
					<season>Dec-17,</season>
					<year>2009</year>
					<fpage>338</fpage>
					<lpage>342</lpage>
					<ext-link ext-link-type="uri" xlink:href="https://pdfs.semanticscholar.org/2319/11ca01bab559aed1 7bd745bb6ae4c97a8d70.pdf">https://pdfs.semanticscholar.org/2319/11ca01bab559aed1 7bd745bb6ae4c97a8d70.pdf</ext-link>
				</element-citation>
			</ref>
			<ref id="B2">
				<label>2</label>
				<mixed-citation>2. Y. Hassan, M. El-Tarhuni&amp; K. Assaleh, “Learning-based spectrum sensing for cognitive radio systems”, J<italic>ournal of Computer Networks and Communications</italic>, vol. 2012, pp. 1-14, 2012. Available: <ext-link ext-link-type="uri" xlink:href="https://www.hindawi.com/journals/jcnc/2012/259824/">https://www.hindawi.com/journals/jcnc/2012/259824/</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Hassan</surname>
							<given-names>Y</given-names>
						</name>
						<name>
							<surname>El-Tarhuni</surname>
							<given-names>M</given-names>
						</name>
						<name>
							<surname>Assaleh</surname>
							<given-names>K</given-names>
						</name>
					</person-group>
					<article-title>Learning-based spectrum sensing for cognitive radio systems</article-title>
					<source>Journal of Computer Networks and Communications</source>
					<volume>2012</volume>
					<fpage>1</fpage>
					<lpage>14</lpage>
					<year>2012</year>
					<ext-link ext-link-type="uri" xlink:href="https://www.hindawi.com/journals/jcnc/2012/259824/">https://www.hindawi.com/journals/jcnc/2012/259824/</ext-link>
				</element-citation>
			</ref>
			<ref id="B3">
				<label>3</label>
				<mixed-citation>3. A. Mate, K. H. Lee &amp; I. T. Lu, “Spectrum sensing based on time covariance matrix using gnu radio and usrp for cognitive radio”, in: 2011 ieee Long Island Systems, Applications and Technology Conference (lisat), Farmingdale, NY, USA, May 6, 2011, pp. 1-6.</mixed-citation>
				<element-citation publication-type="confproc">
					<person-group person-group-type="author">
						<name>
							<surname>Mate</surname>
							<given-names>A</given-names>
						</name>
						<name>
							<surname>Lee</surname>
							<given-names>K. H</given-names>
						</name>
						<name>
							<surname>Lu</surname>
							<given-names>I. T</given-names>
						</name>
					</person-group>
					<source>Spectrum sensing based on time covariance matrix using gnu radio and usrp for cognitive radio</source>
					<conf-name>2011 ieee Long Island Systems, Applications and Technology Conference (lisat)</conf-name>
					<publisher-name>Farmingdale</publisher-name>
					<publisher-loc>NY, USA</publisher-loc>
					<day>06</day>
					<month>05</month>
					<year>2011</year>
					<fpage>1</fpage>
					<lpage>6</lpage>
				</element-citation>
			</ref>
			<ref id="B4">
				<label>4</label>
				<mixed-citation>4. G. Nautiyal&amp; R. Kumar, “Spectrum sensing in cognitive radio using matlab”, <italic>International Journal of Engineering and Advanced Technology (IJEAT)</italic>, vol. 2, no. 5, pp. 529-532, Jun. 2013.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Nautiyal</surname>
							<given-names>G</given-names>
						</name>
						<name>
							<surname>Kumar</surname>
							<given-names>R</given-names>
						</name>
					</person-group>
					<article-title>Spectrum sensing in cognitive radio using matlab</article-title>
					<source>International Journal of Engineering and Advanced Technology (IJEAT)</source>
					<volume>2</volume>
					<issue>5</issue>
					<fpage>529</fpage>
					<lpage>532</lpage>
					<month>06</month>
					<year>2013</year>
				</element-citation>
			</ref>
			<ref id="B5">
				<label>5</label>
				<mixed-citation>5. Z. Xuping&amp; P. Jianguo, “Energy-detection based spectrum sensing for cognitive radio”, in <italic>iet Conference on Wireless, Mobile and Sensor Networks (ccwmsn07)</italic>, Shangai, China, Dec. 12-14, 2007, pp. 944-947.</mixed-citation>
				<element-citation publication-type="confproc">
					<person-group person-group-type="author">
						<name>
							<surname>Xuping</surname>
							<given-names>Z</given-names>
						</name>
						<name>
							<surname>Jianguo</surname>
							<given-names>P</given-names>
						</name>
					</person-group>
					<source>-detection based spectrum sensing for cognitive radio</source>
					<conf-name>iet Conference on Wireless, Mobile and Sensor Networks (ccwmsn07)</conf-name>
					<publisher-loc>Shangai, China</publisher-loc>
					<season>Dec-14,</season>
					<year>2007</year>
					<fpage>944</fpage>
					<lpage>947</lpage>
				</element-citation>
			</ref>
			<ref id="B6">
				<label>6</label>
				<mixed-citation>6. H. Arslan, Cognitive Radio, Software Defined Radio, and Adaptive Wireless Systems (Signals and Communication Technology). New York: Springer-Verlag, 2007.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Arslan</surname>
							<given-names>H</given-names>
						</name>
					</person-group>
					<source>Cognitive Radio, Software Defined Radio, and Adaptive Wireless Systems (Signals and Communication Technology)</source>
					<publisher-loc>New York</publisher-loc>
					<publisher-name>Springer-Verlag</publisher-name>
					<year>2007</year>
				</element-citation>
			</ref>
			<ref id="B7">
				<label>7</label>
				<mixed-citation>7. M. Rahman, A. Haniz, S. Khadka, S., Iswandi, Gahadza, M., Kim, M., ichi Takada, J. “De velopment of spectrum sensing system with gnu radio and usrp to detect emergency radios”, ieice, The Institute of Electronics, Information and Communication Engineers, Sendai, Japan, Technical Report SR2009-57, Oct. 2009.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Rahman</surname>
							<given-names>M</given-names>
						</name>
						<name>
							<surname>Haniz</surname>
							<given-names>A</given-names>
						</name>
						<name>
							<surname>Khadka</surname>
							<given-names>S</given-names>
						</name>
						<name>
							<surname>Iswandi</surname>
							<given-names>S</given-names>
						</name>
						<name>
							<surname>Gahadza</surname>
							<given-names>M</given-names>
						</name>
						<name>
							<surname>Kim</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>ichi Takada</surname>
							<given-names>J</given-names>
						</name>
					</person-group>
					<chapter-title>De velopment of spectrum sensing system with gnu radio and usrp to detect emergency radios</chapter-title>
					<source>ieice, The Institute of Electronics, Information and Communication Engineers</source>
					<publisher-loc>Sendai, Japan</publisher-loc>
					<month>10</month>
					<year>2009</year>
				</element-citation>
			</ref>
			<ref id="B8">
				<label>8</label>
				<mixed-citation>8. A. Fehske, J. Gaeddert&amp; J. Reed, “A new approach to signal classification using spectral correlation and neural networks”, in DySPAN 2005. First IEEE International Symposium on New Frontiers in Dynamic Spectrum Access Networks, Baltimore, USA, Nov. 8-11, 2005. pp. 144-150.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Fehske</surname>
							<given-names>A</given-names>
						</name>
						<name>
							<surname>Gaeddert</surname>
							<given-names>J</given-names>
						</name>
					</person-group>
					<chapter-title>A new approach to signal classification using spectral correlation and neural networks</chapter-title>
					<source>DySPAN 2005. First IEEE International Symposium on New Frontiers in Dynamic Spectrum Access Networks</source>
					<publisher-loc>Baltimore, USA</publisher-loc>
					<season>Nov-11,</season>
					<year>2005</year>
					<fpage>144</fpage>
					<lpage>150</lpage>
				</element-citation>
			</ref>
			<ref id="B9">
				<label>9</label>
				<mixed-citation>9. S. Chaudhari, “Spectrum sensing for cognitive radios: Algorithms, performance, and limita tions”, Ph. D. thesis, Aalto University, Greater Helsinki, Finland, 2012.</mixed-citation>
				<element-citation publication-type="thesis">
					<person-group person-group-type="author">
						<name>
							<surname>Chaudhari</surname>
							<given-names>S</given-names>
						</name>
					</person-group>
					<source>Spectrum sensing for cognitive radios: Algorithms, performance, and limita tions</source>
					<comment content-type="degree">Ph. D</comment>
					<publisher-name>Aalto University</publisher-name>
					<publisher-loc>Greater Helsinki, Finland</publisher-loc>
					<year>2012</year>
				</element-citation>
			</ref>
			<ref id="B10">
				<label>10</label>
				<mixed-citation>10. S. Da, G. Xiaoying, C. Hsiao-hwa &amp; Q. Liang, “Fast cycle frequency domain feature detection for cognitive radio systems”, <italic>Arxiv</italic>, p. 4, Ar. 6, 2009. Available: <ext-link ext-link-type="uri" xlink:href="https://archive.org/details/ arxiv-0903.1183">https://archive.org/details/ arxiv-0903.1183</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Da</surname>
							<given-names>S</given-names>
						</name>
						<name>
							<surname>Xiaoying</surname>
							<given-names>G</given-names>
						</name>
						<name>
							<surname>Hsiao-hwa</surname>
							<given-names>C</given-names>
						</name>
						<name>
							<surname>Liang</surname>
							<given-names>Q</given-names>
						</name>
					</person-group>
					<article-title>Fast cycle frequency domain feature detection for cognitive radio systems</article-title>
					<source>Arxiv</source>
					<fpage>4</fpage>
					<lpage>4</lpage>
					<day>06</day>
					<month>04</month>
					<year>2009</year>
					<ext-link ext-link-type="uri" xlink:href="https://archive.org/details/ arxiv-0903.1183">https://archive.org/details/ arxiv-0903.1183</ext-link>
				</element-citation>
			</ref>
			<ref id="B11">
				<label>11</label>
				<mixed-citation>11. M. Calabro, “A Cooperative Spectrum Sensing Network with Signal Classification Capabili ties”. Ph. D. thesis, Worcester Polytechnic Institute, Worcester, Massachusetts, 2010.</mixed-citation>
				<element-citation publication-type="thesis">
					<person-group person-group-type="author">
						<name>
							<surname>Calabro</surname>
							<given-names>M</given-names>
						</name>
					</person-group>
					<source>A Cooperative Spectrum Sensing Network with Signal Classification Capabili ties</source>
					<comment content-type="degree">Ph. D</comment>
					<publisher-name>Worcester Polytechnic Institute, Worcester</publisher-name>
					<publisher-loc>Massachusetts</publisher-loc>
					<year>2010</year>
				</element-citation>
			</ref>
			<ref id="B12">
				<label>12</label>
				<mixed-citation>12. A. G. Ferrer, E.G. Prieto&amp; D. Peña, “Exploring ica for time series decomposition”, <italic>Working Paper 11-16, Statistics and Econometrics Series 11</italic>, May 2011. Available: <ext-link ext-link-type="uri" xlink:href="http://orff.uc3m.es/ bitstream/handle/10016/11285/ws111611.pdf?sequence=1">http://orff.uc3m.es/ bitstream/handle/10016/11285/ws111611.pdf?sequence=1</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Ferrer</surname>
							<given-names>A. G</given-names>
						</name>
						<name>
							<surname>Prieto</surname>
							<given-names>E.G</given-names>
						</name>
						<name>
							<surname>Peña</surname>
							<given-names>D</given-names>
						</name>
					</person-group>
					<article-title>Exploring ica for time series decomposition</article-title>
					<source>Working Paper 11-16, Statistics and Econometrics Series 11</source>
					<month>05</month>
					<year>2011</year>
					<ext-link ext-link-type="uri" xlink:href="http://orff.uc3m.es/ bitstream/handle/10016/11285/ws111611.pdf?sequence=1">http://orff.uc3m.es/ bitstream/handle/10016/11285/ws111611.pdf?sequence=1</ext-link>
				</element-citation>
			</ref>
			<ref id="B13">
				<label>13</label>
				<mixed-citation>13. L. Molgedey&amp; H. G. Schuster, “Separation of a mixture of independent signals using time delayed correlations”, <italic>Physical Review Letters</italic>, vol. 72, 3634-3637, 1994.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Molgedey</surname>
							<given-names>L</given-names>
						</name>
						<name>
							<surname>Schuster</surname>
							<given-names>H. G</given-names>
						</name>
					</person-group>
					<article-title>Separation of a mixture of independent signals using time delayed correlations</article-title>
					<source>Physical Review Letters</source>
					<volume>72</volume>
					<fpage>3634</fpage>
					<lpage>3637</lpage>
					<year>1994</year>
				</element-citation>
			</ref>
			<ref id="B14">
				<label>14</label>
				<mixed-citation>14. V. Krishnaveni, S. Jayaraman, P. M. Kumar, K. Shivakumar&amp; K. Ramadoss, “Comparison of independent component analysis algorithms for removal of ocular artifacts from electro encephalogram”, <italic>Meas. Sci. Rev. J</italic>, vol. 5, no. 2, pp. 67-78, 2005.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Krishnaveni</surname>
							<given-names>V</given-names>
						</name>
						<name>
							<surname>Jayaraman</surname>
							<given-names>S</given-names>
						</name>
						<name>
							<surname>Kumar</surname>
							<given-names>P. M</given-names>
						</name>
						<name>
							<surname>Shivakumar</surname>
							<given-names>K</given-names>
						</name>
						<name>
							<surname>Ramadoss</surname>
							<given-names>K</given-names>
						</name>
					</person-group>
					<article-title>Comparison of independent component analysis algorithms for removal of ocular artifacts from electro encephalogram</article-title>
					<source>Meas. Sci. Rev. J</source>
					<volume>5</volume>
					<issue>2</issue>
					<fpage>67</fpage>
					<lpage>78</lpage>
					<year>2005</year>
				</element-citation>
			</ref>
			<ref id="B15">
				<label>15</label>
				<mixed-citation>15. Hongli, Sun, Y.: “The study and test of ica algorithms”, in <italic>2005 Proceedings International Conference on Wireless Communications, Networking and Mobile Computing</italic>
 <italic>, vol.</italic> 1, Wuhan, China, Sept. 23-26, 2005, pp. 602-605.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">Hongli<name>
							<surname>Sun</surname>
							<given-names>Y</given-names>
						</name>
					</person-group>
					<chapter-title>The study and test of ica algorithms</chapter-title>
					<source>2005 Proceedings International Conference on Wireless Communications, Networking and Mobile Computing</source>
					<volume>1</volume>
					<publisher-loc>Wuhan, China</publisher-loc>
					<season>Sep-26,</season>
					<year>2005</year>
					<fpage>602</fpage>
					<lpage>605</lpage>
				</element-citation>
			</ref>
			<ref id="B16">
				<label>16</label>
				<mixed-citation>16. T. Kolenda, L. K. Hansen &amp; J. Larsen, “Signal detection using ica: Application to chat room topic spotting”, in <sub>
 <sup>
 <italic>3rd International Conference on Independent Component Analysis and Blind Source Separation, ica’2001</italic>
</sup> 
</sub> , San Diego, USA, Dec. 9-13, 2001, pp. 540-545. Available: <ext-link ext-link-type="uri" xlink:href="http:// cogsys.imm.dtu.dk/publications/2001/kolenda.ica2001.pdf">http:// cogsys.imm.dtu.dk/publications/2001/kolenda.ica2001.pdf</ext-link></mixed-citation>
				<element-citation publication-type="confproc">
					<person-group person-group-type="author">
						<name>
							<surname>Kolenda</surname>
							<given-names>T</given-names>
						</name>
						<name>
							<surname>Hansen</surname>
							<given-names>L. K</given-names>
						</name>
						<name>
							<surname>Larsen</surname>
							<given-names>J</given-names>
						</name>
					</person-group>
					<source>Signal detection using ica: Application to chat room topic spotting</source>
					<conf-name>3rd International Conference on Independent Component Analysis and Blind Source Separation, ica’2001</conf-name>
					<publisher-loc>San Diego, USA</publisher-loc>
					<season>Dec-13,</season>
					<year>2001</year>
					<fpage>540</fpage>
					<lpage>545</lpage>
					<ext-link ext-link-type="uri" xlink:href="http:// cogsys.imm.dtu.dk/publications/2001/kolenda.ica2001.pdf">http:// cogsys.imm.dtu.dk/publications/2001/kolenda.ica2001.pdf</ext-link>
				</element-citation>
			</ref>
			<ref id="B17">
				<label>17</label>
				<mixed-citation>17. H. G. Ma, Q. B. Jiang, Z. Q. Liu, G. Liu &amp; Z. Y. Ma, “A novel blind source separation method for single-channel signal”, <italic>Signal Processing</italic>, vol. 90, no. 12, pp. 3232-3241, 2010.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Ma</surname>
							<given-names>H. G</given-names>
						</name>
						<name>
							<surname>Jiang</surname>
							<given-names>Q. B</given-names>
						</name>
						<name>
							<surname>Liu</surname>
							<given-names>Z. Q</given-names>
						</name>
						<name>
							<surname>Liu</surname>
							<given-names>G</given-names>
						</name>
						<name>
							<surname>Ma</surname>
							<given-names>Z. Y</given-names>
						</name>
					</person-group>
					<article-title>A novel blind source separation method for single-channel signal</article-title>
					<source>Signal Processing</source>
					<volume>90</volume>
					<issue>12</issue>
					<fpage>3232</fpage>
					<lpage>3241</lpage>
					<year>2010</year>
				</element-citation>
			</ref>
			<ref id="B18">
				<label>18</label>
				<mixed-citation>18. S. S. Kalamkar&amp; A. Banerjee, “On the performance of generalized energy detector under noise uncertainty in cognitive radio”, in <italic>National Conference on Communications (ncc)</italic>, Delhi, India, Feb. 15-17, 2013. pp. 1-5</mixed-citation>
				<element-citation publication-type="confproc">
					<person-group person-group-type="author">
						<name>
							<surname>Kalamkar</surname>
							<given-names>S. S</given-names>
						</name>
						<name>
							<surname>Banerjee</surname>
							<given-names>A</given-names>
						</name>
					</person-group>
					<source>On the performance of generalized energy detector under noise uncertainty in cognitive radio</source>
					<conf-name>National Conference on Communications (ncc)</conf-name>
					<publisher-loc>Delhi, India</publisher-loc>
					<season>Feb-17,</season>
					<year>2013</year>
					<fpage>1</fpage>
					<lpage>5</lpage>
				</element-citation>
			</ref>
		</ref-list>
		<fn-group>
			<fn fn-type="other" id="fn1">
				<label>1</label>
				<p>This article is a result from the research” Validation of automatic learning algorithms used in cognitive radio telecommunications systems “, funded by the University of Medellin, during the period 2013-2, 2015-2</p>
			</fn>
		</fn-group>
	</back>
</article>