DIGITAL OSCILLOSCOPE HAVING FRACTIONAL CALCULUS OPERATION AND DISPLAY FUNCTION

Abstract

The present invention provides a system for data mapping and storing in digital three-dimensional oscilloscope, wherein the fixed coefficients, which are calculated according the parameters and settings of a digital oscilloscope, are stored into a fixed coefficient memory CO RAM, the fixed coefficients are outputted to N fractional operation units through N−1 D flip-flop delay units to multiply with the acquired data x(n) and then be accumulated, thus N fractional calculus results are obtained. In this way, N fractional calculus results can be obtained by performing L/N fractional calculus operations. N fractional calculus results are sent to a signal processing and display module, in which they are converted into a display data through a drawing thread, and the display data are sent to LCD for displaying, thus the fractional calculus operation and display of a input signal in a digital oscilloscope is realized.

Claims

1. A digital oscilloscope having fractional calculus operation and display function, comprising: a digital oscilloscope control module; an ADC module, wherein an input analog signal (a signal being measured) x(t) is sampled and quantified to obtain an sampled data ADC(m), m is the sampling point of the sampled data ADC(m); a signal acquisition and storage module, wherein the sampled data ADC(m) is decimated according to selected time base, and then under the control of the digital oscilloscope control module, the sampled data after decimation (acquired data) is stored into the acquisition memory ADC RAM of the signal acquisition and storage module in order based on storage address according to a set pre trigger depth and trigger signal, when the acquisition memory ADC RAM is stored full, which means L acquired data is stored, the acquisition and storage is stopped, the oscilloscope data done signal DSO_DONE outputted by the acquisition memory ADC RAM is turned from 0 to 1, at this moment, the acquired data stored in the acquisition memory ADC RAM can be denoted by x(n), n is storage address, n=0, 1, 2, . . . , L−1; a signal processing and display module; wherein further comprising: a fractional calculus module, which is used for reading out the acquired data x(n) from the acquisition memory ADC RAM to perform a fractional differentiation operation or a fractional integration operation and obtaining a fractional calculus result data y(n), which is sent to the signal processing and display module, in which the acquired data x(n) is turned into a display data and the display data is sent to a LCD for displaying; the fractional calculus module comprising: a fixed coefficient memory CO RAM, which is used for storing the L fixed coefficients c(0), c(1), c(2), . . . , c(L−1) sent from the digital oscilloscope control module, to the j.sup.th fixed coefficient c(j), its value is: $c\left(j\right)=\frac{b_j{L}^{\alpha }}{{\left({\mathrm{XT}}_{\mathrm{base}}\right)}^{\alpha }},j=0,1,2,.Math.,L-1$ where L is the number of the pixels at the horizontal direction of the display area of the digital oscilloscope, X is the number of the divisions at the horizontal direction of the display area of the digital oscilloscope, T.sub.base is the time base of the digital oscilloscope, α is operation order, when α<0, the fractional calculus operation is a fractional integration operation, when α>0, the fractional calculus operation is a fractional differentiation operation, b.sub.j is a fixed parameter, its value is: $b_j=\left(1-\frac{1+\alpha }{j}\right)b_{j-1},b_0=1$ N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1, wherein the output of the fixed coefficient memory CO RAM is connected to the input of the 1.sup.st D flip-flop delay unit D.sub.1, from the 2.sup.nd D flip-flop delay unit D.sub.2 on, the input of a D flip-flop delay unit is connected to the output of its previous D flip-flop delay unit; N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N, wherein each fractional operation unit comprises one multiplier and one accumulator (ACC); to the p.sup.th fractional operation unit THD.sub.p, p=1, 2, . . . , N−1, one (port A) of the two inputs of its multiplier MULTI.sub.p is connected to the acquired data x(n) outputted by the signal acquisition and storage module, another (port B) of the two inputs of its multiplier MULTI.sub.p is connected to the input of the p.sup.th D flip-flop delay unit D.sub.p, the input of its accumulator ACC.sub.p is connected to the output of its multiplier MULTI.sub.p, the output of the accumulator ACC.sub.p is taken as the output of the p.sup.th fractional operation unit THD.sub.p; to the N.sup.th fractional operation unit THD.sub.N, one of the two inputs of its multiplier MULTI.sub.N is connected to the acquired data x(n) outputted by the signal acquisition and storage module, another of the two inputs of its multiplier MULTI.sub.N is connected to the output of the N−1.sup.th D flip-flop delay unit D.sub.N-1, the input of its accumulator ACC.sub.N is connected to the output of its multiplier MULTI.sub.N; a fractional operation results memory FO RAM, which is used for storing the operation results of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N; a fractional operation control module FO CTRL, which is used for controlling the fractional operation, when the oscilloscope data done signal DSO_DONE is turned from 0 to 1 is detected, the fractional operation control module FO CTRL performs the following steps: (1): initialing initializing fractional operation number i to 1, the flag i_done of the i.sup.th fractional calculus operation to 0, fractional operation done signal FO_DONE to 0; (2): resetting resetting the initial values of the N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1 to 0, setting the read address ADDR.sub.adcram of the acquisition memory ADC RAM to N*i−1, resetting the read address ADDR.sub.coram of the fixed coefficient memory CO RAM to 0, resetting the initial values of the accumulators of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N to 0, resetting the write address ADDR.sub.forma of the fractional operation results memory FO RAM to 0; (3): reading data to multiply and accumulate driven by a system clock clk, the fractional operation control module FO CTRL initiates a operation of reading the acquisition memory ADC RAM and a operation of reading the fixed coefficient memory CO RAM at each system clock, the total number of reading is k, k=N*i; to each reading, the read address ADDR.sub.adcram is subtracted by 1, and the read address ADDR.sub.coram is added by 1; the data outputted by the read data port of the acquisition memory ADC RAM in chronological order are acquired data x(k−1), . . . , x(1), x(0), namely, the data sent to the port As of the multipliers of the N fractional operation units are acquired data x(k−1) at 0.sup.th system clock, the data sent to the port As of the multipliers of the N fractional operation units are acquired data x(k−2) at 1.sup.st system clock, the data sent to the port As of the multipliers of the N fractional operation units are acquired data x(k−2) at 2.sup.nd system clock, and so on, the data sent to the port As of the multipliers of the N fractional operation units are acquired data x(0) at (k−1).sup.th system clock; the data outputted by the read data port of the fixed coefficient memory CO RAM in chronological order are the fixed data c(0), c(1), . . . , c(k−1), the fixed data c(0), c(1), . . . , c(k−1) are serially sent to the N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1, namely, at 0.sup.th system clock, the data sent to the port B of the multiplier of the 1.sup.st fractional operation unit THD.sub.1 is the fixed data c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, at 1.sup.st system clock, the data sent to the port Bs of the multipliers of the 1.sup.st, the 2.sup.nd fractional operation units THD.sub.1, THD.sub.2 are respectively the fixed data c(1), c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, at 2.sup.nd system clock, the data sent to the port Bs of the multipliers of the 1.sup.st, the 2.sup.nd and the 3.sup.rd fractional operation units THD.sub.1, THD.sub.2, THD.sub.2 are respectively the fixed data c(2), c(1), c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, and so on, until at the (k−1).sup.th system clock, the data sent to the port Bs of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N are respectively c(k−1), . . . , c(2), c(1), c(0); after each reading, all the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N perform a multiplying operation and an accumulating operation in turn, when the number of the accumulations reaches k, the i.sup.th fractional calculus operation is accomplished, the i.sup.th fractional operation done signal i_done is turned from 0 to 1; (4): storing the results of the i.sup.th fractional calculus operation when the i.sup.th fractional operation done signal i_done is turned from 0 to 1 is detected, the outputs of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N are taken as the results y(k−1)y(k−2), y(k−3), . . . , y(k−N) of the i.sup.th fractional calculus operation and merged into a merged data, the merged data are stored into the fractional operation results memory FO RAM at address ADDR.sub.foram under the same system clock clk; (5) judging whether a frame of acquired data, namely L acquired data are calculated judging whether the current fractional operation number i is less than L/N, if yes, the write address ADDR.sub.foram is added by 1, the current fractional operation number i is added by 1, the i.sup.th fractional operation done signal i_done is set to 1, then returning to step (2); if not, the fractional operation done signal FO_DONE is set to 1, stopping the fractional calculus operation; the digital oscilloscope control module is detecting the fractional operation done signal FO_DONE, when the fractional operation done signal FO_DONE is turned from 0 to 1 is detected, the digital oscilloscope control module performs the following steps: (1): sending the acquired data x(n), n=1, 2, . . . , L−1 stored in the acquisition memory ADC RAM and the fractional calculus result data y(n), n=1, 2, . . . , L−1 stored the fractional operation results memory FO RAM to signal processing and display module, in which the acquired data x(n) and the fractional calculus result data y(n) are respectively converted into a display data through a drawing thread, and the two display data are sent to LCD for displaying; (2): resetting the oscilloscope data done signal DSO_DONE and the fractional operation done signal FO_DONE to 0 to enable the signal acquisition and storage module and the fractional calculus module.

