#### Abstract

Analysis of vertex-varying spectral content of signals on graphs challenges the assumption of vertex invariance and requires the introduction of vertex-frequency representations as a new tool for graph signal analysis. Local smoothness, an important parameter of vertex-varying graph signals, is introduced and defined in this paper. Basic properties of this parameter are given. By using the local smoothness, an ideal vertex-frequency distribution is introduced. The local smoothness estimation is performed based on several forms of the vertex-frequency distributions, including the graph spectrogram, the graph Rihaczek distribution, and a vertex-frequency distribution with reduced interferences. The presented theory is illustrated through numerical examples.

#### 1. Introduction

Graph signal processing is a new and quickly developing field. Many practical signals can be considered as graph signals. The theory and methods for processing the graph signals are introduced and presented in [1–5]. Graph signal processing applications in biomedical systems [6, 7] and analysis of big data [8] provide insight into the graph framework advantages and real-world potential.

In the case of large graphs, we may not be interested in the analysis of the entire graph signal, but rather interested in its local behavior. Signals with varying local vertex behaviors are a class of signals called nonstationary graph signals. One approach to the analysis of nonstationary graph signals is vertex-frequency analysis [7, 9–15], which is a counterpart of time-frequency analysis [16–18] in classic signal processing.

The main representatives of the vertex-frequency representations are local vertex spectrum and its energetic version, graph spectrogram. Window functions are used to localize graph signals in a neighborhood of the considered vertex [9, 12, 15].

Another important class of vertex-frequency representations, called the vertex-frequency energy distributions, were recently introduced in [13, 14]. This class is a counterpart to the class of quadratic time-frequency distributions in classic signal analysis. It has been shown that the graph version of the Rihaczek distribution is of special interest for graph signals since it does not require a localization window. The reduced interference distributions can be derived from the Rihaczek distribution by using appropriate kernel functions. This class of representations, under certain conditions, satisfies marginal properties in both the vertex domain and the spectral domain.

An important concept that is used in classic time domain signal analysis for the description of local signal behavior around a time instant is the instantaneous frequency. The local smoothness is introduced in this paper as an extension of the instantaneous frequency concept to graph signal analysis. The local smoothness is defined by using the graph signal Laplacian matrix. The vertex-frequency representations can be highly concentrated along the local spectral index, corresponding to the local signal smoothness. This property is used to define local smoothness estimators based on the vertex-frequency representations.

After an introduction, we will review the fundamental theory of graph signal processing. This review will include the graph Fourier transform and the global signal smoothness in Section 2. Then, the local signal smoothness will be introduced and its properties derived within Section 3. The vertex-frequency representations, along with their connections to the local signal smoothness, will be presented in Section 4. The theory will be illustrated through a nonstationary graph signal example.

#### 2. Graph Signals

A graph is defined as a set of vertices and a set of edges connecting these vertices. In signal processing, such a structure can be considered as the domain of a signal. The signal values are defined at the graph vertices. The graph Fourier transform (graph spectrum) is defined through the eigenvalue decomposition of the graph Laplacian matrix. Here, we will present a review of the graph spectrum and the global signal smoothness index calculated using the Laplacian matrix [10].

##### 2.1. Graph Signal and Spectrum

A weighted undirected graph with vertices will be considered. The edge weights are nonzero if there is an edge between the vertices and . If there is no edge between the vertices and , the corresponding weight is equal to zero, . The weight matrix is a matrix whose elements are . It is a symmetric matrix (since the underlying graph is undirected), with zeros on the main diagonal.

The definition of the graph Laplacian, using the weight matrix and its elements , is given bywhere is a diagonal matrix, called the degree matrix. Its diagonal elements are obtained from as , while for .

The Laplacian matrix, like any other quadratic matrix, can be written using its eigenvectors and eigenvalues asIn this decomposition, the matrix consists of the matrix eigenvectors, denoted by , as its columns. The diagonal matrix of eigenvalues , , is denoted by . The eigenvectors and eigenvalues of are calculated from . Here we will consider the case with simple eigenvalues, whose multiplicity is one.

