METHOD FOR IMPROVED REPRODUCTION OF DIGITAL SIGNALS

Abstract: The present invention provides a method for digital signal processing, referred to as sample-estimate-hold signal reconstruction, which provides a method of reconstructing output signals more representative of the original input analog signal than the conventional signal processing procedures used in D/A conversion, referred to as sample-hold signal reconstruction. The method uses the history of the analog signal being reconstructed from its digitized form to produce output values which reduce the absolute difference between the integral of the actual signal and the integral of its digital representation over local intervals in the range of the signal. The method is applicable to a wide range of signal processing applications including force-reflection systems, audio and data compression to mention just a few.

Description: METHOD FOR IMPROVED REPRODUCTION OF DIGITAL SIGNALS FIELD OF THE INVENTION The present invention relates to a process for improved reconstruction of analog signals from digitally stored signals. BACKGROUND OF THE INVENTION Digital electronics are rapidly replacing their analog counterparts in many applications of signal processing such as consumer electronics, machine control, robotics, telephony, navigation to mention just a few. Also, computers using digital representations of signals are being used to process signals that may be acquired from an analog source or output to an analog device. Inherent in all these applications is the assumption that the underlying analog signal, such as sound, force, or light intensity, is being represented by a sequence of values corresponding to the amplitude of the signal sampled at fixed intervals. The discreet values are stored in digital format as a sequence of integers. In some applications, such as digital music on a compact disk, an entire signal is stored, but only fetched or examined value by value. In other applications, such as telephony and robotics control, the digital signal is transient and is being processed in real time. In most applications of digital signal processing there are two conversion phases: 1 ) the original analog input signal is converted to a digital sequence by an A/D converter, and eventually, 2) a possibly altered digital sequence is converted to the output analog signal by a D/A converter. The input process is called sampling and is done at a frequency chosen for the application. The reconstructed analog output signal is generated by holding a constant voltage for each interval using the values in the digital sequence in turn. This sample-hold paradigm is inherent in D/A converters used in most digital electronic applications. A basic limitation or drawback of the known methods of reconstruction of an analog signal from its digital representation is that the conversion process forces the output analog function to be a non-continuous piecewise constant function, but does not force the choice of the representing function. An important technological area in which this type of digital processing limitation is problematic is in robotic force control. Specifically, when a robotic device is used to present a human operator with forces that simulate, for example, contact with a wall then it is commonly observed that the virtual wall feels "alive". Rather than the wall acting as a passive barrier, the virtual wall appears to exert a repulsive force so that the return signal coming back from the robot appears to have acquired energy not present in the collision with the wall. This failure to accurately simulate contact also occurs when the robot that presents the human operator with force is connected, via a digital link, with a remote robot that is commanded to interact with real objects and reports the force of interaction over the digital link. In both cases energy appears to have been acquired. The inventors have determined that the error in the energy is proportional to the square of the error in the force. Since the energy is determined by the force-time integral (area under the curve of the force signal) this suggests that finding a piecewise constant function that reduces the area difference with the original continuous signal would have superior behaviour in controlling the robot. It would therefore be advantageous to provide a method of signal processing that gives an analog output signal reconstructed from a digital representation of the original analog input signal such that its integral more closely tracks the integral of the original input signal. SUMMARY OF THE INVENTION The object of the present invention is to provide a method of processing digital signals that gives a reconstructed analog output signal from its digital representation that is more closely representative of the original analog signal than is normally achievable using current signal processing procedures. Another object of the present invention is to provide an improved method for compression of analog signals. Another object of the present invention is to provide a computational method of calculating force reflection in force-reflecting systems. Another object of the present invention is to provide a method of point of view positioning in an interactive visual system. The present invention provides a method for digital signal processing, referred to as sample-estimate-hold signal reconstruction, that is a method of reconstructing output signals more representative of the original input analog signal signal than the conventional signal processing procedures used in D/A conversion, referred to as sample-hold signal reconstruction. The method uses the history of the analog signal being reconstructed from its digitized form to produce output values which reduce the absolute difference between the integral of the actual signal and the integral of its digital representation over local intervals in the range of the signal. In one aspect of the invention there is provided a method of reconstructing an analog signal from a digitized analog signal, comprising; providing sample values y k cfi an analog signal function of time (t)such that y k =F(t 0 +kh)whereh\s the sampling interval, t 0 is the initial time and =0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+ constructing an output analog signal function ^)^., t k ≤ t < t k+i , k = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fF(t)dt- fF(t)dt\ is less than an absolute difference given by | fH(t)dt- fF(t)dt\ locally in intervals defined by a preselected number of previous sample values^; and outputting the reconstructed analog signal function (t) or storing the sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form. In this aspect of the invention each term ^ may be computed as a functional approximation of F(t) in an interval \t k , t Λ+1 ]by calculating an estimation function g k (t) that predicts F(t) in the interval [t k , t A+1 ]by using an effective number of previous values y k ,y k . ,... and integrating the estimation function over the interval [t k , t k+l ] to give j g k (t)dt; i tkti calculating a predicted value to hold using an equation y k =— ( f g k (t)dt); h J calculating a correction term by using an estimation function e k (t) to estimate the integral in the previous interval [t k _ v t and subtracting the predicted area estimation in that interval, the correction being given by Λ =-(f k (t)dt-fg k _ l (t)dt); and adding the correction term to the predicted value to give a value for y k given by y k = -(fg k ( . t)dt + fe k (t)dt-fg k _.(t)dt). In this aspect of the invention the estimation function g k (t) may be equal to p k m (t), which is the polynomial approximation of order m that fits (t) at values y k ,y k . λ ,-,y k . m , and wherein the estimation functioned) equals p k (ή, which is the polynomial approximation of order n that fits (t) at values , _ι, .. ., _ M . In another aspect of the invention there is provided a method for compression of an analog signal. The method comprises a) sampling the analog signal at a first frequency v x or relieving from a storage medium an analog signal previously sampled at frequency V j to produce a digitized signal, filtering the digitized signal to remove frequencies above a preselected frequency, subsampling the filtered signal at a second frequency v 2 wherein v 2 <v, ; b) reconstructing the analog signal from the digitized signal sampled at frequency v 2 , by the steps of providing sample values y k af the analog signal of time (t)such that y k =F(t 0 +kh) where h is the sampling interval, t 0 is the initial time and £=0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+l constructing a sequence { y k } in real time where y k , k=0,l,2, ...μ are computed values of a piecewise constant function of time (t) forming a reconstruction of said signal function such that F(t)=y k , t k ≤ t < t k+l , k = 0,l,2, ...μ , wherein said sequence { y k } is computed by using an effective numerical estimation technique so that an absolute difference given by locally in intervals defined by a preselected number of previous sample values y k is less than an absolute difference given by | fff(t)dt- fF(t)dt\ over the same intervals; and c) outputting the reconstructed analog signal F(t) or storing the sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form. The present invention provides a method of reconstructing an analog signal from a digitized analog signal, comprising; providing sample values y k o an analog signal function of time (t)such that y k =F(t 0 +kh)whereh\s the sampling interval, t 0 is the initial time and £=0,1,2, ... μ ; constructing an output analog signal function^t)^, t k ≤ t < t k+ . , k = 0,1,2, ...μ, by computing a sequence { y k } in real time, wherein each term y k is computed as a functional approximation of F(t) in an interval [t k , t t+1 ] by the steps of calculating an estimation function g k (t) that predicts F(t) in the interval [t k , t^ j ] using an effective number of previous values y k ,y k _ , ... and integrating the estimation function over the interval [t k , t k+l ] to give f g k (t)dt, and calculating a predicted value to hold>' given by y k =— ( [ g ήdt) , h J calculating a correction term by using an estimation function e k (t) to estimate the integral in the previous interval [t k _ v t and subtracting the predicted area estimation in that interval, the correction being given by Λ = f e k (t)dt- f g k _. (t)dt), and Α-l adding the correction term to the predicted value to give a value for y k given by k = - ( / ' g k (i)dt + j " e k (t)dt - f g k (t)dt ) and outputting the reconstructed analog signal function (t) or storing the sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form. In another aspect of the invention there is provided a computational method of calculating force reflection in force-reflecting systems The method comprises the steps of providing sample values y k of an analog signal function of tιme (t)such that y k =F(t 0 +kh) where his the sampling interval, t 0 ιs the initial time and £=0,1, , constructing a sequence { y k } in real time where y k , £=0,1,2, , are computed values of a piecewise constant function of time F(ι) that is a reconstruction of said continuous force function such that y k =F(t 0 +kh) , each term y k said sequence { y k } being computed as a functional approximation of (t) in an interval [t t , t, by I) calculating an estimation function g k (f) that predicts F(t) over the interval [t k , t k+l ] by using sufficient previous values y k ,y k l , and integrating the estimation function over the interval [t k , t k+1 ] to give f g k (t)dt to determine a predicted value to hold II) calculating a correction term by using an estimation function e A (t) to estimate the integral in the previous interval [t k v t and subtracting the predicted area estimation in that interval, the correction being given by iii) adding the correction term to said predicted value to give a value for y k 'hi given by y k = - (fg k (t)dt + fe k (t)dt - fg k _. (t)dt), and 'k 'k-\ 'k-l iv) outputting the reconstructed analog signal function F(t) . The present invention also provides a method of determining a time-dependent three-dimensional affine transformation of a point x, of the form T(x,t) =R(t) x+d(t) where R is an orthogonal matrix having elements which vary according to sample values and d is a translation vector having elements which vary according to sample values, comprising: (a) computing a unit rotation axis vector 7(f)-- r y ( ) and rotation rXi) angle θ(t) from the matrix R(t) ; (b) producing a reconstructed rotation vector axis r(t)-- r y {t) in which each element of the vector r(t), j=x,y,z is reconstructed independently by: computing sample values y k of the function of time such that 3^=7- (t 0 +£/?) where Λ is the sampling interval, t 0 is the initial time and £=0,l,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t<t ktl ; and constructing a function r(t)=y k , t k < t < t k+1 , k = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by is less than an absolute difference given by | fH(t)dt- fr(t)dt\ locally in intervals defined by a preselected number of previous sample values^; (c) producing a reconstructed rotation angle θ(t) by: computing sample values y k of the function of time such that y k =Q(t 0 +kh) where h is the sampling interval, t 0 is the initial time and £=0,l,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+l ; and constructing a function Q(t)=y k , t k ≤ t < t k+x , k = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fθ(t)dt- fθ(t)dt\ is less than an absolute difference given by | ΪH(t)dt- fd(t)dt\ locally in intervals defined by a preselected number of previous sample values^; (d) normalizing said reconstructed rotation axis, and computing a reconstructed matrix R(t)from the normalized reconstructed rotation axis; (e) producing a reconstructed displacement vector in which each element of the vector d(t), j=χ,y,z is reconstructed independently by: computing sample values y k of the function of time such that t =< (t 0 +£?) where /7 is the sampling interval, t 0 is the initial time and £=0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k < t<t k ^ ; and constructing a function d(t)=y k , t k < t< t fe+1 , k = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fd(t)dt- fd(t)dt\ is less than an absolute difference given by | fH(t)dt- ϊdft)dt\ locally in intervals defined by a preselected number of previous sample values^.; and (f) transforming the pointx by f( x, t) = R(t) x +d(t) using the reconstructed R(t) and d(t). BRIEF DESCRIPTION OF THE DRAWINGS The improved process for digital signal processing will now be described, by way of example only, reference being had to the accompanying drawings, in which Figure 1 shows a continuous function and its digital reproduction by a piecewise constant function according to the method of the prior art, Figures 2 and 3 graphically illustrates the sample-estimate-hold method of reconstruction of a digital signal according to the present invention, Figure 4 shows the same continuous function as in Figure 1 but with the reconstructed piecewise constant function being produced using the present sample-estimate-hold method, Figure 5 illustrates a piecewise-hnear polynomial that is approximated by an integral preserving stepwise function, Figure 6 illustrates a piecewise-linear polynomial that is approximated by a stepwise function at a higher sampling rate Figure 7 illustrates the time-displacement behaviour of the virtual- wall impact using sample-and-hold on equal intervals, Figure 8 illustrates the time-displacement behaviour of the virtual- wall impact using sample-estimate-hold on equal intervals, Figure 9 shows the time-displacement behaviour of the virtual-wall using sample-estimate-hold, producing reflected forces 10 times for each force that was determined, Figure 10 illustrates the time-displacement behaviour of the virtual- wall impact using sample-estimate-hold, determining forces 10 times for each force reflection, Figure 11 is a plot of power spectral density of an original music signal, Figure 12 is a plot of the difference in power spectral density of the zero-order-hold representation of the filtered and subsampled music signal, using the prior art method of digital reproduction, and the original music signal of Figure 11 , Figure 13 is a plot of the difference in power spectral density of the music signal reconstructed using the present sample-estimate-hold method after filtering and subsampling and the original music signal of Figure 11 , Figure 14 is a plot of power spectral density of an original voice signal, Figure 15 is a plot of the difference in power spectral density of the zero-order-hold representation of the filtered and subsampled voice signal, using the prior art method of digital reproduction, and the original voice signal of Figure 14, and Figure 16 is a plot of the difference in power spectral density of the voice signal reconstructed using the present sample-estimate-hold method after filtering and subsampling and the original voice signal of Figure 14 DETAILED DESCRIPTION OF THE INVENTION As stated above, from a mathematical perspective, the D/A conversion process involves an underlying continuous function being represented by a piecewise constant function. Figure 1 shows a continuous function and its reproduction with the piecewise constant function, being the prior art sample-hold method. The problem can be stated mathematically as follows: Let y k , £=0,l,2,...μ , be the sample values of some continuous signal function of time F(t) such that y k =F(t 0 +kh) where h is the sampling interval, and t 0 is the initial time. Let H(t) be the piecewise constant function of time such that H(t)=y k , t k < t <t A . +1 , £ = 0,l,2,...μ . H(t) is the output analog function using the sample-hold method. Let F(i) be the piecewise constant function of time such that F(t)=y k , t k ≤ t < t k+1 , k = 0,l,2,...μ . F(t) is the output analog function to be constructed using the sample-estimate-hold method. The problem is to choose the sequence { y k } so that an absolute difference | fF(t)dt- fF(t)dt\ locally over intervals involving a few sample points is smaller than an absolute difference | fH(t)dt- fF(t)dt\ over the same intervals. In order to construct the sequence { y k } each term y k is computed in order that (t) is a functional approximation of F(t) in the interval [t k , using a small finite number of previous values y k ,y k _ λ , ... • The previous values y k ,y k . λ , ... are used to provide an estimation function that predicts the trend of F(t) in the interval [t k , t k+l ] . Let p k '(t) be a polynomial approximation of F(t) using y k and m previous values Then the predicted estimate to hold, y k is found by equating the integral of the estimation function in the interval [t k , t fc+1 ] , '* 1 f p k m (t)dt to the area under the constant value hy k This choice of the predicted estimate makes the integrals of the polynomial estimation function and the constant function equal in the interval [t k , t k+1 ] Provided the polynomial function is a close approximation of F(t) in the interval the error in the integral will be small However, if the estimation is not close over the interval the integral of the polynomial over the interval will not be as close to the integral of F(t) as desirable By using the actual value for y k once it is known, a correction for the overestimation or underestimation in the integral in the current interval can be used to adjust the value of the integral in the next interval This is graphically illustrated by Figures 2 and 3 The correction in the current interval, y k is made by estimating the overestimation or underestimation of the integral in the previous interval using the polynomial /?/(/) that approximates F(t) using y k and n previous values The correction is determined by equating fp k "(t)dt - hy k p _ λ to hy k c Thus, the final estimate for h 1 y k is given by Λ The sample-estimate-hold method is determined by the choice of m and n When m=0 and n=0 , the constructed sequence is y 0 =y 0 , y k = 2y k -y k £ = 1, 2, μ There is no requirement that m and n be equal Consider the case when m=0 and «=1 Then the constructed sequence is $ 0 = O . S> k = (hy k + - (y k _. + y k ) - hy k _ χ )= y k - i^ , A = 1, 2, A simple practical case to treat is when both polynomial approximations are linear, i e m=\ and n=\ Then f p k (t)dt = — h (y k +y k+l ), where J 2 Λ , = 2h(y k -y k .)/h +y k _ { Thus, =1 (y A + 2(y,- ,) + >* ,) = - -^-i The corrector term is Λ C = 0Α ~2 A- I + Λ- 2 ) and the final estimate is y k =y k +y k = 2y k - —y k , + — ' A _ 2 To begin the sequence, set y 0 =y 0 and y. = 2y. -y 0 The sequence { y k } is the proposed sequence of constant values in the piecewise constant functιon (t) which has the property that its integral more accurately reflects the integral in local intervals of the underlying analog signal sampled as the sequence { y k } Figure 4 shows the same continuous function as in Figure 1 but with the reconstructed piecewise constant function being produced using the present sample-estimate-hold method It can be seen that the reconstructed staircase function overlaps the continuous function being above it and below it In the foregoing analysis the numerical estimation of the sequence {y k } is based on using a predictor-corrector technique using polynomial approximations However, the sample-estimate-hold methodology can use any numerical estimation that is designed to reduce the absolute difference between the integral of the original analog signal and the integral of its digital reconstruction in local intervals In the foregoing analysis a polynomial of order m is used in the predictor step to estimate the trend in the signal and a polynomial of order n is used in the corrector step to estimate the error in integration in the previous interval However, it is noted that in a predictor-corrector scheme as disclosed herein there can be a lag in locking on to the trend As the order of the approximation increases more previous samples are needed to produce the reconstructed signal so that there is potentially a trade-off between tracking lag and accuracy Therefore, although higher order polynomials may be used, the inventors have found that for many of the different applications studied a first order estimation for both the predictor and corrector estimation polynomials works well The method uses the history of the analog signal being reconstructed from its digitized form to produce output values which reduce the absolute difference between the integral of the actual signal and the integral of its digital reconstruction in local intervals over the range of the signal The method disclosed herein utilizes the history of the signal being processed from its sampled data to give a better fit to the integral of the underlying continuous function represented by the digital sequence By using knowledge at several previous sample points in the reconstruction, first and second derivative information in the signal is being captured and utilized to provide a more accurate reconstruction From basic calculus, if the sampling interval is chosen sufficiently small, the signal reconstruction using the prior art sample-hold method can be made accurate to within any desired threshold However, for practical applications other considerations such as storage space, bandwidth, or sampling limitations require that the sampling interval be chosen just small enough to provide the desired performance Thus, the sample-estimate-hold method disclosed herein is contemplated by the inventors to be useful in a broad range of applications where it can produce higher fidelity or reduce bandwidth or storage space The sample-estimate-hold method disclosed herein for digital signal processing is very advantageous because it is simple, readily implemented and there is no need of parameter tuning only the order of the prediction and correction functions needs to be selected Non-limiting examples of the efficacy of the digital signal processing method disclosed herein as applied in several different technological applications will now be described EXAMPLE 1 Robotic Force Reflection Control As briefly mentioned in the introduction, an important problem in the field of force-reflecting systems and telerobotics is poor rendering of contact particularly of contact with stiff surfaces There are numerous possible sources of poor performance, including poor contact models, sampling errors, and delays due to computation or data transmission A haptic, or force-reflecting, interface is a device that can sense a human operator's hand motion, and actuate the operator with forces or torques, without significantly impeding the hand motion The term "haptic" refers to both tactile and kinesthetic senses These devices have been used for teleoperation of a remote robot, the characterization of human motor performance, conducting psychophysical experiments on haptics, and the interaction of a human with a virtual environment There are many designs of devices that can provide stimuli for the haptic sense A frequent observation of users is that the interactions feel "alive" or "active" which arises from the system providing more energy to the operator than is provided to the system, rather than feeling like a passive mechanical system The problem of liveliness often becomes very pronounced as the stiffness of the interaction increases, and can ultimately lead to an uncontrollable instability Previous work on this liveliness problem has concentrated primarily on (a) providing more realistic models of contact for a simulated environment, (b) stabilizing effects due to delays in a teleoperation system, see R J Anderson and M W Spong, "Bilateral Control Of Teleoperators With a Time Delay", IEEE Transactions On Automatic Control", Vol 34(5), pp 494-501 , 1989, or (c) analysing effects of sampling errors Colgate et al designed a virtual wall composed of springs and dampers, see J E Colgate and G Schenkel, "Passivity Of A Class Of Sample-Data Systems Application To Haptic Interfaces", Proceedings of the American Control Conference, Baltimore MD 1994, pp 3236-3240. In all cases, the main problem faced was how to obtain passivity in virtual walls. Some researchers have suggested that some criteria have to be defined on the basis of stiffness and damping coefficients in order to make the virtual wall passive, see M. Minsky, M. Ouh-Young and F.P.J.B.O. Steele, "Feeling And Sensing: Issues In Force Display", Computer Graphics, Vol. 24(2), pp. 235-243, 1990.Colgate and Schenkel provide a bound on the sampling frequency which, for a given force-reflecting system and a given virtual wall, will ensure passivity. Their model of system behaviour included zero-order hold in its usual form. The inventors have found that the sample-and-hold method of reconstructing force reflection plays a significant part in the unrealistic presentation of contact information. Force-reflecting systems are computer controlled, so the forces of interaction are acquired at a finite rate. When the forces are reflected to the human operator, the combined operation of a digital-to-analog conversion and an amplifier amounts to a sample-and-hold of the reflected force. Using the technology now available, the effect is a fundamental property of both the discrete domain and of the physical interfaces of a computer to the world. Applying the present sample-estimate-hold method of digital signal processing clarifies some of the issues that arise when implementing a force-reflecting system using a digital computer. The method is numerically based and is completely independent of the stiffness and damping coefficient of the virtual wall thereby allowing the creation of a transparent force rendering of a virtual wall, even of a wall that has zero damping. It is shown herebelow that the way a force signal is sampled is intimately related to the conservation of energy of a dynamical system. The sample-estimate-hold method is applied to the force reflecting control problem for calculating the reflected force while preserving the fundamental characteristic of the force. Adaptations of the methods are provided for the case of determining and reflecting forces at equal rates, and where the determination and reflection take place in integer ratios. The importance of an accurate estimation of impulse in computing the energy transfer in a contact task is evident from an analysis of the dynamics equations governing the system. The dynamics equations of motion are derived from the Lagrangian of a one-degree-of-freedom mechanism but it will be understood that the derivation can be generalized. To derive the dynamics equations, a model is used comprising a one-degree-of-freedom force-reflecting system as a single link of mass m that has spatial coordinate q(t) and velocity coordinate q'(t) . The Lagrangian equation of motion in the presence of a forcing function f(t) is mq "(t) =βt) - c(q '(t)) -p(q(f)) for a damping force c(q '(t)) and a potential force p(q(t)). The kinetic energy of the link at time can be calculated as £(/,) = -m ( ff(t)dt) 2 in the case when there is no damping or potential energy present. This equation implies that, if a one-degree-of-freedom mechanism moves due just to an external forcing function, then errors in estimating the impulse over the time interval magnify the apparent energy quadratically In the more general case of damped motion in a potential field this quadratic term will be altered by damping 'l In a discretely sampled system the integral ff(t)dt is not directly observable but must be reconstructed from observations of f(t) To maintain the principle of conservation of energy, the integral must closely approximate the true integral of the force during the sampling period A particularly poor choice of 'l the integral (βj)dt, and the most commonly used method, is to sample f(t 0 ) and r o hold it during the interval [t 0 ,t,] to produce the estimate ff(t 0 )dt^f(t 0 )(t 1 -t 0 ) A Ό major problem with this approach is any trends in /(t)over the interval are ignored This equation is the usual prior art sample-and-hold method that is used in many force-reflecting systems A major problem is that the error introduced by using this equation to produce the estimate is magnified quadratically in the energy, as indicated in the previously developed equation The present sample-estimate-hold method is adapted to reduce the error introduced into the integral of the signal by the prior-art sample-and-hold method As discussed more generally above, the method is disclosed using first-order polynomials in both the prediction and correction segments to approximate the force-time integral However, the order of the polynomial used to predict and correct the integral may be higher, the prediction and the correction part need not use the same order of polynomial, and non-polynomial models may be used a) Equal-Time Sampling Suppose that force is determined by acquisition from a sensor or generation from a numerical simulation over time, and is then reflected to a human operator. Assume that the sampling intervals of the determined force and the reflected force are fixed and identical, and let the sampling interval be h. The function F(t) s the smooth underlying function of force with time that is to be approximated from the force data obtained at sampling intervals. In the prior art sample-and-hold method, the integral f F(t)dt is approximated by hy k . The present method is to use the previous force readings and the current to predict a new reading for the next instant. If the readings are used to extrapolate to a new reading, then the force-time integral for the upcoming interval can be estimated. Following the general development described earlier, using linear functions for prediction and correction, let y k be the new zero-order-hold force that is to be output at time t = t A . and let the estimated value for time t = t kJr . be y k l : y k+l =2y k -y k _ l . The force-time integral is then approximated as f F(t)dι - —(y k +y k λ ) from which the new force value to be output is found from •i 2 3 ] yζ = — y k - —y k _ v However this equation does not account for the error incurred in the previous time instant by the prediction scheme. If this error is not taken into account the integral of F(t) will be incorrect. Similar to the general procedure described above, a correction term is calculated by estimating the integral from the previous instant to the current instant and subtracting from this area the area of the previous prediction. The graphical illustration of this technique is very similar to that shown in Figures 2 and 3 but with the ordinate axis being Force The correction term is given by y k = —(y k - 2y k _ χ +y k _ 2 ) Addition of the estimation term and the correction term gives =2 Λ -^Λ-ι + - 2 This last equation produces a new value of the force being sampled based on its previous zero-order representations This new value of force, when held over the corresponding sampling interval, approximates the force-time integral locally over intervals of a few sample points more closely than the usual sample-and-hold value would As a result, this new scheme introduces less error in the force-time integral The error in the energy is quadratically less b) Rapid Output Sampling When a force reflecting system is used by a human operator to interact with a virtual environment, it is possible that the controller for the device can produce output forces at a higher rate than the dynamic simulator can calculate the forces In such a situation, the force reflection can be further improved by careful use of the zero-order hold of the output force Essentially, at a given instant the piecewise-linear prediction of the force to be output can be approximated by a series of steps shown in Figure 6 Let h } be the sampling interval for the input of forces (e g from a dynamics simulator or remote robotic sensor), and let 7 o be the sampling interval for output of forces to the motor amplifiers of the force-reflection system Suppose h o and/?, are related by h o = rh t where r is a positive integer As in the case of equal-time sampling where y k is chosen so that it approximates h. i t k . ι - fF(t)dt, with rapid output sampling, the integral f F(t)dt can be approximated by a plurality of stepwise approximations using a boundary condition that the force at time t k+1 being given by the estimated value y k+l . Then set an equation of force for each step S equal to ^ = d - p( r ~ 1 -j),forj = 0,\,...,r-\ , where p is an increment between steps. To maintain the same integral, set y k = __ — S , and substitute the right hand side y=o r ' of the above equation forS into this equation and solve forp to give p =2.0 ^ " ) . (/- I) The value of p can then be found from this formula using the 3 1 equations y k - λ = 2y k -y k _ . and y k = 2y k - —y k _ 1 +—y k - 2 f rom tπe predictor-corrector technique. Then p can be substituted into the right side of the equation for S to derive an expression for S in terms of current and previous determined forces. These output force valuesS are provided to the force-reflecting system at a higher frequency than the acquisition in the relationship in h 0 =rh l . By doing so, the determined forces are extrapolated as a first-order polynomial which is then approximated by m stepwise-constant outputs. c) Rapid Input Sampling When the output rate is limited in order to avoid mechanical resonances or to avoid exceeding the mechanical bandwidth of the force- reflecting system the forces of interaction can be determined at a higher rate than they can be reflected. In this case, the inputs of a time interval can be used to form a better approximation of the smooth function that is assumed to be producing the determined forces. Let h l be the sampling interval for force determination and h o be a sampling interval for force reflection, selecting the interval for force reflection to beh o =rh ι where r is an integer. Let t = t k denote the k th time that a reflected force value is calculated. Between times t = t i _, and t = t k there will be r+ 1 force values determined. These can be used in various ways but to be consistent with the linear-function approach these r + 1 force values are used to construct a least-squares linear polynomial at time £ given by P k (t) = σ k t+a k where σ k and r are found by conventional regression equations. This polynomial A (t) is used to '*;' h estimate the integral / F(f)dt by — (P k (t k ) + P k (t k , λ )) to predict a value to be '* h held. It is also used to estimate the integral f F(t)dt by —(P k (t k _ x ) +P k (t k )), for use in the correction. The , used to predict in the previous interval is also used in determining the correction term. The prediction and correction terms are added together to obtain the output value y k given by ( * ( '* )+ * (/ *- ι))+ ^ • This equatioπ can be simplified, or the quantities in parentheses can be calculated and used directly. Validation of Methods for Robotic Force Reflection Control Validation of the sample-estimate-hold method for force reflection was obtained using a light, low friction five-bar linkage two-degree-of-freedom planar haptic interface. The computer was a PC-compatible microcomputer (Intel 486 processor, 66 MHZ clock). The analog interfaces were provided by commercial boards that provided 12-bit conversion. The Cartesian position of the device handle was found by reading signals from the potentiometers attached to the motor shafts, converting from voltages to radians, and then using forward kinematics to find the position of the device tip. Reflection was performed by converting the desired Cartesian force to motor torques using the inverse Jacobian transformation. The desired torques were converted to desired voltage outputs by previously calibrating the combined sensitivity of the linear amplifiers and the actual motors used. To avoid mechanical resonances and to provide results that may be comparable to results from other devices all tests were conducted at a nominal sampling rate of 400 Hz. Clock jitter was less than 1 %. To ensure that the master mechanical system was passive the human operator was excluded from the tests and instead, the force-reflecting mechanism was tilted by approximately 15 degrees so that a diluted gravitational force acted on the mechanism. Four kinds of tests were conducted on a virtual wall of stiffness 0.35 N-mm. The wall stiffness value was chosen to be the stiffest attainable without reaching saturation of the analog-output stages. The reflection value was calculated as the penetration times the wall stiffness. The first test was a simple implementation of a virtual wall, in which penetration of a virtual wall was reflected immediately. The second test was implementation of the present sample-estimate-hold method designed for equal sampling intervals. The third test was implementation of the sample-estimate-hold method for rapid reflection, in which forces were determined at 400 Hz and reflected at 4000 Hz The final test was implementation of the sample-estimate hold method for rapid determination, in which forces were determined at 4000 Hz and reflected at 400 Hz For each test the displacements and reflection forces were logged and subsequently analysed The area of the hysteresis curve, which is the difference in area from initial impact to maximum penetration and from maximum penetration to final exit, was calculated The hysteresis results are presented in Table I All results are for only the first, largest impact of the handle position with the virtual wall Table I Hysteresis Values From Experiments Case Energy Loss ( mJ) Uncompensated -1 214 Equal Sampling -0 027 Multiple Reflections +0 903 Multiple Acquisitions +0 951 As can be seen from the hysteresis results in Table I, the prior art sample-hold reflection method adds significant energy to the simulated impact of the mechanism with a virtual wall The sample-estimate-hold method for equal sampling adds very little energy, and the other methods with unequal sampling actually absorb energy to effectively dampen the impact In the presence of a potential field (e g gravity) it is expected that the energy provided by an active wall would eventually balance the potential energy of peak rebound, that is, a limit cycle is expected This is confirmed by the displacement-time plots that display exactly this expected phenomenon Figure 7 shows a perpetual-motion bounce of 5 mm magnitude when sample-hold is used while Figure 8 shows that, with sample-estimate-hold, the bounce magnitude is reduced to less than a millimetre Figure 9 gives an experimental demonstration of the passivity of the sample-estimate-hold method for rapid reflection, and Figure 10 demonstrates the passivity for rapid determination The results clearly show that the sample- estimate-hold method significantly reduces the excess energy injected into the system by the prior art sample-and-hold method Stabilization is distinctly improved by increasing the sampling rate in either the force determination stage or the force reflection stage because the sample-estimate-hold methods take considerable advantage of the additional information available In other words the present methods make use of the additional information available and provide better performance than is available by slowing down the overall system to achieve synchronization Tracking in these cases is not perfect but it will be appreciated that it is far better than in using sample-and-hold to reflect contact forces The present method is advantageous over known filtering methods since the latter typically employ several tuning parameters that must be adjusted for each implementation whereas, with the present sample-estimate-hold method, once the order of the method is chosen then all other parameters are completely determined thereby eliminating the need for parameter tuning. A particularly advantageous aspect of this method is the ability to improve performance in a system where the acquisition and reflection rates differ. The sample-estimate-hold methodology has a wide range of potential applications in robotics and may be used wherever zero-order sampling is required. For example, sampling and estimating position or velocity input to a controller, or implementing force control in a traditional manipulator all require the sampling of signals. The inventors also contemplate that the method disclosed herein may be used to examine to provide solutions to problems that arise from time delay. EXAMPLE 2 Reproduction of Digital Acoustic Signals Acoustic signals may be viewed as analog signals in which air pressure is applied against a transducer such as the ear or microphone. A loudspeaker applies a force to a diaphragm which then variably pressurizes the volume of air in the listening area in reconstructing the same forces originally experienced by the microphone. Therefore, energy conservation is an important consideration for improved audio signal reconstruction for the same reasons that use of energy conservation is important in the method for reconstruction of the forces in the robotic force reflection application in Example 1. Digital Music Signal A digital music sample was taken from a compact disk and converted to a monaural signal, 16 bits per sample, with an original sampling rate of 44.1 kHz on a SUN workstation. An algorithm embodying the sample- estimate-hold method was applied to the digital music sample to produce a reconstructed analog music signal. The reconstructed music sample was crisper at higher frequencies and appeared, based on qualitative analysis, to provide a more realistic sound for various musical instruments. The following subsampling test was conducted to show that the music signal reconstructed using the present sample-estimate-hold method better reflects the original audio signal before it was digitized. In the test the assumption is made that the original digital sample using the sample-hold process is the target signal to reproduce. Then a digital signal at of the frequency is selected to represent the digitization of the target signal. It is well known by those skilled in the art that just sampling at the rate will introduce sampling bias. A second order Butterworth filter, with a cut-off frequency of 7 kHz and Hanning windowing, was used to remove the high frequencies, and the subsample at 22.05 kHz was selected. This produced a loss of intermediate frequencies due to the filtering in addition to the truncation of the high frequencies by the sampling. The power spectrum of the original music signal is shown in Figure 11. Figures 12 and 13 give the difference between the power spectra of the filtered then subsampled and the filtered then subsampled then reconstructed music signals, respectively, and the original music signal. Figure 12 shows that the result of the prior art sample-hold reconstruction of the subsampled signal loses about 5 dB in the middle frequencies, but is noticeably poorer when compared to the original signal in the higher frequencies. Figure 13 illustrates the performance of the present sample-estimate-hold process applied to the filtered and subsampled signal wherein the higher frequencies have been lifted to more closely mimic the target signal. Two important conclusions can be drawn from this subsampling experiment. First, using a combination of filtering, subsampling and the sample- estimate-hold reconstruction method disclosed herein, a 2-to-1 compression of the audio signal is achieved with very little degradation in sound quality. This confirms that a 50% improvement in bandwidth is achievable for audio music signals, or conversely substantially twice the sound fidelity is achievable for the same bandwidth. The second conclusion is that using the sample-estimate-hold reconstruction of the original digital music signal is most likely producing a reconstruction that is closer to the original analog music signal. Digital Voice A digital voice sample was recorded on an analog cassette and then sampled using a 16-bit sound board on an IBM compatible PC to obtain a monaural signal with 16 bits per sample at a sampling rate of 16 kHz. Two subsampling tests were conducted. In the first the voice signal is compressed by a process of filtering and subsampling to obtain a 2-to-1 compression and then reconstructed to produce a signal that is much crisper than the filtered and subsampled voice signal and quite close to the uncompressed signal in quality The digital voice signal was filtered with a 3rd-order Butterworth filter with a cutoff frequency of 2 5kHz was used and the signal was subsampled to 8 kHz This produced a loss of intermediate frequencies due to the filtering in addition to the truncation of the high frequencies by the sampling The power spectrum of the original voice signal is shown in Figure 14 Figures 15 and 16 give the difference between the power spectra of the filtered then subsampled, the filtered then subsampled then reconstructed voice signals, respectively, and the original voice signal Figure 15 shows that the result of the prior art sample-hold reconstruction of the subsampled signal loses about 5 dB in the middle frequencies, but is noticeably poorer when compared to the original signal in the higher frequencies Figure 16 illustrates the performance of the present sample- estimate-hold process applied to the filtered and subsampled signal wherein the higher frequencies have been lifted to closely mimic the target signal The second subsampling test achieves a further compression by using a logarithmic amplitude encoding of the filtered and subsampled signal In this test the 16 bit 8kHz signal was compressed to 8 bits using the mu-law representation This achieved a further 2-to-1 compression of the voice signal Then the 8-bιt mu-law representation was decompressed and the sample- estimate-hold reconstruction was used to produce a 16 bit, 8 kHz signal that has a frequency spectrum that is almost identical to the sample-estimate-hold reconstruction from the first test An important conclusion can be drawn from this subsampling experiment. In using a combination of filtering, subsampling and the sample- estimate-hold reconstruction method disclosed herein, a 2-to-1 compression of the voice audio signal is achieved with only a little degradation in sound quality. In applications where transmission bandwidth or storage space is critical, a further 2-to-1 compression can be achieved by using a prior-art logarithmic amplitude compression method in conjunction with the sample-estimate-hold method disclose herein. EXAMPLE 3 Point of View Positioning In Virtual Reality An interactive visual display of virtual environments typically uses a 3D head-tracking system to determine how the virtual scene is to be rendered. A common problem in such systems is that the head tracking exhibits lag. This lag may be significantly reduced using the present sample-estimate-hold process to provide a better reconstruction of the viewer's motion. Rendering a fixed scene from a mobile viewpoint using computer graphics involves a rigid-body transformation of every scene point that is to be drawn. To transform a point x in scene coordinates to viewer coordinates, the transformation is typically expressed in the form T(x ] ,t) =R(t)x +d(t) where R is an orthogonal rotation matrix and d is a translational displacement vector. Applying the present sample-estimate-hold process directly to the orthogonal rotation matrix R is not advisable since the resultant matrix will not in general be orthogonal. It is preferable to compute the unit rotation axis r(ι) and rotation angle θ(t) from the matrixR(t). Given R, define a skew-symmetric matrix B = (R -R τ )/2 whereR τ is a transposed matrix of the matrix R,B is of the form 0 -b z b y B b z 0 -b -b y b 0 From this define a scaled and unit rotation axis and a rotation angle respectively _ (trace(R) - \) so that the transformation T is a rotation about r by angle θ followed by a displacement d. A reconstructed rotation axis (t) is produced using the present sample-estimate-hold method in which each element of the vector r(i) is computed independently. Next, a reconstructed rotation angle θ(t) is produced using the sample-estimate-hold method. The reconstructed rotation axis is normalized then a reconstructed matrix R(t) is computed from the normalized reconstructed rotation axis by computing the reverse transformation - cos(θ))) +cos(θ) (r r (\ - cos(θ))) -r.sin(θ) (r r (1 - cos(θ))) +r sin(θ) R - (r χ r y (\ - cos(θ))) +r z sin(θ) (/^ 2 (1 - cos(θ))) + cos(θ) (r y r : (\ - cos(θ))) -r x sin(θ) (r χ r z ( 1 - cos(θ))) - r y sin(θ) (r z ( 1 - cos(θ))) + r x sin(θ) (r z ( 1 - cos(θ))) + cos(θ) where The efficacy of the sample-estimate-hold process was qualitatively established by rendering a complex scene at 25 Hz using a head-mounted display. The lag normally accompanying the traditional method was significantly reduced and in particular the vertigo-inducing effects were dramatically reduced. An apparatus for point of view positioning in an interactive visual system would include a computer representation of a real or artificial scene, which is computed or relieved from a storage medium, a visual display means, sensors for sensing or determining the orientation and/or position from which the computer-generated visual scene is to be computed, processing means connected to the display means and the sensors. The processor determines a transformation of points in the scene, where each point x t is transformed using a transformation T x p t) = R(t) χ l +d(t) where R is an orthogonal matrix the elements of which vary according to the orientation and position sample values and where d is the translation vector the elements of which vary according tothe orientation and/or position sample values. The processor means includes means for computing the reconstructed rotation matrix R(t)and a reconstructed translation vectord(t) as discussed above so that R(t) and d(t) are digital representations of reconstructed analog signals. The processor renders the transformed points f(x t ,t) to the display means by using the reconstructed analog signals in digital or analog form. The orientation and/or position of an anatomical part of a person or animal may be sensed by a mechanical linkage, acoustic range sensing, sensing of magnetic fields, sensing of electric fields, sensing of reflected light, sensing of emitted light or combinations thereof and in which the orientation and/or position from which the visual scene is to be rendered is computed from the orientation and/or position of the anatomical part The anatomical part may be the head, eye, a hand, a finger, a foot, a toe or any combination thereof of a user Further, the hand and/or foot motions may be sensed by means of steering wheels, joysticks or foot pedals Similarly, an apparatus may be constructed using the signal processing method of the present invention for orienting and/or positioning a remote sensor for display This apparatus includes a remote sensor for remote sensing, an actuator means for mechanically, electrically, magnetically, hydrau cally, or pneumatically changing the orientation and/or position of the remote sensing means, and a sensing means for sensing or determining the desired orientation and/or position of the remote sensing means as sample values A display means connected to the remote sensing means A processing means is connected to the sensing means and the actuator means with the processor operable to determine a transformation of orientation of the remote sensing means at a given time t as a rotation R(t) and/or a position of the remote sensor as a translation vector d(t) , where R is an orthogonal matrix the elements of which vary according to the orientation and position sample values and where d is a translation vector the elements of which vary according to said orientation and position sample values. The processor computes the reconstructed rotation matrix R(t) and a reconstructed translation vector d(t) according to the method of signal reconstruction disclosed herein so that R(t) and d(t) are digital representations of reconstructed analog signals. An actuator is provided for changing the orientation and/or position of the remote sensing means to an orientation specified by the R(t) and/or the position specified by the d(t) by using the reconstructed analog signals in digital or analog form. The remote sensor may a camera, radio detecting and ranging device, laser detecting and ranging device, ultrasound detecting and ranging device, force sensor or tactile sensor to mention just a few. While the sample-estimate-hold method forming the present invention has been described and illustrated with respect to several digital signal processing applications, those skilled in the art will appreciate that the present method can be applied to may other signal processing applications, and as such, they are covered by the claims as defined herebelow.

