CN105659938B - Based on the generalization observing and controlling channel simulation method of the outer discharge technique of high accuracy time delay - Google Patents

Abstract

The invention provides a kind of generalization observing and controlling channel simulation method based on the outer discharge technique of high accuracy time delay, its basic implementation process is as follows: utilize local oscillation signal to input signal S_up(t) carry out quadrature frequency conversion, obtain complex baseband input signal S_upB(t); Sliding-model control is carried out to star ground time delays τ (t); With T_sFor the sampling interval is to complex baseband input signal S_upB(t) sample, obtain the sample sequence of complex baseband input signal, then initial complete cycle of α is carried out to it_inTime delay obtain S_upB1(n); From S_upB1(n) corresponding data in, is selected to carry out filtering to decimal filtering wave by prolonging time device; By sample sequence with T_sFor D/A conversion is carried out at interval, and after reconstruction filtering, obtain the Delay reconstruction signal S of complex baseband signal_uprB(t); According to star ground time delays τ (t), accurately produce signal S_upD(t), and by its orthogonal modulation arriveUpper, to be fixed the Delay reconstruction signal S of local oscillation signal_upL1(t); By the Delay reconstruction signal S of complex baseband signal_uprBAnd the Delay reconstruction signal S of fixing local oscillation signal (t)_upL1(t) following formula is pressed synthetic, to obtain final analog output signal S_upr(t)。<pb pnum="1" />

Description

Based on the generalization observing and controlling channel simulation method of the outer discharge technique of high accuracy time delay

Technical field

The invention belongs to the technical field such as channel simulation, space flight measurement and control, be specifically related to one and prolong based on high accuracyTime outer discharge technique generalization observing and controlling channel simulation method.

Background technology

Along with the develop rapidly of space telemetry and control technology, a large amount of new TT & C architecture are presented and are applied. From existingHave technology, for existing various TT & C architecture, the single machine test method of answering machine, earth station thanMore ripe. But, be not also applicable to the generalization observing and controlling channel simulation equipment of any TT & C architecture at present,Can insert between TT&C Earth Station and TT&C Transponder, for verifying that TT&C system is preferring comprehensively, fullyFunction and dynamic property under operative scenario. Therefore, in the urgent need to finding a kind of unitized observing and controlling channel mouldPlan technology, can when the prior information the unknowns such as TT & C architecture, signal form, signal parameter, accuratelyThe dynamic transmission effect that simulation uplink and downlink signal causes when transmission.

Generalization observing and controlling channel simulation technology, needs to solve the large model of broadband, random waveform measurement and control signal in essenceEnclose, high accuracy, high speed dynamic transmission time delay problem of modelling. For above-mentioned dynamic transmission time delay simulation demand,Conventional method mainly comprises regeneration time delay forwarding method, postpones collimation method, digital RF storage method etc. at present. AgainRaw time delay forwards method according to prior informations such as signal system, signal form, signal parameters, first manages in simulationDevice this locality is born corresponding intermediate frequency or base band measurement and control signal again; Subsequently according to the dynamic scene of need simulations, againSuperposed simulation information on raw signal; Finally, produced final analog output signal by operations such as up-conversions,To realize the simulation of time delay and Doppler frequency shift. Postponing collimation method utilizes microwave numerical-control delay line, microwave optical fiber to prolongLate line, SAW delay line etc. directly carry out time delay to input signal, by the different sections of pilot delay lineCombination realizes the simulation of different delayed time, and the switch speed that controls different delayed time section is realized the mould of change of distanceIntend. Digital RF storage method is sampled to input signal and is stored, and after one section of time delay and conversion, sends out againBe shot out, write the simulation that realizes initial distance with the time reading by control data, write by controllingThe speed entering and read realizes the simulation of change of distance.

But all there is defect separately in said method:

1), when regeneration time delay forwards method and carries out signal regeneration and analog information stack in simulator this locality, all needThe prior informations such as known signal system, signal form, signal parameter. Thus the method is applicable to single moreDynamic transmission time delay simulation under system, single mode of operation, is difficult to realize generalization truly.

2) the time delay simulation context of delay achieved by collimation method directly depends on length, the simulation essence of delay lineDegree depends on the precision of delay line, the switch speed that delay number is dynamically depended in simulation. Therefore, the method difficultyRealizing on a large scale, high accuracy, at a high speed dynamic time delay simulation;