2. A digital oscilloscope having fractional calculus operation and display function of claim 1, wherein: after sending the acquired data x(n) and the fractional calculus result data y(n) to signal processing and display module, the vertical sensitivity corresponding to the acquired data x(n) is multiplied with the number of the vertical divisions of the display area to obtain a display range R1, the vertical sensitivity corresponding to the fractional calculus result data y(n) is multiplied with the number of the vertical divisions of the display area to obtain another display range R2; the acquired data x(n) is normalized as: $d_x\left(n\right)=\left(\frac{x\left(n\right)}{2^v-1}-0.5\right)\times R_1$ the fractional calculus result data y(n) is normalized as: $d_y\left(n\right)=\left(\frac{y\left(n\right)}{\left(2^v-1\right)\times R_2/R_1}-0.5\right)\times R_2$ where ν is the digit number of the ADC module; the data d.sub.x(n) and d.sub.y(n) are sent to a drawing thread, through which the data d.sub.x(n) and d.sub.y(n) are respectively converted into a display data, and the two display data are sent to LCD 4 for displaying.

Description

BRIEF DESCRIPTION OF THE DRAWING

[0048] The above and other objectives, features and advantages of the present invention will be more apparent from the following detailed description taken in conjunction with the accompanying drawings, in which:

[0049] FIG. 1 is a diagram of a hardware architecture of digital implementation of fractional calculus operation in prior art;

