A Robust Image Watermarking Scheme Using Z-Transform, Discrete Wavelet Transform and Bidiagonal Singular Value Decomposition

Watermarking is a widely used solution to the problems of authentication and copyright protection of digital media especially for images,videos,and audio data.Chaos is one of the emerging techniques adopted in image watermarking schemes due to its intrinsic cryptographic properties.This paper proposes a new chaotic hybrid watermarking method combining Discrete Wavelet Transform(DWT),Z-transform(ZT)and Bidiagonal Singular Value Decomposition(BSVD).The original image is decomposed into 3-level DWT,and then,ZT is applied on the HH3 and HL3 sub-bands.The watermark image is encrypted using Arnold Cat Map.BSVD for the watermark and transformed original image were computed,and the watermark was embedded by modifying singular values of the host image with the singular values of the watermark image.Robustness of the proposed scheme was examined using standard test images and assessed against common signal processing and geometric attacks.Experiments indicated that the proposed method is transparent and highly robust.

Z-Transform Based Instantaneous Unit Hydrograph for Hilly Watersheds

Present study emphasizes the applicability of linear theory concept onto hilly watersheds. For this purpose, Z-transform technique was used to derive the instantaneous unit hydrograph (IUH) from the transfer function of autoregressive and moving average (ARMA) type linear difference equation. Parameters of the ARMA type rainfall-runoff process were estimated by least-squares method. The derived IUH from Z-transform (i.e. ARMA-IUH) has been used to compute the hydrologic response i.e. direct runoff hydrograph (DRH). Fur-ther, the superiority of the proposed approach has been tested by comparing the results through the results obtained from the Nash-IUH. Analyzing the results obtained from ARMA-IUH and Nash-IUH for the two hilly watersheds of North Western Himalayas shows the applicability of the linear theory concept even in turbulent flow conditions which are frequently encountered in hilly terrains under similar conditions of flow.

Analysis of Higher Order System with Impulse Exciting Functions in Z-Domain

This paper deals with mathematical modelling of impulse waveforms and impulse switching functions used in electrical engineering. Impulse switching functions are later investigated using direct and inverse z-transformation. The results make possible to present those functions as infinite series expressed in pure numerical, exponential or trigonometric forms. The main advantage of used approach is the possibility to calculate investigated variables directly in any instant of time;dynamic state can be solved with the step of sequences (T/6, T/12) that means very fast. Theoretically derived waveforms are compared with simulation worked-out results as well as results of circuit emulator LT spice which are given in the paper.

Digital implementation of fractional order PID controller and its application

A new discretization scheme is proposed for the design of a fractional order PID controller. In the design of a fractional order controller the interest is mainly focused on the s-domain, but there exists a difficult problem in the s-domain that needs to be solved, i.e. how to calculate fractional derivatives and integrals efficiently and quickly. Our scheme adopts the time domain that is well suited for Z-transform analysis and digital implementation. The main idea of the scheme is based on the definition of Grünwald-Letnicov fractional calculus. In this case some limited terms of the definition are taken so that it is much easier and faster to calculate fractional derivatives and integrals in the time domain or z-domain without loss much of the precision. Its effectiveness is illustrated by discretization of half-order fractional differential and integral operators compared with that of the analytical scheme. An example of designing fractional order digital controllers is included for illustration, in which different fractional order PID controllers are designed for the control of a nonlinear dynamic system containing one of the four different kinds of nonlinear blocks: saturation, deadzone, hysteresis, and relay.

Flexible Time Domain Averaging Technique

Time domain averaging(TDA) is essentially a comb filter,it cannot extract the specified harmonics which may be caused by some faults,such as gear eccentric.Meanwhile,TDA always suffers from period cutting error(PCE) to different extent.Several improved TDA methods have been proposed,however they cannot completely eliminate the waveform reconstruction error caused by PCE.In order to overcome the shortcomings of conventional methods,a flexible time domain averaging(FTDA) technique is established,which adapts to the analyzed signal through adjusting each harmonic of the comb filter.In this technique,the explicit form of FTDA is first constructed by frequency domain sampling.Subsequently,chirp Z-transform(CZT) is employed in the algorithm of FTDA,which can improve the calculating efficiency significantly.Since the signal is reconstructed in the continuous time domain,there is no PCE in the FTDA.To validate the effectiveness of FTDA in the signal de-noising,interpolation and harmonic reconstruction,a simulated multi-components periodic signal that corrupted by noise is processed by FTDA.The simulation results show that the FTDA is capable of recovering the periodic components from the background noise effectively.Moreover,it can improve the signal-to-noise ratio by 7.9 dB compared with conventional ones.Experiments are also carried out on gearbox test rigs with chipped tooth and eccentricity gear,respectively.It is shown that the FTDA can identify the direction and severity of the eccentricity gear,and further enhances the amplitudes of impulses by 35%.The proposed technique not only solves the problem of PCE,but also provides a useful tool for the fault symptom extraction of rotating machinery.

