AU2003288416A1

AU2003288416A1 - Contact sensitive device

Info

Publication number: AU2003288416A1
Application number: AU2003288416A
Authority: AU
Inventors: Nicholas Patrick Roland Hill; Darius Martin Sullivan
Original assignee: New Transducers Ltd
Current assignee: NVF Tech Ltd
Priority date: 2002-12-06
Filing date: 2003-12-03
Publication date: 2004-06-30
Also published as: KR101121891B1; JP2006509301A; GB0228512D0; EP1570414A2; TWI348115B; BR0317049A; CN1720498A; WO2004053781A2; CN100405266C; KR20050088103A; CA2506789A1; AR042333A1; MXPA05006040A; TWI308709B; JP4613068B2; WO2004053781A3; TW200817987A; TW200424916A

Description

WO2004/053781 PCT/GB2003/005268 TITLE: CONTACT SENSITIVE DEVICE 10 DESCRIPTION 15 TECHNICAL FIELD The invention relates to contact sensitive devices. BACKGROUND ART Visual displays often include some form of touch sensitive screen. This is becoming more common with the 20 emergence of the next generation of portable multimedia devices such as palm top computers. The most established technology using waves to detect contact is Surface Acoustic Wave (SAW), which generates high frequency waves on the surface of a glass screen, and their attenuation by the 25 contact of a finger is used to detect the touch location. This technique is "time-of-flight", where the time for the WO2004/053781 PCT/GB2003/005268 2 disturbance to reach one or more sensors is used to detect the location. Such an approach is possible when the medium behaves in a non-dispersive manner i.e. the velocity of the waves does not vary significantly over the frequency range of interest. 5 In W001/48684 and PCT/GB2002/003073 to the present applicant, two contact sensitive devices and methods of using the same are proposed. In both applications, the device comprises a member capable of supporting bending wave vibration and a sensor mounted on the member for measuring 10 bending wave vibration in the member and for transmitting a signal to a processor whereby information relating to a contact made on a surface on the member is calculated from the change in bending wave vibration in the member created by the contact. 15 By bending wave vibration it is meant an excitation, for example by the contact, which imparts some out of plane displacement to the member. Many materials bend, some with pure bending with a perfect square root dispersion relation and some with a mixture of pure and shear bending. The 20 dispersion relation describes the dependence of the in-plane velocity of the waves on the frequency of the waves. Bending waves provide advantages, such as increased robustness and reduced sensitivity to surface scratches, etc. However, bending waves are dispersive i.e. the bending wave 25 velocity, and hence the "time of flight", is dependent on frequency. In general, an impulse contains a broad range of WO2004/053781 PCT/GB2003/005268 3 component frequencies and thus if the impulse travels a short distance, high frequency components will arrive first. In WO01/48684 and PCT/GB2002/003073 a correction to convert the measured bending wave signal to a propagation signal from a 5 non-dispersive wave source may be applied so that techniques used in the fields of radar and sonar may be applied to detect the location of the contact. DISCLOSURE OF INVENTION According to one aspect of the invention, there is 10 provided a contact sensitive device comprising a member capable of supporting bending waves, three sensors mounted on the member for measuring bending wave vibration in the member, whereby each sensor determines a measured bending wave signal and a processor which calculates a location of a contact on 15 the member from the measured bending wave signals, characterised in that the processor calculates a phase angle for each measured bending wave signal, then calculates a phase difference between the phase angles of least two pairs of sensors from which the location of the contact is determined. 20 According to a second aspect of the invention, there is provided a method of determining information relating to a contact on a contact sensitive device comprising the steps of providing a member capable of supporting bending waves and three sensors mounted on the member for measuring bending wave 25 vibration in the member, applying a contact to the member at a location, using each sensor to determine a measured bending WO 2004/053781 PCT/GB2003/005268 4 wave signal and calculating the location of a contact from the measured bending wave signal characterised by calculating a phase angle for each measured bending wave signal, calculating a phase difference between the phase angles of at least two 5 pairs of sensors and determining the location of the contact from the at least two calculated phase differences. The following features may be applied to both the device and the method with the processor being adapted to provide many of the calculations or processing steps of the method. 10 Reflected waves may be suppressed by placing an absorber in contact with the edges of the member. The mechanical impedance of the absorber and member may be selected so as to minimise reflections of bending waves from the edges of the member. In particular, the impedances may be selected such 15 that bending wave energy is strongly absorbed in a frequency band around a chosen frequency co. The impedance of the absorber may be selected to be both resistive and compliant. The impedances may be selected to satisfy the following equation: 20 ZT =-iZ (oo) where ZT is the termination impedance of the absorber and ZB is the mechanical impedance of the edge of the member. The absorber may be made from foamed plastics which may have open or closed cells and may be polyurethane or 25 polyvinylchloride. For example the foam may be a soft PVC WO 2004/053781 PCT/GB2003/005268 5 predominantly closed cell foam such as MIERSTM or a medium to high density, open cell polyurethane foam. Another class of foams that has been found to be suitable are acrylic closed cell foams. These may have a high degree of damping and 5 relatively high stiffness. Such properties are particularly suited to edge termination of stiff, heavy materials such as glass. Examples include 3M serial numbers 4956, 4910, 4950 and 4655. The absorber may extend substantially around the periphery of the member. The absorber may act as a mounting 10 which supports the member in a frame or to another surface. The member may comprise a raised pattern on its surface whereby a contact drawn across the surface provides a variable force to the member to generate bending waves in the member. The pattern may be periodic, or quasi-periodic with a 15 statistically well-defined spatial distribution of undulations. The pattern may be random whereby a contact travelling over the surface of the member generates a random bending wave signal. The random relief pattern may be an anti reflective coating, an anti-glare surface finish or an etched 20 finish such as those which are found on many known transparent panels placed in front of electronic displays. Each measured bending wave signal may be processed by a band-pass filter with a pass-band centred at the chosen frequency w0o and having a bandwidth Aco. The bandwidth Aco of 25 the filter is preferably chosen to address the Doppler effect WO2004/053781 PCT/GB2003/005268 6 whereby a bending wave arrives at a point with a different frequency from its original frequency. Accordingly, the bandwidth preferably obeys the relationship: Aco >> 2k(O o ma 5 where Vm. is the maximum lateral velocity of the contact across the surface, e.g. if the contact is provided by a stylus, va is the maximum velocity that a user is capable of moving the stylus. The phase of each filtered signal may be measured by 10 comparison with a reference signal. The reference signal may have a frequency c0. The measured phase is the average phase difference between the input and reference signals, optimally measured over the interval 2n/Aw. Alternatively, the reference signal may be derived from a filtered signal from a second 15 sensor in which case, the measured phase is the phase difference between two input signals. The phase differences may be calculated at intervals of 2n/A(o which may be intervals of less than 10ms. The reference and input signals may be fed to a phase detector. Output from 20 the phase detectors may be fed through low-pass filters having frequency cut-offs of approximately Ac/2, then through digitisers and finally through a processor to calculate the phase angle 0. The instantaneous phases, O,(t) and Om(t), of two measured 25 bending wave signals may satisfy the phase difference equation WO2004/053781 PCT/GB2003/005268 7 equation: AO,,, = 0t -m = k(o o )Ax, +2mz,, where Axm = x, - Xm (Xm and x, being the distance from the contact location to each sensor labelled m and 1 5 respectively), and k(o) is the wavevector. This equation may be satisfied if the path length difference between two sensors is less than the coherence length of the bandpass filter, which is defined as A COk(co) 2# 00 10 The coherence condition is therefore Ax,,, <<x. If the coherence condition is not satisfied, the above phase equation may not be satisfied. Thus, values of n, and the phase angle difference are required to determine the location of the contact. The shape 15 of the member may be selected to constrain the magnitude of Axj, to values less than half of one wavelength, ie. Ax,.<,/k(wo). In this case, where all possible values of Axim satisfy the condition |Axi<r/k(mo), there is only one value of n~, which is the integer n, satisfying -2ml 20 Alternatively, n may be estimated or inferred in some way. Each phase angle difference in combination with the range of possible values of the integer na may be used to generate a series of path length differences thereby defining a series WO2004/053781 PCT/GB2003/005268 8 of discrete hyperbolic curves on the surface of the member, denoting possible locations of the contact. The location of the contact may be determined by plotting each hyperbola defined by each path length difference and selecting a point 5 at which a large number of the hyperbolae intersect or nearly intersect. This point is likely to be the true location of the contact. Where nlm is unknown, the minimum number of series of hyperbolae required to determine the contact location is three 10 and the likelihood of determining the correct location of the contact is increased by increasing the number of hyperbolae to be plotted. Multiple sensors may be used whereby a phase angle difference may be calculated for each pair of sensors thus generating multiple hyperbolae. In this embodiment, the 15 minimum number of sensors is three. Alternatively, where n is unknown, the measured bending wave signal from each sensor may be divided into two or more discrete frequency bands whereby a phase angle difference may be calculated for each frequency band and for each pair of 20 sensors. Although multiple phase angle differences may be calculated from a single pair of sensors, the phase angle differences at different frequencies are derived from the same path length difference. Thus the minimum number of sensors is three. The dividing of the frequency bands may be achieved by 25 processing the bending wave signals by at least two band-pass filters having different pass-band frequencies. For example, WO 2004/053781 PCT/GB2003/005268 9 using two band-pass filters having frequencies oo+ms and coo-@, the phase angle differences AOa, AOb from two sensors may be defined as AO. = k(co +co )Ax + 2ml 5 AO = k(o -m )Ax+2 2Ib where Ax is a single path-length difference defined by the contact and the position of the sensors. Therefore the values of na and nb may be selected so that the measured phase angle differences infer similar values of 10 the path-length difference. There may be only one combination of values (na, nb) for which this is possible. In this case the true value of the path-length difference may be determined. The correct combination (na, nb) may be determined as the combination of values that minimise the expression: 15A0, - 27m, A0 b

-

2 mnb I5 k(co +o) k(coo-co) The path length difference may then be estimated as: Ax_1 AO,-2im. A b-2mlb 2(mo+m,) k(mo-o)5 Where this process is repeated with two pairs of sensors, two path length differences may be determined, which in turn 20 may be used to determine the location of the contact. Alternatively, where na is unknown, an initial determination of the location of the contact may be made using the methods taught in W001/48684 and PCT/GB2002/003073 (as WO2004/053781 PCT/GB2003/005268 10 summarised in Figure 11). Thereafter it may be assumed that the contact moves more slowly than the bending waves and hence that the phase angle differences change by small increments over the timescale At. Accordingly, each value of n may be 5 chosen to minimise the change in path length difference. Measured phase angle differences may contain random errors which may result in the selection of incorrect values of n. This error may be alleviated by evaluating the likelihood of successive sequences of n, for example by a 10 state-space estimator such as the well known Kalman filter. The sequence having the maximum measure of likelihood is selected. The state-space estimator provides an estimate of the internal state of a system of which noisy measurements are 15 made. A necessary input to the state-space estimator is a statistical description of the evolution of the system state. An example of such a state is the set of coordinates that describes the position and velocity of an object in contact with the member. It is widely known that the Kalman filter and 20 other state-space estimators may provide a measure of likelihood that a sequence of observed, noisy measurements are consistent with the model of the system state. A state-space estimator may therefore be employed to take a sequence of a pair of path-length differences (say Ax 12 and 25 Ax34) taken at different times (say tL, t2, t 3 ,...), to estimate WO 2004/053781 PCT/GB2003/005268 11 the system-state, i.e. the position and velocity of the contact, at those times. Moreover, the overall likelihood of those values of path-length difference being consistent with the model of the system may be evaluated. 5 Where the sequence of path-length differences are obtained from a sequence of phase-angle differences and a set of integers (n = n(t) , n ( t 2 ), n ( t 3 ),...), the measure of likelihood generated by the state-space estimator may be used to infer the likelihood that the correct values of n were 10 chosen. It follows that a method for choosing the correct sequence of integers, n, is to find the sequence for which the state-space estimator gives the maximum measure of likelihood. As mentioned above, the state space estimator uses some statistical description of the evolution of the system state. 15 A suitable model for the motion of the contact may be a simple random walk. Alternatively, the model may employ a detailed statistical description of how the user moves the stylus or finger. One example is a statistical description of how the user moves a pen while writing text or individual characters. 20 The processor may further be adapted to include in the determination procedure any available information about where the contact can be expected. For example, if the member is an input device for a graphical user interface where the user is presented with a choice of 'buttons' to press, it may be 25 useful to assume that any contact on the member occurs within the discrete areas corresponding to the buttons.

WO 2004/053781 PCT/GB2003/005268 12 Alternatively, a map of the probability at which a contact is likely to occur and which is based on the expected behaviour of the user may be used. The device may comprise a software application with a graphical user interface (GUI) 5 which interacts with the operating system by means of an application program interface (API) in which the API is adapted to generate the probability map. The probability map may be based on the location, size, and frequency of use of objects presented by the graphical user interface. The 10 probability map may also be based on information about the relative likelihood of the various GUI elements being activated. The following characteristics may apply to all embodiments of the invention. The device may comprise means 15 for recording measured bending wave signals from the or each sensor over time as the contact moves across the member. The information relating to the contact may be calculated in a central processor. The sensors may be mounted at or spaced from an edge of the member. The sensors may be in the form of 20 sensing transducers which may convert bending wave vibration into an analogue input signal. The member may be in the form of a plate or panel. The member may be transparent or alternatively non-transparent, for example having a printed pattern. The member may have 25 uniform thickness. Alternatively, the member may have a more complex shape, for example a curved surface and/or variable WO2004/053781 PCT/GB2003/005268 13 thickness. The device may be a purely passive sensor with the bending wave vibration and hence the measured bending wave signals being generated by an initial impact or by frictional 5 movement of the contact. The contact may be in the form of a touch from a finger or from a stylus which may be in the form of a hand-held pen. The movement of a stylus on the member may generate a continuous signal which is affected by the location, pressure and speed of the stylus on the member. The 10 stylus may have a flexible tip, e.g. of rubber, which generates bending waves in the member by applying a variable force thereto. The variable force may be provided by tip which alternatively adheres to or slips across a surface of the member. As the tip moves across of the member a tensile 15 force may be created which at a certain threshold, causes any adhesion between the tip and the member to break, thus allowing the tip to slip across the surface. The bending waves may have frequency components in the ultrasonic region (>20 kHz). 20 The member may also be an acoustic radiator and an emitting transducer may be mounted to the member to excite bending wave vibration in the member to generate an acoustic output. The frequency band of the audio signal of the transducer preferably differs from and does not overlap the 25 frequency band of the measurements from the sensors. The audio signal may thus be filtered, for example, the audio band WO 2004/053781 PCT/GB2003/005268 14 may be limited to frequencies below 20kHz, and the vibration measurements may be limited to frequencies above 20kHz. A sensor may have dual functionality and act as the emitting transducer. 5 The or each emitting transducer or sensor may be a bender transducer which is bonded directly to the member, for example a piezoelectric transducer. Alternatively, the or each emitting transducer or sensor may be an inertial transducer which is coupled to the member at a single point. The inertial 10 transducer may be either electrodynamic or piezoelectric. A contact sensitive device according to the invention may be included in a mobile phone, a laptop or a personal data assistant. For example, the keypad conventionally fitted to a mobile phone may be replaced by a continuous moulding which is 15 touch sensitive according to the present invention. In a laptop, the touchpad which functions as a mouse controller may be replaced by a continuous moulding which is a contact sensitive device according to the invention. Alternatively, the contact sensitive device may be a display screen, e.g. a 20 liquid crystal display screen comprising liquid crystals which may be used to excite or sense bending waves. The display screen may present information relating to the contact. BRIEF DESCRIPTION OF DRAWINGS The invention is diagrammatically illustrated, by way of 25 example, in the accompanying drawings, in which: Figure 1 is a schematic plan view of a touch sensitive WO 2004/053781 PCT/GB2003/005268 15 device according to one aspect of the invention; Figure 2 is a schematic perspective view of the device of Figure 1; Figure 3 is a schematic side view of a one-dimensional 5 beam; Figure 4a is a graph showing the amplitude of the reflection coefficient against frequency (Hz), the amplitude is unitless since it is a ratio; Figure 4b is a graph showing the phase (in radians) of 10 the reflection coefficient against frequency (Hz); Figures 5a and 5b are schematic perspective views of alternative touch sensitive devices; Figure 6 is a flowchart of a method of finding the location of a contact according to the invention; 15 Figure 7a is a schematic block diagram of apparatus used for calculating phase angles; Figure 7b is a schematic block diagram of apparatus used with that of Figure 7a; Figures 8a to 8d are plan views of apparatus according to 20 the invention showing the hyperbolae of path length differences; Figure 9 is a schematic block diagram of alternative apparatus used for calculating phase angles; Figure 10 is a flow chart showing an alternative method 25 of calculating the location of the contact; Figure 11 is a flow chart showing a method of calculating WO2004/053781 PCT/GB2003/005268 16 the location of the contact using the dispersion corrected correlation function; Figure 11a is a graph of dispersion corrected correlation function against time, and 5 Figure 12a is a schematic block diagram of a contact sensitive device which also operates as a loudspeaker, and Figure 12b is a method of separating audio signal and measured bending wave signal in the device of Figure 12a. DETAILED DESCRIPTION 10 Figure 1 shows a contact sensitive device 10 comprising a transparent touch sensitive plate 12 mounted in front of a display device 14. The display device 14 may be in the form of a television, a computer screen or other visual display device. A stylus 18 in the form of a pen is used for writing 15 text 20 or other matter on the touch sensitive plate 12. The transparent touch sensitive plate 12 is a member, e.g. an acoustic device, capable of supporting bending wave vibration. As shown in Figure 2, four sensors 16 for measuring bending wave vibration in the plate 12 are mounted 20 on the underside thereof. The sensors 16 are in the form of piezoelectric vibration sensors and are mounted one at each corner of the plate 12. At least one of the sensors 16 may also act as an emitting transducer for exciting bending wave vibration in the plate. In this way, the device may act as a 25 combined loudspeaker and contact sensitive device. A mounting 22 made of foamed plastics is attached to the WO2004/053781 PCT/GB2003/005268 17 underside of and extends substantially around the periphery of the plate 12. The mounting 22 has adhesive surfaces whereby the member may be securely attached to any surface. The mechanical impedance of the mounting and plate are selected so 5 as to minimise reflections of bending waves from the plate edges. The relationship between mechanical impedance of the mounting and the plate may be approximated by considering the one dimensional model shown in Figure 3. The model comprises a 10 waveguide 34 in the form of a beam which terminates at an edge mounting 36 having a termination impedance. An incident wave 38 travelling down the waveguide 34 is reflected by the mounting 36 to form a reflected wave 40. The incident and reflected waves are plane waves travelling in the direction 15 perpendicular to the edge. Assuming the mounting 36 satisfies the following boundary conditions: (i) the termination impedance only couples into the lateral velocity, i.e. it does not provide any torque resistance; whereby the bending moment is equal to zero at the edge and 20 (ii) The ratio of the lateral shear force and the velocity at the edge is equal to the terminal impedance; the reflection coefficient at the mounting is given by: R () = - ZT/ZB(m)-i Z- /ZR (O) + 1 where ZT is the termination impedance of the mounting and 25 ZB is the mechanical impedance of the end of the waveguide, WO 2004/053781 PCT/GB2003/005268 18 given by Bk (o)(1) 2a where k(o) , is the wavevector which may be expressed in terms of the bending stiffness, B, and mass per unit area, 9, 5 of the panel, k1/ B Thus, the reflection coefficient is determined by the ratio of the impedances at the end of the waveguide and the mounting. Furthermore, the impedance of the waveguide is 10 proportional to the square root of frequency and is both real and reactive in equal weights (i.e. n/4 phase angle). Accordingly, the reflection coefficient is likely to be strongly frequency dependent. The reflection coefficient vanishes, i.e. bending wave 15 energy is strongly absorbed in a frequency band around mo, if the following condition is satisfied: Zr = -iZB CO) Thus, the termination impedance of the mounting must have both real and imaginary components, or, equivalently, the 20 mounting should be both resistive and compliant. The plate may be, for example, 1 mm thick polycarbonate sheet which has mass per unit area, p=1.196 kg M -2 and bending stiffness, B=0.38 Nm. The equations above can be used WO 2004/053781 PCT/GB2003/005268 19 to calculate the impedances of the plate and absorber required to strongly absorb bending wave energy around the chosen angular frequency so = 2n(900 Hz). The impedance, per unit width for a 1mm beam 5 approximation of the plate is -2 ZB(w0) = (1+i) 33.8 N s m -2 The properties of the absorber which provide the desired absorption are thus: Resistance per unit width, 10 Re(ZT )= Im[Z (coO )]=33.8Nsm - 2 Stiffness per unit width, -ilm(Z)coo = Re[ZB(co)]co = 1.91x105 Nm - 2 . The reflection coefficient is a unitless complex number. Figures 4a and 4b are graphs showing the amplitude and phase 15 of the reflection coefficient R(w) varying with frequency. The amplitude of the reflection coefficient is zero and its phase is reversed for oo approximately equal to 900Hz. In Figures 5a and 5b, the plate 12 has uniform surface roughness in the form of a raised surface pattern 28,29. The 20 stylus 18 is drawn across the surface along a path 30 and as it crosses a raised part or line of the pattern it generates bending waves 32 in the member. Thus contact from the stylus 18 provides a source of bending wave vibration in the member. In Figure 5a, the surface pattern 28 is a periodic pattern of 25 raised crossed lines and in Figure 5b, the surface pattern 29 WO2004/053781 PCT/GB2003/005268 20 is a random relief pattern. In the embodiments of Figures 2, 5a and 5b, as the contact moves over the rough surface of the member, bending waves radiate isotropically in the member from the point of 5 contact. The displacement of the member at a distance, x, from the point of contact is related to the displacement at the point of contact by a transfer function, H(c;x). At distances larger than the wavelength, A=27/k(m), the transfer function can be approximated as, 10 H(c;x) A ex where A is a constant and k(o), is the wavevector defined previously. Although H(a>;x) strictly only applies to bending waves on an infinite plate, since the mounting is strongly absorb bending wave vibrations, the relationship is satisfied. 15 The transfer function shows that where a source of bending waves emits a purely sinusoidal frequency with angular frequency, cOo, the phase difference A012 between displacements at two locations which are at distances, x, and x 2 , from the point of contact for the source is: 20 exp(iA12)= exp[ik(o 0 Xxl- x 2 )] This implies the following relationship between the phase angle difference, the path length difference Ax=(x,-x) and an integer n 12 . A0 2 = 0 -0 2 = k(mo)Axz 12 + 2n 2 WO 2004/053781 PCT/GB2003/005268 21 Figure 6 shows the steps in the method for using this equation to determine the contact location: a) Measure a bending wave signal with each sensor to give measured bending wave signals W.(t) and W(t), 5 b) Calculate the phase angles 0,(t) and 0,(t) of the measured bending wave signals Wi(t) and W.(t), c) Calculate the difference between the two phase angles 0,(t) and O (t) d) Calculate the location of the contact from 10 k(coo)Axy = AO,, - 2 zn Figure 7a shows a schematic block diagram of a device for calculating the phase angle 6j of a bending wave signal W (t) measured by one of the sensors. The signal Wj (t) is a random signal and is thus uncorrelated over long time scales. The 15 signal is first amplified by an amplifier 42 and then processed by an analogue band-pass filter 44 with a pass-band centred at co and a bandwidth of Aco. A moving source of bending waves may demonstrate the Doppler effect, whereby a bending wave which has frequency m 0 20 and is emitted by a source moving at velocity v towards a point on a member arrives at that point with a different frequency defined by co-k(o 0 j,. The maximum angular frequency shift between bending waves at two different points on the member is therefore 2k( 0 x) , where v, is the maximum WO 2004/053781 PCT/GB2003/005268 22 velocity of the moving source. If the angular frequency shift becomes larger than the width of the band pass filter, the phase difference equation above does not hold. Accordingly, the bandwidth Aco of the filter 44 is set to be greater than 5 this maximum frequency shift and thus obeys the relationship: Aco >> 2k(o)mV After processing by the filter 44, the resulting filtered signal W'j(t) is an amplitude and phase modulated carrier with frequency a and is defined by: 10 i(t) = A, (t)sin kot +8 0t where A-(t) and O (t) are the amplitude and phase of the signal. Both fluctuate over a timescale At determined by the bandwidth of the filter, namely At=2x/Am. The maximum frequency at which independent phase angle measurements may be 1 15 taken from the output of the bandpass filter is Since a touch sensor typically provides an updated measurement of the contact position every 10ms, the condition for the minimum frequency of positional measurement is At<10ms. The filtered signal W'j(t) is then passed simultaneously 20 to two analogue phase detectors 46. Such detectors are well known in the art, for example, see p644 of "The Art of Electronics" by Horowitz and Hill. Reference signals each having frequency coo but a phase difference of n/2 are also fed to the two phase detectors. The outputs of the phase WO2004/053781 PCT/GB2003/005268 23 detectors are passed through low-pass filters 48 each having frequency cut-offs of approximately Aco/2. The outputs of the low-pass filters are proportional to cos(89) and sin(O8) respectively. These outputs are then digitised by digitisers 5 50 and processed by processor 52 to give the phase angle O. Figure 7b shows how the reference signals used in Figure 7a may be generated. A second bending wave signal Wi(t) is measured at a second sensor. The signal is fed through an amplifier 42 and analogue band-pass filter 44 to generate a 10 filtered signal W'j(t). The filtered signal W'%(t) forms the reference signal which is fed directly to one phase detector 46. The filtered signal is also fed to the second phase detector 46 via a device which shifts its phase by n/2. The phase shifted signal is used as the reference signal to the 15 second phase detector 46. Figures 8a to 8d show how the phase angle differences and hence the path length differences may be used to calculate the location of the contact. The equation in step (d) of Figure 6 defines a hyperbolic curve which can be overlaid on 20 the plate 12. Figure Ba shows the three hyperbolic curves 26 which are generated using three different values of nim and the calculated phase angle difference for a pair of sensors 16 mounted one on each end of the short sides of the plate 12. Similarly Figures 8b and 8c show the hyperbolic curves 26 25 which are generated by the phase angle difference and WO2004/053781 PCT/GB2003/005268 24 different values of nlm for two other pairs of sensors. Figure Sd shows all the hyperbolic curves created by the sensors. The contact location 24 is the point of intersection of three hyperbolic curves, one from each pair of sensors. 5 From the contact location 24, the correct value of njm may be inferred. A method of inferring n is implemented using the embodiment shown in Figure 9. The bending wave signal W(t) measured by each sensor is simultaneously processed by two 10 band-pass filters 48,54. Two phase angles, one for each filter, are calculated, for example as described in Figure 7. The filters 48, 54 have slightly different pass-band frequencies whereby two phase angle differences, one for each pass-band frequency, are provided by each pair of sensors. 15 The phase angle differences AOa, Ab from the sensors may be defined as AO, = k(wo + c)Ax+2-na AOb = k(o -ms)Ax +2nb where Ax is a single path-length difference defined by 20 the contact and the position of the sensors. The correct combination (na, nb) may be determined as the combination of values that minimise the expression: AO,-2nna AOb2-2nZnb k(coo +pC) k(co - ) The path length difference may then be estimated as: WO 2004/053781 PCT/GB2003/005268 25 AX1 A 0-2m, AOb- 2 moib 2 k(co+co,) k(mo-co) Another pair of sensors may then be used to determine a second path length difference. Each path length difference defines a hyperbolic curve on the panel. The intersection 5 point of these two hyperbolic curves is the location of the contact. As in Figures 8a to 8d, the hyperbolae are plotted and the point at which the largest number of hyperbolae intersect is likely to be the true location of the contact. Figure 10 shows an alternative method for calculating the 10 location of the contact from the equation above, namely i. Measure a pair of bending wave signals W(t) and W(t), one signal being measured by a sensor; ii. Calculate the dispersion corrected correlation function of the two signals using the method described in Figures 15 11 and lla; iii. Calculate the initial position of the contact using the dispersion corrected correlation function, as described in Figures 11 and 1la; iv. Remeasure bending wave signals WH(t) and W(t); 20 v. Calculate the phase angle of each signal - for example as described in Figures 7a and 7b; vi. Calculate the difference between the phase angles; vii. Select the value of nim which minimises the change in the path length difference; WO 2004/053781 PCT/GB2003/005268 26 viii.Plot the hyperbola defined by k(eo)Ax, = AOU - 2my ix. Repeat steps (iv) to (viii), remeasuring the bending wave signals at regular intervals At, for example At=2t/Ao. 5 At step (viii) a minimum of two hyperbolae from different pairs of sensors are required to determine the position of the contact. Therefore the entire process must be performed simultaneously for at least two pairs of sensors. Thus a minimum number of two phase angle differences must be 10 determined. Two phase angle differences may be generated by using two sensors and splitting the signal into two frequency bands as described in Figure 9. Alternatively, multiple sensors may be used so that multiple phase angle differences may be calculated using different pairs of sensors. 15 Figure 11 shows a method of calculating the dispersion corrected correlation function to reveal the difference in path length between the contact location and the sensors. The method set out below summarises the information in PCT/GB2002/003073. The method comprises the following steps: 20 (a) Measure two bending wave signals W(t) and W(t); (b) Calculate the Fourier transform of W,(t) and W 2 (t) to arrive at (o) and ,W 2 (m) and hence the intermediate function W,(w) W2(); where W() is the complex conjugate Fourier transform, t represents time o is 2f where f is frequency.

WO 2004/053781 PCT/GB2003/005268 27 (c) Calculate a second intermediate function M(w) which is a function of '(CA)W2(w) (d) and (e) at the same time as performing steps (a) to (c), the frequency stretching operation f(co)=v(p/B) Yj is 5 calculated using the predetermined panel dispersion relation k = (p/B) XVJ,;. (f) M(o0) and f()v(pB)4 are combined to arrive at the dispersion corrected correlation function: G(t)- M[f()]exp(ot)do; and 27cm 10 (g) the dispersion corrected correlation function is plotted against time with a peak occurring at time t 12 as shown in Figure 11a; (h) Ax2 is calculated from tz 2 ; Ax 12 is the path-length difference between the path lengths x, and x 2 from the first 15 and second sensors to the contact. (i) Ax 12 defines a hyperbolae which may be plotted as in Figure 7 to calculate the location of the contact. As with the method of Figure 10, a minimum of two hyperbolae are required to determine the location of the 20 contact. Thus the ways of generating more hyperbolae discussed above apply equally to this method. The second intermediate function M(o) may simply be W,() W*() which gives a standard dispersion corrected WO 2004/053781 PCT/GB2003/005268 28 correlation function. Alternatively, M(co) may be selected from the following functions which all yield phase equivalent functions to the standard dispersion corrected correlation function: 5 a) M(co) = ()W^ ,) b) M(o) = W, (CO) W 2 (CO) ri, (CO) P ; (CO) c) M(co) = W(m) () (co)'2() where <o(x) is a real valued function d) M(co) = () (o)y(c) where V(co) is a real valued 10 function Alternatively, M(co) may be the function D(0) which is the Fourier transformation of the correlation function D(t): D(t) = W,(t + t')W (t')dt' The steps are calculate D(t); calculate D(0) and apply a 15 frequency stretching operation to arrive at the dispersion corrected correlation function: G(t)=-L fbL)[f(co)]exp(ict)dmo. Alternatively, at step (f) the following dispersion corrected correlation function may be calculated: w G(t)= F[f(co) *f (co)] q 12 [f(co)] exp(icot)dco 2we 20 where WO 2004/053781 PCT/GB2003/005268 29 Al(c)= .()W j(w)exp[- ik(co)Axj II where {j()} and W2(JCO) are the Fourier transformation and complex conjugate Fourier transformation of two measured bending wave signals {WIj(t)} and W2(t) and {Ax} is the path 5 length difference. A sensor may act as both the first and second sensor whereby the dispersion corrected correlation function is an autocorrelation function. The autocorrelation function may be calculated applying the same steps for the dispersion 10 corrected correlation function using W(t)= W 2 (t) Figure 12a shows a contact sensitive device which also operates as a loudspeaker. Figure 12b shows a method for partitioning the audio signal and measured signal into two distinct frequency bands so that the contribution of the audio 15 signal to the processed measured signal is suppressed. The device comprises a member 106 in which bending waves are generated by an emitting transducer or actuator 108 and the contact. The emitting transducer applies an audio signal to the member 106 to generate an acoustic output. Before being 20 applied to the member, the audio signal is filtered by a low pass filter 112 which, as shown in Figure 12b, removes the audio signal above a threshold frequency f 0 . As shown in Figure 12b, the contact generates a signal which has a power output which is substantially constant over WO2004/053781 PCT/GB2003/005268 30 a large frequency band. The signal from the contact and the audio signal sum to give a combined signal which is passed through a high pass filter 114 to remove the signal above the threshold frequency f 0 . The filtered signal is then passed to 5 a digitiser 116 and onto a processor 118.