Graph signal samples, , , are sensed/defined at each graph vertex . These signal samples can be written in vector form as an vector:

The graph discrete Fourier transform (GDFT) of a signal is defined by [10]The coefficients in the GDFT are calculated as the projections of the considered graph signal to the eigenvectors

The inverse graph discrete Fourier transform (IGDFT) follows from the property that holds for the Laplacian matrix eigenvectors, where is an identity matrix. The IGDFT relation is , with

The GDFT concept can be extended to the directed graphs. The cases of repeated eigenvalues can also easily be included in the analysis [19–21].

##### 2.2. Global Graph Signal Smoothness

In classic signal analysis, when the signal domain is time, the signal smoothness can be defined through a second-order difference . Since classic time domain signal processing can be considered as graph signal processing on a circular graph, the second-order difference can be written as , where is the Laplacian of the circular graph. The signal smoothness can be measured as cumulative energy of the signal changes . It can also be calculated as . In matrix notation, we get the quadratic form . This approach can be extended to general (non-circular) graphs.

From the Laplacian eigendecomposition, we haveorsince for an eigenvector holds . For an arbitrary eigenvector and the corresponding eigenvalue , here we omitted index for notation simplicity. The quadratic form of an eigenvector is equal to the corresponding eigenvalue. This quadratic form can be used as a measure of the signal smoothness. We can write the quadratic form as Since , the last relation can also be rewritten asThe sum of the previous two relations producesObviously, a small value of corresponds to slow eigenvector variations , within the neighboring/connected vertices. This means that the eigenvectors calculated with a small represent a low-pass (slow-varying) part of the graph signal.

Since the eigenvalues of the Laplacian matrix are equal to the quadratic form , they are nonnegative. It is known that at least one eigenvalue of the Laplacian is zero. The corresponding eigenvector is constant, i.e., maximally smooth signal.

The graph signal smoothness is defined, in general, using a full analogy with (8). By normalizing the quadratic form,With the signal energy, the smoothness index definition is obtained as

An example of the time domain signals and graph signals, with various values of the global smoothness , is presented in Figure 1. It is obvious that small values correspond to the smooth (slow-varying) signals and that large values of indicate fast-varying signals.

**(a) Constant signal**

**(b) Slow varying signal**

**(c) Fast varying signal**

Now consider the signal whose form is given by a weighted sum of the eigenvectors, The global smoothness of this signal is It is obvious that , where and .

The smoothness of graph signals is used in graph topology learning [22], vertex ordering, and graph clustering [23]. Since , corresponding to , is constant, the vertex ordering can be done using the next smoothest eigenvector (called the Fiedler vector). The vertices are ordered according to the indices of the sorted values. Regions with similar values can be used for the graph clustering.

#### 3. Local Graph Signal Smoothness

The local graph signal smoothness is introduced next. Its properties are analyzed in the second part of this section.

##### 3.1. Local Graph Signal Smoothness Definition

Assume the simplest case, when the analyzed signal is proportional to the th Laplacian eigenvector,In a full analogy to classic spectral analysis, we can say that the signal of this form is a monocomponent signal, since its spectrum has only one nonzero coefficient at the th position. We can define the spectral index of this component, or its smoothness index, asIt is equal to the corresponding eigenvalue.

The smoothness index can be related to the frequency in the time domain signal analysis [10]. The classic Fourier analysis may be obtained as a special case of the GDFT on a circular undirected graph. For this graph, the eigenvectors are periodic functions with frequencies . The smoothness index is obtained from asWe can conclude that the eigenvalue corresponds to the squared classic signal analysis frequency . If continuous-time is considered, instead of discrete-time, or the case with a small is considered in the discrete-time domain (18), we would get

