CN105468860A - Method for implementing digital signal processing algorithms based on Praat script - Google Patents

Method for implementing digital signal processing algorithms based on Praat script Download PDF

Info

Publication number
CN105468860A
CN105468860A CN201510870167.8A CN201510870167A CN105468860A CN 105468860 A CN105468860 A CN 105468860A CN 201510870167 A CN201510870167 A CN 201510870167A CN 105468860 A CN105468860 A CN 105468860A
Authority
CN
China
Prior art keywords
praat
script
digital signal
signal processing
program
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.)
Pending
Application number
CN201510870167.8A
Other languages
Chinese (zh)
Inventor
库尔班吾布力
贾振红
汪烈军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xinjiang University
Original Assignee
Xinjiang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xinjiang University filed Critical Xinjiang University
Priority to CN201510870167.8A priority Critical patent/CN105468860A/en
Publication of CN105468860A publication Critical patent/CN105468860A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2117/00Details relating to the type or aim of the circuit design
    • G06F2117/08HW-SW co-design, e.g. HW-SW partitioning

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

The invention relates to the technical field of signal and information processing and particularly relates to a method for implementing digital signal processing algorithms based on a Praat script. The method utilizes a voice processing software Praat environment to implement various digital signal processing algorithms by writing a Praat script program. All algorithm implementation (simulation) results can be observed in modes of figures, images and sound; a simulation implementation (program running) speed is high; a running result is visual and vivid; an audio-visual effect can be achieved by combining algorithm simulation; and the method is a high-efficiency implementing method for various digital signal processing algorithms.

Description

