\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Detection and Parameter Extraction of Low Probability of Intercept Frequency Hopping Signals using the Spectrogram and the Reassigned Spectrogram},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{2015-03-15 (revised: 15 March 2015)}
\def\TheID{\makeatother }
\def\TheDate{2015-03-15}
\title{Detection and Parameter Extraction of Low Probability of Intercept Frequency Hopping Signals using the Spectrogram and the Reassigned Spectrogram}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Daniel L.  Stevens}

             \author[2]{Stephanie A.  Schuckers}

             \affil[1]{  Air Force Research Laboratory, Rome, NY}

\renewcommand\Authands{ and }

\date{\small \em Received: 10 February 2015 Accepted: 1 March 2015 Published: 15 March 2015}

\maketitle


\begin{abstract}
        


Digital intercept receivers are currently moving away from Fourier-based analysis and towards classical time-frequency analysis techniques, such as the Wigner-Ville distribution, Choi- Williams distribution, spectrogram, and scalogram, for the purpose of analyzing low probability of intercept radar signals (e.g. triangular modulated frequency modulated continuous wave and frequency shift keying). Although these classical time-frequency techniques are an improvement over the Fourier-based analysis, they still suffer from a lack of readability, due to cross-term interference, and a mediocre performance in low SNR environments. This lack of readability may lead to inaccurate detection and parameter extraction of these radar signals. In this paper, the use of the Hough transform, because of its ability to suppress cross-term interference, separate signals from cross-terms, and perform well in the presence of noise, is proposed as an improved signal analysis technique. With these qualities, the Hough transform has the potential to produce better readability and consequently, more accurate signal detection and parameter extraction metrics.

\end{abstract}