For the time domain signals with a time-varying spectrum, the concept of instantaneous frequency is introduced. Several approaches to the instantaneous frequency exist [16–18]. In general, for a signal with varying frequency, we can define instantaneous frequency by considering the signal behavior in the vicinity of the considered time instant . If the signal form at the instant and its small neighborhood is close to the form of a sinusoidal signal with frequency , then we can say that the instantaneous frequency of the considered signal, at the considered time instant , is equal to . In this case, the frequency can be estimated by using a few samples around the considered time instant [24]. Another method to find the instantaneous frequency is to approximate the signal by a second-order polynomial around ,If we compare this signal with a sinusoidal signal expansion at the instant , for a small , we can conclude that, for , the sinusoidal signal that fits the signal defined by (19), around the considered time instant , has the frequency such thatNow we can conclude that the instantaneous frequency of the considered signal at a time instant is . If , we can use the ratio of higher-order derivatives in order to obtain the signal’s instantaneous frequency (assuming that ).

The discrete-time definition of the squared instantaneous frequency iswhere is the second-order difference of the considered signal.

In the previous section, we show that the second-order difference of a time domain signal corresponds to the elements of , where is the Laplacian of a circular graph. An example of a signal with time-varying smoothness is presented in Figure 2. In the first part, , the signal is slow-varying (with a small local smoothness), then a fast varying part of the signal follows, and in the last part, the signal is moderately smooth.

In analogy with (22), we will introduce the local smoothness for a signal defined on an arbitrary graph as

We have assumed that .

##### 3.2. Properties of the Local Smoothness

Some of the properties of the local smoothness are described next.(1)Consider a monocomponent signal Its local smoothness is vertex independent. This smoothness is equal to the global smoothness since In the time domain signal analysis, this property means that the instantaneous frequency of a sinusoidal signal is equal to its frequency.(2)Assume a piecewise monocomponent signal where are subsets of the vertices such that for , and each vertex belongs to a subset . Within each subset, the considered signal is proportional to the eigenvector . For each interior vertex , i.e., a vertex whose neighborhood lies in the same set , the local smoothness is An example of a piecewise monocomponent graph signal is presented in Figure 3. Three subsets of vertices , , and are considered. They are marked in colors in Figure 3. The component spectral indices are , , and . For subset , the boundary vertices are 1, 4, 6, 17, and 22. For subset , the boundary vertices are 23, 24, 29, and 34. For subset , the boundary vertices are 35, 39, 53, 62, and 64. All other vertices are interior vertices. The local smoothness of the piecewise monocomponent graph signal from Figure 3 is calculated and presented in Figure 4. The obtained results are exact for each interior vertex (presented with dots in Figure 4). For the boundary vertices, the results are not exact since we include samples from all neighboring vertices in the local smoothness calculation. Some of them are outside the considered set . The results for the boundary vertices are indicated by the cross marks.(3)An ideal vertex-frequency distribution can be defined as It has been assumed that the local smoothness is rounded to the nearest eigenvalue. For the graph and the signal presented in Figure 3, the ideal vertex-frequency distribution is shown in Figure 5. This distribution can be used as a local smoothness estimator since, for each vertex , the maximum of is positioned at . The index of the eigenvalue that corresponds to the local smoothness is obtained as

#### 4. Vertex-Frequency Representations

The energy vertex-frequency distributions follow the concept of the time-frequency energy distributions in classic signal analysis. The estimation of the local smoothness can be obtained by using the vertex-frequency representations that localize the graph signal energy on the local smoothness. Here we will present the vertex-frequency energy distribution, a reduced interference vertex-frequency distribution, and the graph signal spectrogram, as the tools for local smoothness estimation.

##### 4.1. Energy Vertex-Frequency Distributions

The energy of a signal is commonly defined as The signal can be written as , where is the GDFT of the signal. The signal energy is now where the distribution of the signal energy in the vertex-frequency domain isThis distribution corresponds to the Rihaczek distribution in classic time-frequency analysis.

A vertex-frequency distribution satisfies the marginal properties, if The marginal properties state that the signal power can be obtained by a summation of over and that the squared signal spectrum can be obtained by a summation of over .

We will show that the vertex-frequency distribution defined by (36) satisfies the local smoothness property (32)since is the inverse GDFT of . In a matrix form, it is equal to

