CN104133404B - A kind of signal processing method and device - Google Patents

Abstract

This application provides a kind of signal processing method, including: according to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；Analog signal is carried out A/D conversion, obtains data signal；According to FIR algorithm and fft algorithm, data signal is carried out spectrum analysis, obtain the first signal；According to presetting correction algorithm, frequency and the amplitude correction of the first signal are processed, obtain secondary signal；According to predesigned compensation rule, the amplitude compensation to secondary signal, obtain the 3rd signal and export to calculating equipment.Use the method, after data signal being carried out spectrum analysis according to FIR algorithm and fft algorithm, carry out frequency and the correcting process of amplitude, amplitude is compensated.It is modified signal processing, has cut down the spectrum leakage phenomenon of fft algorithm, alleviated the fence effect of fft algorithm simultaneously；And signal is carried out amplitude compensation, compensate for the logical fluctuation of FIR band and the problem of intermediate zone decay, improve the precision of signal processing results.

Description

A kind of signal processing method and device

Technical field

The present invention relates to signal transacting field, in particular, relate to a kind of signal processing method and dress Put.

Background technology

Along with development and the further investigation of control algolithm of large scale integrated circuit, Digital Signal Processing with Its precision is high, reliability is high, flexibility is big, be prone to the advantages such as large-scale integrated, all obtains in every field To extensively application.

Digital Signal Processing is that Serial No. is processed by the method using numerical computations, translates the signals into Become certain form suited the requirements.As, data signal is filtered limiting other frequency bands or filter and makes an uproar Sound and interference, or they are separated with other signals, signal is carried out spectrum analysis to understand signal Frequency spectrum composition, and then signal is identified.