[0050] FIG. 2 is a diagram of a digital oscilloscope in prior art;

[0051] FIG. 3 is a diagram of a digital oscilloscope having fractional calculus operation and display function in accordance with the present invention;

[0052] FIG. 4 is a diagram of the fractional operation module shown in FIG. 3;

[0053] FIG. 5 is a waveform comparison graph of a acquired data and the corresponding fractional operation result data, where U1 is the oscillogram of the acquired data, U2 is the oscillogram of corresponding fractional operation result data.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

[0054] Hereinafter, preferred embodiments of the present invention will be described with reference to the accompanying drawings. It should be noted that the similar modules are designated by similar reference numerals although they are illustrated in different drawings. Also, in the following description, a detailed description of known functions and configurations incorporated herein will be omitted when it may obscure the subject matter of the present invention.

[0055] The fractional calculus operation in the present invention is realized according to the definition of G-L Equation. Let x(t) as an analog signal of continuous time, G-L fractional calculus operation is defined as follows:

[00003]$\begin{array}{cc}{D}^{\alpha }x\left(t\right)=\underset{h.fwdarw.0}{\lim }{h}^{-\alpha }\underset{j=0}{\overset{\infty }{.Math.}}{\left(-1\right)}^j\left(\frac{n}{j}\right)x\left(t-\mathrm{jh}\right)& \left(1\right)\end{array}$

[0056] where D stands for fractional operator, α is operation order, when α<0, the operation of fractional calculus is a fractional integration operation, when α>0, the operation of fractional calculus is a fractional differentiation operation, h is operation step, (n/h) satisfies the following condition:

[00004]$\left(\frac{n}{j}\right)=C_n^j=\frac{n!}{j!\left(n-j\right)!}$

[0057] So, when t.sub.1−t.sub.0=+∞, equation (1) needs to perform infinite calculations, that is to say, large amount of data are needed to store, which will consume large amount of storage resources. But the storage resources generally are limited, and large amount of data calculations need more operation time. Therefore, fractional calculus operation can be performed only to the signal of finite length. According to finite length calculation principle, equation (1) can be reformulated as follows:

[00005]$\begin{array}{cc}{ }_{t_0}{D}_{t_1}^{\alpha }x\left(t\right)=\underset{h.fwdarw.0}{\lim }{h}^{-\alpha }\underset{j=0}{\overset{\frac{t_1-t_0}{h}}{.Math.}}{\left(-1\right)}^j\left(\frac{n}{j}\right)x\left(t-\mathrm{jh}\right)& \left(2\right)\end{array}$

[0058] where t.sub.0 and t.sub.1 respectively are the starting time and the ending time. When fractional calculus operation is implemented digitally, the signal is discretized, equation (2) can be further reformulated as follows:

[00006]$\begin{array}{cc}\begin{array}{c}{ }_{t_0}^{}{D}_{t_1}^{\alpha }x\left[n\right]=\underset{h.fwdarw.0}{\lim }\frac{1}{h^{\alpha }}\underset{j=0}{\overset{.Math.\frac{t_1-t_0}{h}.Math.}{.Math.}}{b}_{j}x\left[n-\mathrm{jh}\right]\\ b_j=\left(1-\frac{1+\alpha }{j}\right)b_{j-1},b_0=1\end{array}& \left(3\right)\end{array}$

[0059] where ┌ ┐ stands for round up.

[0060] A digital oscilloscope has X divisions and L pixels at the horizontal direction of the display area, its time base is T.sub.base, so the sampling interval T.sub.S_LCD of waveform display is (X*T.sub.base)/L. In the present invention, combining with the display theory of digital oscilloscope, the sampling interval T.sub.S_LCD of waveform display is taken as the operation step h. To the displaying of a frame of waveform, there has:

[00007]$.Math.\frac{t_1-t_0}{h}.Math.\ge L$

[0061] When the operation order a and the time base T.sub.base of digital oscilloscope is fixed, the parameters b.sub.j and 1/h.sup.α are fixed, then there has:

[00008]$c\left(j\right)=\frac{b_j}{h^{\alpha }},j=0,1,2,3,.Math.,L-1$

[0062] Equation (3) can be rewritten as follows:

[00009]$\begin{array}{cc}\begin{array}{c}y\left[n\right]={.Math.}_{j=0}^{L-1}c\left[j\right]x\left[n-j\right]\\ \mathrm{where}c\left(j\right)=b_j{T_{S\_\mathrm{LCD}}}^{-\alpha }=\frac{b_j{L}^{\alpha }}{{\left({\mathrm{XT}}_{\mathrm{base}}\right)}^{\alpha }},j=0,1,2,.Math.,L-1.\end{array}& \left(4\right)\end{array}$

[0063] Then we can realize a fractional calculus operation of fixed L points in a digital oscilloscope by performing multiplication and accumulation to an acquired data.

[0064] FIG. 3 is a diagram of a digital oscilloscope having fractional calculus operation and display function in accordance with the present invention.

[0065] In one embodiment, as shown in FIG. 3, a digital oscilloscope having fractional calculus operation and display function in accordance with the present invention comprises an ADC module 1, a signal acquisition and storage module 2, a signal processing and display module 3, a LCD 4, a digital oscilloscope control module 5 and a fractional calculus module 6, where the fractional calculus module 6 is the innovation of the present invention.

[0066] An input analog signal (a signal being measured or a signal under measurement) x(t) is sampled and quantified by the ADC module 1 to obtain an sampled data ADC(m), m is the sampling point of the sampled data ADC(m). In the signal acquisition and storage module 2, the sampled data ADC(m) is decimated according to selected time base, and then under the control of the digital oscilloscope control module 5, the sampled data after decimation (acquired data) is stored into the acquisition memory ADC RAM 202 of the signal acquisition and storage module 2 in order based on storage address according to a set pre trigger depth and trigger signal. When the acquisition memory ADC RAM 202 is stored full, which means L acquired data is stored, the acquisition and storage is stopped, the oscilloscope data done signal DSO_DONE outputted by the acquisition memory ADC RAM is turned from 0 to 1. At this moment, the acquired data stored in the acquisition memory ADC RAM can be denoted by x(n), n is storage address, n=0, 1, 2, . . . , L−1. The signal acquisition and storage of digital oscilloscope belongs to the prior art and is not detailed here.

[0067] In the present invention, a fractional calculus module 6 is specially innovated, which is used for reading out the acquired data x(n) from the acquisition memory ADC RAM 202 to perform a fractional differentiation operation or a fractional integration operation and obtaining a fractional calculus result data y(n), which is sent to the signal processing and display module 3, in which the acquired data x(n) is turned into a display data and the display data is sent to the LCD 4 for displaying.

[0068] In the embodiment, as shown in FIG. 4, the fractional calculus module 6 comprises a fixed coefficient memory CO RAM 601, N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1 602, N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603, a fractional operation results memory FO RAM 604 and a fractional operation control module FO CTRL 605.

[0069] The fixed coefficient memory CO RAM 601 is used for storing the L fixed coefficients c(0), c(1), c(2), . . . , c(L−1) sent from the digital oscilloscope control module 5. To the j.sup.th fixed coefficient c(j), its value is:

[00010]$\begin{array}{cc}c\left(j\right)=\frac{b_j{L}^{\alpha }}{{\left({\mathrm{XT}}_{\mathrm{base}}\right)}^{\alpha }},j=0,1,2,.Math.,L-1& \left(5\right)\end{array}$