For the vertex-frequency distribution defined by (36), the local smoothness bandwidth (33) may be written in terms of , , and , since corresponds to the elements of .

*Example*. The distribution of the graph signal from Figure 3 is illustrated in Figure 6. The marginal properties (sums over and over ) are presented below and right of the distribution image. Both marginal properties are satisfied, as expected. It is important to note that this distribution does not use a localization window. From the vertex-frequency representation, we can identify the signal components and the cross-terms. The cross-terms, well known in classic time-frequency analysis, are produced by mixing the signal components in the calculation of the distribution values . The third signal component of the signal analyzed in Figure 6 exists at vertices only, and the distribution is nonzero for lower vertex indices at . Also, there is no signal component at , but is obviously not equal to zero.

##### 4.2. Vertex-Frequency Distributions with Reduced Interference

In order to reduce the cross-terms interferences and to preserve the marginal properties, a general class of reduced interference time-frequency distributions is extended to the graph signals [14]. The frequency domain definition of the reduced interference energy distribution iswhere is a kernel function. For , the graph Rihaczek distribution (36) follows. The exponential kernel, a counterpart to the Choi-Williams kernel in classic time-frequency analysis, is defined aswherefor and .

The reduced interference vertex-frequency distribution is presented in Figure 7. Here we have used the exponential kernel. It is notable that the cross-terms are reduced as compared to Figure 6, while the marginal properties are preserved in this case.

Now we will consider a general case and review the conditions that the distribution kernel should satisfy in order to preserve the marginal properties.

A sum of all values should be equal to the signal energyThis relation is satisfied if

The vertex marginal property of the distribution is satisfied if sinceNote that the eigenvectors are orthonormal, producing

Moreover, if this condition is satisfied, then the vertex moment property holds

The frequency marginal property holds if A sum of over the vertex index is

If the frequency marginal property holds, then the frequency moment property holds as well,

The local smoothness property (32) of is satisfied if This can be written asThe local smoothness property is satisfied if

The reduced interference distributions can be used as estimators of the local smoothness. The local smoothness is estimated as the eigenvalue that corresponds to the position of the maximum in , for a considered vertex ,

The reduced interference distribution , along with the marginal properties, is presented in Figure 7. Performance of the distribution as a local smoothness estimator will be illustrated through an example at the end of this section.

##### 4.3. Vertex-Frequency Spectrogram

In the classic time-frequency analysis, the short-time Fourier transform and the spectrogram are well-developed tools for analysis of nonstationary signals. Their extension to the graph signals leads to the vertex-frequency spectrogram. It can be calculated as the spectrum of a signal multiplied by an appropriate localization window function The window function should localize the signal content around the vertex . In general, it is vertex-dependent, in contrast to the classic time domain spectrogram, where commonly the same window (with a shift in time) is used.

In a special case, when , the localized vertex spectrum is equal to the standard spectrum for each ; that is, no vertex localization is performed. The second special case is a maximally localized windowThe localized vertex spectrum, in this case, is equal to the signal , for each , and we do not have any spectral resolution.

The spectrogram of a graph signal is defined as

The vertex marginal of the spectrogram iswhere Parseval’s theorem is used. It is obvious that the vertex marginal property is not satisfied for a general localization window . Only for a very specific case when , the vertex marginal is equal to the signal energy .

A summation over and should produce the total signal energy. For the vertex spectrogram, we getIf the localization windows are such that holds for all , then the vertex spectrogram is energy unbiased

