IIR filters analysis
1. Introduction
Analysis methods based on determining system performance specifications by step response, as well as indirect methods: polezero plot, magnitude response and integral analysis methods are applied in automatic control theory for performance estimation of linear systems [1,2,3]. However, in many cases the mentioned methods result in crude performance estimation of a linear system (filter) operation. Furthermore, direct methods of linear system performance specifications (settling time, accuracy, overshoot etc.) characterization require a huge amount of calculations being performed.
Specification or estimation of signal processing performance criteria are usual tasks in frequency filter analysis. In some cases it is considered to be enough to examine a filter behavior at average statistical parameters of a useful signal and its disturbance. In other cases, for instance for robust filters, it is rather more complicated – one needs to determine the limit of variables for signal processing performance specifications at any possible input signal parameters variation.
The author offers to use the filter analysis methods, developed by him on the basis of spectral representations of the Laplace transform, to solve efficiently the problem of signal processing performance specifications determination by frequency filters at different variations of input signal parameter [4,5,6]. The mentioned methods are based on using consistent mathematical models for input signals and filter impulse characteristics by means of a set of continuous/discrete semiinfinite or finite damped oscillatory components. Similar models can be applied for simple semiinfinite harmonic and aperiodic signals or filter impulse characteristics, compound signals of any form, including signals with composite envelopes, as well as pulse signals (radio and video pulses).
The application of signal/filter frequency and frequencytime representations, based on Laplace transform, allowed developing simple and effective direct methods for performance analysis of signal processing by analog and digital filters.
To simplify the task of analog and digital filter signal processing performance analysis the author offers two methods for performance expressanalysis of signal processing by frequency filters using filter frequency responses based on Laplace transform: frequency and frequencytime analysis methods [7].
The frequency method from an indirect analysis method for signal processing, in fact, has transformed into a direct analysis method by means of Laplace transform spectral representations. This method is the most effective in cases, where only two main performance specifications: signal processing speed and accuracy – are required to be evaluated.
The frequencytime analysis method is being applied in cases, where there is a need to evaluate signal processing speed and accuracy, as well as the history of transient processes in a filter, for instance, to control oscillation of transient process in a filter. It is suggested to perform the analysis by using frequency responses based on filter transfer function, dependent on time, in sections of input signal complex frequencies.
In case of FIR filter an effective estimation of signal processing performance specifications can be carried out by using filter frequency response 3D analysis based on Laplace transform in sections of input signal complex frequencies, considering their change. To evaluate signal processing performance specifications for IIR filters one will need along with sections of 3D filter frequency responses to use sections of 3D signal spectrum on filter impulse response complex frequencies [5].
The issues about the application of the analog and digital filter analysis methods, developed by the author for signal processing performance analysis by frequency filters, are considered below.
2. IIR filters analysis
2.1. Signal processing performance analysis by analog IIR filters
Let us consider signal processing performance analysis by IIR filters at semiinfinite input signals on the basis of analysis methods based on the Laplace transform spectral representations.
Three methods of frequency filter analysis are suggested by the author for the timeandfrequency representations positions of signals and linear systems in coordinates of complex frequency [5,6]. Let us consider the first two methods for signal processing performance analysis by frequency filters.
The mathematical description for the generalized input signal and IIR filter at time and frequency domains for the first (item2) and the second (item 3) analysis methods, mathematical expressions for calculating forced and free components of a filter reaction by the first method (items 4,5), components of a filter reaction by the second method (item6) are given in the Table 1.
1.  Input signal  
2.  Filter  
3.  Time dependent transfer function 