Claims

1. A contact sensitive device comprising a member capable of 5 supporting bending waves, three sensors mounted on the member for measuring bending wave vibration in the member, whereby each sensor determines a measured bending wave signal and a processor which calculates a location of a contact on the member from the measured bending wave signals, characterised 10 in that the processor calculates a phase angle for each measured bending wave signal and a phase difference between the phase angles of least two pairs of sensors so that at least two phase differences are calculated from which the location of the contact is determined. 15

2. A contact sensitive device according to claim 1 comprising an absorber at the edges of the member whereby reflected waves are suppressed.

3. A contact sensitive device according to claim 2, wherein the mechanical impedance of the absorber and the member are 20 selected so as to minimise reflections of bending waves from the edges of the member.

4. A contact sensitive device according to claim 3, wherein the impedances are selected so that bending wave energy is strongly absorbed in a frequency band around a chosen 25 frequency wo. WO 2004/053781 PCT/GB2003/005268 32

5. A contact sensitive device according to claim 4, wherein the impedances are selected to satisfy the following equation: Z r = - i Z (oo) where ZT is the termination impedance of the absorber and 5 ZB is the mechanical impedance of the edge of the member.

6. A contact sensitive device according to claim 4 or claim 5, comprising a band-pass filter for filtering each measured bending wave signal, the filter having a pass-band centred at the chosen frequency ao and a bandwidth of Aw. 10