确定采样系统采样点间响应的因数采样法

计算机控制系统的被控对象的输出往往是连续的,采样点的输出响应不能完全反映系统的输出性能,所以研究采样控制系统采样点间响应的方法尤为重要。本文利用文献[1]的因数采样法概念,给出用于求取开环、闭环采样控制系统采样点间响应的因数采样法的推导方法及计算公式,在计算采样系统采样点间响应具有实际应用价值。

Numerical Computations of Nonlocal Schrodinger Equations on the Real Line

The numerical computation of nonlocal Schrödinger equations (SEs) on the whole real axis is considered. Based on the artifcial boundary method, we frst derive the exact artifcial nonrefecting boundary conditions. For the numerical implementation, we employ the quadrature scheme proposed in Tian and Du (SIAM J Numer Anal 51:3458-3482, 2013) to discretize the nonlocal operator, and apply the z-transform to the discrete nonlocal system in an exterior domain, and derive an exact solution expression for the discrete system. This solution expression is referred to our exact nonrefecting boundary condition and leads us to reformulate the original infnite discrete system into an equivalent fnite discrete system. Meanwhile, the trapezoidal quadrature rule is introduced to discretize the contour integral involved in exact boundary conditions. Numerical examples are fnally provided to demonstrate the efectiveness of our approach.

On Sample Size Determination When Comparing Two Independent Spearman or Kendall Coefficients

One of the most commonly used statistical methods is bivariate correlation analysis. However, it is usually the case that little or no attention is given to power and sample size considerations when planning a study in which correlation will be the primary analysis. In fact, when we reviewed studies published in clinical research journals in 2014, we found that none of the 111 articles that presented results of correlation analyses included a sample size justification. It is sometimes of interest to compare two correlation coefficients between independent groups. For example, one may wish to compare diabetics and non-diabetics in terms of the correlation of systolic blood pressure with age. Tools for performing power and sample size calculations for the comparison of two independent Pearson correlation coefficients are widely available;however, we were unable to identify any easily accessible tools for power and sample size calculations when comparing two independent Spearman rank correlation coefficients or two independent Kendall coefficients of concordance. In this article, we provide formulas and charts that can be used to calculate the sample size that is needed when testing the hypothesis that two independent Spearman or Kendall coefficients are equal.

Roll angular rate extraction based on modified spline-kernelled chirplet transform

The roll angular rate is much crucial for the guidance and control of the projectile.Yet the high-speed rotation of the projectile brings severe challenges to the direct measurement of the roll angular rate.Nevertheless,the radial magnetometer signal is modulated by the high-speed rotation,thus the roll angular rate can be achieved by extracting the instantaneous frequency of the radial magnetometer signal.The objective of this study is to find out a precise instantaneous frequency extraction method to obtain an accurate roll angular rate.To reach this goal,a modified spline-kernelled chirplet transform(MSCT)algorithm is proposed in this paper.Due to the nonlinear frequency modulation characteristics of the radial magnetometer signal,the existing time-frequency analysis methods in literature cannot obtain an excellent energy concentration in the time-frequency plane,thereby leading to a terrible instantaneous frequency extraction accuracy.However,the MSCT can overcome the problem of bad energy concentration by replacing the short-time Fourier transform operator with the Chirp Z-transform operator based on the original spline-kernelled chirplet transform.The introduction of Chirp Z-transform can improve the construction accuracy of the transform kernel.Since the construction accuracy of the transform kernel determines the concentration of time-frequency distribution,the MSCT can obtain a more precise instantaneous frequency.The performance of the MSCT was evaluated by a series of numerical simulations,high-speed turntable experiments,and real flight tests.The evaluation results show that the MSCT has an excellent ability to process the nonlinear frequency modulation signal,and can accurately extract the roll angular rate for the high spinning projectiles.

Description and reconstruction of one-dimensional photonic crystal by digital signal processing theory

A method of describing one-dimensional photonic crystals （1DPCs） based on Z-domain digital signal processing theory is presented. The analytical expression of the target band gap spectrum in the digital domain is obtained by the autocorrelation of its impulse response. The feasibility of this method is verified by reconstructing two simple 1DPC structures with a target photonic band gap obtained by the traditional transfer matrix method. This method provides an effective approach to function-guided designs of interference-based band gap structures for photonic applications.

Digital synthesis of programmable photonic integrated circuits