In prior art, FIR (Finite impact Response, have limit for length's impulse response) algorithm and FFT (Fast Fourier Transform, FFT) algorithm is the main method of spectrum analysis.As Traditional signal processing flow shown in Fig. 1, primary signal is via signal condition, A/D After (Analog-to-digital, analog/digital) conversion, deliver to digital signal processing unit and carry out frequency spectrum and divide Analysis, then calculates.Wherein, this spectrum analysis link is made up of FIR and FFT.

Owing to low order FIR exists the logical big and problem of intermediate zone decay that fluctuates of band, cause the letter of different frequency point Number attenuation coefficient difference, can cause the signal of critical frequency to judge by accident.And there is impact when application in fft algorithm The spectrum leakage problem of measurement accuracy.

Therefore, use method of the prior art when carrying out signal transacting, easily occur that signal transacting is tied The problem that really precision is relatively low.

Summary of the invention

In view of this, the invention provides a kind of signal processing method and device, solve in prior art The problem that signal processing results precision is relatively low.

For achieving the above object, the present invention provides following technical scheme:

A kind of signal processing method, including:

According to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；

Described analog signal is simulated/numeral A/D conversion, obtains data signal；

According to having limit for length's impulse response FIR algorithm and FFT fft algorithm to described numeral letter Number carry out spectrum analysis, obtain the first signal；

It is modified processing to the frequency of described first signal and amplitude according to presetting correction algorithm, obtains the Binary signal；

According to predesigned compensation rule, the amplitude compensation to described secondary signal, obtain the 3rd signal and export To calculating equipment.

Above-mentioned method, it is preferred that described correction algorithm includes FFT interpolation algorithm, described foundation is preset Correction algorithm is modified process specifically includes frequency and the amplitude of described first signal:

The discrete spectrum of described secondary signal is calculated according to the FFT interpolation algorithm preset；

Select the spectral line of predetermined number in described discrete spectrum；

Calculate the spectral line amplitude that each spectral line is corresponding；

It is calculated correction according to the window function preset and described spectral line amplitude and default fitting condition Information, and according to described update information, frequency and the amplitude of described first signal are modified, obtain Binary signal.

Above-mentioned method, it is preferred that described according to predesigned compensation rule, the amplitude to described secondary signal Compensation specifically includes:

Compensate from default FIR filter and table is searched the compensated information corresponding with the frequency of described secondary signal；

The amplitude of described secondary signal tabled look-up compensation according to described compensated information, obtain the 3rd signal.

Above-mentioned method, it is preferred that in the described discrete spectrum of described selection, the spectral line of predetermined number includes:

Determine the spectral line of amplitude maximum in described discrete spectral line；

Select remaining spectral line in the described predetermined number adjacent with described maximum spectral line.

Above-mentioned method, it is preferred that window function that described foundation is preset and described spectral line amplitude, and in advance If fitting condition, frequency and amplitude to described first signal be modified including:

The first function and the second function is constructed respectively according to described amplitude；

According to described default window function, the fitting function of MATLAB, respectively to described first function and Second Function Fitting, obtains the actual frequency value function of described first signal and the multinomial of actual magnitude function Formula fitting formula, and utilize described fitting polynomial formulas that frequency and the amplitude of described first signal are carried out Correcting process.

Above-mentioned method, it is preferred that described window function uses Blackman window function.

Above-mentioned method, it is preferred that described determine that in described discrete spectral line, the spectral line of amplitude maximum includes:

It is as the criterion with described predeterminated position, described discrete spectral line is used the sequence of bubble sort method, obtains amplitude Maximum spectral line.

Above-mentioned method, it is preferred that described predetermined number includes: 2,3 or 5.

A kind of signal processing apparatus, including:

Conditioning module, for nursing one's health primary signal according to the conditioning condition preset, obtains simulation letter Number；

Modular converter, for described analog signal is carried out A/D conversion, obtains data signal；

Spectrum analysis module, for carrying out frequency spectrum according to FIR algorithm and fft algorithm to described data signal Analyze, obtain the first signal；

Correcting module, for repairing frequency and the amplitude of described first signal according to default correction algorithm Just process, obtain secondary signal；

Compensating module, for according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtaining 3rd signal also exports to calculating equipment.

Above-mentioned device, it is preferred that described correction algorithm includes FFT interpolation algorithm, described correcting module Specifically for:

The discrete spectrum of described secondary signal is calculated according to the FFT interpolation algorithm preset；

Select the spectral line of predetermined number in described discrete spectrum；

Calculate the spectral line amplitude that each spectral line is corresponding；

It is calculated correction according to the window function preset and described spectral line amplitude and default fitting condition Information, and according to described update information, frequency and the amplitude of described first signal are modified, obtain Binary signal.

Understand via above-mentioned technical scheme, compared with prior art, the invention provides at a kind of signal Reason method, including: according to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal； Described analog signal is carried out A/D conversion, obtains data signal；Foundation FIR algorithm and fft algorithm are to institute State data signal and carry out spectrum analysis, obtain the first signal；According to presetting correction algorithm to described first letter Number frequency and amplitude be modified process, obtain secondary signal；According to predesigned compensation rule, to described The amplitude compensation of secondary signal, obtains the 3rd signal and exports to calculating equipment.Use the method, at root Carry out data signal after spectrum analysis completes, proceeding frequency and amplitude according to FIR algorithm and fft algorithm Correcting process, and amplitude is compensated.It is modified signal processing, has cut down fft algorithm Spectrum leakage phenomenon, alleviates the fence effect of fft algorithm simultaneously；And signal is carried out amplitude compensation, more Mend the logical fluctuation of FIR band and the problem of intermediate zone decay, improve the precision of signal processing results.

Accompanying drawing explanation

In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to reality Execute the required accompanying drawing used in example or description of the prior art to be briefly described, it should be apparent that below, Accompanying drawing in description is only embodiments of the invention, for those of ordinary skill in the art, not On the premise of paying creative work, it is also possible to obtain other accompanying drawing according to the accompanying drawing provided.

Fig. 1 is traditional signal processing flow figure；

The flow chart of a kind of signal processing method embodiment 1 that Fig. 2 provides for the application；

The flow chart of a kind of signal processing method embodiment 2 that Fig. 3 provides for the application；

The idiographic flow of step S307 in a kind of signal processing method embodiment 2 that Fig. 4 provides for the application Figure；

The flow chart of a kind of signal processing method embodiment 3 that Fig. 5 provides for the application；

The response diagram of example in a kind of signal processing method embodiment 3 that Fig. 6 provides for the application；

The frequency of example and amplitude attenuation in a kind of signal processing method embodiment 3 that Fig. 7 provides for the application Graph of a relation between coefficient；

The structural representation of a kind of signal processing apparatus embodiment 1 that Fig. 8 provides for the application.

Detailed description of the invention

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out Clearly and completely describe, it is clear that described embodiment is only a part of embodiment of the present invention, and It is not all, of embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art are not doing Go out the every other embodiment obtained under creative work premise, broadly fall into the scope of protection of the invention.

The flow chart of a kind of signal processing method that the embodiment of the present application 1 as shown in Figure 2 provides, including:

Step S201: according to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；

Wherein, this conditioning condition can be the condition of follow-up A/D transformational structure input signal.

Concrete, this conditioning process includes: power amplification, filtering, linear compensation, isolates and/or protects Protecting, the process of conditioning can be set according to actual conditions, and the present embodiment does not repeats.

Step S202: described analog signal is carried out A/D conversion, obtains data signal；

Wherein, owing to follow-up FIR algorithm and fft algorithm need according to digital format signal, and pass through Obtain after conditioning for analog format signal, therefore, this analog signal is carried out A/D conversion, makes this letter Number be converted to analog signal.

Step S203: described data signal is carried out spectrum analysis according to FIR algorithm and fft algorithm, To the first signal；

Wherein, successively data signal is carried out spectrum analysis according to FIR algorithm and fft algorithm, obtain the One signal.

Concrete, this spectrum analysis process is: uses fft algorithm to realize the conversion to data signal, adopts The filtering to data signal is realized with FIR algorithm.

But, use fft algorithm, this conversion process easily occur spectrum leakage phenomenon；And use FIR There is the logical big and problem of intermediate zone decay that fluctuates of band in the filtering of algorithm.

Therefore, in follow-up step S204-205, so that the problems referred to above are made up.

Step S204: frequency and the amplitude of described first signal are modified place according to presetting correction algorithm Reason, obtains secondary signal；

Wherein, before carrying out signal transacting, preset correction algorithm.

Concrete, according to this correction algorithm, frequency and amplitude to the first signal are modified respectively, with Overcome the problem that fft algorithm carries out the spectrum leakage occurred during digital conversion.

Concrete makeover process illustrates in follow-up embodiment, and the present embodiment does not repeats.

Step S205: according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtain the 3rd Signal also exports to calculating equipment.

Wherein, before carrying out signal transacting, predesigned compensation rule.

Concrete, according to this offset rule, the secondary signal through revising is carried out amplitude compensation, with gram What clothes employing FIR algorithm filtering caused carries logical fluctuation greatly and the problem of intermediate zone decay.

Wherein, the 3rd signal exports to follow-up calculating equipment, to carry out follow-up calculating process.

During it should be noted that the application carries out spectrum analysis to data signal, calculating for FIR of employing Method and fft algorithm, in actual enforcement, when carrying out spectrum analysis, it is possible to only with FIR algorithm and FFT One in algorithm, accordingly, when use for above-mentioned two kinds of algorithms time, perform step S204 and 205, The signal completing spectrum analysis carries out the correction of frequency, amplitude, and amplitude compensation；When use only When being FIR algorithm, only carrying out step S205, the signal completing spectrum analysis carries out amplitude compensation；When During the simply fft algorithm used, only carry out step S204, the signal that spectrum analysis is completed carry out frequency, The correction of amplitude.

Certainly, in reality is implemented, when carrying out spectrum analysis, also can be only with FIR algorithm and FFT Other algorithms outside algorithm, when the result of spectrum analysis also has the logical fluctuation of band big and the problem of intermediate zone decay Time, then perform step S205；When the result of spectrum analysis also has spectrum leakage phenomenon, then perform step S204。

It should be noted that in actual enforcement, the method in the application may migrate to embedded software or DSP (digital signal processor, digital signal processor), inputs this embedded software by primary signal Or DSP realizes the method.

To sum up, present embodiments provide a kind of signal processing method, including: according to the conditioning condition preset Primary signal is nursed one's health, obtains analog signal；Described analog signal is carried out A/D conversion, is counted Word signal；According to FIR algorithm and fft algorithm, described data signal is carried out spectrum analysis, obtain the first letter Number；It is modified processing to the frequency of described first signal and amplitude according to presetting correction algorithm, obtains the Binary signal；According to predesigned compensation rule, the amplitude compensation to described secondary signal, obtain the 3rd signal also Output is to calculating equipment.Use the method, according to FIR algorithm and fft algorithm, data signal is being carried out frequently After analysis of spectrum completes, proceed the correcting process of frequency and amplitude, and amplitude is compensated.Right Signal is modified processing, and has cut down the spectrum leakage phenomenon of fft algorithm, has alleviated fft algorithm simultaneously Fence effect；And signal is carried out amplitude compensation, compensate for the logical fluctuation of FIR band and the problem of intermediate zone decay, Improve the precision of signal processing results.

In implementing, this correction algorithm can use FFT interpolation algorithm.

Such as Fig. 3, it is shown that the flow chart of a kind of signal processing method embodiment 2 that provides of the application, including:

Step S301: according to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；

Step S302: described analog signal is carried out A/D conversion, obtains data signal；

Step S303: described data signal is carried out spectrum analysis according to FIR algorithm and fft algorithm, To the first signal；

Wherein, step S301-303 is consistent with step S201-203 in embodiment 1, and the present embodiment does not does Repeat.

Step S304: calculate the discrete spectrum of described first signal according to the FFT interpolation algorithm preset；

Concrete, according to the first sample frequency preset, the first signal is sampled, obtain discrete signal； Utilize Fourier transformation and the window function preset, the time-domain signal of this discrete signal be converted to frequency-region signal, Obtain the continuous frequency spectrum of this first signal；This continuous frequency spectrum is continued sampling and obtains some discrete spectral lines, This sample frequency is relevant to the first sample frequency.

Such as, this first signal is frequency f₀, amplitude be A, phase angle beEndless sine wave signal x (t), The expression formula of this first signal is:

With sample frequency f_sAfter x (t) is sampled, obtain the discrete signal of following form:

The forms of time and space of the window function used is w (n), and its continuous frequency spectrum is W (2 π f), i.e.

$W\left(2\mathrm{\π}f\right)=\underset{n=-\∝}{\overset{n=+\∝}{\Σ}}w\left(n\right)e^{-j2\mathrm{\π}fn}---\left(3\right)$

Utilizing Fourier transformation, will be brill red with signal and for frequency domain, the continuous frequency spectrum obtaining x (t) signal be

Ignore spectral peak secondary lobe impact at negative frequency, obtain

To this link spectral with sample frequencySampling, wherein, N is data length, obtains

The formula (6) that this sampling obtains is some discrete spectral lines.

Step S305: select the spectral line of predetermined number in described discrete spectrum；

Concrete, comprise determining that the spectral line of amplitude maximum in described discrete spectral line；Select and described maximum Remaining spectral line in the described predetermined number that spectral line is adjacent.

Wherein determine the process of the spectral line of amplitude maximum in described discrete spectral line particularly as follows: with described default position It is set to standard, described discrete spectral line is used the sequence of bubble sort method, obtains the spectral line of amplitude maximum.

Wherein, described predetermined number includes: 2,3 or 5.

Wherein, when for 2, the spectral line chosen includes: near the spectral line of amplitude maximum and this spectral line Secondary big spectral line；When for 3, the spectral line chosen includes: phase before and after the spectral line of amplitude maximum and this spectral line 2 adjacent spectral lines；When for 5, the spectral line chosen includes: before the spectral line of amplitude maximum and this spectral line 4 the most adjacent spectral lines, the most each two spectral lines.

The present embodiment illustrates as a example by choosing 3 spectral lines.

Ifk₀Possible position integer, its spectral line greatest around is k₁, k₁Left side spectral line is k₂、 Right side spectral line is k₃, selection spectral line is k₁、k₂And k₃。

Step S306: calculate the spectral line amplitude that each spectral line is corresponding；

Wherein, calculate the offset distance distance values of this spectral line and predeterminated position according to the spectral line selected in step S305, And it is calculated, according to the discrete formula that obtains of sampling, the amplitude that the spectral line of selection is corresponding.

Such as, selecting spectral line is k₁、k₂And k₃, the amplitude of each bar spectral line is respectively Introduce parameter alpha=k₀-k₁。

Due to k₁For k₀Spectral line greatest around, so the span of α is [-0.5,0.5], has reacted k₀Partially From k₁Distance, obtain k₁=k₀-α、k₂=k₀-α-1、k₂=k₀-α+1, substitutes into formula (6) and obtains

$y_1=\stackrel{\‾}{X}\left(k_1\mathrm{\Δ}f\right)=\frac{A}{2}W\left(2\mathrm{\π}\right(\frac{k_1\frac{f_s}{N}-f_0}{f_s}\left)\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(k_1-k_0)}{N}\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α})}{N}\right)---\left(7\right)$