7. A contact sensitive device according to claim 6, wherein the bandwidth Aco of the filter obeys the relationship: Aco >> 2k(m 0 )hm where vmn is the maximum lateral velocity of the contact.

8. A contact sensitive device according to any one of claims 15 2 to 7, wherein the absorber is made from foamed plastics.

9. A contact sensitive device according to any one of the preceding claims, wherein the member comprises a raised pattern on its surface whereby a contact drawn across the surface provides a force to the member to generate bending 20 waves in the member.

10. A contact sensitive device according to claim 9, wherein the pattern is random whereby a contact travelling over the surface of the member generates a random bending wave signal.

11. A contact sensitive device according to claim 10, wherein 25 the pattern is formed from an anti-reflective coating, an WO2004/053781 PCT/GB2003/005268 33 anti-glare surface finish or an etched finish.

12. A contact sensitive device according to any one of the preceding claims, comprising at least two band-pass filters which have different pass-band frequencies and which 5 simultaneously process the bending wave signals measured by a pair of sensors whereby a phase angle difference for each pass-band frequency is provided by a pair of sensors.

13. A contact sensitive device according to any one of the preceding claims, comprising four sensors on the member. 10

14. A contact sensitive device according to any one of the preceding claims, comprising means for determining the initial location of the contact using the dispersion corrected correlation function of pairs of measured bending wave signals and means for determining subsequent locations of the contact 15 using the phase angle difference between pairs of measured bending wave signals.