4.  Forced components  
5.  Free components  
6.  Filter reaction 
The operation of the real part extraction on the right side of the expression in the items 1,2,3 for
The first method is a complex amplitude method generalization for definition of forced and free components for filter reaction at semiinfinite or finite input signals [6]. The advantages of this method are related to simple algebraic operations, which are used for determining the parameters of linear system reaction (filter, linear circuit) components to an input action described by a set of semiinfinite or finite damped oscillatory components. To analyze a filter it is needed to use simple algebraic operations and operate a set of complex amplitudes and frequencies of forced and free filter reaction components. In this case, there are simple relations between complex amplitudes of output signal forced components and complex amplitudes of an input signal (item 4 Table 1), between complex amplitudes of output signal free components and complex amplitudes of a filter impulse function (item5).
The timeandfrequency approach in the second analysis method applies to a filter transfer function, i.e. time dependent transfer function of the filter is used [6,8]. In that case, instead of two sets of filter reaction components only one of them may be used.
Analysis methods given in the Table 1 enable to reduce effectively the computational costs when performing a filter analysis by using simple algebraic operations to determine the forced and free components of a filter reaction to an input action as a set of damped oscillatory components. Therefore, the considered analysis methods for linear systems (filters) can be effectively applied for performance analysis of signal processing by frequency filters.
Let us consider a simple example of performance analysis of signal processing by a highpass secondorder filter relating to signal processing task in power system protection and automation devices [9,10]. The filter is used to extract a sinusoidal component of commercial frequency and eliminate disturbance as a free component of transient processes in a control object. In this case, a change of a useful signal initial phase is acceptable.
All the initial data and dependencies which are necessary for the analysis are represented in the Table 2. IIR filter parameters are specified, the mathematical description of an input signal with specified sizes of changing for useful signal and disturbance parameters affecting their spectrum is given in the Table 2 as well.
The impulse function of highpass secondorder filter contains a delta function of Dirac which is used for determining complex amplitudes of forced components when defining
The analysis results should ensure the following performance criteria of signal processing by a filter:
a filter settling time should be less than 30 ms at 5% acceptable total error of signal processing at any value of disturbance parameters within the specified range,
an acceptable error at frequency deviation of useful signal from the nominal value of 50Hz within the range ±5 Hz should not be more than 5%,
an acceptable overshoot should not be more than 10%.
As it follows from the Table 1, simple algebraic operations are applied to determine complex amplitudes, as well as forced and free components of a filter output signal.
When using Mathematica, Mapple, Matlab, Mathcad and other stateofart mathematical software for determining forced and free components of an output signal it is necessary to specify only complex amplitude vectors of an input signal and a filter impulse function, as well as complex frequency vectors of an input signal and a filter. In this case all the necessary calculations, related to a filter analysis, would be carried out automatically. If it is needed to determine complex amplitudes of a filter impulse function at specified transfer function the readymade formulas may be used [6], which can be easily applied in the mathematical software mentioned above.
1.  Input signal  
2.  Highpass filter  
3.  Forced components of complex amplitudes  
4.  Forced components  
5.  Complex amplitudes of free components  
6.  Free components  
7.  Error 
All the examples in the present chapter are given using the mathematical software Mathcad. Mathcad was chosen due to pragmatic considerations related to assuring the maximum visibility of the examples for filter analysis, as in Mathcad mathematical expressions are given in the form, closest to universally accepted mathematical notation [11,12].
An example of a filter computation using Mathcad at the specified filter parameters and the following input signal parameters:
For determining or estimating performance specifications of signal processing by the investigated IIR filter one need either to improve the software (Figure 1) or to reduce the amount of calculations by simplifying the analysis task. Let us consider the second option first.
The easiest operation is to define the error level in signal processing by a filter at frequency deviation of useful sinusoidal signal within the range ±5 Hz from the nominal value of 50 Hz. This error, as it is known, may be determine by an average amplitudefrequency response of a filter. In this case the value of a filter amplitudefrequency response in the areas of frequency
A filter settling time can be defined by total damping of a free component
A damping time of disturbance free component
Variables
A presice settling time
Let us consider an estimation of total damping
Dependencies
In case of
An overshoot level in a filter can be determined by the program improvement on the Mathcad example, shown on the Figure 1.
The performance analysis results for signal processing of the investigated IIR filter are represented in the Table 3. Performance specifications determined by using one traditional method  by step response of a secondorder lowpass filter are given in the Table 3 as well. In this case the description for the lowpass filter was obtained on the basis of the investigated secondorder lowpass filter by applying a well known frequency transformation [14].
1.  Settling time, s  0,0172  0,0275 
2.  Maximum error in the steadystate mode, %  0  5 
3.  Maximum overshoot level, %  1,17  11,81 
4.  Additional error at a frequency variation of an useful signal    5 
As it follows from the Table 3, in the considered example there are some substantial differences in the performance specifications, gained by the traditional method and on the basis of their direct determination.
2.2. IIR filters analysis at dissemination of timeandfrequency approach to transfer function filter
As it follows from the Table 1 (item 3 and item 6), estimation of functioning performance for filters can be performed at dissemination of timeandfrequency approach to filter transfer function, in other words by using time dependent transfer function of a filter [5,6].
In that case, comparing to the first method where two groups of filter reaction components are using – forced and free components, only the first group of components is using. The information about the transient process is containing in the time dependent complex amplitudes of a filter reaction
The necessary dependences for determination of output signal components of a filter are represented in the Table 1, an example for highpass filter analysis at the specified parameters of an useful signal is given on the bottom part of the fig.1.
The advantage of the considered method is connected to determination envelopes for every component of a filter output signal, based on which the total envelope of a filter output signal and a variation law of initial phase can be defined. This information can be effectively used for performance analysis of signal processing by frequency filters. For instance, when determining an overshoot level (oscillativity) of a transient process in a filter.
2.3. IIR filter expressanalysis method
It follows from the Table 1, that quality indexes estimation for filter operation can be carried out on the basis of interim calculation results –
The expressanalysis methods for filters, including performance analysis of signal processing, were developed based on investigation of 3D and 4D frequency responses [7]. It is enough to consider the sections
The expressanalysis methods mentioned above can be effectively applied for FIR filters as well, the detailed explanation will be given further in the present chapter.
2.4. Digital IIR filters analysis
Under the definition of digital filters in the chapter discrete filters are ment. In many cases it is justified, for instance, in cases of using microcontrollers or digital signal processor with high digit capacity and especially for microprocessors with support for floatingpoint operations [15,16].
When using discrete filters their analysis has a lot of similarities with the analysis of analog filtersprototypes. There is a small difference only when it comes to transition from images to originals. The main expressions for determining components of a digital filter output signal when injecting on the filter input a signal as a set of discrete semiinfinite damped oscillatory components are given in the Table 4.
An example for digital filter analysis as a continuation of the example of the analog filterprototype analysis (fig. 1) is represented on the fig. 3. The mathematical description for the digital filter was obtained by the method of invariant impulse responses at the discrete sampling step
1.  Input signal  
2.  Filter  
3.  Time dependent transfer function  
4.  Forced components  
5.  Free components  
6.  Filter reaction 
1.  Input signal  
2.  Filter  
3.  Time dependent transfer function  
4.  Forced components  
5.  Free components  
6  Filter reaction 
2.5. IIR filters analysis at finite signals
Performance analysis of processing finite signals by IIR filters as a set of damped oscillatory components with finite duration may be performed on the basis of dependencies for IIR filters at semiinfinite signals [5,6].
All the needed expressions were obtained on the basis of the expressions from the Table 1 using a time shift and the principle of additivity.
Let us consider the IIR filter analysis at compound input signals as a set of sequentially adjacent finite signals [17]. The calculation of IIR filter reaction for this case is represented in the Table 5.
In this case every component of an input signal in a general way has a different shift and a different duration, where
The expressions represented in the Table 5 can be significantly simplified, if the IIR filter analysis at a finite signal, for instance, at injection a finite signal on its input from
An example of filter calculation, analogous to the example on the fig.1, but using an input signal with finite duration is given on the fig.4.
The dependences given in the present section can be effectively used for not only analysis of passing through IIR filter one or another finite signal, but also for performance analysis of signal processing by filters.
For digital IIR filter analysis at injection finite signals on their inputs analogous mathematical operations are applied. All the necessary dependences may be obtained on the basis of the formulas from the Table 4 [13].
3. FIR filters performance analysis
3.1. Particularities of the analysis
Mathematical description for FIR filters can be obtained on the basis of the IIR filter description (Table 1) by using twice as many of filter impulse function components.
1.  Input signal  
2.  FIR filter  
3.  Forced components  
4.  Free components  
5.  Filter reaction  
6.  Time dependent transfer function  
7.  Filter reaction 
The additional components have the same set of complex frequencies and differ by time shift and values of complex amplitudes in a way to ensure the finitude of a filter impulse characteristic [5,6]. According to this approach to mathematical description, IIR filters are special cases of FIR filters.
The inputoutput dependences for FIR filters can be obtained on the basis of analogous dependences of FIR filters by using time shift and principle of additivity operations [6].
Comparing to IIR filters, FIR filters have finite duration of transient processes of their own, which are defined by a filter length. This to a certain extent simplifies performance analysis for signal processing of this type of filters, especially when using the suggested expressanalysis methods for signal processing performance by FIR filters.
3.2. Analog FIR filters
Basic expressions for FIR filter analysis at injection on the filter input a set of semiinfinite damped oscillatory components are given in the Table 6.
Due to characteristics of FIR filters, among forced and free components in an filter output signal there is the third group of components
1.  Input signal  
2.  FIR filter: impulse function, time dependent transfer function 