[0070] where L is the number of the pixels at the horizontal direction of the display area of the digital oscilloscope, X is the number of the divisions at the horizontal direction of the display area of the digital oscilloscope, T.sub.base is the time base of the digital oscilloscope, α is operation order, when α<0, the fractional calculus operation is a fractional integration operation, when α>0, the fractional calculus operation is a fractional differentiation operation, b.sub.j is a fixed parameter, its value is:

[00011]$\begin{array}{cc}{b}_j=\left(1-\frac{1+\alpha }{j}\right)b_{j-1},b_0=1& \left(6\right)\end{array}$

[0071] The output of the fixed coefficient memory CO RAM 601 is connected to the input of the 1.sup.st D flip-flop delay unit D.sub.1, from the 2.sup.nd D flip-flop delay unit D.sub.2 on, the input of a D flip-flop delay unit is connected to the output of its previous D flip-flop delay unit.

[0072] Each fractional operation unit comprises one multiplier and one accumulator (ACC). To the p.sup.th fractional operation unit THD.sub.p, p=1, 2, . . . , N−1, one (port A) of the two inputs of its multiplier MULTI.sub.p is connected to the acquired data x(n) outputted by the signal acquisition and storage module 2, another (port B) of the two inputs of its multiplier MULTI.sub.p is connected to the input of the p.sup.th D flip-flop delay unit D.sub.p, the input of its accumulator ACC.sub.p is connected to the output of its multiplier MULTI.sub.p, the output of the accumulator ACC.sub.p is taken as the output of the p.sup.th fractional operation unit THD.sub.p; to the N.sup.th fractional operation unit THD.sub.N, one of the two inputs of its multiplier MULTI.sub.N is connected to the acquired data x(n) outputted by the signal acquisition and storage module 2, another of the two inputs of its multiplier MULTI.sub.N is connected to the output of the N−1.sup.th D flip-flop delay unit D.sub.N-1, the input of its accumulator ACC.sub.N is connected to the output of its multiplier MULTI.sub.N.

[0073] The fractional operation results memory FO RAM 604 is used for storing the operation results of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N.

[0074] The fractional operation control module FO CTRL 605 is used for controlling the fractional operation. When the oscilloscope data done signal DSO_DONE is turned from 0 to 1 is detected, the fractional operation control module FO CTRL performs the following steps:

[0075] (1): Initialing

[0076] Initializing fractional operation number i to 1, the flag i_done of the i.sup.th fractional calculus operation to 0, fractional operation done signal FO_DONE to 0.

[0077] (2): Resetting

[0078] Resetting the initial values of the N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1 602 to 0, setting the read address ADDR.sub.adcram of the acquisition memory ADC RAM 202 to N*i−1, resetting the read address ADDR.sub.coram of the fixed coefficient memory CO RAM 601 to 0, resetting the initial values of the accumulators of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 to 0, resetting the write address ADDR.sub.foram of the fractional operation results memory FO RAM 602 to 0.

[0079] (3): Reading data to multiply and accumulate

[0080] Driven by a system clock clk, the fractional operation control module FO CTRL 605 initiates a operation of reading the acquisition memory ADC RAM 202 and a operation of reading the fixed coefficient memory CO RAM 601 at each system clock, the total number of reading is k, k=N*i; to each reading, the read address ADDR.sub.adcram is subtracted by 1, and the read address ADDR.sub.coram is added by 1.

[0081] The data outputted by the read data port of the acquisition memory ADC RAM 202 in chronological order are acquired data x(k−1), . . . , x(1), x(0), namely, the data sent to the port As of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 are acquired data x(k−1) at 0.sup.th system clock, the data sent to the port As of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 are acquired data x(k−2) at 1.sup.st system clock, the data sent to the port As of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 are acquired data x(k−2) at 2.sup.nd system clock, and so on, the data sent to the port As of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 are acquired data x(0) at (k−1).sup.th system clock.

[0082] The data outputted by the read data port of the fixed coefficient memory CO RAM 601 in chronological order are the fixed data c(0), c(1), . . . , c(k−1), the fixed data c(0), c(1), . . . , c(k−1) are serially sent to the N−1 D flip-flop delay units D.sub.1, D.sub.2, D.sub.3, . . . , D.sub.N-1 602, namely, at 0.sup.th system clock, the data sent to the port B of the multiplier of the 1.sup.st fractional operation unit THD.sub.1 is the fixed data c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, at 1.sup.st system clock, the data sent to the port Bs of the multipliers of the 1.sup.st, the 2.sup.nd fractional operation units THD.sub.1, THD.sub.2 are respectively the fixed data c(1), c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, at 2.sup.nd system clock, the data sent to the port Bs of the multipliers of the 1.sup.st, the 2.sup.nd and the 3.sup.rd fractional operation units THD.sub.1, THD.sub.2, THD.sub.2 are respectively the fixed data c(2), c(1), c(0), the data sent to the port Bs of the multipliers of the rest of the fractional operation units are 0, and so on, until at the (k−1).sup.th system clock, the data sent to the port Bs of the multipliers of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N are respectively c(k−1), . . . , c(2), c(1), c(0).