3) the time delay simulation precision achieved by radio frequency storage method depends primarily on the speed of read-write clock, forRealize the simulation of high-precision time delay, the read-write clock of high speed must be possessed. Therefore, high-precision at someThe application scenario of degree, due to the restriction of high clock rate, often cannot meet engineering demand.

Summary of the invention

The object of the invention is the defect in order to overcome prior art, in order to solve TT & C architecture, signal shapeIn the situation of the prior information such as formula, signal parameter the unknown, realize on a large scale, high accuracy, pass dynamically at a high speedDefeated time delay problem of modelling, proposes a kind of generalization observing and controlling channel simulation based on the outer discharge technique of high accuracy time delayMethod.

The inventive method is achieved through the following technical solutions:

Based on a generalization observing and controlling channel simulation method for the outer discharge technique of high accuracy time delay, it was implemented substantiallyJourney is as follows:

Step 1: utilize frequency for f_upLLocal oscillation signal to input signal S_up(t) carry out quadrature frequency conversion, obtainComplex baseband input signal S_upB(t)；

Step 2: to star ground time delays τ (t) with T_sCarry out sliding-model control for interval, obtain the star of discretizationGround time delays τ (n), and n=0,1 ..., N_s, according to formula (4);

$\begin{array}{c}\mathrm{\&tau;}\left(n\right)=\mathrm{\&tau;}\left(0\right)+{\mathrm{\&tau;}}_v\left(n\right)\\ =\&lsqb;{\mathrm{\&alpha;}}_{in}+{\mathrm{\&alpha;}}_d+{\mathrm{\&alpha;}}_v\left(n\right)\&rsqb;T_s\\ =\&lsqb;{\mathrm{\&alpha;}}_{in}+{\mathrm{\&alpha;}}_{vd}\left(n\right)\&rsqb;T_s\end{array}---\left(4\right)$

Wherein, τ (0) represents initial time delay, τ_v(n) time delay changing, α is represented_inFor the number of cycles of initial time delayNumber, α_dFor the Fractional number of initial time delay, α_v(n) for changing the number of cycles of time delay,α_vd(n)＝α_d+α_v(n)。

Step 3: with T_sFor the sampling interval is to complex baseband input signal S_upB(t) sample, obtain complex base band defeatedEnter the sample sequence { S of signal_upB(n)，n＝0，1，...，N_sAnd storage, then to S_upB(n) initial complete cycle of α is carried out_inTime delay obtain S_upB1(n)；

Step 4: according to α_vd(n), from S_upB1(n) corresponding data in, is selected to carry out filtering to decimal filtering wave by prolonging time device,S is obtained to realize the time delay of initial Fractional time delay and period of change_uprB(n)；

Step 5: by { S_uprB(n)，n＝0，1，...，N_sWith T_sFor D/A conversion is carried out at interval, and after reconstruction filteringObtain the Delay reconstruction signal S of complex baseband signal_uprB(t)；

Step 6: in carry out step 3, according to star ground time delays τ (t), accurately produces signal S_upD(t)，And its orthogonal modulation is arrived(t) upper, to be fixed the Delay reconstruction signal S of local oscillation signal_upL1(t)

$\begin{array}{c}{S}_{upD}\left(t\right)=\exp \&lsqb;-j2{\mathrm{\&pi;f}}_{upL}\mathrm{\&tau;}\left(t\right)\&rsqb;\\ S_{upL1}\left(t\right)=S_{upD}\left(t\right)\&times;S_{upL}^{*}\left(t\right)\\ =\exp \{j2{\mathrm{\&pi;f}}_{upL}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\}\end{array}---\left(18\right)$

Wherein,For local oscillation signal S_upL(t) complex conjugate;

Step 7: by the Delay reconstruction signal S of complex baseband signal_uprBAnd the Delay reconstruction letter of fixing local oscillation signal (t)Number S_upL1(t) following formula is pressed synthetic, to obtain final analog output signal S_upr(t)。

Further, the detailed process of step 4 of the present invention is:

Step by step 1: calculate one group of parallel fixed coefficient c according to (11) formula and (12) formula_m(n)；

$h_2\left(n\right)=\underset{k\&NotEqual;n}{\overset{N}{\underset{k=0}{\mathrm{\&Pi;}}}}\frac{d-k}{n-k},n=0,1,\mathrm{2...}N---\left(11\right)$

$h_2\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}{c}_m\left(n\right)d^m,n=0,1,\mathrm{2...}N---\left(12\right)$

Step by step 2: by α_vd(n) integer part and fractional part is decomposed into, namely

Step by step 3: by c_m(n), d (n) and D_in(n) formula (13) is brought into, the S calculating_uprB(n), namely

$S_{uprB}\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}\&lsqb;\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}{c}_m\left(k\right)S_{upB1}(n-D_{in}(n)-k)\&rsqb;d^m\left(n\right)---\left(15\right).$

