CN106023973B - A kind of working method of electronic percussion instrument - Google Patents
A kind of working method of electronic percussion instrument Download PDFInfo
- Publication number
- CN106023973B CN106023973B CN201610317076.6A CN201610317076A CN106023973B CN 106023973 B CN106023973 B CN 106023973B CN 201610317076 A CN201610317076 A CN 201610317076A CN 106023973 B CN106023973 B CN 106023973B
- Authority
- CN
- China
- Prior art keywords
- signal
- pressure
- value
- pressure sensor
- working method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10H—ELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
- G10H3/00—Instruments in which the tones are generated by electromechanical means
- G10H3/12—Instruments in which the tones are generated by electromechanical means using mechanical resonant generators, e.g. strings or percussive instruments, the tones of which are picked up by electromechanical transducers, the electrical signals being further manipulated or amplified and subsequently converted to sound by a loudspeaker or equivalent instrument
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10H—ELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
- G10H2250/00—Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
- G10H2250/131—Mathematical functions for musical analysis, processing, synthesis or composition
- G10H2250/215—Transforms, i.e. mathematical transforms into domains appropriate for musical signal processing, coding or compression
- G10H2250/235—Fourier transform; Discrete Fourier Transform [DFT]; Fast Fourier Transform [FFT]
Abstract
The present invention relates to instrument devices fields, disclose a kind of working method of electronic percussion instrument.Pass through the working method of the electronic percussion instrument, it can be according to the offset of the X/Y axis on the striking speed of dynamic spatial location and pressure measurement user and impact key, and corresponding signal spectrum figure is obtained by Fourier space transformation, finally under the Mappings of Helmholtz resonator, sound of the band frequency between 20-50Hz is issued.This working method can be applied to the scenes such as composition creation, instrument playing practice, and will play a great role in music formative education and field of music therapy.
Description
Technical field
The present invention relates to instrument devices fields, and in particular, to a kind of working method of electronic percussion instrument.
Background technique
Traditional instrument is more demanding to player, and can only play out a kind of sound of single musical instrument.People need into
The academic program of row many years can be only achieved the level for playing certain musical instrument freely.But common Amateur Music fan is come
It says, musical instrument that is easy to operate and can playing out a variety of music is only the best solution tool for meeting its demand.
Summary of the invention
Aiming at the problem that above-mentioned current traditional instrument equipment, the present invention provides a kind of work sides of electronic percussion instrument
Method, can be according to the offset of the X/Y axis on the striking speed of dynamic spatial location and pressure measurement user and impact key, and passes through
Fourier space transformation obtains corresponding signal spectrum figure, finally under the Mappings of Helmholtz resonator, issues one section
Sound of the frequency between 20-50Hz.This working method can be applied to the scenes such as composition creation, instrument playing practice, and will be
It plays a great role in music formative education and field of music therapy.
The technical solution adopted by the present invention provides a kind of working method of electronic percussion instrument, which is characterized in that including
Following steps: S101. calculates the dynamic spatial location P of impact point using a station method topology, while being surveyed using pressure sensing
Calculate to amount technology cycle the dynamic pressure value F of impact point；S102. to the dynamic spatial location P and the dynamic pressure
Value F carries out Fourier transform, obtains with the periodic signal f (t) of triangle Fourier space formal expansion；S103. to the week
Phase signal f (t) carries out spectrum analysis processing, obtains the first signal spectrum figure of periodic signal f (t)；S104. to first letter
Number spectrogram carries out discrete sampling processing, obtains the second signal spectrogram as made of limited sine-wave superimposed；S105. basis
Helmholtz resonator principle carries out voice mapping processing to the second signal spectrogram, obtains frequency between 20-50Hz's
Low-frequency sine electric signal；S106. the low-frequency sine electric signal is sent to loudspeaker, carries out voice broadcasting.
Specifically, in the step of the dynamic pressure value F for periodically calculating impact point using pressure sensing measuring technique
Include the following steps: that S201. periodically acquires the induced voltage signal that multiple pressure sensors generate in rapid；S202. basis
Following formula calculates pressure suffered by each pressure sensor and value FSumWith centre-of-pressure position (X0, Y0):
In formula, n is the sum of pressure sensor, X0For the abscissa of centre-of-pressure position, Y0For the vertical of centre-of-pressure position
Coordinate, U0nFor the inductive voltage value of n-th of pressure sensor, U'0nFor the operating voltage value of n-th of pressure sensor, matrix K
For the inductive voltage value of each pressure sensor and operating voltage value are mapped as the pressure and value FSumWith the Center of Pressure
Position (X0, Y0) constant coefficient mapping matrix, kfnFor the matrix element for being located at the first row and the n-th column in the constant coefficient mapping matrix
Element, kxnFor the matrix element for being located at the second row and the n-th column in the constant coefficient mapping matrix, kynSquare is mapped for the constant coefficient
It is located at the matrix element of the third line and the n-th column in battle array；S203. by the centre-of-pressure position (X0, Y0) empty with the dynamic of impact point
Between position P corresponded to, and according to the pressure and value FSumWith the centre-of-pressure position (X0, Y0) calculate impact point dynamic
Pressure value F.
Further specifically, the expression formula of the constant coefficient mapping matrix K are as follows:
In formula, n is the sum of pressure sensor, UnFor the reference voltage of n-th of pressure sensor, RnfFor n-th of pressure
The feedback resistance of sensor, kknFor the characteristic curve of n-th of pressure sensor, pnxFor the horizontal seat of n-th of pressure sensor
Mark, pnyFor the ordinate of n-th of pressure sensor.
Specifically, the triangle Fourier expansion form of the periodic signal f (t) is as follows:
In formula, ω0=2 π t are fundamental frequency, m ω0For m subfrequency, AmFor the corresponding Fourier triangle of m subharmonic
Shape fourier coefficient, θmFor m subharmonic first phase.
Specifically, the step S103 includes the following steps: that S301. carries out at adding window frame by frame the periodic signal f (t)
Reason, obtains several signal frames；S302. corresponding FFT spectrum is calculated frame by frame；S303. pseudo-colours is carried out to all FFT spectrums
Mapping, obtains the first signal spectrum figure of the periodic signal f (t).
Further specifically, using the framing windowed function of Matlab system: hamming in the step S301
(N), hanning (N), blackman (N) and bartlett (N) carry out windowing process frame by frame to the periodic signal f (t),
In, N is frame length；In the step S302, use the FFT function of Matlab system: fft (S) and plural number take conjugate function:
Conj (Z) calculates corresponding FFT spectrum, wherein S is the frame signal after adding window, and Z is the calculated result of fft (S) function；In
In the step S303, all FFT spectrums are carried out using the False color mapping function colormap (MAP) of Matlab system
False color mapping, wherein MAP is False color mapping vector.
Further specifically, further including following steps after the step S303: S304. uses the frequency of Matlab system
Spectrogram explicit function: imagese (t, f, C) shows the first signal spectrum figure, wherein t is time coordinate, and f is frequency seat
Mark, C are the colored level value from amplitude spectrum after the mapping of pseudo- color.
Specifically, including the following steps: that S401. obtains the period letter using period point count in the step S104
The cycle T of number f (t)；S402. discrete sampling processing is carried out to the periodic signal f (t) using the cycle T, obtained corresponding
Discrete sampling value sequence f (Δ ti), (i=0,1 ..., N0)；S403. the discrete sampling value sequence f (Δ ti) is carried out
Judgement is read, and is clipped the part less than a signal period from tail portion, is obtained data f ' (the Δ t* containing the integer signal period
I), (i=0,1 ..., N0-1)；S404. simultaneous solution obtains the parameter of each secondary frequency point as follows, obtain by limited sine-wave superimposed and
At second signal spectrogram:
In formula, Cm、A0、AmAnd BmThe fourier series expression of periodic function is corresponded to for the second signal spectrogram
Formula:
In parameter.
Specifically, including the following steps: S501. according to interframe frequency matching criterior to described in the step S105
Each secondary sine wave in binary signal spectrogram carries out peak value matching, obtains matched sine wave parameter；S502. to it is matched just
String wave parameter carries out interpolation processing；S503. by the sine wave parameter of interpolated processing and preset Helmholtz resonator scale
Voice mapping is carried out, low-frequency sine electric signal of the frequency between 20-50Hz is obtained.
To sum up, using the working method of electronic percussion instrument provided by the present invention, have the following beneficial effects: that (1) is adopted
With the dynamic spatial location of station method topology measurement impact point, the need of dynamic spatial location real―time precision measurment can be met
It asks, can operate in insertion type equipment；(2) pressure measurement sensing technology is used, multiple diaphragm pressures is periodically acquired and passes
The information of sensor, can accurate calculation go out the pressure value of tested impact point；(3) it has created a kind of by position and pressure sensing parameter exhibition
It opens as triangular form Fourier space, and signal spectrum, which divides, to be realized to Fourier space processing result using based on Matlab system
Analysis and display, make corresponding spectrogram, belong to multi-field intersection innovation achievement；(4) according to the design of Helmholtz resonator
Principle, will measure resulting sine-wave superimposed and carry out mapping processing can obtain the sine wave signal of one section of specific frequency, and frequency is situated between
Between 20-50Hz, after loudspeaker plays, due to the uniqueness of its frequency dimension, the music of generation can give human psychological and Sheng
Pleasant sensation and euphoria in reason have therapeutic effect medically.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the flow chart of work methods of electronic percussion instrument provided by the invention.
Fig. 2 is that a station method measures mathematical model in the working method of electronic percussion instrument provided by the invention.
Fig. 3 is to realize signal spectral analysis with Matlab system in the working method of electronic percussion instrument provided by the invention
With the signal dynamics spectrogram of a voice of display methods.
Specific embodiment
Hereinafter reference will be made to the drawings, describes the work of electronic percussion instrument provided by the invention in detail by way of example
Method.It should be noted that the descriptions of the manners of these embodiments are used to facilitate the understanding of the present invention, but composition pair
Restriction of the invention.
Various technologies described herein can be used for but be not limited to instrument devices field, can be also used for other similar necks
Domain.
The terms "and/or", only a kind of incidence relation for describing affiliated partner, indicates that there may be three kinds of passes
System, for example, A and/or B, can indicate: individualism A, individualism B exist simultaneously tri- kinds of situations of A and B, the terms
"/and " it is to describe another affiliated partner relationship, indicate may exist two kinds of relationships, for example, A/ and B, can indicate: individually depositing
In A, two kinds of situations of individualism A and B, in addition, character "/" herein, typicallying represent forward-backward correlation object is a kind of "or" pass
System.
Embodiment one
Fig. 1 shows the flow chart of work methods of electronic percussion instrument provided by the invention, and Fig. 2 shows offers of the present invention
Electronic percussion instrument working method in a station method measure mathematical model, Fig. 3 shows electronic percussion provided by the invention
The signal dynamics frequency spectrum of a voice of signal spectral analysis and display methods is realized in the working method of device with Matlab system
Figure.The working method of the electronic percussion instrument provided in this embodiment includes the following steps.
S101. the dynamic spatial location P of impact point is calculated using a station method topology, while being measured using pressure sensing
Calculate to technology cycle the dynamic pressure value F of impact point.
In step s101, one station method topology of the use calculates the method for the dynamic spatial location P of impact point such as
Under: as shown in Fig. 2, setting P point as the dynamic spatial location (coordinate points of i.e. tested impact key) of impact point, pass through angle sensor
It is α with the vertical direction angle for measuring coordinate system central point O that device measurement, which obtains point P, and horizontal direction angle is β, passes through distance and passes
It is d that sensor measurement, which obtains O point to the distance of P point, so as to which the spatial position coordinate of point P is calculated according to system model
Value.
In step s101, specifically, periodically calculating the dynamic of impact point using pressure sensing measuring technique described
Include the following steps: that S201. periodically acquires the induced voltage letter that multiple pressure sensors generate in the step of state pressure value F
Number；S202. pressure suffered by each pressure sensor and value F are calculated according to the following formulaSumWith centre-of-pressure position (X0, Y0):
In formula, n is the sum of pressure sensor, X0For the abscissa of centre-of-pressure position, Y0For the vertical of centre-of-pressure position
Coordinate, U0nFor the inductive voltage value of n-th of pressure sensor, U'0nFor the operating voltage value of n-th of pressure sensor, matrix K
For the inductive voltage value of each pressure sensor and operating voltage value are mapped as the pressure and value FSumWith the Center of Pressure
Position (X0, Y0) constant coefficient mapping matrix, kfnFor the matrix element for being located at the first row and the n-th column in the constant coefficient mapping matrix
Element, kxnFor the matrix element for being located at the second row and the n-th column in the constant coefficient mapping matrix, kynSquare is mapped for the constant coefficient
It is located at the matrix element of the third line and the n-th column in battle array；S203. by the centre-of-pressure position (X0, Y0) empty with the dynamic of impact point
Between position P corresponded to, and according to the pressure and value FSumWith the centre-of-pressure position (X0, Y0) calculate impact point dynamic
Pressure value F.
The pressure sensor can be, but not limited to as thin-film pressure sensor, by taking hand dish as an example, in the strike of hand dish
Position is disposed with 8 diaphragm pressure sensors, and each thin-film pressure sensor has high sensitivity, can pat in user
When, by generate induced voltage signal record in real time user striking speed, impact position on X/Y axis faint offset and
Pressure Pa on the position of power.The expression formula of the constant coefficient mapping matrix K are as follows:
In formula, n is the sum of pressure sensor, UnFor the reference voltage of n-th of pressure sensor, RnfFor n-th of pressure
The feedback resistance of sensor, kknFor the characteristic curve of n-th of pressure sensor, pnxFor the horizontal seat of n-th of pressure sensor
Mark, pnyFor the ordinate of n-th of pressure sensor.In step s 201, it can be, but not limited to acquire 8 using 10ms as the period
The induced voltage signal of the generation of diaphragm pressure sensor,
Then according to step S202 and step S203 can accurate calculation go out the dynamic pressure value of impact point.
S102. Fourier transform is carried out to the dynamic spatial location P and dynamic pressure value F, obtained with triangle
The periodic signal f (t) of Fourier space formal expansion.
In step s 102, the triangle Fourier expansion form of the periodic signal f (t) is as follows:
In formula, ω0=2 π t are fundamental frequency, m ω0For m subfrequency, AmFor the corresponding Fourier triangle of m subharmonic
Shape fourier coefficient, θmFor m subharmonic first phase.
S103. spectrum analysis processing is carried out to the periodic signal f (t), obtains the first signal frequency of periodic signal f (t)
Spectrogram.
In step S104, specifically, the step S103 includes the following steps: S301. to the periodic signal f (t)
Windowing process frame by frame is carried out, several signal frames are obtained；S302. corresponding FFT (Fast Fourier is calculated frame by frame
Transformation, fast Fourier transform) frequency spectrum；S303. False color mapping is carried out to all FFT spectrums, obtained described
The first signal spectrum figure of periodic signal f (t).Further specifically, can be, but not limited to based on Matlab system (MATrix
LABoratory, matrix labotstory are a business mathematics software produced by The MathWorks company, the U.S., for calculating
The advanced techniques computational language and interactive environment that method exploitation, data visualization, data analysis and numerical value calculate；In addition to matrix
Outside the common functions such as operation, drafting function/data image, MATLAB can also be used to create user interface and with call other languages
The program that speech (including C, C++ and FORTRAN) is write) signal spectral analysis is realized to Fourier space processing, it makes accurately
First signal spectrum figure: in the step S301, using the framing windowed function of Matlab system: hamming (N),
Hanning (N), blackman (N) and bartlett (N) carry out windowing process frame by frame to the periodic signal f (t), wherein N
For frame length；In the step S302, use the FFT function of Matlab system: fft (S) and plural number take conjugate function: conj
(Z) FFT spectrum of the periodic signal f (t) is calculated, wherein S is the frame signal after adding window, and Z is the calculating of fft (S) function
As a result；In the step S303, using Matlab system False color mapping function colormap (MAP) to FFT spectrum into
Row False color mapping, wherein MAP is False color mapping vector.MAP is one 64 × 3 matrix, first, second, and third
The value of column respectively indicates red, green and blue saturation degrees.For example, MAP=[0 0 0] indicates the mapping of a black；MAP
=[1 1 1] indicate a pure white mapping；MAP=[1 0 0] indicates a pure red mapping.
Further specifically, further including following steps after the step S303: S304. uses the frequency of Matlab system
Spectrogram explicit function: imagese (t, f, C) shows the first signal spectrum figure, wherein t is time coordinate, and f is frequency seat
Mark, C are the colored level value from amplitude spectrum after the mapping of pseudo- color.C is the colour from amplitude spectrum after False color mapping
Level value.If the framing number of whole section of voice is M, then t is the row vector of M dimension, its value is corresponding to every frame starting point
Time serial number, the effective sampling point of the frequency domain of frequency spectrum are N/2, so f is the row vector of a N/2, each element corresponds to frequency spectrum width
Angle value.Correspondingly, C is the matrix of a M × (N/2) dimension.
As shown in Figure 3 one realizes the signal of a voice of signal spectral analysis and display methods with Matlab system
Dynamic spectrum figure.When calling spgram () function, partial parameters used are set as: Winsiz=256, Shift=64, Base
=0, Coltype=1 (MAP=hot).Descending be mapped to of spectrum value is shown by becoming clear to dark pseudo-colours in this figure
Show, it is possible thereby to be clearly observed the resonance peak structure of voiced segments.
S104. discrete sampling processing is carried out to the first signal spectrum figure, obtains being formed by limited sine-wave superimposed
Second signal spectrogram.
In the step S104, specifically, including the following steps: that S401. obtains the period using period point count
The cycle T of signal f (t)；S402. discrete sampling processing is carried out to the periodic signal f (t) using the cycle T, obtains and corresponds to
Discrete sampling value sequence f (Δ ti), (i=0,1 ..., N0)；S403. to the discrete sampling value sequence f (Δ ti) into
Row judgement is read, and is clipped the part less than a signal period from tail portion, is obtained the data f ' (Δ containing the integer signal period
T*i), (i=0,1 ..., N0-1)；S404. simultaneous solution obtains the parameter of each secondary frequency point as follows, obtains by limited sine-wave superimposed
Made of second signal spectrogram:
In formula, Cm、A0、AmAnd BmThe fourier series expression of periodic function is corresponded to for the second signal spectrogram
Formula:
In parameter.
S105. voice mapping processing is carried out to the second signal spectrogram according to Helmholtz resonator principle, obtained
Low-frequency sine electric signal of the frequency between 20-50Hz.
In the step S105, specifically, including the following steps: S501. according to interframe frequency matching criterior to described
Each secondary sine wave in second signal spectrogram carries out peak value matching, obtains matched sine wave parameter；S502. to matched
Sine wave parameter carries out interpolation processing；S503. by the sine wave parameter of interpolated processing and preset Helmholtz resonator sound
Rank carries out voice mapping, obtains low-frequency sine electric signal of the frequency between 20-50Hz.
In step S501, by carrying out the matching, the sine wave freuqency that may be matched, amplitude and phase ginseng
Number, but due to the time variation of these parameters, will directly be will lead to these parameter synthesis voices discontinuous between frame and frame
Property, the quality of synthetic speech signal can be reduced.Therefore, it is necessary to take measures and make smoothed between frame between frames
It crosses, a kind of common method is exactly to carry out interpolation to sine wave signal parameter, and specific method can be, but not limited to are as follows: to
Matched phase and amplitude carries out cube interpolation and an interpolation, the speech signal parameter after then acquiring interpolation respectively.In addition,
Still illustrate the parameter preset of the Helmholtz resonator scale by taking hand dish as an example: being equipped with a basis in the bottom centre of hand dish
Hole " Gu " of Helmholtz resonator principle design, is then equipped with 7 lines, around central point basis line " Ding " at an upper portion thereof
It is arranged from down to height with zigzag, and each line and the Helmholtz resonator of main body is enable to empathize, concurrently set basis
Sound is D2, and center " Ding " is D3, and five ranks are double bass A3 in sound ring, the D4 and A4 spent in scales including its eight down,
" Gu " is D5.The sine wave parameter of interpolated processing and preset Helmholtz resonator scale are finally subjected to voice mapping,
Obtain low-frequency sine electric signal of the frequency between 20-50Hz.
S106. the low-frequency sine electric signal is sent to loudspeaker, carries out voice broadcasting.
According to above step, a kind of work of electronic percussion instrument based on position and pressure sensing technology can be finally realized
Method can be applied to the scenes such as composition creation, instrument playing practice.
The working method of the electronic percussion instrument provided in this embodiment has the following technical effect that (1) using a station
Method topology measures the dynamic spatial location of impact point, can meet the needs of dynamic spatial location real―time precision measurment, can transport
For in insertion type equipment；(2) pressure measurement sensing technology is used, the letter of multiple diaphragm pressure sensors is periodically acquired
Breath, can accurate calculation go out the pressure value of tested impact point；(3) having created a kind of by position and pressure sensing parametric evolving is triangle
Type Fourier space, and to the analysis of Fourier space processing result realization signal spectrum and shown using based on Matlab system
Show, make corresponding spectrogram, belongs to multi-field intersection innovation achievement；It (4), will according to the design principle of Helmholtz resonator
The sine wave signal of one section of specific frequency can be obtained by measuring resulting sine-wave superimposed and carrying out mapping processing, and frequency is between 20-
Between 50Hz, after loudspeaker plays, due to the uniqueness of its frequency dimension, the music of generation can to human psychological and physiologically
Pleasant sensation and euphoria have therapeutic effect medically.
As described above, the present invention can be realized preferably.For a person skilled in the art, religion according to the present invention
It leads, the working method for designing various forms of electronic percussion instruments does not need creative labor.The present invention is not being departed from
Principle and spirit in the case where these embodiments are changed, modify, replace, integrating and modification still falls within guarantor of the invention
It protects in range.
Claims (9)
1. a kind of working method of electronic percussion instrument, which comprises the steps of:
S101. the dynamic spatial location P of impact point is calculated using a station method topology, while using pressure sensing measuring technique
Periodically calculate the dynamic pressure value F of impact point, wherein the dynamic for calculating impact point using a station method topology is empty
Between position P method it is as follows: set P point as the dynamic spatial location of impact point, point P obtained by angular transducer measurement and is surveyed
The vertical direction angle for measuring coordinate system central point O is α and horizontal direction angle is β, obtains O point extremely by range sensor measurement
The distance of P point is d, and the spatial position coordinate value of point P is finally calculated according to system model；
S102. Fourier transform is carried out to the dynamic spatial location P and dynamic pressure value F, obtained in triangle Fu
The periodic signal f (t) of leaf progression form expansion；
S103. spectrum analysis processing is carried out to the periodic signal f (t), obtains the first signal spectrum figure of periodic signal f (t)；
S104. discrete sampling processing is carried out to the first signal spectrum figure, obtains as made of limited sine-wave superimposed the
Binary signal spectrogram；
S105. voice mapping processing is carried out to the second signal spectrogram according to Helmholtz resonator principle, obtains frequency
Between the low-frequency sine electric signal of 20-50Hz；
S106. the low-frequency sine electric signal is sent to loudspeaker, carries out voice broadcasting.
2. a kind of working method of electronic percussion instrument as described in claim 1, which is characterized in that passed described using pressure
Include the following steps: in the step of calculating to sensed quantity technology cycle the dynamic pressure value F of impact point
S201. the induced voltage signal that multiple pressure sensors generate periodically is acquired；
S202. pressure suffered by each pressure sensor and value F are calculated according to the following formulaSumWith centre-of-pressure position (X0, Y0):
In formula, n is the sum of pressure sensor, X0For the abscissa of centre-of-pressure position, Y0For the vertical seat of centre-of-pressure position
Mark, U0nFor the inductive voltage value of n-th of pressure sensor, U'0nFor the operating voltage value of n-th of pressure sensor, matrix K is
The inductive voltage value of each pressure sensor and operating voltage value are mapped as the pressure and value FSumWith the Center of Pressure position
Set (X0, Y0) constant coefficient mapping matrix, kfnFor the matrix element for being located at the first row and the n-th column in the constant coefficient mapping matrix
Element, kxnFor the matrix element for being located at the second row and the n-th column in the constant coefficient mapping matrix, kynSquare is mapped for the constant coefficient
It is located at the matrix element of the third line and the n-th column in battle array；
The expression formula of the constant coefficient mapping matrix K are as follows:
In formula, n is the sum of pressure sensor, UnFor the reference voltage of n-th of pressure sensor, RnfFor n-th of pressure sensing
The feedback resistance of device, kknFor the characteristic curve of n-th of pressure sensor, pnxFor the abscissa of n-th of pressure sensor,
pnyFor the ordinate of n-th of pressure sensor；
S203. by the centre-of-pressure position (X0, Y0) carried out with the dynamic spatial location P of impact point it is corresponding, and according to the pressure
Power and value FSumWith the centre-of-pressure position (X0, Y0) calculate impact point dynamic pressure value F.
3. a kind of working method of electronic percussion instrument as described in claim 1, which is characterized in that the periodic signal f (t)
Triangle Fourier expansion form it is as follows:
In formula, ω0=2 π t are fundamental frequency, m ω0For m subfrequency, AmFor corresponding Fourier triangle Fu of m subharmonic
Vertical leaf system number, θmFor m subharmonic first phase.
4. a kind of working method of electronic percussion instrument as described in claim 1, which is characterized in that the step S103 includes
Following steps:
S301. windowing process frame by frame is carried out to the periodic signal f (t), obtains several signal frames；
S302. corresponding FFT spectrum is calculated frame by frame；
S303. False color mapping is carried out to all FFT spectrums, obtains the first signal spectrum figure of the periodic signal f (t).
5. a kind of working method of electronic percussion instrument as claimed in claim 4, which is characterized in that include the following:
In the step S301, the framing windowed function of Matlab system: hamming (N) is used, hanning (N),
Blackman (N) and bartlett (N) carries out windowing process frame by frame to the periodic signal f (t), wherein N is frame length；
In the step S302, use the FFT function of Matlab system: fft (S) and plural number take conjugate function: conj (Z) meter
Calculate corresponding FFT spectrum, wherein S is the frame signal after adding window, and Z is the calculated result of fft (S) function；
In the step S303, using the False color mapping function colormap (MAP) of Matlab system to all FFT frequencies
Spectrum carries out False color mapping, wherein MAP is False color mapping vector.
6. a kind of working method of electronic percussion instrument as claimed in claim 5, which is characterized in that the step S303 it
After further include following steps:
S304. use the spectrogram explicit function of Matlab system: imagese (t, f, C) shows the first signal spectrum figure,
Wherein, t is time coordinate, and f is frequency coordinate, and C is the colored level value from amplitude spectrum after the mapping of pseudo- color.
7. a kind of working method of electronic percussion instrument as described in claim 1, which is characterized in that in the step S104
Include the following steps:
S401. the cycle T of the periodic signal f (t) is obtained using period point count；
S402. discrete sampling processing is carried out to the periodic signal f (t) using the cycle T, obtains corresponding dis-crete sample values
Sequence f (Δ ti), (i=0,1 ..., N0)；
S403. judgement reading is carried out to the discrete sampling value sequence f (Δ ti), clips less than one signal period from tail portion
Part, obtain the data f ' (Δ t*i) containing the integer signal period, (i=0,1 ..., N0-1)；
S404. simultaneous solution obtains the parameter of each secondary frequency point as follows, obtains the second signal frequency spectrum as made of limited sine-wave superimposed
Figure:
In formula, Cm、A0、AmAnd BmThe fourier series expression formula of periodic function is corresponded to for the second signal spectrogram:
In parameter, m is harmonic wave
Number.
8. a kind of working method of electronic percussion instrument as described in claim 1, which is characterized in that in the step S105
Include the following steps:
S501. peak value matching is carried out to each secondary sine wave in the second signal spectrogram according to interframe frequency matching criterior,
Obtain matched sine wave parameter；
S502. interpolation processing is carried out to matched sine wave parameter；
S503. the sine wave parameter of interpolated processing and preset Helmholtz resonator scale are subjected to voice mapping, obtained
Low-frequency sine electric signal of the frequency between 20-50Hz.
9. a kind of working method of electronic percussion instrument as claimed in claim 2, which is characterized in that the pressure sensor is
Thin-film pressure sensor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610317076.6A CN106023973B (en) | 2016-05-12 | 2016-05-12 | A kind of working method of electronic percussion instrument |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610317076.6A CN106023973B (en) | 2016-05-12 | 2016-05-12 | A kind of working method of electronic percussion instrument |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106023973A CN106023973A (en) | 2016-10-12 |
| CN106023973B true CN106023973B (en) | 2019-10-29 |
Family
ID=57099822
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610317076.6A Expired - Fee Related CN106023973B (en) | 2016-05-12 | 2016-05-12 | A kind of working method of electronic percussion instrument |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106023973B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109360545A (en) * | 2018-11-14 | 2019-02-19 | 重庆市星贯众文化艺术传播有限公司 | Electronic hand dish |
| TWI741675B (en) * | 2020-07-13 | 2021-10-01 | 高頓科技有限公司 | Human-machine sensor input component and human-machine sensor input system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6111181A (en) * | 1997-05-05 | 2000-08-29 | Texas Instruments Incorporated | Synthesis of percussion musical instrument sounds |
| CN1420261A (en) * | 2001-09-19 | 2003-05-28 | 通用电气公司 | System and method for suppressing pressure wave using correction signal |
| JP2007248551A (en) * | 2006-03-14 | 2007-09-27 | Yamaha Corp | Waveform data producing method, waveform data producing device, program, and waveform memory producing method |
| CN101465121A (en) * | 2009-01-14 | 2009-06-24 | 苏州瀚瑞微电子有限公司 | Method for implementing touch virtual electronic organ |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8178769B2 (en) * | 2010-07-08 | 2012-05-15 | Mark David Steele | Universal drum pedal instrument mounting stand |
-
2016
- 2016-05-12 CN CN201610317076.6A patent/CN106023973B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6111181A (en) * | 1997-05-05 | 2000-08-29 | Texas Instruments Incorporated | Synthesis of percussion musical instrument sounds |
| CN1420261A (en) * | 2001-09-19 | 2003-05-28 | 通用电气公司 | System and method for suppressing pressure wave using correction signal |
| JP2007248551A (en) * | 2006-03-14 | 2007-09-27 | Yamaha Corp | Waveform data producing method, waveform data producing device, program, and waveform memory producing method |
| CN101465121A (en) * | 2009-01-14 | 2009-06-24 | 苏州瀚瑞微电子有限公司 | Method for implementing touch virtual electronic organ |
Non-Patent Citations (1)
| Title |
|---|
| 利用傅里叶级数分析决定乐音音色的因素;赵纪青 等;《湖北广播电视大学学报》;20081130;第28卷(第11期);第074页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106023973A (en) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105976800A (en) | Electronic device, method for recognizing playing of string instrument in electronic device | |
| CN101853650B (en) | Music processing and output system and method thereof | |
| CN102568453B (en) | Performance apparatus and electronic musical instrument | |
| CN105741639B (en) | A kind of micro- sense palm musical instrument for simulating bowstring kind musical instrument | |
| US11437006B2 (en) | Systems and methods for music simulation via motion sensing | |
| CN106023973B (en) | A kind of working method of electronic percussion instrument | |
| Ouni et al. | VisArtico: a visualization tool for articulatory data | |
| JP2013504096A (en) | A watch with an integrated semitone gradation | |
| CN102568455A (en) | Performance apparatus and electronic musical instrument | |
| CN105810185A (en) | Multifunctional digital musical instrument | |
| CN107566969A (en) | A kind of enclosed environment internal low-frequency Reconstruction of Sound Field method | |
| EP2133865A2 (en) | Sound synthesizer | |
| GB2559135A (en) | Transducer apparatus for an edge-blown aerophone and an edge-blown aerophone having the transducer apparatus | |
| CN105807907B (en) | Body-sensing symphony performance system and method | |
| Ferguson et al. | Seeing sound: Real-time sound visualisation in visual feedback loops used for training musicians | |
| US10341795B2 (en) | Log complex color for visual pattern recognition of total sound | |
| CN205486954U (en) | Little sense palm musical instrument of simulation bowstring class musical instrument | |
| CN107705774A (en) | A kind of electronic musical instrument analog machine | |
| Mickiewicz | Particle image velocimetry and proper orthogonal decomposition applied to aerodynamic sound source region visualization in organ flue pipe | |
| WO2017150963A1 (en) | A color coded music notation system | |
| Barrass | Sculpting a sound space with information properties | |
| KR101754983B1 (en) | Digital Multifunctional Musical Instrument | |
| McLachlan | A spatial theory of rhythmic resolution | |
| Park et al. | Discrete isomorphic completeness and a unified isomorphic layout format | |
| Huovinen et al. | Bird's‐Eye Views of the Musical Surface: Methods for Systematic Pitch‐Class Set Analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191029 Termination date: 20200512 |