3.  Filter reaction 
Basic expressions for FIR filter analysis at signal injection on a filter input as a set of damped oscillatory components with finite duration are represented in the Table 7. The calculation of a filter reaction is given in the Table 7 using only the second analysis method for the case when duration of all the components of a filter impulse function is equal
Let us consider an analysis example of FIR filters which are used in one of the most perspective intelligent electronic devices (IED)  Phasor Measurement Units (PMU) [19].
A brief description of a basic algorithm for PMU signal processing on the example of an analog systemprototype is given in the item 1 of the Table 8 [20].
An input signal of intelligent electronic devices is represented by a set of complex amplitudes
An average FIR filter
The considered filters should have low sensitivity to a change of damping coefficient
1.  Block scheme of algorithm 

2.  Signal description 
Semiinfinite signal Signal as a set of finite components 
3.  Input signal of a filter 
At semiinfinite input signal device 
4.  Algorithm  
5.  Average FIR filter 
FIR filter analysis is performed at input signal of a device as a set of semiinfinite or finite damped oscillatory components according to the algebraic expressions from the Table 7 and the Table 8.
An example of FIR filter analysis using Mathcad at compound input signals as a set of sequentially adjacent finite signals is given on the fig.6.
Each of five sets of finite components represented comply with a particular power system regime. The first set is related to a normal regime of a power system and is represented by sections of sinusoidal component of industrial frequency and higher harmonics. The second set of finite signals corresponds to an accidental regime and contains of finite sinusoidal and exponential component, the third set shows nocurrent regime, the fourth set is connected to a change of PMU input signals envelope due to automation operation, the fifth set represents a swings regime.
A filter processes a complex signal
A plot for PMU input signal
It appears from the plot that the investigated filter has a feature which is connected to the absence of overshoot (oscillation) of transient process in the filter in its traditional definition, in other words at a stepwise growth of a signal. and presence of overshoot at a stepwise reduction of signal.
In fact, due to particular characteristics of impulse function, from the traditional point of view, an overshoot is absent in the considered FIR filter, as at the end of transient process in the filter signal behavior is close to aperiodic process. In this case oscillation is noticeable on the initial stage of transient process.
Analysis of FIR filter and signal processing algorithm in total is carried out only at fixed values of input signal parametres in the example. With some further ordinary improvement, analog to the example on the fig.1, it is possible to determine performance specifications for FIR filter signal processing at any variation of useful signal and disturbance parameters. However, it is much easier to use the specific expressanalysis methods.
3.3. Expressanalysis methods for signal processing performance
Let us consider an example of expressanalysis of signal processing performance of a FIR filter, which mathematical description is given in the Table 7. In addition, let us consider performance analysis for signal processing of a filter with the impulse function
To check the adherence to the mentioned conditions, it is enough to consider amplitudefrequency response of the filters in complex frequency coordinates.
The amplitudefrequency response for the first and second filters is given in the Fig.8. The multiplier
To check the compliance with the requirements for filters to the specified quality indexes of signal processing for the specified input signal, it would be enough to consider two sections of 3D amplitudefrequency response: in the sections
The filters 1 and 2 in the section of 3D frequency response
In the case of a need along with estimating the speed and accuracy of signal processing to appraise the history of transient processes in a filter, for instance, to control oscillation of transient process in a filter, it is eligible to conduct the analysis using amplitudefrequency response based on the filter transfer function, which depends on time.
It is necessary for analysis to depress a damped oscillatory component with the complex frequency
The graph presented describes graphically the fact, that the second filter does not ensure the suppressing of an input signal component with the complex frequency
3.3. Digital FIR filters
As for IIR filters (Table 1 and 4), inputoutput dependencies for digital (discrete) FIR filter may be obtained by discretization of expressions for analog FIR filter and transition from the Laplace transform to Z transform [13].
In case of using the mathematical software Matlab, Mathcad etc. for digital FIR filter analysis it is enough to get the data about complex amplitudes and input signal frequencies
In some cases only mathematical description of analog filterprototypeis specified. In case of FIR filters, different methods of transition from analog filterprototype description to digital filter description, for instance, method of differential equation discretization, method of invariant impulseresponses, bilinear transformation are applied [14]. The mentioned methods expand also for the case of FIR filters in the author’s paper [13].
The transfer function of analog filterprototype with finite impulse response according to the Table 6 can be described in the following way
The expression for a digital filter transfer function
where
All the mentioned constants, except
1.  Method of differential equation discretization  
2.  Method of invariant impulse responses  
3.  Bilinear transformation 
The additional indexes of the constants in the Table 9, corresponding to order numbers of transition methods, are given for the constants, which do not coincide with the analog filterprototype parameters.
4. Performance expressanalysis at modulated signal
4.1. Mathematical description of signals
The signals as a set of semiinfinite or finite damped oscillatory components were considered above. Compound signals of different forms, including compound periodical and quasiperiodic signals, nonstationary signals and signals with compound envelopes can be synthesized on the basis of the collection of components mentioned above [5,6]. The mentioned models also make it possible to describe the majority of impulse signals, which are widely applicable in radio engineering (radio pulse and video pulse).
The analysis methods considered above can be applied for that kind of signals.
However, the more general case is obviously more interesting, when signals with compound dependencies of envelopes and signal total phase are applied. Semiinfinite or finite signals with compound envelopes, the most frequently used signals in radio engineering, can be described by the following model
or in general case it would be
For the considered signal models, for instance, of the signal (1),their finite duration can be specified by the finite duration of a complex amplitude
Using the model (1), signals with amplitude, phase and frequency modulation which are commonly used in radio engineering signals can be described. Among with the similar models input signals in sound processing systems, automation devices of power systems at electromechanical transients process can be described.
It is known, that using the signal (1), in case when amplitude, frequency or initial phase are time functions, it is impossible to define the law of amplitude variation and a filter initial phase by values of amplitudefrequency response and phasefrequency response on an input signal frequency [21,22]. Thus, the analysis methods considered above cannot be applied directly.
There are exceptional cases when at some variation laws
However, signals (1) and (2) can be decomposed into components of a set of damped oscillatory components using the Prony's method and its modifications [23]. Leaving behind the issue concerning decomposition into “real” components by using the mentioned methods it is important to note that even general signal approximation (1) and (2) by Prony's method for calculations enables to apply the considered filter analysis methods.
Another option for solution of the filter analysis at signal type (1) is connected to modification of the suggested analysis methods. Modification of the second analysis method is represented below.
4.2. Filter analysis
Let us consider IIR filter analysis at input signals (1), as well as for special cases, when amplitude
To develop the necessary dependences one can use expressions, obtained for IIR filter at injection on its input a set of finite damped oscillatory components. Let us decompose the signal into time steps, during which signal time dependent parameters are mostly constant. In this case time step
In case of an even step of signal decomposition one will obtain the following expression to determine a filter ouput signal
where
As before, only algebraic operations with complex amplitudes, frequencies and values of time dependent transfer function on complex frequency of an input signal are applied for the filter analysis at the input signals type (1).
The same approach can be used for digital filters as well – one should replace continuous time
The filter analysis on the basis of the expression (3) is approximate and can be considered as numerical method. To determine explicit dependencies one need to perform passage to the limit
The required dependency of inputoutput can be defined also by using the convolution integral by substitution of the expression for the input signal (1). Performing discretization with the passage to the limit, the following inputoutput dependency for an analog filterprototype can be obtained
where
In case of amplitude, phase modulation or their combination, the expression for inputoutput would significantly simplify
For IIR filter the following expression for a derivative of time dependent transfer function takes place
In case of
It follows from the comparison of the inputoutput dependencies (4) and (5), that in the second case complex amplitude of IIR filter output signal on the input injection as a damped oscillatory component is determined by multiplication of a signal complex amplitude by the value of time dependent transfer function on the input signal frequency
The dependencies enable to perform filter analysis, including performance analysis of signal processing, at input signals with a composite form. Solving problems of this kind is relevamt not only for radio engineering and communication systems, but also for other industries.
For the example of considered in the sections 3.2 and 3.3 analysis of frequency filters, which are used in PMU, according to a new version of IEEE C37.118.2 standard it is necessary to perform testing for the mentioned devices not only at stepwise change of amplitude and initial phase of sinusoidal component of input signal with commercial frequency, but also at PMU input signals at electromechanical transient process in a power system. In the mentioned regimes of power system operation in the controlled currents and voltage envelopes and total phases are time functions.
For FIR filter analysis at input actions (1) the input  output expression can be obtained on the basis of the dependences given in the Table 6. Suchlike dependences for IIR filters and FIR filters can be obtained based on the dependences given before for the first analysis method.
5. The application of the analysis methods
Let us consider possible areas of application for the suggested performance analysis methods for signal processing by frequency filters, which are used in intelligent electronic devices of power systems, in automation devices, in radio engineering and communication systems, as well as in other fields of engineering where digital signal processing is commonly applied.
The prospectives of power system development in nearest future are related to technology Smart Grid implementation and the application of automatic control and regulation systems of a new generation. Power system control improvement involves a wide application of fast action IED based on synchronized measurement of current and voltage phasors of a fundamental harmonic on the basis of IEEE C37.1182011 and IEC 61850905 standards.
Up to date IED should ensure performance signal processing in conditions of an intense electromagnetic and electromechanical transients process. Mathematical model of an input signal IED in a normal and an accidental regimes of power systems in some cases can be represented by a set of semiinfinite or finite damped oscillatory components, in other cases – by analogous models, in which complex amplitudes and frequencies of mentioned components are time functions.
Most of IED power systems should ensure performance of signal processing at any possible combination of input signal parameters. The suggested analysis methods enable to solve effectively problems of determination of performance specifications for frequency filters, used in IED power systems. The examples for the performance analysis of signal processing by frequency filters and algorithms of signal processing for IED of power systems based on the phasor measurement technology are considered in the sections 3 and 4 of the present chapter. The example for IIR filter analysis for general devices of relay protection and automation is given in the section 2.
The analysis methods for analog IIR filters considered in the chapter can be applied also for the linear circuit analysis. It is important to note that for absolute majority of microprocessor control systems and measuring systems the information sources have analogue nature. In this case for analysis of controlled objects equivalent circuits based on linear electric circuits are used. The illustrative examples: equivalent circuit of power systems, power plants, electriс grids and powersupply systems. Application of the suggested analysis methods enables a consistent approach for the regime analysis of the controlled object operation, analysis of analog filtersprototype and digital filters [5,6].
The suggested analysis methods may be applied for performance analysis of signal processing by frequency filters in uptodate measuring devices, automation devices, radio devices, in communications systems, sound processing systems and other devices, where digital signal processing is commonly used.
Advantages and particularities of the suggested analysis methods are related to uniform analysis methods for analog filterprototypes (IIR filters and FIR filters), as well as for digital filters.
The majority of pulse signals (radio and video pulses), which are commonly used in radio engineering, can be described by a set of finite damped oscillatory components. For performance analysis of pulse signal processing by IIR filters and FIR filters the analysis methods considered in sections 2 and 3 of the chapter can be applied.
The author suggests to use the methods, considered in the section 4 of the present chapter for performance analysis of modulated signal processing by filters.
The synthesis methods mentioned above can be also effectively applied for typical signal filtering problems, including problems of a useful signal extraction against the white noise.
In this case, the white noise realizations can be described as a set of timeshifted fast damping exponents of different digits. The initial values and appearance time of the mentioned exponential components are random variables, which variation law ensures the white noise to have specified spectral characteristics.
6. Performance analysis of signal processing as a step of filter synthesis
Guaranteeing the necessary quality of signal processing is utterly important at frequency filter synthesis. The application of the considered above approaches for the signal type (1) and the filter impulse characteristics (2) enables to reject the traditional approach, related to formulation the requirements to a filter amplitudefrequency response in different fields (band pass, stop band, transition region). In the case of filter synthesis, it is enough to lay down the requirements to the filter frequency characteristics on the basis of Laplace transform on complex frequencies of an input signal with the allowance for their change. Thus, according to the approach described above, the signal processing performance analysis is one of the steps of filter synthesis.
This approach enables to formalize the task of filter synthesis and to gain optimal solutions in combination with methods of multicriterion optimization with limitations[5]. Using methods of multicriterion optimization, by limitations values of filter frequency responses are ment and in some cases values of input signal spectrum based on the Laplace transform spectral representations on complex frequencies of impulse function of a filter and input signal [5,20].
Synthesized in that way filters ensure the specified performance specifications for signal processing at any possible variation of input signal parameters.
Among with frequency filter synthesis this approach gives an opportunity to perform time window synthesis for horttime Fourier transform, as well as synthesis of father and mother wavelets for cases when an input signal can be represented by a set of semiinfinite or finite damped oscillatory components [5,20].
7. Conclusion
The effective methods of performance analysis for signal processing based on the Laplace transform spectral representations and the uniformity of mathematical models of input signals as well as the filter impulse functions as a sets of continuous/discrete semiinfinite or finite damped oscillatory components were developed for a direct determination of IIR filters and FIR filters performance specification. input signals as well as the filter impulse functions as a set of continuous/discrete semiinfinite or finite damped oscillatory components were developed. Simple semiinfinite harmonic, aperiodic signals, compound signals and impulse characteristics of any form can be synthesized on the basis of components set mentioned above, including signal composite envelopes, as well as pulse signals (radio and video pulses). The uniformity of the mathematical description of the signals and filters enables, on one hand, allows to employ both a consistent and a compact form of their characterization in the configuration of a set of complex amplitudes, complex frequencies and time parameters. On the other hand, it simplifies significantly performance analysis for signal processing by analog or digital filters at any possible useful signal and disturbance parameters variation by reducing significantly the amount of calculations.
The analysis methods can be used in case of mathematical models as well  where complex amplitudes and/or complex frequencies are time functions.
To simplify the task of the analysis, two methods are suggested for performance expressanalysis of signal processing by frequency filters using filter frequency responses based on Laplace transform: frequency and frequencytime analysis methods.
The application of the suggested methods for performance analysis of signal processing as one of the steps of filters synthesis enables to automatize filter design for filters with low sensitivity to a signal parameter change within the specified range.
References
 1.