Beneficial effect:

The first, the generalization observing and controlling channel simulation method that the present invention is based on the outer discharge technique of high accuracy time delay is not paid close attention toThe prior informations such as TT & C architecture, signal form, signal parameter, by the signal profit to broadband, random waveformCarrying out direct time delay by the method that initial complete cycle counts and decimal filtering wave by prolonging time device combines forwards to realize dynamicThe simulation of state transmission delay. Thus the method is a kind of generalization analogy method truly.

The second, the maximum delay simulation context achieved by the present invention, depends primarily on the storage of sampled dataThe degree of depth; Achieved time delay precision, depends primarily on the time delay essence of counting precision, decimal filtering wave by prolonging time deviceDegree; Achieved time delay simulation dynamically, depends primarily on the renewal speed of decimal filtering wave by prolonging time device delay parameterThe generation precision of degree and local oscillator Delay reconstruction signal. Thus, only need to select jumbo memory device, high-precisionThe decimal filtering wave by prolonging time device, high-order DDS technology etc. of degree, can realize on a large scale, high accuracy, at a high speed dynamicallyTime delay simulation.

Three, the present invention does not pay close attention to the prior informations such as TT & C architecture, signal form, signal parameter, can realizeOn a large scale, high accuracy, dynamic transmission delay simulation at a high speed.

Accompanying drawing explanation

Fig. 1 is system block diagram of the present invention.

Detailed description of the invention

Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.

Generalization observing and controlling channel simulation method based on the outer discharge technique of high accuracy time delay may be summarized as follows: first,Utilize fixing local oscillator to carry out quadrature demodulation to the radio-frequency input signals of broadband, random waveform, obtain broadband, appointThe complex baseband signal of meaning waveform; Subsequently, complex baseband signal is sampled, and according to star ground time delays τ (t)Adopt to count the method combining with decimal filtering wave by prolonging time device initial complete cycle discrete complex baseband signal is carried outHigh accuracy time delay, by the Delay reconstruction signal that obtains complex baseband signal after D/A conversion reconstruction filtering; WithTime, according to star ground time delays τ (t), produce the Delay reconstruction signal of fixing local oscillation signal; Finally, by complex radicalWith the Delay reconstruction signal Delay reconstruction signal syntheses corresponding with fixing local oscillation signal, obtain final simulation defeatedGo out signal.

Consider institute of the present invention extracting method, its up-downgoing principle of simulation is identical, therefore, below above behavior exampleThe implementation process of the method is described.

Transmit as S if earth station is up_upAnd suppose that it is signal system, signal form, signal parameter (t)Etc. the broadband rf signal of prior information the unknown, t star ground time delays is τ (t), then according to real thingReason process, arrives the up output signal S of answering machine_upr(t) according to formula (1);

S_upr(t)＝S_up[t-τ(t)](1)

Institute of the present invention extracting method target is just the S that transmits of satellite receiver_up(t), prolong according to the star ground timeTime τ (t), accurate analog produces and arrives the up input signal S of answering machine_upr(t)。

The inventive method is achieved through the following technical solutions:

Based on the generalization channel simulation method of the outer discharge technique of high accuracy time delay, its basic implementation process is as follows:When the transmission delay to up channel is simulated, now S_up(t) up input signal is represented, when to descendingWhen the transmission delay of channel is simulated, now S_up(t) descending input signal is represented.

Step 1: utilize frequency for f_upLLocal oscillation signal to input signal S_up(t) carry out quadrature frequency conversion, obtainComplex baseband input signal S_upB(t)；

S_upL(t)＝exp(-j2πf_upLt)(2)

S_upB(t)＝S_up(t)×exp(-j2π_upLt)(3)