$y_2=\stackrel{\‾}{X}\left(k_2\mathrm{\Δ}f\right)=\frac{A}{2}W\left(2\mathrm{\π}\right(\frac{k_2\frac{f_s}{N}-f_0}{f_s}\left)\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(k_2-k_0)}{N}\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α}-1)}{N}\right)---\left(8\right)$

$y_3=\stackrel{\‾}{X}\left(k_3\mathrm{\Δ}f\right)=\frac{A}{2}W\left(2\mathrm{\π}\right(\frac{k_3\frac{f_s}{N}-f_0}{f_s}\left)\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(k_3-k_0)}{N}\right)=\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α}+1)}{N}\right)---\left(9\right)$

Step S307: according to the window function preset and described spectral line amplitude and the fitting condition meter preset Calculation obtains update information, and repaiies frequency and the amplitude of described first signal according to described update information Just, secondary signal is obtained；

Wherein, according to above-mentioned spectral line and the amplitude of spectral line, combine the window function preset, be calculated One group with pending the with parameter group corresponding to signal, this parameter group is carried out with the fitting formula preset Matching, obtains the update information of amplitude, frequency, according to this update information to the amplitude of the first signal and frequency Rate is modified, and obtains secondary signal.

As shown in Figure 4, the particular flow sheet of this step S307, including:

Step S401: construct the first function and the second function respectively according to described amplitude；

Example in continuation, according to spectral line k₁、k₂And k₃Amplitude constructed fuction 1 and function 2:

Function 1:

$\frac{y_2+y_1}{y_3+y_1}=\frac{\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α}-1)}{N}\right)+\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α})}{N}\right)}{\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α}+1)}{N}\right)+\frac{A}{2}W\left(\frac{2\mathrm{\π}(-\mathrm{\α})}{N}\right)}=g\left(\mathrm{\α}\right)---\left(10\right)$

Function 2:

$y_1+y_2+y_3=\frac{A}{2}\left(W\right(\frac{2\mathrm{\π}(-\mathrm{\α})}{N})+W(\frac{2\mathrm{\π}(-\mathrm{\α}-1)}{N})+W(\frac{2\mathrm{\π}(-\mathrm{\α}+1)}{N}\left)\right)=\frac{A}{2}\mathrm{\ρ}\left(\mathrm{\α}\right)---\left(11\right)$

Step S402: according to described default window function, the fitting function of MATLAB, respectively to described First function and the second Function Fitting, obtain actual frequency value function and the actual magnitude of described first signal The fitting polynomial formulas of function, and utilize the described fitting polynomial formulas frequency to described first signal It is modified processing with amplitude.