A kind of implementation method of the digital signal processing algorithm based on Praat script
Technical field
The present invention relates to Signal and Information Processing technical field, particularly relate to a kind of implementation method of the digital signal processing algorithm based on Praat script.
Background technology
Digital signal processing (DSP) is one of field very important in Science and engineering, and it develops rapidly in 50 years in the past.Now, DSP is commonly used.Because it is in communication, space exploration, consumer electronics product, robot, medicine, instrument and meter, military project, widespread use in many key areas such as automobile and seismology.DSP has made user remove noise signal, accelerates communication speed, and stores more data, and provide many advantages in its simulation process.
In DSP theoretical system, there is the expression formula of many complexity, formula, algorithm, waveform, frequency spectrum etc.Within some time, understand these is very difficult.Up to the present, there is the method using traditional programming language to realize in DSP algorithm realization method, as: Fortran, Pascal, C, C++, Mathematical, Java etc.Also the method using graphical programming environment is had, e.g., Matlab, Simulink, LabVIEW, SystemView, and MacromediaFlash etc.In addition, also have the implementation method using software and hardware to combine, e.g., CodeComposerStudio (CCS), Filed-ProgrammableGatearray (FPGA) etc.
In these multiple implementation methods, use the simulated effect of traditional programming language implementation method, imagery expression effect undesirable; Use graphical programming environment method can only with image simulation result; And the difficulty of the implementation method combined with software and hardware is large, simulation result is single.
Summary of the invention
The object of the invention is to solve the problem, provide the many algorithms that a kind of method by Praat Script Programming realizes digital signal processing, their simulated effects can show with figure, image and sound.
The present invention is achieved through the following technical solutions:
Based on an implementation method for the digital signal processing algorithm of Praat script, comprise the following steps:
1) first, start Praat software, enter Praat environment, its major interfaces " Praat target " is opened;
2) under " control " menu of Praat main window, " new praat script " option is selected, the perform region of writing praat shell script can be entered, or under this menu, select " opening praat script ", the shell script preserved can be opened and edit;
3) 2) described in perform region write the algorithm of digital signal processing according to praat script syntax rule realize program, and preserve the program of program or the editor write by the hold-over command in the File menu of " script window ";
4) 2) described in perform region " RUN " menu in the " RUN " order selected can Run Script program, obtain the Realization of Simulation result of digital signal processing algorithm.
Further, described major interfaces is main window, and described perform region is window.
Compared with prior art, the invention has the beneficial effects as follows: this method utilizes speech processing software Praat environment, realized the many algorithms of digital signal processing by the method for writing Praat script journey.Algorithm realization (emulation) result all can be observed by the mode of figure, image and sound, the Realization of Simulation (program operation) speed is fast, operation result visual pattern, can the vision of combination algorithm the Realization of Simulation and auditory effect, be a kind of efficient implementation method of digital signal processing many algorithms.
Accompanying drawing explanation
Fig. 1 is structured flowchart of the present invention;
Fig. 2 is the surface chart starting Praat and open its script window;
Fig. 3 is for writing Praat shell script surface chart;
Fig. 4 is for opening voice document surface chart;
Fig. 5 is the figure of the new audio files that algorithm realization produces;
Fig. 6 is an original signal waveforms;
Fig. 7 is another original signal waveforms;
Fig. 8 is the oscillogram of above two (Fig. 6 and Fig. 7) signal plus;
Fig. 9 is the figure of the drawing Praat shell script of complex tone;
Figure 10 is the oscillogram of complex tone;
Figure 11 is power spectrum density and the autocorrelation function graph of white noise;
Figure 12 is the window of the shell script that white noise produces;
Figure 13 be white noise and spectrogram;
Figure 14 is the part shell script figure providing signal addition algorithm;
Figure 15 is the result figure of two original signals and their addition using " in real time " pattern to carry out;
Figure 16 is the partial code figure of the program that realizes of the logical and high-pass filtering of the low pass of signal, band;
Figure 17 is the later Praat target window figure of three kinds of filtering;
Figure 18 is original signal oscillogram;
Figure 19 is the later signal waveforms of low-pass filtering;
Figure 20 is the later signal waveforms of bandpass filtering;
Figure 21 is the later signal waveforms of high-pass filtering.
In figure: 1 starts Praat software, 2 select " new praat script " option to write praat shell script under " control " menu, 3 select " opening praat script " under " control " menu, edit the shell script preserved, 4 preserve the program of program or the editor write by the hold-over command in the File menu of " script window ", 5 " reading " menus clicking Praat main window, and select " reading from file " to read literary composition and filesselected, the " RUN " order of the " RUN " menu of 6 script windows can Run Script program Run Script program.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and Examples, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
Embodiment 1
Refer to Fig. 1-Fig. 8, the present invention is the many algorithms utilizing Praat script to realize digital signal processing.Provide the associated picture of each link in the operating process of Praat shell script below.After starting Praat software, its major interfaces (main window) " Praat target " is opened, and this window has " control ", " newly ", and four menus such as " reading " and " writing ", as shown in Fig. 2 (a).Then, under " control " menu of Praat main window select " new praat script " option, can enter write praat shell script perform region (window) as shown in Fig. 2 (b).In above Fig. 2 (b), its script window of Praat has " file ", " editor ", four menus such as " searching " and " RUN ".Below menu bar is exactly the perform region of compile script program.Realizing program and preserving of the algorithm of digital signal processing is write according to praat script syntax rule in this region.During save routine, the hold-over command in the File menu of " script window " is utilized to preserve the program of program or the editor write.As shown in Figure 3.In above Fig. 3, only provide a part for this program, the row that in program, symbol " # " starts is the comment section of program.After writing program, click " reading " menu of Praat main window, and select " reading from file " just can file reading, and choose audio files by mouse-click, as shown in Figure 4.In above Fig. 4, in the perform region on the left side, show the file be opened.After mouse clicks these files, background document becomes blueness, illustrates selected.Finally, can Run Script program in the " RUN " order of the " RUN " menu utilizing script window, obtain the Realization of Simulation result of digital signal processing algorithm, namely, new audio files is produced, as shown in Figure 5, when general, voice document length is chosen for and is no more than 10 seconds.In Figure 5, after select File, the Edit button on the right of clicking is observed and is realized result.Two signals originally and their be added after signal as shown in Fig. 6-Fig. 8.In fig. 8, the first half is the waveform being added later signal time domain, waveform can see signal duration, the parameters such as amplitude.Lower Half is the waveform being added later signal spectrum (after Fourier transform), can obtain signal frequency composition, fundamental frequency, resonance peak, the parameters such as energy from this waveform.In addition, not only intuitively arrive the waveform after process, the effect of process can be heard further, the understanding to this algorithm can be deepened.Therefore, by the Realization of Simulation visual result that the present invention's (realizing digital signal processing algorithm with Praat shell script) obtains, image, not only obtains the signal waveform after algorithm process, can compare with original signal, but also the effect of the signal after process can be heard.This method, by conjunction with vision, auditory effect, reaches the effect of comprehensive simulating.
Embodiment 2
The complex tone plot analysis of Praat shell script
The plot of the here matrix that has been this code establishing, but there is no label that is inner and outside, and be the sine wave that 4 row 3 arrange, nx=3, ny=4, i=0tonx, j=0tony its to generate formula be 1/2*sin (2*pi* (' nx'* (' j'-1)+' i') * x), its cycle is (1, 1)=1, (1, 2)=2, (1, 3)=3, (2, 1)=4, (2, 2)=5, (2, 3)=6, (3, 1)=7, (3, 2)=8, (3, 3)=9, (4, 1)=10, (4, 2)=11, (4, 3)=12, also can represent by matrix, duration is all 0.5 second.It realizes shell script as shown in Figure 9, and program operation result as shown in Figure 10, as can be seen from Figure 10, realizes the multiple sinusoidal signal that (drawing) meets above condition simultaneously.
Embodiment 3
The generation of white noise
White noise is a kind of random signal in common communication system with smooth power spectrum density.DSP process introduces the white noise of several types and continuous time, the white noise signal of infinite bandwidth.The power spectrum density of white noise is defined as follows:
P n ( &omega; ) = n 0 2 , ( - &infin; < &omega; < + &infin; ) - - - ( 1 )
Wherein, n0 is a constant.Its autocorrelation function is expressed as:
R n ( &tau; ) = 1 2 &pi; &Integral; - &infin; + &infin; ( n 0 2 ) e j &omega; &tau; d &omega; = n 0 2 &delta; ( &tau; ) - - - ( 2 ) ;
The power spectrum density of white noise and the figure of autocorrelation function are as represented in Figure 11.
If Figure 12 is to the window of shell script showing white noise, this order code comprises annotation (wherein with # beginning), space, variable, with circulation (for circulation) etc., the duration of this signal arranges 1 second with " positive duration (secondary) 1.0 ".The operation result of this program as shown in figure 13.
Embodiment 4
Signal addition algorithm
Signal addition is a rudimentary algorithm in signal transacting.The addition of two signals is:
z(t)=x(t)+y(t)(3)
Wherein, x (t) and y (t) are two signals, z (t) be they and.Under special circumstances, one of them signal is constant (or X (t)=C, wherein C is constant).Write the additive operation that a Praat shell script realizes two signals.As provided the part shell script of signal addition algorithm in Figure 14.
It should be noted that select addition scheme to be very important before working procedure.Because it has two kinds of patterns, one of them is " pointwise value " pattern, and it is suitable for the addition of discrete signal, and another kind is " in real time " pattern, for the addition in different time domains and sampling rate signal.Figure 15 is the result be added with their with two original signals that " in real time " pattern is carried out.
In fig. 12, A, B are respectively two original signals, and C is that they are added later waveform.
Embodiment 5
Multiple filtering
Utilize the method not only can realize all kinds of (low pass, band lead to, are with resistance and high pass) filtering of digital signal, and the multiple filtering of signal can be realized by a shell script.Figure 16 is a part for the program that realizes of the logical and high-pass filtering of the low pass of signal, band.
After this program is run, the waveform of the signal after Praat target window produces these three kinds of filtering respectively.As shown in figure 17.
In Figure 14 window, topmost (selected) signal is former (before filtering) signal, and three signals below it are respectively the later signal waveform of low-pass filtering, bandpass filtering and high-pass filtering, these signals utilize the Edit button on the right to find out.Signal and the later signal waveform of high-pass filtering after the later signal of original signal, low-pass filtering, bandpass filtering, as shown in Figure 18-Figure 21.
From Figure 18-Figure 21, not only can find out original signal and the later signal waveform of three kinds of filtering, the relation of frequency spectrum and difference, but also their filter effect can be heard, realize this three kinds of filtering algorithms visually.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.

Claims (2)

1. based on an implementation method for the digital signal processing algorithm of Praat script, it is characterized in that, comprise the following steps:
1) first, start Praat software, enter Praat environment, its major interfaces " Praat target " is opened;
2) under " control " menu of Praat main window, " new praat script " option is selected, the perform region of writing praat shell script can be entered, or under this menu, select " opening praat script ", the shell script preserved can be opened and edit;
3) 2) described in perform region write digital signal processing algorithm according to praat script syntax rule realize program, and preserve the program of program or the editor write by the hold-over command in the File menu of " script window ";
4) 2) described in perform region " RUN " menu in the " RUN " order selected can Run Script program, obtain the Realization of Simulation result of digital signal processing algorithm.
2. the implementation method of a kind of digital signal processing algorithm based on Praat script according to claim 1, it is characterized in that: described major interfaces is main window, described perform region is window.
CN201510870167.8A 2015-12-01 2015-12-01 Method for implementing digital signal processing algorithms based on Praat script Pending CN105468860A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510870167.8A CN105468860A (en) 2015-12-01 2015-12-01 Method for implementing digital signal processing algorithms based on Praat script

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510870167.8A CN105468860A (en) 2015-12-01 2015-12-01 Method for implementing digital signal processing algorithms based on Praat script