Step 2: to star ground time delays τ (t) with T_sCarry out sliding-model control for interval, obtain the star of discretizationGround time delays τ (n), and n=0,1 ..., N_s, and decompose by formula (4);

$\begin{array}{c}\mathrm{\&tau;}\left(n\right)=\mathrm{\&tau;}\left(0\right)+{\mathrm{\&tau;}}_v\left(n\right)\\ =\&lsqb;{\mathrm{\&alpha;}}_{in}+{\mathrm{\&alpha;}}_d+{\mathrm{\&alpha;}}_v\left(n\right)\&rsqb;T_s\\ =\&lsqb;{\mathrm{\&alpha;}}_{in}+{\mathrm{\&alpha;}}_{vd}\left(n\right)\&rsqb;T_s\end{array}---\left(4\right)$

Wherein, τ (0) represents initial time delay, τ_v(n) time delay changing, α is represented_inFor the number of cycles of initial time delayNumber, α_dFor the Fractional number of initial time delay, α_v(n) for changing the number of cycles of time delay. α_vd(n) be initialFractional time delay number and the number of cycles sum that changes time delay, the i.e. α of time delay_vd(n)＝α_d+α_v(n)。

Wherein, intRepresent downward floor operation.

Step 3: with T_sFor the sampling interval is to complex baseband input signal S_upB(t) sample, obtain complex base band defeatedEnter the sample sequence { S of signal_upB(n)，n＝0，1，...，N_sAnd storage. Meanwhile, start to count from storing the moment,Every storage sampled point, count value adds 1, to be counted to α_inThe sample sequence of storage exports by time, with logicalCross storage to realize S_upB(n) initial complete cycle of α is carried out_inTime delay, obtain S_upB1(n)；

S_upB1(n)＝S_upB(n-α_in)(6)

Step 4: according to α_vd(n), from S_upB1(n) corresponding data in, is selected to carry out filtering to decimal filtering wave by prolonging time device,To realize the time delay of initial Fractional time delay and period of change, obtain { S_uprB(n)，n＝0，1，...，N_s}

$\begin{array}{c}{S}_{uprB}\left(n\right)=S_{upB1}\&lsqb;n-{\mathrm{\&alpha;}}_{vd}\left(n\right)\&rsqb;\\ =S_{upB}\&lsqb;n-{\mathrm{\&alpha;}}_{in}-{\mathrm{\&alpha;}}_{vd}\left(n\right)\&rsqb;\end{array}---\left(7\right)$

This step principle is as follows:

For easy analysis, first suppose α_vd(n)=D, D is fixing non-integer-period time delay. For realizing oneFixing non-integer-period time delay, can by S_upB1(n) realized by having the delayed time system of the shown frequency response of (8) formula

H_id(e^jω)＝exp(-jωD)(8)

If non-integer-period time delay D is decomposed into integer part D further_inWith fractional part d, namely

Then (8) formula can be equivalent to the cascade of following two delayed time systems

$\begin{array}{c}{H}_{id}\left(e^{j\mathrm{\&omega;}}\right)=H_1\left(e^{j\mathrm{\&omega;}}\right)\&times;H_2\left(e^{j\mathrm{\&omega;}}\right)\\ =\exp (-{\mathrm{j\&omega;D}}_{in})\&times;\exp (-j\mathrm{\&omega;}d)\end{array}---\left(10\right)$

Wherein,H can be realized by time delay complete cycle₂(e^jω)＝e^-jωdCorresponding decimal can be designedFiltering wave by prolonging time device is realized.

The present invention utilizes Lagrange type decimal filtering wave by prolonging time to approach delayed time system H₂(e^jω), its impact ringsAnswer h₂(n) following form is there is

$h_2\left(n\right)=\underset{k\&NotEqual;n}{\overset{N}{\underset{k=0}{\mathrm{\&Pi;}}}}\frac{d-k}{n-k},n=0,1,\mathrm{2...}N---\left(11\right)$

Wherein, N represents the length of Lagrange type decimal filtering wave by prolonging time device;

Lagrange type decimal filtering wave by prolonging time device is e at low frequency place^-jωdAccurately approximate. Therefore, for broadbandLimited baseband signal, Lagrange type decimal filtering wave by prolonging time device is a good tool that realizes decimal time delay.

But the shock response of Lagrange type decimal filtering wave by prolonging time device is the function of Fractional time delay d.Like this, for the decimal time delay of continuous variable, need calculate according to different Fractional time delays correspondingShock response. For the ease of dynamically updating of Fractional time delay, the present invention utilizes the setting P of decimal time delay dEach coefficient of rank approximation by polynomi-als, namely

$h_2\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}{c}_m\left(n\right)d^m,n=0,1,\mathrm{2...}N---\left(12\right)$