Wherein, in the present embodiment, this window function uses Blackman window function.

OrderSimulate the inverse function of g (x), i.e.

α=g^-1(β)=a₀+a₁β+a₂β²+...+a_nβⁿ (12)

Wherein, a₀、a₁、a₂...a_nCan any real number of value.

Concrete available MATLAB (MATrix LABoratory, matrix labotstory) associated polynomials The inverse function of interpolation instrument matching g (x).

It follows that the actual frequency values of the first signal is

f₀=k₀Δ f=(k₁+α)Δf (13)

In like manner, polynomial interopolation is utilized to simulate the inverse function of ρ (x) in formula (11), i.e.

A=N^-1(y₁+y₂+y₃)ρ^-1(α)=N^-1(y₁+y₂+y₃)(b₀+b₁α+b₂α²+...+b_nαⁿ) (14)

Wherein, b₀、b₁、b₂...b_nCan any real number of value.

Wherein, the time-domain expression of Blackman window function is

$w\left(n\right)=0.42-0.5cos\left(\frac{2\mathrm{\π}n}{N}\right)+0.08\left(\frac{2\mathrm{\π}}{N}2n\right)---\left(15\right)$

Owing to the frequency of the first processed signal is consistent, then α can be calculated, and β also can be calculated, In like manner,Also can be calculated.

Change sample frequency, such as f=f+0.001, repeats in above-mentioned steps and samples x (t), with The calculating process needed after and, finally draws one group of α, β, γ value.

The fitting function utilizing MATLAB uses polyfit, xx=polyfitg^-1(β,α,3)、 Yy=polyfitg^-1(α,γ,3)。

Polynomial fit function p=polyfitg^-1(x, y, n) in, the top step number of n representative polynomial, x, y are Will the data of matching, p is the multinomial coefficient p=[p that matching obtains₁,p₂,...,p_n,p_n+1], and meet Y=p₁xⁿ+p₂x^n-1+...+p_nx+p_n+1。

In conjunction with above-mentioned fitting function, obtain the fitting polynomial formulas (matching of frequency values of formula (12) Amount)

α=1.9458-2.6332 β+0.7089 β²-0.0221β³ (16)

In conjunction with above-mentioned fitting function, obtain the fitting polynomial formulas (the matching amount of amplitude) of formula (14)

A=N^-1(y₁+y₂+y₃)(2.1770+0.0006α+0.5631α²-0.0105α³) (17)

According to the matching amount in formula (16) and (17), to amplitude A in the first signal and frequency f₀Enter Row compensates, and obtains secondary signal.

It should be noted that the present embodiment have employed the mode of approximation by polynomi-als, FFT correction algorithm number Little according to amount of calculation, simplify calculating process.

Further, during compensating, the method that have employed sampling, reduce the expense of system.

Step S308: according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtain the 3rd Signal also exports to calculating equipment.