\keywords{radar detection, hough transform, low probability of intercept.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\let\tabcellsep& 	 	 		 
\section[{Introduction a) Frequency hopping techniques}]{Introduction a) Frequency hopping techniques}\par
low probability of intercept (LPI) radar that uses frequency hopping techniques changes the transmitting frequency in time over a wide bandwidth in order to prevent an intercept receiver from intercepting the waveform. The frequency slots used are chosen from a frequency hopping sequence, and it is this unknown sequence that gives the radar the advantage over the intercept receiver in terms of processing gain. The frequency sequence appears random to the intercept receiver, and so the possibility of it following the changes in frequency is remote  {\ref [PAC09]}. This prevents a jammer from reactively jamming the transmitted frequency [ADA04]. Frequency hopping radar performance depends only slightly on the code used, given that certain properties are met. This allows for a larger variety of codes, making it more difficult to intercept 1 Though classical time-frequency analysis techniques are a great improvement over Fourier analysis techniques, they may suffer from poor timefrequency localization, as described above. This may result in degraded readability of time-frequency representations, potentially leading to inaccurate LPI radar signal detection and parameter extraction metrics, and as such, can lead to decisions based on inaccurate information. 
\section[{c) Reassignment method}]{c) Reassignment method}\par
A promising avenue for overcoming this deficiency is the utilization of the reassignment method. The reassignment method, which can be applied to most energy distributions [HIP00], has, in theory, a perfectly localized distribution for chirps, tones and impulses [BOA03], making it a good candidate for the analysis of certain LPI radar signals, such as the triangular modulated frequency modulated continuous wave (FMCW) (which can be viewed as back-to-back chirps) and the frequency shift keying (FSK) (which can be viewed as tones). 
\section[{d) Spectrogram and reassigned spectrogram}]{d) Spectrogram and reassigned spectrogram}\par
The spectrogram is defined as the magnitude squared of the STFT[BOA03], [HIP00], [HLA92], [MIT01], [PAC09]. For non-stationary signals, the STFT is usually in the form of the spectrogram [GRI08].\par
The STFT of a signal ??(??)is given in equation 2.5 as:?? ?? (??, ð??"ð??"; ?) = ? ??(??)? +? ?? (????)?? ??? 2??ð??"ð??"?? ???? (2.5)\par
Where ?(??) is a short time analysis window localized around ?? = 0 and ð??"ð??" = 0. Because multiplication by the relatively short window ?(?? ? ??) effectively suppresses the signal outside a neighborhood around the analysis point ?? = ??, the STFT is a 'local' spectrum of the signal ??(??) around ??. Think of the window ?(??) as sliding along the signal ??(??) and for each shift ?(?? ? ??) we compute the usual Fourier transform of the product function ??(??)?(?? ? ??) . The observation window allows localization of the spectrum in time, but also smears the spectrum in frequency in accordance with the uncertainty principle, leading to a trade-off between time resolution and frequency resolution. In general, if the window is short, the time resolution is good, but the frequency resolution is poor, and if the window is long, the frequency resolution is good, but the time resolution is poor.\par
The STFT was the first tool devised for analyzing a signal in both time and frequency simultaneously. For analysis of human speech, the main method was, and still is, the STFT. In general, the STFT is still the most widely used method for studying non-stationary signals  {\ref [COH95]}.\par
The spectrogram (the squared modulus of the STFT) is given by equation 2.6 as:?? ?? (??, ð??"ð??") = ?? ??(??) +? ?? ?(?? ? ??)?? ??? 2??ð??"ð??"?? ????? 2 (2.6)\par
The spectrogram is a real-valued and nonnegative distribution. Since the window h of the STFT is assumed of unit energy, the spectrogram satisfies the global energy distribution property. Thus we can interpret the spectrogram as a measure of the energy of the signal contained in the time-frequency domain centered on the point (t, f) and whose shape is independent of this localization.\par
Here are some properties of the spectrogram: 1) time and frequency covariance -the spectrogram preserves time and frequency shifts, thus the spectrogram is an element of the class of quadratic time-frequency distributions that are covariant by translation in time and in frequency (i.e. Cohen's class); 2) time-frequency resolution -the time-frequency resolution of the spectrogram is limited exactly as it is for the STFT; there is a trade-off between time resolution and frequency resolution. This poor resolution is the main drawback of this representation; 3) interference structure -as it is a quadratic (or bilinear) representation, the spectrogram of the sum of two signals is not the sum of the two spectrograms (quadratic superposition principle); there is a crossspectrogram part and a real part. Thus, as for every quadratic distribution, the spectrogram presents interference terms; however, those interference terms are restricted to those regions of the time-frequency plane where the signals overlap. Thus if the signal components are sufficiently distant so that their spectrograms do not overlap significantly, then the interference term will nearly be identically zero[COH95], [HLA92], [ISI96].\par
The original idea of reassignment was introduced in an attempt to improve the spectrogram [OZD03]. As with any other bilinear energy distribution, the spectrogram is faced with an unavoidable trade-off between the reduction of misleading interference terms and a sharp localization of the signal components. 
\section[{Global Journal of Researches in Engineering ( ) Volume XX Issue IV Version I 13}]{Global Journal of Researches in Engineering ( ) Volume XX Issue IV Version I 13}\par
Year 2020F © 2020 Global Journals\par
We can define the spectrogram as a twodimensional convolution of the WVD of the signal by the WVD of the analysis window, as in equation 2.9:?? ?? (??, ð??"ð??"; ?) = ? ?? ?? +? ?? (??, ??)?? ? (?? ? ??, ð??"ð??" ? ??)???? ???? (2.9)\par
Therefore, the distribution reduces the interference terms of the signal's WVD, but at the expense of time and frequency localization. However, a closer look at equation 2.9 shows that ?? ? (?? ? ??, ð??"ð??" ? ??) delimits a time-frequency domain at the vicinity of the (??, ð??"ð??") point, inside which a weighted average of the signal's WVD values is performed. The key point of the reassignment principle is that these values have no reason to be symmetrically distributed around (??, ð??"ð??") , which is the geometrical center of this domain. Therefore, their average should not be assigned at this point, but rather at the center of gravity of this domain, which is much more representative of the local energy distribution of the signal [BOA03]. Reasoning with a mechanical analogy, the local energy distribution ?? ? (?? ? ??, ð??"ð??" ? ??)?? ?? (??, ??) (as a function of ?? and ??) can be considered as a mass distribution, and it is much more accurate to assign the total mass (i.e. the spectrogram value) to the center of gravity of the domain rather than to its geometrical center. Another way to look at it is this: the total mass of an object is assigned to its geometrical center, an arbitrary point which except in the very specific case of a homogeneous distribution, has no reason to suit the actual distribution. A much more meaningful choice is to assign the total mass of an object, as well as the spectrogram value, to the center of gravity of their respective distribution [BOA03]. This is exactly how the reassignment method proceeds: it moves each value of the spectrogram computed at any point (??, ð??"ð??") to another point (?? ,ð??"ð??" ?) which is the center of gravity of the signal energy distribution around (??, ð??"ð??") (see equations 2.10 and 2.11) [LIX08]:?? (??; ??, ð??"ð??") = ? ?? ?? ? (?????,ð??"ð??"??? )?? ?? (??,?? )???? ???? +? ?? ? ?? ? (?????,ð??"ð??"??? )?? ?? (??,?? )???? ???? +? ?? (2.10) ð??"ð??" ?(??; ??, ð??"ð??") = ? ?? ?? ? (?? ? ??, ð??"ð??" ? ??)?? ?? (??, ??)???? ???? +? ?? ? ?? ? (?? ? ??, ð??"ð??" ? ??)?? ?? (??, ??)???? ???? +? ?? (2.11)\par
and thus leads to a reassigned spectrogram (equation (2.12)), whose value at any point (?? ? , ð??"ð??" ? ) is the sum of all the spectrogram values reassigned to this point:?? ?? (??) (?? ? , ð??"ð??" ? ; ?) = ? ?? ?? +? ?? (??, ð??"ð??"; ?)????? ? ? ?? (??; ??, ð??"ð??")??? ?ð??"ð??" ? ? ð??"ð??" ?(??; ??, ð??"ð??")? ???? ??ð??"ð??" (2.12)\par
One of the most interesting properties of this new distribution is that it also uses the phase information of the STFT, and not only its squared modulus as in the spectrogram. It uses this information from the phase spectrum to sharpen the amplitude estimates in time and frequency. This can be seen from the following expressions of the reassignment operators: (??). This leads to an efficient implementation for the reassigned spectrogram without explicitly computing the partial derivatives of phase. The reassigned spectrogram may thus be computed by using 3 STFTs, each having a different window (the window function h; the same window with a weighted time ramp t*h; the derivative of the window function h with respect to time (dh/dt)). Reassigned spectrograms are therefore very easy to implement, and do not require a drastic increase in computational complexity.?? (??; ??, ð??"ð??") = ? ??? ?? (\par
One of the most important properties of the reassignment method is that the application of the reassignment process to any distribution of Cohen's class theoretically yields perfectly localized distributions for chirp signals, frequency tones, and impulses, since the WVD does so also. As mentioned earlier, this is one of the reasons that the reassignment method can be chosen as a signal process analysis tool for analyzing LPI radar waveforms such as triangular modulated FMCW waveforms (which can be viewed as back-toback chirps) and FSK waveforms (which can be viewed as frequency tones).\par
The reassignment method provides readability improvement. The components are much better localized and very concentrated. In order to rectify the classical time-frequency analysis deficiency of poor time-frequency localization, there needs to be a method that produces more concentrated distributions, which the reassignment method does. This squeezing quality of the reassignment method lead to improved readability -which leads to more accurate metrics extractedwhich in turn, creates a more informed and safer intercept receiver environment. 
\section[{II.}]{II.} 
\section[{Methodology}]{Methodology}\par
The methodologies detailed in this section describe the processes involved in obtaining and comparing metrics between the classical time-frequency analysis technique of the spectrogram vs. the reassigned spectrogram, for the detection and characterization of low probability of intercept frequency hopping radar signals.\par
The tools used for this testing were: MATLAB (version 8.3), Signal Processing Toolbox (version 6.21), Wavelet Toolbox (version 4.13), Image Processing Toolbox (version 9.0), Time-Frequency Toolbox (version 1.0) (http://tftb.nongnu.org/).\par
All testing was accomplished on a desktop computer (Dell Precision T1700; Processor -Intel Xeon CPU E3-1226 v3 3.30GHz; Installed RAM -32.0GB; System type -64-bit operating system, x64-based processor).\par
Testing was performed for the 4 component frequency hopping waveform, whose parameters were chosen for academic validation of signal processing techniques. Due to computer processing resources they were not meant to represent real-world values. The number of samples for each test was chosen to be 512, which seemed to be the optimum size for the desktop computer. Testing was performed at three different SNR levels: 10dB, 0dB, and the lowest SNR at which the signal could be detected. The noise added was white Gaussian noise, which best reflects the thermal noise present in the IF section of an intercept receiver [PAC09]. Kaiser windowing was used, when windowing was applicable. 100 runs were performed for each test, for statistical purposes. The plots included in this paper were done at a threshold of 5\% of the maximum intensity and were linear scale (not dB) of analytic (complex) signals; the color bar represented intensity. The signal processing tools used for each task were the spectrogram and the reassigned spectrogram.\par
The 4 component frequency hopping signal (prevalent in the LPI arena [AMS09]) had the following parameters: sampling frequency=5KHz; carrier frequencies=1KHz, 1.75KHz, 0.75KHz, 1.25KHz; modulation bandwidth=1KHz; modulation period=.025sec.\par
After each particular run of each test, metrics were extracted from the time-frequency representation. The different metrics extracted were as follows:\par
Plot (processing) time: Time required for plot to be displayed.\par
Percent detection: Percent of time signal was detectedsignal was declared a detection if any portion of each of the signal components (4 signal components for frequency hopping) exceeded a set threshold (a certain percentage of the maximum intensity of the timefrequency representation).\par
Threshold percentages were determined based on visual detections of low SNR signals (lowest SNR at which the signal could be visually detected in the timefrequency representation) (see Figure \hyperref[fig_1]{1}). For visually detected low SNR plots (like this one), the percent of max intensity for the peak z-value of each of the signal components was noted (here 98\%, 78\%, 75\%, 63\%), and the lowest of these 4 values was recorded (63\%). Ten test runs were performed for both timefrequency analysis tools (spectrogram and reassigned spectrogram) for this waveform. The average of these recorded low values was determined and then assigned as the threshold for that particular time-frequency analysis tool. Note -the threshold value assigned for the spectrogram was 60\%. For percent detection determination, these threshold values were included in the time-frequency plot algorithms so that the thresholds could be applied automatically during the plotting process. From the threshold plot, the signal was declared a detection if any portion of each of the signal components was visible (see Figure \hyperref[fig_3]{2}).  Modulation bandwidth: Distance from highest frequency value of signal (at a threshold of 20\% maximum intensity) to lowest frequency value of signal (at same threshold) in Y-direction (frequency). 
\section[{Global Journal of Researches in Engineering}]{Global Journal of Researches in Engineering}\par
The threshold percentage was determined based on manual measurement of the modulation bandwidth of the signal in the time-frequency representation. This was accomplished for ten test runs of each time-frequency analysis tool (spectrogram and reassigned spectrogram), for the 4 component frequency hopping waveform. During each manual measurement, the max intensity of the high and low measuring points was recorded. The average of the max intensity values for these test runs was 20\%. This was adopted as the threshold value, and is representative of what is obtained when performing manual measurements. This 20\% threshold was also implemented for determining the modulation period and the time-frequency localization (both are described below).\par
For modulation bandwidth determination, the 20\% threshold value was included in the time-frequency plot algorithms so that the threshold could be applied automatically during the plotting process. From the threshold plot, the modulation bandwidth was manually measured (see Figure  {\ref 4}). From this threshold plot, the modulation period was measured manually from the left side of the signal (left red arrow) to the right side of the signal (right red arrow) in the x-direction (time). This was done for all 4 signal components, and the average value was determined.\par
Time-frequency localization: From Figure \hyperref[fig_6]{6}, the time-frequency localization is a manual measurement (at a threshold of 20\% maximum intensity) of the 'thickness' (in the y-direction) of the center of each of the 4 frequency hopping signal components, and then the average of the 4 values are determined. The average frequency 'thickness' is then converted to: percent of the entire y-axis. Lowest detectable SNR: The lowest SNR level at which at least a portion of each of the signal components exceeded the set threshold listed in the percent detection section above.\par
For lowest detectable SNR determination, these threshold values were included in the time-frequency plot algorithms so that the thresholds could be applied automatically during the plotting process. From the threshold plot, the signal was declared a detection if any portion of each of the signal components was visible.\par
The lowest SNR level for which the signal was declared a detection is the lowest detectable SNR (see Figure  {\ref 7}).\par
The data from all 100 runs for each test was used to produce the actual, error, and percent error for each of these metrics listed above.\par
The metrics from the spectrogram were then compared to the metrics from the reassigned spectrogram. By and large, the reassigned spectrogram outperformed the spectrogram, as will be shown in the results section.  \hyperref[tab_1]{1}, the reassigned spectrogram outperformed the spectrogram in average percent error: carrier frequency (0.74\% vs. 0.93\%), modulation bandwidth (10.82\% vs. 25.70\%), modulation period (9.30\% vs. 11.84\%), and time-frequency localization (ydirection) (4.05\% vs. 9.09\%);and in average: percent detection (86.84\% vs. 67.24\%), and lowest detectable SNR (-3.5db vs. -2.7db), while the spectrogram outperformed the reassigned spectrogram in average plot time (4.72s vs. 7.62s).\par
Figure  {\ref 8} shows comparative plots of the spectrogram vs. the reassigned spectrogram (4 component frequency hopping) at SNRs of 10dB (top), 0dB (middle), and -3dB (bottom).  IV. 
\section[{Discussion}]{Discussion}\par
This section will elaborate on the results from the previous section.\par
From Table \hyperref[tab_1]{1}, the performance of the spectrogram and the reassigned spectrogram will be summarized, including strengths, weaknesses, and generic scenarios in which each particular signal analysis tool might be used.\par
The spectrogram outperformed the reassigned spectrogram in average plot time (4.72s vs 7.62s). However, the spectrogram was outperformed by the reassigned spectrogram in every other category. The spectrogram's extreme reduction of cross-term interference is grounds for its good plot time, but at the expense of signal localization (i.e. it produces a 'thicker' signal (as is seen in Figure  {\ref 8}) -due to the trade-off between cross-term interference and signal localization). This poor signal localization ('thicker' signal), coupled with the reassigned spectrogram's 'squeezing' quality, can account for the spectrogram being outperformed by the reassigned spectrogram in the areas of: average percent error of modulation bandwidth, modulation period, time-frequency localization (y-direction), lowest detectable SNR, and percent detection. Note that average percent detection and lowest detectable SNR are both based on visual detection in the time-frequency representation. Figure  {\ref 8} clearly shows that the signals in the reassigned spectrogram plots are more readable than those in the spectrogram plots, which accounts for the reassigned spectrogram's better average percent detection and lowest detectable SNR. The spectrogram might be used in a scenario where a short plot time is necessary, but where accurate parameters are not as vital. Such a scenario might be a 'quick and dirty' check to see if a signal is present, without accurate extraction of its parameters. The reassigned spectrogram might be used in a scenario where you need accurate parameters, in a low SNR environment, in a quick time frame.\par
V. 
\section[{Conclusions}]{Conclusions}\par
Digital intercept receivers, whose main job is to detect and extract parameters from low probability of intercept radar signals, are currently moving away from Fourier-based analysis and towards classical timefrequency analysis techniques, such as the spectrogram, for the purpose of analyzing low probability of intercept radar signals. Based on the research performed for this paper (the novel direct comparison of the spectrogram versus the reassigned spectrogram for the signal analysis of low probability of intercept frequency hopping radar signals) it was shown that the reassigned spectrogram by-and-large outperformed the spectrogram in analyzing these low probability of intercept radar signals -for reasons brought out in the discussion section above. More accurate characterization metrics could well translate into saved equipment and lives.\par
Future plans include analysis of additional low probability of intercept radar waveforms, using additional time-frequency analysis and reassignment method techniques.\begin{figure}[htbp]
\noindent\textbf{2020}\includegraphics[]{image-2.png}
\caption{\label{fig_0}Year 2020 F©}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1}\includegraphics[]{image-3.png}
\caption{\label{fig_1}Figure 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-4.png}
\caption{\label{fig_2}(}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2}\includegraphics[]{image-5.png}
\caption{\label{fig_3}Figure 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3}\includegraphics[]{image-6.png}
\caption{\label{fig_4}Figure 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{45}\includegraphics[]{image-7.png}
\caption{\label{fig_5}Figure 4 :Figure 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-8.png}
\caption{\label{fig_6}Figure 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-9.png}
\caption{\label{fig_7}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{78}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Figure 7 :Figure 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.85\textwidth}}
presents the overall test metrics for the\\
signal processing analysis techniques used in this\\
testing (spectrogram versusre assigned spectrogram).\end{longtable} \par
 
\caption{\label{tab_1}Table 1}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.522787610619469\textwidth}P{0.1654867256637168\textwidth}P{0.16172566371681416\textwidth}}
carrier frequency\tabcellsep 0.93\%\tabcellsep 0.74\%\\
modulation bandwidth\tabcellsep 25.70\%\tabcellsep 10.82\%\\
modulation period\tabcellsep 11.84\%\tabcellsep 9.30\%\\
time-frequency localization-y\tabcellsep 9.09\%\tabcellsep 4.05\%\\
percent detection\tabcellsep 67.24\%\tabcellsep 86.84\%\\
lowest detectable snr\tabcellsep -2.7db\tabcellsep -3.5db\\
plottime\tabcellsep 4.72s\tabcellsep 7.62s\\
From Table\tabcellsep \tabcellsep \end{longtable} \par
 