Wherein, c_m(n) be polynomial Coefficients of Approximation.

Like this, S_uprB(n) can by S_upB1(n) acquisition is calculated as follows

$\begin{array}{c}{S}_{uprB}\left(n\right)=S_{upB1}(n-D)\\ =S_{upB1}(n-D_{in})*h_2\left(n\right)\\ =\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}{h}_2\left(k\right)S_{upB1}(n-D_{in}-k)\\ =\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}\&lsqb;\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}{c}_m\left(k\right)d^m\&rsqb;S_{upB1}(n-D_{in}-k)\\ =\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}\&lsqb;\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}{c}_m\left(k\right)S_{upB1}(n-D_{in}-k)\&rsqb;d^m\end{array}---\left(13\right)$

Wherein, * represents convolution algorithm;

Known according to the form of (13) formula, for S_upB1(n) non-integer time delay D is carried out, need according to D'sInteger part D_inFrom S_upB1(n) in, select data and passed through one group of that walk abreast, fixed coefficient c_m(n) FIRWave filter, the fractional part d of the output recycling D of wave filter is weighted summation.

Like this, this step can by implementing step by step as follows:

Step by step 1: calculate one group of parallel fixed coefficient c according to (11) formula and (12) formula_m(n)；

Step by step 2: by α_vd(n) integer part and fractional part is decomposed into, namely

Wherein, D_in(n) at S_upB1(n) in, select for calculating required N+1 point data, d (n) is for when calculatingBe weighted.

Step by step 3: by described d (n) and D_in(n) formula (13) is brought into, the S calculating_uprB(n), namely

$S_{uprB}\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}\&lsqb;\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}{c}_m\left(k\right)S_{upB1}(n-D_{in}(n)-k)\&rsqb;d^m\left(n\right)---\left(15\right)$

Repeat step by step 2-step by step 3, to obtain n=0,1 ..., N_s{ the S in each moment_uprB(n)，n＝0，1，...，N_s}。

Step 5: by { S_uprB(n)，n＝0，1，...，N_sWith T_sFor D/A conversion is carried out at interval, and after reconstruction filteringObtain the Delay reconstruction signal S of complex baseband signal_uprB(t)；

$\begin{array}{c}{S}_{uprB}\left(t\right)=S_{upB}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\\ =S_{up}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\exp \{-j2{\mathrm{\&pi;f}}_{upL}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\}\end{array}---\left(16\right)$

Step 6: in carry out step 3, according to star ground time delays τ (t), accurately produces signal S_upD(t)

S_upD(t)＝exp[-j2πf_upLτ(t)](17)

And its orthogonal modulation is arrivedUpper, to be fixed the Delay reconstruction signal S of local oscillation signal_upL1(t)

$\begin{array}{c}{S}_{upL1}\left(t\right)=S_{upD}\left(t\right)\&times;S_{upL}^{*}\left(t\right)\\ =\exp \{j2{\mathrm{\&pi;f}}_{upL}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\}\end{array}---\left(18\right)$

Wherein,For local oscillation signal S_upL(t) complex conjugate.

Step 7: by the Delay reconstruction signal S of complex baseband signal_uprBAnd the Delay reconstruction letter of fixing local oscillation signal (t)Number S_upL1(t) following formula is pressed synthetic, to obtain final analog output signal S_upr(t)

$\begin{array}{c}{S}_{upr}\left(t\right)=S_{uprB}\left(t\right)\&times;S_{upL1}\left(t\right)\\ =S_{up}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\end{array}---\left(19\right)$

Through step 1-step 7, the present invention is in prior informations such as TT & C architecture, signal form, signal parametersIn unknown situation, (simulation of down channel is former to achieve transmission delay simulation to TT&C system up channelManage identical).

In sum, these are only preferred embodiment of the present invention, be not intended to limit guarantor of the present inventionProtect scope. Within the spirit and principles in the present invention all, any amendment of doing, be equal to replacement, improvement etc.,Within protection scope of the present invention all should be included in.

Claims (2)

1. the generalization observing and controlling channel simulation method based on the outer discharge technique of high accuracy time delay, is characterized in that,Its basic implementation process is as follows:

Step 1: utilize frequency for f_upLLocal oscillation signal to input signal S_up(t) carry out quadrature frequency conversion, obtainComplex baseband input signal S_upB(t)；

Step 2: to star ground time delays τ (t) with T_sCarry out sliding-model control for interval, obtain the star of discretizationGround time delays τ (n), and n=0,1 ..., N_s}；