Wherein, step S308 is consistent with step S205 in embodiment 1, and the present embodiment does not repeats.

It should be noted that window function in the present embodiment use for Blackman window function, certainly, It is not limited to this, in actual enforcement, also can use rectangular window function, Hanning window function or hamming window Function, is not explained in the present embodiment.

To sum up, a kind of information processing method that the present embodiment provides, use FFT interpolation algorithm to the first letter Number amplitude and frequency be modified, obtain secondary signal.Due to this FFT interpolation algorithm, and above walk The rapid fft algorithm that described data signal carries out spectrum analysis employing is corresponding, it is possible to have targetedly Problem to the spectrum leakage using fft algorithm to occur when carrying out spectrum analysis.

In implementing, this can use FIR filter compensation table to compensate the amplitude of signal.

Such as Fig. 5, it is shown that the flow chart of a kind of signal processing method embodiment 3 that provides of the application, including:

Step S501: according to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；

Step S502: described analog signal is carried out A/D conversion, obtains data signal；

Step S503: described data signal is carried out spectrum analysis according to FIR algorithm and fft algorithm, To the first signal；

Step S504: frequency and the amplitude of described first signal are modified place according to presetting correction algorithm Reason, obtains secondary signal；

Wherein, step S501-504 is consistent with step S201-204 in embodiment 1, and the present embodiment does not does Repeat.

Step S505: compensate lookup table from default FIR filter corresponding with the frequency of described secondary signal Compensated information；

Wherein, this FIR filter compensation table comprises the corresponding relation of frequency and attenuation coefficient.

Concrete, according to the frequency of this secondary signal, from this FIR filter compensation table, search the benefit of correspondence Repay information, i.e. attenuation coefficient.

It should be noted that this compensating coefficient table can be unified being arranged in calculating equipment, it is possible to for It is preset in other storage devices, the present embodiment does not limits.

Step S506: the amplitude of described secondary signal tabled look-up compensation according to described compensated information, To the 3rd signal and export to calculating equipment.

Wherein, according to the compensated information found and attenuation coefficient, the amplitude of secondary signal is compensated, To make up, use FIR algorithm the to carry out logical fluctuation of band that bandpass filtering brings is big and intermediate zone decays.

Compensation table is described below sets up process.

With sample frequency as 500Hz, exponent number is 128,2～is described in detail as a example by 12Hz bandpass filtering.

MATLAB software is utilized to generate one group of filter factor according to above-mentioned requirements；MATLAB is utilized to draw Go out the frequency response chart that category coefficient is corresponding, this response diagram such as Fig. 6；Fig. 6 medium frequency is responded transverse axis number It is converted into actual frequency according to by normalized frequency, and only focuses on intermediate zone and the logical interior frequency component of band, Longitudinal axis data are converted into attenuation coefficient, obtain the frequency shown in Fig. 7 and the relation between amplitude attenuation coefficient Figure；Observation Fig. 7 understands, and 0Hz component decays to 0.5 times of original component, and 7Hz component decays to former 1 times of beginning component, 14Hz component decays to 0.25 times of original component.Obvious, this filter passband hinders Band fluctuation is big, if not processing, it is difficult to meet requirement of engineering.Said frequencies response is changing into form Form, so that completing to compare frequency values after FFT calculating terminates to carry out the correction of amplitude.

To sum up, a kind of information processing method that the present embodiment provides, use FIR filter to compensate table to signal Amplitude compensate, obtain the 3rd signal.This FIR filter compensation table is default form, when carrying out During compensation, only need to make a look up from this compensation table, data meter according to the frequency of signal to be compensated Calculation amount is little, and this FIR filter compensates table, with described data signal is carried out frequency spectrum the dividing of preceding step The FIR algorithm that analysis uses is corresponding, it is possible to have targetedly to when using FIR algorithm to carry out spectrum analysis The logical big and problem of intermediate zone decay that fluctuates of band occurred.

A kind of signal processing method embodiment provided with above-mentioned the application is corresponding, and the application also provides for A kind of signal processing apparatus embodiment.

Such as Fig. 8, it is shown that the structural representation of a kind of signal processing apparatus embodiment 1 that the application provides, Including: conditioning module 801, modular converter 802, spectrum analysis module 803, correcting module 804 and benefit Repay module 805；

Conditioning module 801, for nursing one's health primary signal according to the conditioning condition preset, is simulated Signal；

Wherein, the conditioning condition of this conditioning module 801 can be follow-up A/D transformational structure input signal Condition.

Concrete, this conditioning process includes: power amplification, filtering, linear compensation, isolates and/or protects Protecting, the process of conditioning can be set according to actual conditions, and the present embodiment does not repeats.