Programmable photonic waveguide meshes can be programmed into many different circuit topologies and thereby provide a variety of functions.Due to the complexity of the signal routing in a general mesh,a particular synthesis algorithm often only accounts for a specific function with a specific cell configuration.In this paper,we try to synthesize the programmable waveguide mesh to support multiple configurations with a more general digital signal processing platform.To show the feasibility of this technique,photonic waveguide meshes in different configurations(square,triangular and hexagonal meshes)are designed to realize optical signal interleaving with arbitrary duty cycles.The digital signal processing(DSP)approach offers an effective pathway for the establishment of a general design platform for the software-defined programmable photonic integrated circuits.The use of well-developed DSP techniques and algorithms establishes a link between optical and electrical signals and makes it convenient to realize the computer-aided design of optics–electronics hybrid systems.

Log-Time Sampling of Signals: Zeta Transform

We introduce a general framework for the log-time sampling of continuous-time signals. We define the zeta transform based on the log-time sampling scheme, where the signal x(t)is sampled at time instants tn=Tlogn,n=1,2,....The zeta transform of the log-time sampled signals can be described by a linear combination of Riemann zeta function, which firmly joins the log-time sampling process to the number theory. The instantaneous sampling frequency of the log-sampled signal equals fn=n/T,n=1,2,..., i.e. it increases linearly with the sampling number. We describe the properties of the log-sampled signals and discuss several applications in nonuniform sampling schemes.

Enhanced Frequency Resolution in Data Analysis

We present a numerical study of the resolution power of Padé Approximations to the Z-transform, compared to the Fourier transform. As signals are represented as isolated poles of the Padé Approximant to the Z-transform, resolution depends on the relative position of signal poles in the complex plane i.e. not only the difference in frequency (separation in angular position) but also the difference in the decay constant (separation in radial position) contributes to the resolution. The frequency resolution increase reported by other authors is therefore an upper limit: in the case of signals with different decay rates frequency resolution can be further increased.

Frequency Domain Convolution of Rational Transfer Functions

The convolution of two rational transfer functions is also rational, but a formula for the convolution has never been derived. This paper introduces a formula for the convolution of two rational functions in the frequency domain by two new methods. The first method involves a partial fraction expansion of the rational transfer functions where the problem gets reduced to the sum of the convolution of the partial fractions of the two functions, each of which can be solved by a new formula. Since the calculation of the roots of a high-order polynomial can be very time-consuming, we also demonstrate new methods for performing the convolution without calculating these roots or undergoing partial fraction expansion. The convolution of two rational Laplace transform denominators can be calculated using their resultant, while that of the two rational Z-transform transfer functions can be found using Newton’s identities. The numerator can be easily found by multiplying the numerator with the initial values of the power series of the result.

Continuous Forcing Spectra of Even Polygonal Chains

Let G be a graph that admits a perfect matching M.A forcing set S for a perfect matching M is a subset of M such that it is contained in no other perfect matchings of G.The cardinality of a forcing set of M with the smallest size is called the forcing number of M,denoted by f(G,M).The forcing spectrum of G is defined as:Spec(G)={f(G,M)|M is a perfect matching of G}.In this paper,by applying the Ztransformation graph(resonance graph)we show that for any polyomino with perfect matchings and any even polygonal chain,their forcing spectra are integral intervals.Further we obtain some sharp bounds on maximum and minimum forcing numbers of hexagonal chains with given number of kinks.Forcing spectra of two extremal chains are determined.

THE ENUMERATION OF SEVERAL CLASSES OF HEXAGONAL SYSTEMS

In [1], we introduced the concept of z-transformation graphs of perfect matchings of hexagonal systems and showed that the z-transformation graph of perfect matchings of a hexagonal system has at most two vertices of degree one. In this paper, we enumerate the hexagonal systems whose z-transformation graphs have a vertex of degree one. In particular, such hexagonal systems with various symmetries are also enumerated.

Oscillations of a delay differential equation with piecewise constant arguments

1 Main results Consider the differential-differenee equaionx’(t)+px(t-1)+qx([t-1]) =0, (1)where p,q∈(0,∞) and[] denotes the greatst-integer function. Recently the oscillations of eq. (1) have been discussed and several very interesting re-sults have been established. However, up to date there exists no literature on

The Global Behavior of Finite Difference-Spatial Spectral Collocation Methods for a Partial Integro-differential Equation with a Weakly Singular Kernel

The z-transform is introduced to analyze a full discretization method fora partial integro-differential equation (PIDE) with a weakly singular kernel. In thismethod, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule isconsidered. The global stability and convergence properties of complete discretizationare derived and numerical experiments are reported.