Voronov A. A. 1985 Basic principles of automatic control theory: Special linear and nonlinear systems.  2.
Nise N. S. 2004 Control System Engineering.  3.
Ogata K. 1997 Modern Control Engineering.  4.
Mokeev A. V. 2007 Spectral expansion in coordinates of complex frequency application to analysis and synthesis filters.  5.
Mokeev A. V. 2011 Application of spectral representations in coordinates of complex frequency for the digital filter analysis and synthesis.  6.
Mokeev A. V. 2008 Signal processing in intellectual electronic devices of electric power systems. 3  7.
Mokeev A. V. 2011 Quality analysis of signal processing using digital filters.  8.
Kharkevich A. A. 1960 Spectra and Analysis.  9.
Vanin V. K. Pavlov G. M. 1991 Relay Protection of Computer Components  10.
Prévé C. 2006 Protection of Electrical Networks  11.
Maxfield B. 2006 Engineering with MathCad: Using MathCad to Create and Organize Your Engineering Calculations.  12.
Pritchard P. 2011 Mathcad: A Tool for Engineering Problem Solving.  13.
Mokeev A. V. 2008 Signal processing in intellectual electronic devices of electric power systems. 4  14.
Smith S. W. 2002 Digital Signal Processing: A Practical Guide for Engineers and Scientists  15.
Ifeachor E. C. Jervis B. W. 2002 Digital Signal Processing: A Practical Approach.  16.
Lyons R. G. 2004 Understanding Digital Signal Processing.  17.
Mokeev A. V. 2009 Analysis of digital filters used for preprocessing of signals of relay protection.  18.
Mokeev A. V. 2009 Frequency filters analysis on the basis of features of signal spectral representations in complex frequency coordinates.  19.
Hector J. Altuve Ferrer. H. J. A. Edmund O. Schweitzer E. O. 2010 Modern Solutions for Protection, Control, and Monitoring of Electric Power Systems  20.
Mokeev A. V. 2011 Signal processing algorithms for intelligent electronic devices using phasor measurement technology.  21.
Fink L. M. 1984 Signals, Interferences and Errors.  22.
Gonorovsky I. S. 1981 Radio circuits and signals.  23.
Marple S. L. 1987 Digital Spectral Analysis with Application.