[0083] After each reading, all the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 perform a multiplying operation and an accumulating operation in turn, when the number of the accumulations reaches k, the i.sup.th fractional calculus operation is accomplished, the i.sup.th fractional operation done signal i_done is turned from 0 to 1;

[0084] (4): Storing the results of the i.sup.th fractional calculus operation

[0085] When the i.sup.th fractional operation done signal i_done is turned from 0 to 1 is detected, the outputs of the N fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, . . . , THD.sub.N 603 are taken as the results y(k−1)y(k−2), y(k−3), . . . , y(k−N) of the i.sup.th fractional calculus operation and merged into a merged data, the merged data are stored into the fractional operation results memory FO RAM 604 at address ADDR.sub.foram under the same system clock clk;

[0086] (5) Judging whether a frame of acquired data, namely L acquired data are calculated

[0087] Judging whether the current fractional operation number i is less than L/N, if yes, the write address ADDR.sub.foram is added by 1, the current fractional operation number i is added by 1, the i.sup.th fractional operation done signal i_done is set to 1, then returning to step (2); if not, the fractional operation done signal FO_DONE is set to 1, stopping the fractional calculus operation.

[0088] The results of the i.sup.th fractional calculus operation (Ni system clock) of the N fractional operation units are as follows:

TABLE-US-00001 TABLE 1 Fractional System clock operation unit clk.sub.0 clk.sub.1 . . . clk.sub.Ni − 2 clk.sub.Ni − 1 output THD.sub.1 c(0)*x(Ni − 1) c(1)*x(Ni − 2) . . . c(Ni − 2)*x(1) c(Ni − 1)*x(0) y(Ni − 1) THD.sub.2 0*x(Ni − 1) c(0)*x(Ni − 2) . . . c(Ni − 3)*x(1) c(Ni − 2)*x(0) y(Ni − 2) THD.sub.3 0*x(Ni − 1) 0*x(Ni − 2) . . . c(Ni − 4)*x(1) c(Ni − 3)*x(0) y(Ni − 3) . . . . . . . . . . . . . . . . . . . . . THD.sub.N 0*x(Ni − 1) 0*x(Ni − 2) . . . c(Ni − N + 1)*x(1) c(Ni − N)*x(0) y(Ni − N)

[0089] The digital oscilloscope control module 5 is detecting the fractional operation done signal FO_DONE, when the fractional operation done signal FO_DONE is turned from 0 to 1 is detected, the digital oscilloscope control module performs the following steps:

[0090] (1): Sending the acquired data x(n), n=1, 2, . . . , L−1 stored in the acquisition memory ADC RAM 202 and the fractional calculus result data y(n), n=1, 2, . . . , L−1 stored the fractional operation results memory FO RAM 604 to signal processing and display module 3, in which the acquired data x(n) and the fractional calculus result data y(n) are respectively converted into a display data through a drawing thread, and the two display data are sent to LCD 4 for displaying.

[0091] In the embodiment, after sending the acquired data x(n) and the fractional calculus result data y(n) to signal processing and display module 3, the vertical sensitivity corresponding to the acquired data x(n) is multiplied with the number of the vertical divisions of the display area to obtain a display range R1, the vertical sensitivity corresponding to the fractional calculus result data y(n) is multiplied with the number of the vertical divisions of the display area to obtain another display range R2. The range of the acquired data x(n) is [0,2.sup.v−1], where ν is the digit number of the ADC module, the data range needed by the drawing thread is:

[00012]$\left[-\frac{R_1}{2},\frac{R_1}{2}\right]$

[0092] Therefore, it is needed to perform a normalization of data displaying to obtain the data d.sub.x(n), d.sub.y(n) which can be recognized by the drawing thread. The normalization of the acquired data x(n) is:

[00013]$d_x\left(n\right)=\left(\frac{x\left(n\right)}{2^v-1}-0.5\right)\times R_1$

[0093] The normalization of the fractional calculus result data y(n) is:

[00014]$d_y\left(n\right)=\left(\frac{y\left(n\right)}{\left(2^v-1\right)\times R_2/R_1}-0.5\right)\times R_2$

[0094] The data d.sub.x(n) and d.sub.y(n) are sent to a drawing thread, through which the data d.sub.x(n) and d.sub.y(n) are respectively converted into a display data, and the two display data are sent to LCD 4 for displaying.

[0095] (2): Resetting the oscilloscope data done signal DSO_DONE and the fractional operation done signal FO_DONE to 0 to enable the signal acquisition and storage module and the fractional calculus module.

[0096] The digital oscilloscope having fractional calculus operation and display function in accordance with present invention operates as follows:

[0097] Step S1: entering the display interface of fractional calculus operation through the interface operation of user;

[0098] Step S2: setting the time base T.sub.base and the operation order α of the digital oscilloscope though the interface operation of user;

[0099] Step S3: the digital oscilloscope control module obtains the time base T.sub.base and the operation order α, and calculates the L fixed coefficients c(1), c(1), c(2), . . . , c(L−1) according to equation (5);

[0100] Step S4: storing the L fixed coefficients c(0), c(1), c(2), . . . , c(L−1) into the fixed coefficient memory CO RAM 601;

[0101] Step S5: obtaining the instructions of acquisition (including but not limited to trigger level, input sensitivity, pretrigger depth, coupling mode, trigger channel, time base) through the interface operation of user, and generating corresponding command to set the digital oscilloscope, setting the oscilloscope data done signal DSO_DONE and the fractional operation done signal FO_DONE to 0;

[0102] Step S6: enabling signal acquisition function and fractional calculus operation function, namely enabling the signal acquisition and storage module and the fractional calculus module;

[0103] Step S7: input analog signal (the signal being measured or the signal under measurement) x(t) is sampled and quantified by the ADC module 1 to obtain sampled data ADC(m), m is sampling point; in the signal acquisition and storage module 2, the sampled data ADC(m) is split into two sampled data, one sampled data is sent to digital comparator module 203, and a trigger signal is generated through trigger signal generating module according to the output of digital comparator module 203, the other sampled data is sent to data flow control module 201 in which the sampled data is decimated according to a decimation coefficient D, the sampled data after decimation (acquired data) is sent to acquisition memory ADC RAM 202, where acquisition memory ADC RAM 202 is a dual port RAM, the decimation coefficient D is determined by selected time base. Before arriving at the pre trigger depth, digital oscilloscope control module 5 only performs write operation on acquisition memory ADC RAM 202; after arriving at the pre trigger depth and before generating trigger signal trig, digital oscilloscope control module 5 performs write operation and read operation simultaneously, and the speeds of write and read operations are the same; after generating trigger signal trig, digital oscilloscope control module 5 only performs write operation on acquisition memory ADC RAM 202, when acquisition memory ADC RAM 202 is stored full, which means a data acquisition is complete, oscilloscope data done signal DSO_DONE is turned from 0 to 1.

[0104] Step S8: in the fractional calculus module 6, when the oscilloscope data done signal DSO_DONE is turned from 0 to 1 is detected by the fractional operation control module FO CTRL 605, the fractional operation control module FO CTRL performs initialing, resetting, reading data to multiply and accumulate, storing the results of the i.sup.th fractional operation, judging and returning, until a frame of acquired data, namely L acquired data are calculated.

[0105] Step S9: in digital oscilloscope control module 5, when the fractional operation done signal FO_DONE is turned from 0 to 1 is detected, it sends the fractional calculus result data y(n), n=1, 2, . . . , L−1 to signal processing and display module 3, in which the fractional calculus result data y(n) is converted into a display data through a drawing thread, and the display data are sent to LCD 4 for displaying. More details are as follows:

[0106] After sending the acquired data x(n) and the fractional calculus result data y(n) to signal processing and display module 3, the vertical sensitivity corresponding to the acquired data x(n) is multiplied with the number of the vertical divisions of the display area to obtain a display range R1, the vertical sensitivity corresponding to the fractional calculus result data yn) is multiplied with the number of the vertical divisions of the display area to obtain another display range R2. Then the acquired data x(n) is normalized as data d.sub.x(n):

[00015]$d_x\left(n\right)=\left(\frac{x\left(n\right)}{2^v-1}-0.5\right)\times R_1$

[0107] The fractional calculus result data y(n) is normalized as data d.sub.y(n):

[00016]$d_y\left(n\right)=\left(\frac{y\left(n\right)}{\left(2^v-1\right)\times R_2/R_1}-0.5\right)\times R_2$