Modular converter 802, for described analog signal is carried out A/D conversion, obtains data signal；

Wherein, owing to follow-up FIR algorithm and fft algorithm need according to digital format signal, and pass through Obtain after conditioning for analog format signal, therefore, modular converter 802 carries out A/D to this analog signal Conversion, makes this signal be converted to analog signal.

Spectrum analysis module 803, for carrying out described data signal according to FIR algorithm and fft algorithm Spectrum analysis, obtains the first signal；

Wherein, data signal is carried out successively by spectrum analysis module 803 according to FIR algorithm and fft algorithm Spectrum analysis, obtains the first signal.

Concrete, this spectrum analysis module 803 spectrum analysis process is: use fft algorithm to realize logarithm The conversion of word signal, uses FIR algorithm to realize the filtering to data signal.

But, use fft algorithm, this conversion process easily occur spectrum leakage phenomenon；And use FIR There is the logical big and problem of intermediate zone decay that fluctuates of band in the filtering of algorithm.

Therefore, at follow-up correcting module 804 and compensating module 805, so that the problems referred to above are made up.

Correcting module 804, for carrying out frequency and the amplitude of described first signal according to default correction algorithm Correcting process, obtains secondary signal；

Wherein, before carrying out signal transacting, correcting module 804 is preset correction algorithm.

Concrete, correcting module 804 is according to this correction algorithm, frequency and the amplitude difference to the first signal It is modified, with the problem overcoming spectrum leakage that fft algorithm occurs when carrying out digital conversion.

Wherein, the 3rd signal exports to follow-up calculating equipment, to carry out follow-up calculating process.

Wherein, described correction algorithm includes FFT interpolation algorithm, this correcting module 804 specifically for: depend on The discrete spectrum of described secondary signal is calculated according to default FFT interpolation algorithm；Select described discrete spectrum The spectral line of middle predetermined number；Calculate the offset distance distance values of each spectral line and predeterminated position, and described skew The spectral line amplitude that distance value is corresponding；According to the window function preset and described spectral line amplitude and the plan preset Conjunction condition is calculated update information, and according to described update information to the frequency of described first signal and width Value is modified, and obtains secondary signal.

Concrete makeover process illustrates in above-mentioned embodiment of the method, and the present embodiment does not repeats.

Compensating module 805, for according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtaining To the 3rd signal and export to calculating equipment.

Wherein, before carrying out signal transacting, to compensating module 805 predesigned compensation rule.

Concrete, compensating module 805, according to this offset rule, carries out width to the secondary signal through revising Value complement is repaid, to overcome the logical fluctuation of the band using FIR algorithm filtering to cause big and the problem of intermediate zone decay.

During it should be noted that the application carries out spectrum analysis to data signal, calculating for FIR of employing Method and fft algorithm, in actual enforcement, when carrying out spectrum analysis, it is possible to only with FIR algorithm and FFT One in algorithm, accordingly, when use for above-mentioned two kinds of algorithms time, trigger correcting module 804 With compensating module 805, the signal completing spectrum analysis carries out the correction of frequency, amplitude, and amplitude is mended Repay；When the simply FIR algorithm used, only trigger compensating module 805, the signal that spectrum analysis is completed Carry out amplitude compensation；When the simply fft algorithm used, only trigger correcting module 804, frequency spectrum is divided The signal analysed carries out the correction of frequency, amplitude.

Certainly, in reality is implemented, when carrying out spectrum analysis, also can be only with FIR algorithm and FFT Other algorithms outside algorithm, when the result of spectrum analysis also has the logical fluctuation of band big and the problem of intermediate zone decay Time, then trigger compensating module 805；When the result of spectrum analysis also has spectrum leakage phenomenon, then trigger and repair Positive module 804.

Implement it should be noted that actual, the method in the application may migrate to embedded software or DSP, inputs primary signal in this embedded software or DSP and realizes the method.

To sum up, present embodiments provide a kind of signal processing apparatus, including: conditioning module, for foundation Primary signal is nursed one's health by the conditioning condition preset, and obtains analog signal；Modular converter, for by institute State analog signal and carry out A/D conversion, obtain data signal；Spectrum analysis module, for calculating according to FIR Method and fft algorithm carry out spectrum analysis to described data signal, obtain the first signal；Correcting module, uses In being modified processing to frequency and the amplitude of described first signal according to default correction algorithm, obtain second Signal；Compensating module, for according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtaining To the 3rd signal and export to calculating equipment..Use this device, according to FIR algorithm and fft algorithm Carry out data signal after spectrum analysis completes, proceeding the correcting process of frequency and amplitude, and right Amplitude compensates.It is modified signal processing, has cut down the spectrum leakage phenomenon of fft algorithm, with Time alleviate the fence effect of fft algorithm；And signal is carried out amplitude compensation, compensate for the logical ripple of FIR band The problem that dynamic and intermediate zone decays, improves the precision of signal processing results.

