1. A method for diagnostics of at least one internal combustion engine, said engine includes: at least two cylinders; a crankshaft; a crankshaft rotation angle sensor (CRAS) measuring rotation angles of said engine and possessing predetermined technical characteristics, including a pulse response characteristic; a vibro-sensor measuring vibrations of said engine, said vibro-sensor possesses predetermined technical characteristics, including a pulse response characteristic; a logic-mark sensor marking a predetermined phase of operation of a predetermined cylinder, selected from said at least two cylinders; said method includes preprocessing a number of signals received substantially from said engine; said number of signals consist of: a CRAS signal ( 1 ) produced by said CRAS, a vibro-sensor signal ( 2 ) produced by said vibro-sensor, and a logic-mark signal ( 3 ) produced by said logic-mark sensor; said method comprising the steps of: a1) receiving said signals ( 1 ), ( 2 ), and ( 3 ) accordingly from said CRAS, said vibro-sensor, and said logic-mark sensor; b1) at each said signal ( 3 ), dividing the signal ( 2 ) into a plurality of idling cycles, wherein each said idling cycle corresponds to a rotation angle of the crankshaft of 720 degrees; c1) verifying exclusion of idling cycles with varying parameters from said plurality of idling cycles; d1) a step for determining an actual position of said crankshaft based on a stochastic filtration and quasicontinuous representation of the signal ( 1 ); e1) a step for determining an actual position of said crankshaft based on the time-varying Kalman filter or a nonlinear stochastic filter; f1) a step for determining the Riesz generator basis function based on said technical characteristics of the CRAS and approximation of said pulse response characteristic; g1) a step for generating a signal ( 4 ) being a quasicontinuous representation of the signal ( 1 ); a2) a step for determining the Riesz basis generator function, based on said technical characteristics of the vibro-sensor and on said pulse response characteristic of the vibro-sensor; b2) a step for generating a signal ( 5 ) being a quasicontinuous representation of the signal ( 2 ); a3) a step for obtaining a signal ( 6 ) and a signal ( 8 ), both being a secondary discreeting of the signal ( 5 ) with a replacement of argument; b3) a step for obtaining a signal ( 7 ), being a secondary discreeting of a reciprocal function, reciprocal to the signal ( 4 ); a4) a step for obtaining a discrete wavelet transform of the signal ( 6 ), said discrete wavelet transform is represented by a signal ( 9 ); b4) a step for obtaining a continuous wavelet transform of the signal ( 8 ), said continuous wavelet transform is represented by a signal ( 11 ); c4) a step for obtaining a discrete wavelet transform of a function, reciprocal to the signal ( 4 ), and producing output data including a number of two-dimensional arrays of predetermined coefficients, said two-dimensional arrays include pair elements, said two-dimensional arrays are represented by a signal ( 10 ); a5) a step for dividing said two-dimensional arrays relatively to a predetermined index in the discrete transforms of steps (a4) and (c4), and relatively to a predetermined index in the continuous transform of step (b4); and either b5) a step for computations of a one-dimensional empirical probability distribution function for each element of said two-dimensional arrays, said one-dimensional empirical probability distribution functions are represented by signals ( 12 ) and ( 13 ); or c5) a step for a computation of a two-dimensional empirical probability distribution function for each said pair element of said two-dimensional arrays, said two-dimensional empirical probability distribution functions are represented by a signal ( 14 ); wherein the step (c5) is provided as an alternative to the step (b5).