[0108] Lastly, the data d.sub.x(n) and d.sub.y(n) are sent to a drawing thread, through which the data d.sub.x(n) and d.sub.y(n) are respectively converted into a display data, and the two display data are sent to LCD 4 for displaying.

[0109] And resetting the oscilloscope data done signal DSO_DONE and the fractional operation done signal FO_DONE to 0 to enable the signal acquisition and storage module and the fractional calculus module.

Example

[0110] In the present example, the input analog signal is a sine wave with the frequency of 500 KHz and the amplitude of 1V, the number of the pixels at the horizontal direction of the display area of the digital oscilloscope is 1000 (namely L=1000), the number of the divisions at the horizontal direction of the display area of the digital oscilloscope is 10 (namely X=10), the number of the pixels at the vertical direction of the display area of the digital oscilloscope is 800 (namely L=1000), the number of the divisions at the vertical direction of the display area of the digital oscilloscope is 8. In a specific implementation, the data acquisition and the fractional calculus operation are realized by a FPGA (model: XC7K160T-2FFG676I), the model of the ADC is AD9628, the number of the channels of the ADC is 2, the sampling rate of the ADC is configured at 100 MSPS through the FPGA.

[0111] In the present example, the time base T.sub.base=1us/div, and the operation order α=0.5, namely the operation of fractional calculus operation is a fractional differentiation operation.

[0112] In the present example, the number of the fractional operation units is 4, namely N=4. The running states of the first round of operations of the 4 fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, THD.sub.4 are shown in Table 2:

TABLE-US-00002 TABLE 2 Fractional operation System clock unit clk.sub.0 clk.sub.1 clk.sub.2 clk.sub.3 output THD.sub.1 c(0)*x(3) c(1)*x(2) c(2)*x(1) c(3)*x(0) y(3) THD.sub.2 0*x(3) c(0)*x(2) c(1)*x(1) c(2)*x(0) y(2) THD.sub.3 0*x(3) 0*x(2) c(0)*x(1) c(1)*x(0) y(1) THD.sub.4 0*x(3) 0*x(2) 0*x(1) c(0)*x(0) y(0)

[0113] The running states of the second round of operations of the 4 fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, THD.sub.4 are shown in Table 3:

TABLE-US-00003 TABLE 3 Fractional System clock operation unit clk.sub.0 clk.sub.1 . . . clk.sub.6 clk.sub.7 output THD.sub.1 c(0)*x(7) c(1)*x(6) . . . c(6)*x(1) c(7)*x(0) y(7) THD.sub.2 0*x(7) c(0)*x(6) . . . c(5)*x(1) c(6)*x(0) y(6) THD.sub.3 0*x(7) 0*x(6) . . . c(4)*x(1) c(5)*x(0) y(5) THD.sub.4 0*x(7) 0*x(6) . . . c(3)*x(1) c(4)*x(0) y(4)

[0114] The running states of the 250.sup.th round of operations of the 4 fractional operation units THD.sub.1, THD.sub.2, THD.sub.3, THD.sub.4 are shown in Table 4:

TABLE-US-00004 TABLE 4 Fractional System clock operation unit clk.sub.0 clk.sub.1 . . . clk.sub.998 clk.sub.999 output THD.sub.1 c(0)*x(999) c(1)*x(998) . . . c(998)*x(1) c(999)*x(0) y(999) THD.sub.2 0*x(999) c(0)*x(998) . . . c(997)*x(1) c(998)*x(0) y(998) THD.sub.3 0*x(999) 0*x(998) . . . c(996)*x(1) c(997)*x(0) y(997) THD.sub.4 0*x(999) 0*x(998) . . . c(995)*x(1) c(996)*x(0) y(996)

[0115] In the present example, after a fractional calculus operation is completed, the waveforms of the acquired data x(n) and the fractional calculus result data y(n) are shown in FIG. 5, where U1 is the waveform of the acquired data x(n), its vertical sensitivity is 0.5V/div, vertical range is 5V, U2 is the waveform of the fractional calculus result data y(n), its vertical sensitivity is 0.2 mV/div, vertical range is 2 mV. As can be seen from FIG. 5, for the input signal is a sine waveform, the fractional calculus result is also a sine waveform, and its phase is advanced than that of the input signal, which shows that the fractional calculus operation is realized.

[0116] While illustrative embodiments of the invention have been described above, it is, of course, understand that various modifications will be apparent to those of ordinary skill in the art. Such modifications are within the spirit and scope of the invention, which is limited and defined only by the appended claims.