\caption{\label{tab_2}Table 1 :}\end{figure}
 			\footnote{© 2020 Global Journals} 		 		\backmatter  			  				\begin{bibitemlist}{1}
\bibitem[Xia and Chen (1999)]{b11}\label{b11} 	 		‘A Quantitative SNR Analysis for the Pseudo Wigner-Ville Distribution’.  		 			X Xia 		,  		 			V Chen 		.  	 	 		\textit{IEEE Transactions on Signal Processing}  		October, 1999. 47  (10)  p. .  	 
\bibitem[Wei et al. (2003)]{b10}\label{b10} 	 		‘Analysis of Multicomponent LFM Signals Using Time-Frequency and The Gray-Scale Inverse Hough Transform’.  		 			G Wei 		,  		 			S Wu 		,  		 			E Mao 		.  	 	 		\textit{IEEE Workshop on Statistical Signal Processing}  		September 28 -October 1, 2003. p. .  	 
\bibitem[Conference Record of the Twenty-Eighth Asilomar Conference on Signals, Systems and Computers ()]{b8}\label{b8} 	 		\textit{Conference Record of the Twenty-Eighth Asilomar Conference on Signals, Systems and Computers},  				1994. p. .  	 
\bibitem[Pac09] Pace ()]{b6}\label{b6} 	 		\textit{Detecting and Classifying Low Probability of Intercept Radar},  		 			[ Pac09] Pace 		,  		 			P 		.  		2009. Norwood, MA: Artech House.  	 
\bibitem[Papandreou et al.]{b7}\label{b7} 	 		\textit{Detection and Estimation of Generalized Chirps Using Time-Frequency Representations},  		 			A Papandreou 		,  		 			G F Boudreaux-Bartels 		,  		 			S Kay 		.  		 	 
\bibitem[Mitra ()]{b4}\label{b4} 	 		\textit{Digital Signal Processing, A Computer-Based Approach, Second Edition},  		 			S Mitra 		.  		2001. Boston, MA: McGraw-Hill.  	 
\bibitem[Rangayyan and Krishnan ()]{b9}\label{b9} 	 		‘Feature Identification in the Time-Frequency Plane by Using the Hough-Radon Transform’.  		 			R Rangayyan 		,  		 			S Krishnan 		.  	 	 		\textit{Pattern Recognition}  		2001. 34 p. .  	 
\bibitem[Li and Bi ()]{b2}\label{b2} 	 		 			X Li 		,  		 			G Bi 		.  		\textit{A New Reassigned Time-Frequency Representation. 16 th European Signal Processing Conference},  				 (Lausanne, Switzerland)  		August 25-29, 2008. p. .  	 
\bibitem[Hlawatsch and Boudreaux-Bartels (1992)]{b0}\label{b0} 	 		‘Linear and Quadratic Time-Frequency Signal Representations’.  		 			F Hlawatsch 		,  		 			G F Boudreaux-Bartels 		.  	 	 		\textit{IEEE Signal Processing Mag}  		April 1992. 9  (2)  p. .  	 
\bibitem[Li and Xiao (2003)]{b3}\label{b3} 	 		‘Recursive Filtering Radon-Ambiguity Transform Algorithm for Detecting Multi-LFM Signals’.  		 			Y Li 		,  		 			X Xiao 		.  	 	 		\textit{Journal of Electronics (China)}  		May 2003. 20  (3)  p. .  	 
\bibitem[Ozd03] Ozdemir (2003)]{b5}\label{b5} 	 		\textit{Time-Frequency Component Analyzer. Dissertation},  		 			[ Ozd03] Ozdemir 		,  		 			A 		.  		Sept. 2003. Ankara, Turkey.  		 			Bilkent University 		 	 
\bibitem[Auger et al. ()]{b1}\label{b1} 	 		\textit{Time-Frequency Toolbox Users Manual},  		 			F Auger 		,  		 			P Flandrin 		,  		 			P Goncalves 		,  		 			O Lemoine 		.  		1996.  		 			Centre National de la Recherche Scientifique and Rice University 		 	 
\end{bibitemlist}
 			 		 	 
\end{document}