Publications (1)

Publication Number Publication Date
CN105468860A true CN105468860A (en) 2016-04-06

Family

ID=55606555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510870167.8A Pending CN105468860A (en) 2015-12-01 2015-12-01 Method for implementing digital signal processing algorithms based on Praat script

Country Status (1)

Country Link
CN (1) CN105468860A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7599821B1 (en) * 2003-04-30 2009-10-06 Cadence Design Systems, Inc. Mixed language simulation cycle for analog mixed-signal circuit design and simulation
CN101727902A (en) * 2008-10-29 2010-06-09 中国科学院自动化研究所 Method for estimating tone
CN104795077A (en) * 2015-03-17 2015-07-22 北京航空航天大学 Voice annotation quality consistency detection method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7599821B1 (en) * 2003-04-30 2009-10-06 Cadence Design Systems, Inc. Mixed language simulation cycle for analog mixed-signal circuit design and simulation
CN101727902A (en) * 2008-10-29 2010-06-09 中国科学院自动化研究所 Method for estimating tone
CN104795077A (en) * 2015-03-17 2015-07-22 北京航空航天大学 Voice annotation quality consistency detection method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
熊子瑜: ""Praat语音分析软件操作手册"", 《新浪博客》 *

Similar Documents

Publication Publication Date Title
Feichtinger Advances in Gabor analysis
Callegari et al. Spectral properties of chaos-based FM signals: Theory and simulation results
CN101483797B (en) Head-related transfer function generation method and apparatus for earphone acoustic system
Chen et al. Warped variational mode decomposition with application to vibration signals of varying-speed rotating machineries
JP2021167977A (en) Voice signal processing method, voice signal processing device, electronic apparatus and storage medium
US11393443B2 (en) Apparatuses and methods for creating noise environment noisy data and eliminating noise
Sandoval et al. Theory of the hilbert spectrum
CN106249282B (en) A kind of frequency domain seismic channel set generation method suitable for AVAF invertings
CN114089806A (en) Waveform sequence creating apparatus, method, electronic device, and storage medium
Bailey et al. A digital matched filter for reverse time chaos
CN105468860A (en) Method for implementing digital signal processing algorithms based on Praat script
CN102543091B (en) System and method for generating simulation sound effect
CN114446316B (en) Audio separation method, training method, device and equipment of audio separation model
Siohan et al. Derivation of extended Gaussian functions based on the Zak transform
US8462984B2 (en) Data pattern recognition and separation engine
CN105093325B (en) It is a kind of quantitative to carry frequency method
Koutsomichalis Mapping and visualization with SuperCollider
CN114019236A (en) Power grid harmonic single-channel aliasing target signal detection method and device
CA3117942A1 (en) Systems and methods of signal analysis and data transfer using spectrogram construction and inversion
CN108418672A (en) A kind of voice communication encryption method based on inverse time chaos matched filter
Liu et al. Automatic radar waveform recognition based on neural network
Muda et al. New Optimized Adaptive Time Series IDS Classifier Algorithm: Beyond Deep Learning
Dobre et al. TIC-TAC based live acoustic watermarking with improved forgery detection performances
Helffrich et al. The phase-locked loop as a tool for signal analysis
CN115798504A (en) Single-channel voice separation method, system, equipment and storage medium

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160406

WD01 Invention patent application deemed withdrawn after publication