The localization windows could be defined in the spectral domain using a generalized graph convolution [9] or in the vertex domain using the vertex neighborhood [12].(i)The localization window, defined in the spectral domain, is equal to [9] where is a window basic function defined in the spectral domain, for example, as where is the amplitude of the window and is a constant that defines the width of the window. The vertex-frequency spectrogram calculated with localization windows defined in the spectral domain is presented in Figure 8. It is obvious that the marginal properties are not satisfied in this case.(ii)The localization window can be defined in the vertex domain. The window function value depends on the distance between the vertices and where is a form of the basic window that corresponds to the classic signal processing window form. Here we will review a method for obtaining the localization window functions, at each vertex, in matrix form [10]. The vertices whose distance is follow from matrix . This matrix is equal to the graph adjacency matrix. The matrix , with elements , is obtained from the weighting matrix The vertices whose distance is follow from the matrix We have used the following notation: for the logical (Boolean) matrix product, for the Hadamard product (element-by-element multiplication), and for the matrix whose all elements are equal to . The elements of the matrix are equal to 1 if the distance between vertices and is 2, and 0 otherwise. Matrix gives the information about all vertices that are connected with walks of the length and a lower walk. The element-by-element multiplication by the matrix removes the vertices connected with walks of length , while the multiplication by removes the diagonal elements. When , a recursive relation for the matrix follows. It will give the information about the vertices at a distance The matrix for the graph localization windows is formed as The graph signal weighted by the localization window is calculated by using the previous matrix as An example of the vertex-frequency spectrogram calculated by the vertex domain localization window is presented in Figure 9.

The relation that would connect the vertex domain spectrogram (56) and a general vertex-frequency distribution (40) is very complex. In order to establish this relation, the general vertex-frequency distribution should be rewritten using a kernel function in the vertex-vertex shift domain. This form is dual to (40)where is the vertex-vertex shift kernel. The conditions for the frequency marginal and the vertex marginal with a vertex-vertex shift kernel are

The kernel that corresponds to the vertex domain spectrogram (56) is For localization windows defined in the spectral domain, the kernel function can be written as This kernel cannot satisfy both marginal properties. The unbiased energy condition is

Here, the local smoothness property cannot be satisfied. However, since the graph spectrogram is concentrated along the local smoothness, we can still use the maximum-based estimator of the local smoothness

#### 5. Numerical Results

Now we will calculate the local signal smoothness by using the Laplacian applied to the graph signal (23) as well as the maximum positions of the Rihaczek distribution, the reduced interference distribution, and the local vertex spectrum with the spectral domain and the vertex domain windows. We will use the signal presented in Figure 3. The spectral domain window, defined by (62) and (63), is calculated with and . The vertex domain window is calculated according to (68), with and , i.e., , , , and .

The obtained results are presented in Figure 10 and in the first row of Table 1. For all considered estimators, the local smoothness is estimated at the vertices where significant signal sample values are detected, . The theoretical value of the local smoothness is presented by a line and the estimated values are presented by dots. The number of outliers (indicated by NO in Table 1) is the number of vertices where the estimated smoothness is not equal to the theoretical one. The mean squared error (MSE) of the local smoothness estimations is also given in Table 1, for each considered case.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

Next, we will consider a noisy signal. The signal is corrupted by a Gaussian noise added to the signal samples. Signal-to-noise ratio (SNR) is varied from db to 50dB. The number of outliers and the MSE are given in Table 1. We can conclude that the direct method of local smoothness estimation using the Laplacian is very sensitive to the noise, while the vertex-frequency based estimations are robust to the noise. For the local vertex spectrum, the number of outliers is slightly increased for the high SNR cases. This increase is caused by the signal smoothing. Note that the Rihaczek distribution and the reduced interference distribution provide better estimation for a high SNR, while for the SNR below 0dB, the local vertex spectrum, calculated with either the spectral or the vertex domain window, results in fewer outliers and a lower MSE.

The ratio of the concentration measures is calculated using the -norm of the Rihaczek distribution and the -norm of the reduced interference distributions [25]. It varied from in the non-noisy case to for the highest noise level.

#### 6. Conclusion

In this paper, the local smoothness of graph signals is introduced and analyzed. Methods for local smoothness estimation, based on the signal Laplacian and the vertex-frequency representations, are given and applied to examples with graph signals. It has been shown that the local smoothness is a counterpart of the instantaneous frequency in classic signal analysis and can be estimated using vertex-frequency distributions. Finally, the vertex-frequency energy distributions, including a reduced interference distribution, and the local vertex spectrogram with two windowing techniques, are considered as the local smoothness estimators.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.