15. A contact sensitive device according to any of the preceding claims, wherein the phase angle determining means comprises a phase detector. 20

16. A contact sensitive device according to claim 15, wherein the processor comprises a low-pass filter and a digitiser for determining the phase angles.

17. A contact sensitive device according to any one of the preceding claims, wherein the member is an acoustic radiator 25 and an emitting transducer is mounted to the member to excite bending wave vibration in the member to generate an acoustic WO 2004/053781 PCT/GB2003/005268 34 output.

18. A contact sensitive device according to claim 17, comprising means ensuring that the acoustic output and measured bending wave signals are in discrete frequency bands. 5

19. A contact sensitive device according to any one of the preceding claims, wherein the member is transparent.

20. A method of determining information relating to a contact on a contact sensitive device comprising the steps of providing a member capable of supporting bending waves and 10 three sensors mounted on the member for measuring bending wave vibration in the member, applying a contact to the member at a location, using each sensor to determine a measured bending wave signal and calculating the location of a contact from the measured bending wave signal characterised by calculating a 15 phase angle for each measured bending wave signal, calculating a phase difference between the phase angles of at least two pairs of sensors and determining the location of the contact from the at least two calculated phase differences.

21. A method according to claim 20, comprising suppressing 20 reflected waves by placing an absorber at the edges of the member.