τ(n)＝τ(0)+τ_v(n)

＝[α_in+α_d+α_v(n)]T_s

＝[α_in+α_vd(n)]T_s

Wherein, τ (0) represents initial time delay, τ_v(n) time delay changing, α is represented_inFor the number of cycles of initial time delayNumber, α_dFor the Fractional number of initial time delay, α_v(n) for changing the number of cycles of time delay,α_vd(n)＝α_d+α_v(n)；

Step 3: with T_sFor the sampling interval is to complex baseband input signal S_upB(t) sample, obtain complex base band defeatedEnter the sample sequence { S of signal_upB(n)，n＝0，1，...，N_s, then to S_upB(n) initial complete cycle of α is carried out_inTime delayObtain S_upB1(n)；

Step 4: according to α_vd(n), from S_upB1(n) corresponding data in, is selected to carry out filtering to decimal filtering wave by prolonging time device,To realize the time delay of initial Fractional time delay and period of change, obtain S_uprB(n)；

Step 5: by { S_uprB(n)，n＝0，1，...，N_sWith T_sFor D/A conversion is carried out at interval, and after reconstruction filteringObtain the Delay reconstruction signal S of complex baseband signal_uprB(t)；

Step 6: in carry out step 3, according to star ground time delays τ (t), accurately produces signal S_upD(t)，And its orthogonal modulation is arrivedUpper, to be fixed the Delay reconstruction signal S of local oscillation signal_upL1(t)；

S_upD(t)＝exp[-j2πf_upLτ(t)]

$\begin{array}{c}{S}_{upL1}\left(t\right)=S_{upD}\left(t\right)\&times;S_{upL}^{*}\left(t\right)\\ =\exp \{j2{\mathrm{\&pi;f}}_{upL}\&lsqb;t-\mathrm{\&tau;}\left(t\right)\&rsqb;\}\end{array}---\left(18\right)$

Wherein,For local oscillation signal S_upL(t) complex conjugate;

Step 7: by the Delay reconstruction signal S of complex baseband signal_uprBAnd the Delay reconstruction letter of fixing local oscillation signal (t)Number S_upL1(t) following formula is pressed synthetic, to obtain final analog output signal S_upr(t)；

S_upr(t)＝S_uprB(t)×S_upL1(t)

＝S_up[t-τ(t)]。

2. according to claim 1 based on the high accuracy time delay generalization observing and controlling channel simulation side of discharge technique outwardMethod, is characterized in that, the detailed process of described step 4 is:

Step by step 1: calculate one group of parallel fixed coefficient c according to (11) formula and (12) formula_m(n)；

$h_2\left(n\right)=\underset{k\&NotEqual;n}{\overset{N}{\underset{k=0}{\mathrm{\&Pi;}}}}\frac{d-k}{n-k},n=0,1,\mathrm{2...}N---\left(11\right)$

$h_2\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}{c}_m\left(n\right)d^m,n=0,1,\mathrm{2...}N---\left(12\right)$

N represents the length of Lagrange type decimal filtering wave by prolonging time device; D represents Fractional time delay, and P represents multinomialFormula matching exponent number;

Step by step 2: by α_vd(n) integer part and fractional part is decomposed into, namely

Step by step 3: by c_in(n), d (n) and D_in(n) formula (13) is brought into, the S calculating_uprB(n), namely

$S_{uprB}\left(n\right)=\underset{m=0}{\overset{P}{\mathrm{\&Sigma;}}}\&lsqb;\underset{k=0}{\overset{N}{\mathrm{\&Sigma;}}}{c}_m\left(k\right)S_{upB1}(n-D_{in}(n)-k)\&rsqb;d^m\left(n\right)---\left(13\right).$