Claims: THEREFORE WHAT IS CLAIMED IS: 1. A method of reconstructing an analog signal from a digitized analog signal, comprising; providing sample values y k an analog signal function of time (t)such that the sampling interval, t 0 is the initial time and £=0,1,2, ... μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+1 ; constructing an output analog signal function P(t)=y k , t k ≤ t < t k ^ , £ = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by I fP(t)dt- fF(t)dt\ is less than an absolute difference given by I fH(t)dt- fF(t)dt\ locally in intervals defined by a preselected number of previous sample values^; and outputting said reconstructed analog signal function (t) or storing said sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form. 2. The method according to claim 1 wherein each term y k is computed as a functional approximation of F(t) in an interval [t k , t k+1 ] by calculating an estimation function g k (t) that predicts F(t) in the interval [t k , t k+l ] by using an effective number of previous values y k ,y k . x , - and '*.ι integrating said estimation function over the interval [t k , t k+l ] to give f g k (t)dt ι ' w { calculating a predicted value to hold using an equation y k =— ( f g k (t)dt), h J calculating a correction term by using an estimation function e^(t) to estimate the integral in the previous interval [t k v t and subtracting the predicted area estimation in that interval, said correction being given by Λ = f e k (t)dt-jg k l (t)dt), and 'k 1 adding said correction term to said predicted value to give a value for y k given by 3 The method according to claim 2 wherein the estimation function g k (t) equals p k "(t) , which is the polynomial approximation of order m that fιts (t) at values y k ,y k ,y k m , and wherein the estimation function e k (t) equals p k (t) , which is the polynomial approximation of order n that fιts£(t) at values , ,y k n 4 The method according to claim 3 wherein m = 0 and n= 1 and then 5 The method according to claim 3 wherein m = 1 and n = 1 and then Λ = 2 Λ --Λ I + -Λ 2 6. The method according to claim 5 wherein the step of providing sample values said signal function comprises providing said sample values at a rate faster than the rate at which said reconstructed signal function is outputted, letting h ι be the sampling interval for signal sampling and h o be a sampling interval for signal reconstruction, selecting the interval for signal reconstruction to beh o =rh ι where r is an integer, letting t=t k denote the k th time that a reconstructed signal value is calculated, constructing an estimation polynomial at time£ given by P k (t) = σ k t+a k by calculating σ k and a k using effective regression equations, using said polynomial P k (t) to estimate the h integral f F(t)dt by —(P k {t k ) +P k ( t k ^)) to predict a value to be held and to J .z h estimate the integral f F(t)dt by — iP k { k. .) +P k (t k )), and using the polynomialP^JJ) from the previous interval to calculate a correction term, and adding the prediction and correction terms together to compute the output value y k given by * ( * ( * )+ * ( '™ ))+[ i^ ^ 7. The method according to claim 5 wherein the step of providing sample values y k of said signal function comprises providing said sample values at a rate slower than the rate at which said reconstructed signal function is outputted, letting h ι be a sampling interval for signal reconstruction and h o be a sampling interval for signal reconstruction, selecting the interval for signal reconstruction to beh o =rh t where r is an integer, approximating the integral f F(t)dt by a plurality of stepwise approximations using a boundary condition that the reconstruction at time t k+χ being given by y k , setting an equation of reconstruction for each step S equal to S =y k+l - p(r - 1 -/), / = 0,...,r-l , where p (y k . -y ) is an increment between steps given by p =2.0 — — , and wherein the step (r - l) of outputting said reconstructed analog signal includes outputting said reconstructed signal values S or storing said reconstructed signal values in said digital storage medium. 8. The method according to claim 1 wherein said sample values { y k } of the analog signal function of time F(t) are digitized from an audio analog signal so that said reconstructed signal is an audio signal. 9. The method according to claim 2 wherein said sample values { y k } of the analog signal function of time F(f) are digitized from an audio analog signal so that said reconstructed signal is an audio signal. 10. The method according to claim 3 wherein said sample values { y k } of the analog signal function of time F(t) are digitized from an audio analog signal so that said reconstructed signal is an audio signal. 11. The method according to claim 10 wherein said analog signal function is music, said digitized audio signals are stored in a storage medium and said reconstructed audio signal is music. 12. The method according to claim 11 wherein said storage medium include but are not limited to compact discs, digital audio tapes, magnetic discs, and magnetic optical discs. 13. The method according to claim 10 wherein said analog signal function is audio voice and said reconstructed audio signal is a voice signal. 14. A method for compression of an analog signal, comprising; a) sampling said analog signal at a first frequency v, or relieving from a storage medium an analog signal previously sampled at frequency v } to produce a digitized signal, filtering said digitized signal to remove frequencies above a preselected frequency, subsampling the filtered signal at a second frequency v 2 wherein v 2 <v p b) reconstructing said analog signal from said digitized signal sampled at frequency v 2 , by the steps of providing sample values y k of said analog signal of time (t)such that the sampling interval, t 0 is the initial time and £=0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k < t <t k constructing a sequence { y k } in real time where y k , £=0,1,2, ...μ are computed values of a piecewise constant function of time (t) forming a reconstruction of said signal function such that F t)=y k , t k < t< t k+l , £ = 0,1,2, ...μ, wherein said sequence { y k } is computed by using an effective numerical estimation technique so that an absolute difference given by I fP(t)dt- fF(t)dt\ locally in intervals defined by a preselected number of previous sample values y k is less than an absolute difference given by | ΪH(t)dt- fF(t)dt\ over the same intervals; and c) outputting said reconstructed analog signal P(t) or storing said sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form. 15. The method according to claim 14 wherein said analog signal is an audio signal. 16. The method according to claim 15 wherein each term y k is computed as a functional approximation of F(t) in an interval [t k , t k+ϊ ] by calculating an estimation function g k (t) that predicts F(t) in the interval [t k , t k+l ] by using an effective number of previous values y k ,y k _ , .- and l l ti ti b i i i f ti () t estimate the integral in the previous interval [t k _ t and subtracting the predicted area estimation in that interval, said correction being given by ad ng sai correcton term to said predicted value to give a value for y k given by <k>ι y k = -(fg k (t)dt + Je k (t)dt-fg k _.'t)dt) . 17. The method according to claim 16 wherein the estimation function g k (t) equals p"'(t), which is the polynomial approximation of order m that fits (t) at values y k ,y k . ■■■,y k . m . and wherein the estimation function e^t) equals p k (t), which is the polynomial approximation of order n that fits (t) at va\uesy k ,y k _ γ ,...,y k _ n . 18. The method according to claim 17 wherein m = 0 and n = 1 and then ~ 3 1 y k = ^y k --y k - l ■ 19. The method according to claim 17 wherein m= 1 and n= 1 and 20. The method according to claim 16 wherein said audio signal is music, said first frequency v 2 is about 44.1 kHz and said second frequency v 2 is about 22.05 kHz. 21. The method according to claim 17 wherein said audio signal is music, said first frequency V j is about 44.1 kHz and said second frequency v 2 is about 22.05 kHz. 22. The method according to claim 18 wherein said audio signal is music, said first frequency V j is about 44.1 kHz and said second frequency v 2 is about 22.05 kHz. 23. The method according to claim 19 wherein said audio signal is music, said first frequency v, is about 44.1 kHz and said second frequency v 2 is about 22.05 kHz. 24. The method according to claim 16 wherein said audio signal is voice, said first frequency . is about 16 kHz and said second frequency v 2 is about 8 kHz. 25. The method according to claim 17 wherein said audio signal is voice, said first frequency V j is about 16 kHz and said second frequency v 2 is about 8 kHz 26 The method according to claim 18 wherein said audio signal is voice, said first frequency v 2 is about 16 kHz and said second frequency v 2 is about 8 kHz 27 A method of reconstructing an analog signal from a digitized analog signal, comprising, providing sample values y k of an analog signal function of time (t)such that y k =F(t 0 +kh) where ? is the sampling interval, t 0 is the initial time and *=0,1,2, μ , constructing an output analog signal function (t)=y k , t k ≤ t < t k j , £ = 0,1,2, μ , by computing a sequence { y k } in real time, wherein each term >~ A is computed as a functional approximation of F(t) in an interval [t A , y the steps of calculating an estimation function g k (t) that predicts F(t) in the interval [t A , t Aτl ] using an effective number of previous values y k ,y k and integrating sa d estimation function over the interval [t k , t A+1 ] to give f g k (t)dt, and ' k calculating a predicted value to ho\dy k given by y k =— ( f g k (t)dt), h J calculating a correction term by using an estimation function e k (t) to estimate the integral in the previous interval [t k v t and subtracting the predicted area estimation in that interval, said correction being given by adding said correction term to said predicted value to give a value for y k given by 'JW h h Λ = ^ (fg k ( )Λ + fe k (t)dt - fg k . (t)dt), h '* ι h l outputting said reconstructed analog signal functιon£(t) or storing said sequence { y k } in a digital storage medium for subsequent relieving and outputting in analog form 28 The method according to claim 27 wherein the estimation function g k (t) equals p k m (t) , which is the polynomial approximation of order m that fιts (t) at values y k ,y k p ,y k „,, and wherein the estimation function e k (t) equals /? A "(t) , which is the polynomial approximation of order n that fιts (t) at values 29 The method according to claim 2 wherein said analog signal F(t) is a force signal acquired from a sensor means or a numerical simulation over time and the reconstructed signal functιon (t) ιs reflected to a human operator, and wherein h is the sampling interval for said force and the reconstructed signal function reflected to said operator 30 The method according to claim 2 wherein said analog signal function F(t) is a force signal acquired from a sensor means or a numerical simulation over time and the reconstructed signal functιon (t) ιs reflected to a human operator, and wherein the step of providing sample values y k of said analog signal function comprises providing said sample values at a rate faster than the rate at which said reconstructed signal function is outputted, letting ? ; be the sampling interval for force acquisition and h o be a sampling interval for force reflection, selecting the interval for force reflection to beh o =rh ι where r is an integer, letting t=t A denote the k th time that a reflected force value is calculated, constructing an estimation polynomial at tιme£ given by P k (t) = o k t +a by calculating σ k and a k using effective regression equations, using said to predict +P k (t k )), and using the . t) from the previous interval to calculate a correction term, and adding the prediction and correction terms together to compute the output value y k given by y W*) + ^A + .)) + [ ( ^*-i) + ^*))-^-i('i.-i) + ^-i(0)] • 31. The method according to claim 30 wherein said force signal is acquired from a sensor means or a numerical simulation over time. 32. The method according to claim 2 wherein said analog signal function F(t) is a force signal and the reconstructed signal function P(t) is reflected to a human operator, and wherein the step of providing sample values y k of said signal function comprises providing said sample values at a rate slower than the rate at which said reconstructed signal function is outputted, letting h t be a sampling interval force acquisition and h o be a sampling interval for output of forces to an activator forming part of a force reflection system, selecting the interval for force reflection to beh o =rh ι where / • is an integer, approximating the integral f F(t)dt by a plurality of stepwise approximations using a boundary condition that the force at time t k ] being given by y k+l , setting an equation of force for each step S equal to S =y k+l - p(r - 1 ~j),j = 0,...,r-\ (y k , -y.) ,where p is an increment between steps given by p = 2.0 , and wherein (r - \) the step of outputting said reconstructed signal function includes outputting said output force values S to said force reflecting system. 33. The method according to claim 32 wherein said force signal is acquired from a sensor means or a numerical simulation over time. 34. A method of determining a time-dependent three-dimensional affine transformation of a point x, of the form T(x,t) =R(t)x+d(t) where R is an orthogonal matrix having elements which vary according to sample values and d is a translation vector having elements which vary according to sample values, comprising: (a) computing a unit rotation axis vector r(t)= r y (t) and rotation angle θ(t) from said matrix R(t) ; (b) producing a reconstructed rotation vector axis r(t)-- r y (t) in which each element of the vector r(t), j=χ,y,z is reconstructed independently by: computing sample values y k of the function of time such that the sampling interval, t 0 is the initial time and £=0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k < t <t k+i ; and constructing a function r j t)=y k , t k < t < t k+l , k = 0,1,2, ...μ, by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fr(t)dt- fr(f)dt\ is less than an absolute difference given by | [H(ήdt- fr(t)dt\ locally in intervals defined by a preselected number of previous sample values^; (c) producing a reconstructed rotation angle θ(t) by: computing sample values y k of the function of time such that ^ A =θ(t 0 +£ ?) where/? is the sampling interval, t 0 is the initial time and £=0,l,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+1 ; and constructing a function Q(t)=y k , t k ≤ t< t k+i , £ = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fβ(t)dt- fθ(t)dt\ is less than an absolute difference given by locally in intervals defined by a preselected number of previous sample valuesjμ A ; (d) normalizing said reconstructed rotation axis, and computing a reconstructed matrix R(t)from the normalized reconstructed rotation axis; d x (i) (e) producing a reconstructed displacement vector d(t)-- d y (t) d (t) in which each element of the vector d(t), j=χ,y,z is reconstructed independently by: computing sample values y k of the function of time such that y k =d(t 0 +kh) where/? is the sampling interval, t 0 is the initial time and £=0,1,2, ...μ , defining a piecewise constant function of time given by H(t) =y k , t k ≤ t <t k+l ; and constructing a function d(t)=y k , t k < t < t k+1 , £ = 0,1,2, ...μ , by computing a sequence { y k } in real time using an effective numerical estimation technique so that an absolute difference given by | fd(t)dt- fd(t)dt\ is less than an absolute difference given by | fH(t)dt- fd(t)dt\ locally in intervals defined by a preselected number of previous sample values_y A , and (f) transforming said pointx by f(χ,t) = R(t) x +d(t) using said reconstructed R(t) and d(t) 35 The method according to claim 34 wherein the step of computing said rotation axis T and rotation angle θ(t) comprises the steps of defining a skew-symmetric matrix B = (R ~R τ )/2 where/? τ is a transposed matrix of the matrix R, writing B in a form 0 -b. b y B = b. 0 -b -b b 0 defining a scaled and unit rotation axis and a rotation angle respectively as (trace(R) - \) b = cos(θ)= so that said transformation T is a rotation about r by angle θ followed by a displacement d 36 The method according to claim 35 wherein the step of computing a reconstructed matrix R from a reconstructed rotation axis / τ (t)and a reconstructed rotation angle θ(t) comprises computing a reverse transformation given by (r x (\ - cos(θ))) + cos(θ) (r χ r y (l - cos(θ))) -r z sin(θ) (r χ f z (\ - cos(θ))) + /^ sin(θ) R = (;y ,(l - cos(θ))) +r z sin(θ) (r y 2 (\ - cos(θ))) +cos(θ) (r z (l - cos(θ))) -r x sin(θ) (f ,( 1 - cos(θ))) - r ; sin(θ) (r γ r z (\ - cos(θ))) + r χ sin(θ) ( ?( 1 - cos(θ))) + cos(θ) where 37. A method of point of view positioning in an interactive visual system, comprising: (a) providing a computer-generated visual scene of a real or artificial scene, which is computed or relieved from a storage medium; (b) sensing orientation and/or position sample values from which said computer-generated visual scene is to be computed; (c) determining a transformation of points in said scene, where each point x is transformed using a transformation T(x t) = R(t) x ~ ι +d(t) where R is an orthogonal matrix the elements of which vary according to said orientation and position sample values and where d is a translation vector the elements of which vary according to said orientation and/or position sample values; (e) computing a reconstructed rotation matrix R(t)and a reconstructed translation vector d(t) according to the method of claim 34 so that R(t) and d(t) are digital representations of reconstructed analog signals; and (f) rendering the transformed points f(x ι ,t) to a display means by using said reconstructed analog signals in digital or analog form. 38. A method for point of view positioning in an interactive visual system according to claim 37 wherein an orientation and/or position of an anatomical part of a person or animal is sensed by a mechanical linkage, acoustic range sensing, sensing of magnetic fields, sensing of electric fields, sensing of reflected light, sensing of emitted light or combinations thereof and in which the orientation and/or position from which the visual scene is to be rendered is computed from said orientation and/or position of said anatomical part. 39. A method of point of view positioning in an interactive visual system according to claim 38, in which the anatomical part is a head, an eye, a hand, a finger, a foot, a toe or any combination thereof of a user. 40. A method according to claim 38 or 39 wherein hand and/or foot motions are sensed by means of steering wheels, joysticks or foot pedals. 41. An apparatus for orienting and/or positioning a remote sensor for display, comprising (a) sensing means for remote sensing; (b) actuator means for mechanically, electrically, magnetically, hydraulically, or pneumatically changing the orientation and/or position of said remote sensing means; (c) sensing or determining a desired orientation and/or position of said remote sensing means as sample values; (d) display means connected to said remote sensing means; (e) a processor means operable to determine a transformation of orientation of said remote sensing means at a given time t as a rotation R(t) and/or a position of said remote sensing means as a translation vector d(t) , where R is an orthogonal matrix the elements of which vary according to said orientation and position sample values and where d is a translation vector the elements of which vary according to said orientation and position sample values; (f) said processor means operable to compute a reconstructed rotation matrix R(t) and a reconstructed translation vector d(t) according to the method of claim 34 so that R(t) and d(t) are digital representations of reconstructed analog signals; and (g) actuator means connected to said processing means for changing the orientation and/or position of said remote sensing means to an orientation specified by said R(t) and/or the position specified by said d(t) by using said reconstructed analog signals in digital or analog form. 42. An apparatus according to claim 41 wherein the remote sensor is a camera, radio detecting and ranging device, laser detecting and ranging device, ultrasound detecting and ranging device, force sensor or tactile sensor. 43. An apparatus according to claims 41 or 42 wherein orientation and/or position of an anatomical part of a person or animal is sensed by a mechanical linkage, acoustic range sensing, sensing of magnetic fields, sensing of electric fields, sensing of reflected light, sensing of emitted light or combinations thereof and in which the orientation and/or position of the remote sensor is computed from said orientation and/or position of said anatomical part. 44. An apparatus according to claim 43 wherein the anatomical part is a head, an eye, a hand, a finger, a foot, a toe or any combination thereof of a user. 45. An apparatus according to claim 43 wherein hand and/or foot motions are sensed by means of steering wheels, joysticks or foot pedals. 46. An apparatus according to claims 41 , 42, 43, 44 or 45 wherein the means of orienting the remote sensor is a robotic linkage, a mobile vehicle, an aircraft, a spacecraft, a submersible vehicle, magnetic levitation or electric levitation.

Inventors: ELLIS RANDY E, JENKINS MICHAEL A, SARKAR NILANJAN, ELLIS, RANDY, E., JENKINS, MICHAEL, A., SARKAR, NILANJAN
Applicants: ELLIS RANDY E, JENKINS MICHAEL A, SARKAR NILANJAN, ELLIS, RANDY, E., JENKINS, MICHAEL, A., SARKAR, NILANJAN
Application year: 1997
Main classification: H03M1
IPC codes: G11B20, H03M1