22. A method according to claim 21, comprising selecting the mechanical impedances of the absorber and the member so as to minimise reflections of bending waves from the edges of the 25 member.

23. A method according to claim 22, comprising selecting the WO2004/053781 PCT/GB2003/005268 35 impedances so that bending wave energy is strongly absorbed in a frequency band around a chosen frequency co,.

24. A method according to claim 23, comprising selecting the impedances to satisfy the following equation 5 Z = - i ZB(Co) where ZT is the impedance of the absorber and ZB is the impedance of the edge of the member.

25. A method according to claim 23 or claim 24, comprising filtering each measured bending wave signal by a band-pass 10 filter having a pass-band centred at the chosen frequency co, and a bandwidth of A.

26. A method according to any one of claims 20 to 25, comprising applying the phase difference equation: AOIM = 0 - M = k(co)Aim +2nn7m 15 to determine the location of the contact, where O' is the phase angle of a measured bending wave signal, xi is the distance from the contact location to each sensor, Axim=xl-xm is the path length difference of two sensors, k(o) is the wavevector and nlm is an unknown integer. 20

27. A method according to claim 26, comprising selecting the member to constrain the magnitude of Axam to values less than one half of a wavelength so that njm is determined from 2AO., - 2met < .

28. A method according to claim 26, comprising determining an WO2004/053781 PCT/GB2003/005268 36 initial location of the contact using the dispersion corrected correlation function of a pair of measured bending wave signals and selecting a value of nlm which minimises change in the path length difference. 5

29. A method according to claim 26, comprising selecting a series of values of nim, combining the series of values with each phase angle difference to define a series of path length differences, plotting the series of graphs of the path length differences, and inferring the true value of nim from a point 10 at which a large number of the graphs intersect.

30. A method according to any one of claims 20 to 29, comprising calculating multiple phase angle differences from the pairs of phase angles, plotting a graph of each path length difference and selecting a point at which a large 15 number of the hyperbolae intersect to be the location of the contact.

31. A method according to any one of claims 20 to 30, comprising dividing the measured bending wave signals from each sensor into at least two discrete frequency bands and 20 calculating a phase angle difference for a pair of sensors for each frequency band.