Claims (9)

1. a signal processing method, it is characterised in that including:

According to the conditioning condition preset, primary signal is nursed one's health, obtain analog signal；

Described analog signal is simulated/numeral A/D conversion, obtains data signal；

According to having limit for length's impulse response FIR algorithm and FFT fft algorithm to described numeral letter Number carry out spectrum analysis, obtain the first signal；

It is modified processing to the frequency of described first signal and amplitude according to presetting correction algorithm, obtains the Binary signal；

According to predesigned compensation rule, the amplitude compensation to described secondary signal, obtain the 3rd signal and export To calculating equipment；

Wherein, described according to predesigned compensation rule, the amplitude compensation of described secondary signal is specifically included:

Compensate from default FIR filter and table is searched the compensated information corresponding with the frequency of described secondary signal；

The amplitude of described secondary signal tabled look-up compensation according to described compensated information, obtain the 3rd signal.

Method the most according to claim 1, it is characterised in that described correction algorithm includes that FFT inserts Value-based algorithm, described foundation is preset correction algorithm and is modified processing to frequency and the amplitude of described first signal Specifically include:

The discrete spectrum of described secondary signal is calculated according to the FFT interpolation algorithm preset；

Select the spectral line of predetermined number in described discrete spectrum；

Calculate the spectral line amplitude that each spectral line is corresponding；

It is calculated correction according to the window function preset and described spectral line amplitude and default fitting condition Information, and according to described update information, frequency and the amplitude of described first signal are modified, obtain Binary signal.

Method the most according to claim 2, it is characterised in that in the described discrete spectrum of described selection The spectral line of predetermined number includes:

Determine the spectral line of amplitude maximum in described discrete spectral line；

Select remaining spectral line in the described predetermined number adjacent with described maximum spectral line.

Method the most according to claim 2, it is characterised in that described foundation preset window function and Described spectral line amplitude, and the fitting condition preset, frequency and amplitude to described first signal are repaiied Just include:

The first function and the second function is constructed respectively according to described amplitude；

According to described default window function, the fitting function of matrix labotstory MATLAB, respectively to described First function and the second Function Fitting, obtain actual frequency value function and the actual magnitude of described first signal The fitting polynomial formulas of function, and utilize the described fitting polynomial formulas frequency to described first signal It is modified processing with amplitude.

Method the most according to claim 2, it is characterised in that described window function uses Blacknam Window function.

Method the most according to claim 3, it is characterised in that described determine in described discrete spectral line The spectral line of amplitude maximum includes:

It is as the criterion with described predeterminated position, described discrete spectral line is used the sequence of bubble sort method, obtains amplitude Maximum spectral line.

Method the most according to claim 2, it is characterised in that described predetermined number includes: 2, 3 or 5.

8. a signal processing apparatus, it is characterised in that including:

Conditioning module, for nursing one's health primary signal according to the conditioning condition preset, obtains simulation letter Number；

Modular converter, for described analog signal is carried out A/D conversion, obtains data signal；

Spectrum analysis module, for carrying out frequency spectrum according to FIR algorithm and fft algorithm to described data signal Analyze, obtain the first signal；

Correcting module, for repairing frequency and the amplitude of described first signal according to default correction algorithm Just process, obtain secondary signal；

Compensating module, for according to predesigned compensation rule, the amplitude compensation to described secondary signal, obtaining 3rd signal also exports to calculating equipment；

Wherein, described compensating module specifically for:

Compensate from default FIR filter and table is searched the compensated information corresponding with the frequency of described secondary signal；

The amplitude of described secondary signal tabled look-up compensation according to described compensated information, obtain the 3rd signal.

Device the most according to claim 8, it is characterised in that described correction algorithm includes that FFT inserts Value-based algorithm, described correcting module specifically for:

The discrete spectrum of described secondary signal is calculated according to the FFT interpolation algorithm preset；

Select the spectral line of predetermined number in described discrete spectrum；

Calculate the spectral line amplitude that each spectral line is corresponding；

It is calculated correction according to the window function preset and described spectral line amplitude and default fitting condition Information, and according to described update information, frequency and the amplitude of described first signal are modified, obtain Binary signal.