CN103226637B - Optimizing method for state space of switch current wavelet filter circuit - Google Patents
Optimizing method for state space of switch current wavelet filter circuit Download PDFInfo
- Publication number
- CN103226637B CN103226637B CN201310146139.2A CN201310146139A CN103226637B CN 103226637 B CN103226637 B CN 103226637B CN 201310146139 A CN201310146139 A CN 201310146139A CN 103226637 B CN103226637 B CN 103226637B
- Authority
- CN
- China
- Prior art keywords
- cascade
- sub
- wavelet filter
- impulse response
- peak value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The invention provides an optimizing method for the state space of a switch current wavelet filter circuit. The optimizing method comprises the following steps: (1) decomposing the analog wavelet filter transfer function into a normalized cascade; (2) evaluating the impulse response peak of each subsystem of the normalizated cascade, and computing the optimal coefficients; (3) optimizing analog wavelet filter transfer function: each subsystem of the normalizated cascade multiplies the optimize coefficient, so as to optimizate the analoge wavelet filter state space, and enable size of the impulse response peak of the optimized subsystem of the normalizated cascade to be in consistency, and enable the dynamic range of the switch current circuit of each subsystem of the normalizated cascade to be the same. The method provided by the utility model is capable of unifying each subsystem and realizing offset current of the circuit, so as to achieve the purpose of lowering circuit power consumption, improving anti-jamming capability, and simplifying circuit design difficulty. The method is suitable for integrating wavelet transfer function through the switch current circuit, that is, design wavelet filter circuit by adopting switch current technology.
Description
Technical field
The invention belongs to Electronics and Information Engineering field, relate to a kind of switching current wavelet filter state space optimization method.
Background technology
Wavelet transformation, with its good time-frequency local characteristics, is widely used in non-stationary and transient signal process, has now become one of the most effective mathematical tool of each engineering field signal transacting.Wavelet transformation also can use hardware implementing with software simulating.Use hardware implementing wavelet transformation, particularly realize with mimic channel, owing to having low in energy consumption, fireballing advantage, be subject to more and more attention.Mimic channel realizes the linear combination that wavelet transformation can regard the yardstick wavelet filter different with displacement as, so the design of wavelet filter is the basis realizing wavelet transformation.Wavelet filter bank has two crucial steps, and one is approaching of wavelet filter, and two is circuit design of wavelet filter.At present, no matter be that time domain is approached or frequency domain approaches, be all finally to obtain wavelet filter frequency domain transfer function H (s).If frequency domain approaches, directly frequency domain transfer function H (s) can be obtained, if time domain is approached, time domain approximating function h (t) frequency domain transfer function H (s) need be obtained by Laplace transform.Circuit design wavelet filter is adopted after obtaining frequency domain transfer function H (s), the mimic channel kind that can be used to design wavelet filter is a lot, and what report was more at present is have low-voltage, low-power consumption, the Switched-Current Circuit technology compatible with standard digital CMOS industry.As the patent of invention " a kind of method for designing wavelet filter " that the patent No. is 2011102989344, it adopts Maclaurin series to approach wavelet frequency domain function, then adopts the biquadratic joint based on switching current differentiator to design the wavelet filter of cascaded structure for elementary cell; The patent No. is the patent of invention " switching current self-adaptation chaos differential evolution method for designing wavelet filter " of 2011100545508, it adopts self-adaptation chaos differential evolution optimization algorithm to carry out wavelet filter time domain and approaches, then adopt Laplace transform to obtain frequency domain transfer function H (s) of wavelet filter, then adopt the biquadratic joint based on switched-current integrator to carry out filter synthesis; The patent No. is the patent of invention " a kind of simulation continuous wavelet transform circuit based on switched current technique " of 2010102503863, adopt based on switched-current integrator single order joint and biquadratic save carry out filter synthesis.
Above-mentioned method for designing does not all consider the state space problem of switching current wavelet filter circuit, peak value and the dynamic range of each economize on electricity road output of wavelet filter circuit cascade are inconsistent, if in order to reduce circuit power consumption, the bias current on each economize on electricity road will be set to different value according to output peak value and dynamic range, cause the bias current disunity on each economize on electricity road, increase difficulty to circuit design; If in order to the convenience designed, the bias current on each economize on electricity road can be unified, but this be will with peak value and dynamic range maximum one economize on electricity road for foundation, unified bias current is higher than the peak value of the dynamical output electric current on this economize on electricity road, and this can increase the power consumption of designed wavelet filter circuit undoubtedly.
Summary of the invention
In order to overcome in existing analog current wavelet filter bank above shortcomings and meet low-power consumption, highly reliable practical application needs, the invention provides a kind of analog wavelet filter status space optimization method, method is simple, highly versatile, be applicable to low pressure, low-power consumption and high frequency, at a high speed, the engineer applied of highly reliable analog wavelet filter circuit design and reality.
The optimization in switching current wavelet filter circuit state space is the optimization of analog wavelet filter transfer function, only has the transport function of an acquisition state space optimum, just can design the circuit of state space optimum.
Technical scheme of the present invention is: switching current wavelet filter circuit state space optimization method, comprises the following steps: (1) normalization cascade decomposition analog wavelet filter transfer function; (2) normalization cascade each sub-system impulse response peak value is asked, calculation optimization coefficient; (3) Optimized Simulated wavelet filter transport function: each sub-system of normalization cascade is multiplied by optimized coefficients, make analog wavelet filter status spatial optimization, namely make the impulse response peak value of each sub-system of cascade in the same size, then make the dynamic range of the Switched-Current Circuit realizing each sub-system identical.
In step (1), described normalization cascade decomposition analog wavelet filter transfer function refers to and wavelet filter transfer function H (s) is decomposed into normalized cascade form, namely
, wherein k represents coefficient; N represents the joint number of cascade; I is subscript, represents the numbering of cascade joint number; H
is () represents the transport function of each sub-system of normalized cascade before optimizing, it is characterized in that the highest term coefficient of molecule and denominator polynomials is 1; ∏ connects to take advantage of negative sign, realizes the cascade of analog wavelet wave filter.
In step (2), described in ask normalization cascade each sub-system impulse response peak value, refer to the transfer function H asking subsystem
is the peak value of () impulse response, can solve in numerical analysis, also can be solved by Matlab software.
Described calculation optimization coefficient refers to that the coefficient k adopting following calculating formula normalization cascade decomposition analog wavelet filter transfer function H (s) to be obtained resolves into independently optimized coefficients
k i
（1），
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
k i represent optimized coefficients, go in the i-th sub-system transport function being assigned to cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of negative sign; Optimized coefficients
k i relational expression k=is met with coefficient k
.
In step (3), described Optimized Simulated wavelet filter transport function: each sub-system of normalization cascade is multiplied by optimized coefficients, and concrete methods of realizing is: the optimized coefficients calculated according to (1) formula
k i the normalization cascade decomposition expression formula of amendment analog wavelet filter transfer function, namely adopts following formula to revise the cascade decomposition expression formula of analog wavelet filter transfer function
（2），
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
represent the transport function of removing each sub-system of normalized cascade after optimizing, wherein
k irepresent that i-th saves the optimized coefficients value be assigned to, H
is () represents the transport function of each sub-system of normalized cascade before optimizing.
In step (3), described in make the impulse response peak value of each sub-system of cascade after optimizing in the same size, concrete methods of realizing is: according to the impulse response peak value size calculating formula of each sub-system of cascade after optimizing
（3），
Substitute into (3) formula by (1) formula,
（4），
(3) formula is with in (4) formula,
a p represent the impulse response peak value of each sub-system of cascade after optimizing; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
k i represent optimized coefficients;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of negative sign; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains;
The described dynamic range of the Switched-Current Circuit realizing each sub-system that then makes is identical, and the current signal dynamic range expression formula namely realizing the Switched-Current Circuit of each sub-system is
（5），
In formula, I
p represent the current signal realizing the Switched-Current Circuit of each sub-system; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of negative sign; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains.
Principle of the present invention is: normalization cascade decomposition analog wavelet filter transfer function can produce a coefficient, distributes this coefficient by optimizing thus reaches beneficial effect of the present invention; Optimize the principle of this coefficient of distribution as shown in calculating formula (1), analysis meter formula (1) is known, if certain saves normalized subsystem impulse response peak value before optimizing
a i comparatively large, then from scale-up factor
kin be assigned to a less optimized coefficients
k i , thus its impulse response peak value is reduced, compress its dynamic range simultaneously; If impulse response peak value
a i less, be then assigned to larger optimized coefficients
k i , thus its impulse response peak value is expanded, expand the current signal dynamic range of its realizing circuit simultaneously; Thus obtain the optimized cascade form of analog wavelet filter transfer function state space, and then the Switched-Current Circuit realizing analog wavelet wave filter is made to have optimum state space.
The switching current wavelet filter circuit state space optimization method of the present invention, by the coefficient that reasonable distribution normalization cascade decomposition analog wavelet filter transfer function produces, make the impulse response peak value of each sub-system of analog wavelet wave filter cascade consistent, when wavelet filter comprehensive by Switched-Current Circuit, the state space of the Switched-Current Circuit of each sub-system of cascade reaches optimization, namely the current signal dynamic range on each economize on electricity road is identical, thus reach the bias current of unified each sub-system realizing circuit, reduce circuit power consumption, improve antijamming capability, simplify the object of circuit design difficulty.The present invention is suitable for further by the comprehensive wavelet filter transport function of Switched-Current Circuit, namely adopts Switched-Current Circuit technical design wavelet filter circuit.
Accompanying drawing explanation
Fig. 1 is implementation step figure of the present invention;
Fig. 2 (a) is the impulse response oscillogram before six rank Gauss wavelet wave filter first segment subsystems are optimized;
Fig. 2 (b) is the impulse response oscillogram after six rank Gauss wavelet wave filter first segment subsystems are optimized;
Fig. 3 (a) is the impulse response oscillogram before six rank Gauss wavelet wave filter second section subsystems are optimized;
Fig. 3 (b) is the impulse response oscillogram after six rank Gauss wavelet wave filter second section subsystems are optimized;
Fig. 4 (a) is the impulse response oscillogram before six rank Gauss wavelet wave filter the 3rd sub-system are optimized;
Fig. 4 (b) is the impulse response oscillogram after six rank Gauss wavelet wave filter the 3rd sub-system are optimized;
Fig. 5 is each sub-system impulse response oscillogram before six rank Gauss wavelet filter optimizations;
Fig. 6 is each sub-system impulse response oscillogram after six rank Gauss wavelet filter optimizations;
Fig. 7 is system impulse response wave shape figure before six rank Gauss wavelet filter optimizations;
Fig. 8 is system impulse response wave shape figure after six rank Gauss wavelet filter optimizations.
Embodiment
Below in conjunction with drawings and Examples, the invention will be further described.
With reference to Fig. 1, the present invention includes following steps: (1) normalization cascade decomposition analog wavelet filter transfer function; (2) normalization cascade each sub-system impulse response peak value is asked, calculation optimization coefficient; (3) each sub-system of normalization cascade is multiplied by optimized coefficients, make analog wavelet filter status spatial optimization, namely make the impulse response peak value of each sub-system of cascade in the same size, then make the state space of the Switched-Current Circuit realizing each sub-system consistent, namely the dynamic range of current signal is identical.
In step (1), described normalization cascade decomposition analog wavelet filter transfer function refers to and wavelet filter transfer function H (s) is decomposed into normalized cascade form, namely
, wherein k represents coefficient; N represents the joint number of cascade; I is subscript, represents the numbering of cascade joint number; H
is () represents the transport function of each sub-system of normalized cascade before optimizing, it is characterized in that the highest term coefficient of molecule and denominator polynomials is 1; ∏ connects to take advantage of negative sign, realizes the cascade of analog wavelet wave filter.
In step (2), described in ask normalization cascade each sub-system impulse response peak value, refer to the transfer function H asking subsystem
is the peak value of () impulse response, can solve in numerical analysis, also can be solved by Matlab software.
Described calculation optimization coefficient refers to that the coefficient k adopting following calculating formula normalization cascade decomposition analog wavelet filter transfer function H (s) to be obtained resolves into independently optimized coefficients
k i
（1），
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
k i represent optimized coefficients, go in the i-th sub-system transport function being assigned to cascade;
a i represent the peak value of the impulse response of the i-th sub-system transport function; ∏ connects to take advantage of negative sign; Optimized coefficients
k i relational expression k=is met with coefficient k
.
In step (3), described Optimized Simulated wavelet filter transport function: each sub-system of normalization cascade is multiplied by optimized coefficients, and concrete methods of realizing is: the optimized coefficients calculated according to (1) formula
k i the cascade decomposition expression formula of amendment analog wavelet filter transfer function, namely adopts following formula to revise the normalization cascade decomposition expression formula of analog wavelet filter transfer function
（2），
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade; H
is () represents the transport function of each sub-system of normalized cascade before optimizing;
represent the transport function of removing each sub-system of normalized cascade after optimizing, wherein
k irepresent that i-th saves the optimized coefficients value be assigned to.
In step (3), described in make the impulse response peak value of each sub-system of cascade after optimizing in the same size, concrete methods of realizing is: according to the impulse response peak value size calculating formula of each sub-system of cascade after optimizing
（3），
Substitute into (3) formula by (1) formula,
（4），
(3) formula is with in (4) formula,
a p represent the impulse response peak value of each sub-system of cascade after optimizing; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
k i represent optimized coefficients;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of negative sign; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains;
The described dynamic range of the Switched-Current Circuit realizing each sub-system that then makes is identical, and the current signal dynamic range expression formula namely realizing the Switched-Current Circuit of each sub-system is
（5），
In formula, I
p represent the current signal realizing the Switched-Current Circuit of each sub-system; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of negative sign; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains.
The inventive method is suitable for any switching current wavelet filter circuit state space optimization.In order to the convenience illustrated, be described for six rank Gauss wavelet filter circuit state space optimizations below.
The patent No. be 2011100545508 patent of invention disclose one " switching current self-adaptation chaos differential evolution method for designing wavelet filter ", it adopts self-adaptation chaos differential evolution optimization algorithm to carry out wavelet filter time domain and approaches, and then adopts Laplace transform to obtain frequency domain transfer function H (s) of six rank Gauss wavelet wave filters as follows
（6），
Comprise three sub-system transport functions in formula, use
,
with
represent, be respectively
（7），
（8），
（9），
Three sub-system
,
with
impulse response waveform respectively as shown in Fig. 2 (a), 3 (a) He 4 (a).As seen from the figure, first segment subsystem
impulse response peak value appeared at for 0 moment, peak value size is 0.0696, and dynamic range is [-0.038 ~ 0.0696]; Second section subsystem
impulse response peak value appear at 0.5s, peak value size is 19.61, and dynamic range is [-4.99 ~ 19.61]; 3rd sub-system
impulse response peak value appear at 1.17s, peak value size is 7.52, and dynamic range is [-3.49 ~ 7.52].As can be seen here, based on the so simple cascade decomposition patten's design analog wavelet filter circuit of (6) formula, then will there is very large difference in the output peak value of its each knot circuit and dynamic range, therefore the bias current of each joint realizing circuit is arranged, and there are the following problems: if in order to reduce circuit power consumption, the bias current on each economize on electricity road will be set to different value according to output peak value and dynamic range, cause the bias current disunity on each economize on electricity road, increase difficulty to circuit design; If in order to the convenience designed, the bias current on each economize on electricity road can be unified, but this be will with peak value and dynamic range maximum one economize on electricity road for foundation, here be second section peak value maximum (being 19.61), unified bias current is higher than the peak value (such as bias current is 20 to 25) of the dynamical output electric current on this economize on electricity road, and this can increase the power consumption of designed wavelet filter circuit undoubtedly.
Frequency domain transfer function below for above-mentioned six rank Gauss wavelet wavelet filters carries out state space optimization.
1st step: wavelet filter transfer function H (s) is decomposed in normalization, is decomposed into normalized cascade form by wavelet filter transfer function H (s) shown in formula (3),
（10），
Obtain coefficient k=0.0696, normalized three sub-system transport function H
1(s), H
2(s) and H
3(s) represent, be respectively into
（11），
（12），
（13），
2nd step: calculation optimization coefficient.Use A
1, A
2and A
3represent the normalized three sub-system transport function H before optimizing respectively
1(s), H
2(s) and H
3s the impulse response peak value of (), carries out numerical analysis and calculates: A
1=1, A
2=19.61, A
3=7.52, and the joint number N=3 of six rank Gauss wavelet wave filter cascades, according to (1) formula, three optimized coefficients k can be calculated
1=2.1733; k
2=0.1108; k
3=0.2890.
3rd step: Optimization of Wavelet filter transfer function.Three the optimized coefficients k tried to achieve by the 2nd step
1, k
2and k
3be multiplied by the normalized three sub-system transfer function H that the 1st step is tried to achieve respectively
1(s), H
2(s) and H
3s (), obtains and goes normalized three sub-system transport functions after state space optimization, use respectively
,
with
represent:
（14），
（15），
（16），
Normalized three sub-system transport functions are gone after being optimized by cascade
,
with
, six rank Gauss wavelet filter transfer functions of state space optimization can be obtained:
（17），
Theoretical analysis: six rank Gauss wavelet wave filters after state space optimization, its three sub-system transport functions comprised
,
with
impulse response peak value size be that 2.713(calculates according to (4) formula), the circuit dynamic range then realizing the Switched-Current Circuit of three sub-system is [-2.713,2.713] (obtaining according to (5) formula).
Simulation analysis: the impulse response waveform of three sub-system after optimization is respectively as shown in Fig. 2 (b), 3 (b) He 4 (b).As seen from the figure, first segment subsystem transport function
impulse response peak value appeared at for 0 moment, peak value size is 2.713, and dynamic range is [-1.178 ~ 2.713]; Second section subsystem transport function
impulse response peak value appear at 0.5s, peak value size is 2.713, and dynamic range is [-0.553 ~ 2.713]; 3rd sub-system transport function
impulse response peak value appear at 1.17s, peak value size is 2.713, and dynamic range is [-1.009 ~ 2.713].As can be seen here, (17) formula after state space optimization, the output peak value of its each sub-system reaches unanimously (peak value size is 2.713), and dynamic range is limited to [-2.713,2.713], the dynamic range of first segment subsystem changes to [-1.178 ~ 2.713] after optimization (belonging to [-2.713,2.713]) from [-0.038 ~ 0.0696] before optimizing; The dynamic range of second section subsystem changes to [-0.553 ~ 2.713] after optimization (belonging to [-2.713,2.713]) from [-4.99 ~ 19.61] before optimizing; The dynamic range of the 3rd sub-system changes to [-1.009 ~ 2.713] after optimization (belonging to [-2.713,2.713]) from [-3.49 ~ 7.52] before optimizing.After optimizing, the realizing circuit of each sub-system can according to identical circuit dynamic range [-2.713,2.713] carry out bias current setting, as bias current is set to 3 to 5, reduce about 85% than the unified bias current (being 20 to 25) before optimizing.
In order to the contrast situation before and after display state space optimization more directly perceived, by the impulse response waveform drawing of each sub-system of six rank Gauss wavelet wave filters in same figure.
The impulse response waveform of six rank Gauss wavelet wave filters each sub-system before state space optimization as shown in Figure 5.
The impulse response waveform of six rank Gauss wavelet wave filters each sub-system after state space optimization as shown in Figure 6.
Comparison diagram 5 and Fig. 6 known, the each sub-system of six rank Gauss wavelet wave filters before optimization, there is significant differences in state space, the peak value of the 3rd sub-system and dynamic range are obviously little than the second subsystem, the peak value of first segment subsystem and dynamic range little to cannot distinguish in the drawings (concrete numerical value can see Fig. 2 (a), 3 (a) and 4 (a)); After optimization of the present invention, the state space of each sub-system is basically identical, the peak value consistent (peak value size is 2.713) of three joints, and dynamic range is limited to [-2.713,2.713].
In order to check whether state space optimization affects analog wavelet wave filter external characteristic, i.e. the input and output characteristic of analog wavelet filter system, spy draws out the impulse response oscillogram of six rank Gauss wavelet overall filter systems.
The impulse response waveform of six rank Gauss wavelet wave filters total system before state space optimization as shown in Figure 7.
Six rank Gauss wavelet wave filters after state space optimization total system impulse response waveform as shown in Figure 8.
Comparison diagram 7 and Fig. 8 known, the each sub-system of six rank Gauss wavelet wave filters before optimization after, the impulse response waveform of total system does not have difference, namely state space optimization does not affect the input and output characteristic of wavelet filter total system, just change internal system structure, optimize internal system state space.
Comprehensive above-mentioned simulation result and analysis show, the present invention effectively can optimize the state space of little wave filter, namely the state space of each sub-system of wavelet filter cascade is made to reach optimization, when with circuit synthesis wavelet filter, the impulse response peak value of the realizing circuit of each sub-system of cascade reaches consistent, circuit dynamic range is identical, thus reaches the bias current unifying each sub-system realizing circuit, the object reducing circuit power consumption, improve antijamming capability, simplify circuit design difficulty.The present invention is suitable for further by the comprehensive wavelet filter transport function of Switched-Current Circuit, namely adopts switched current technique design wavelet filter circuit.
Claims (4)
1. switching current wavelet filter circuit state space optimization method, is characterized in that: comprise the following steps: (1) normalization cascade decomposition analog wavelet filter transfer function; (2) normalization cascade each sub-system impulse response peak value is asked, calculation optimization coefficient; (3) Optimized Simulated wavelet filter transport function: each sub-system of normalization cascade is multiplied by optimized coefficients, make analog wavelet filter status spatial optimization, namely make the impulse response peak value of each sub-system of cascade after optimizing in the same size, then make the dynamic range of the Switched-Current Circuit realizing each sub-system identical;
In step (2), described calculation optimization coefficient refers to that the coefficient k adopting following calculating formula normalization cascade decomposition analog wavelet filter transfer function H (s) to be obtained resolves into independently optimized coefficients
k i :
，
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
k i represent optimized coefficients, go in the i-th sub-system transport function being assigned to cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of symbol; Optimized coefficients
k i relational expression k=is met with coefficient k
.
2. switching current wavelet filter circuit state space optimization method according to claim 1, it is characterized in that: in step (1), described normalization cascade decomposition analog wavelet filter transfer function refers to and wavelet filter transfer function H (s) is decomposed into normalized cascade form, namely
, wherein k represents coefficient, and k>0; N represents the joint number of cascade; I is subscript, represents the numbering of cascade joint number; H
is () represents the transport function of each sub-system of normalized cascade before optimizing, it is characterized in that the highest term coefficient of molecule and denominator polynomials is 1; ∏ connects to take advantage of symbol, realizes the cascade of analog wavelet wave filter.
3. switching current wavelet filter circuit state space optimization method according to claim 1, it is characterized in that: in step (3), described Optimized Simulated wavelet filter transport function: each sub-system of normalization cascade is multiplied by optimized coefficients, and concrete methods of realizing is: according to the optimized coefficients calculated
k i the cascade decomposition expression formula of amendment analog wavelet filter transfer function, namely adopts following formula to revise the cascade decomposition expression formula of analog wavelet filter transfer function
，
In formula, i is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
represent the transport function of removing each sub-system of normalized cascade after optimizing, wherein
k irepresent that i-th saves the optimized coefficients value be assigned to, H
is () represents the transport function of each sub-system of normalized cascade before optimizing.
4. switching current wavelet filter circuit state space optimization method according to claim 1, it is characterized in that: in step (3), the described impulse response peak value of each sub-system of cascade after optimizing that makes is in the same size, and particular content is: the impulse response peak value size calculating formula of each sub-system of the cascade after optimization is
，
In formula,
a p represent the impulse response peak value of each sub-system of cascade after optimizing; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of symbol; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains;
The described dynamic range of the Switched-Current Circuit realizing each sub-system that then makes is identical, and the current signal dynamic range expression formula namely realizing the Switched-Current Circuit of each sub-system is
，
In formula, I
p represent the current signal realizing the Switched-Current Circuit of each sub-system; I is subscript, represents the numbering of cascade joint number; N represents the joint number of cascade;
a i represent the peak value of the impulse response of the transport function of each sub-system of normalized cascade before optimizing; ∏ connects to take advantage of symbol; K represents the coefficient that normalization cascade decomposition analog wavelet filter transfer function H (s) obtains.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310146139.2A CN103226637B (en) | 2013-04-25 | 2013-04-25 | Optimizing method for state space of switch current wavelet filter circuit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310146139.2A CN103226637B (en) | 2013-04-25 | 2013-04-25 | Optimizing method for state space of switch current wavelet filter circuit |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103226637A CN103226637A (en) | 2013-07-31 |
| CN103226637B true CN103226637B (en) | 2015-07-08 |
Family
ID=48837082
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310146139.2A Active CN103226637B (en) | 2013-04-25 | 2013-04-25 | Optimizing method for state space of switch current wavelet filter circuit |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103226637B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176214A (en) * | 2011-03-08 | 2011-09-07 | 长沙河野电气科技有限公司 | Method for designing wavelet filter by means of switching current and adaptive chaotic differential evolution |
| US8069200B2 (en) * | 2005-04-28 | 2011-11-29 | Qsigma, Inc. | Apparatus and method for implementing floating point additive and shift operations |
-
2013
- 2013-04-25 CN CN201310146139.2A patent/CN103226637B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8069200B2 (en) * | 2005-04-28 | 2011-11-29 | Qsigma, Inc. | Apparatus and method for implementing floating point additive and shift operations |
| CN102176214A (en) * | 2011-03-08 | 2011-09-07 | 长沙河野电气科技有限公司 | Method for designing wavelet filter by means of switching current and adaptive chaotic differential evolution |
Non-Patent Citations (1)
| Title |
|---|
| 连续小波变换的对数模拟滤波器实现;黄姣英等;《电子测量与仪器学报》;20090615;第23卷(第6期);第52-56页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103226637A (en) | 2013-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104866665A (en) | Hybrid simulation method including power electronic equipment based on interface equivalence and interaction | |
| Chen et al. | Three-phase power-flow by direct ZBR method for unbalanced radial distribution systems | |
| CN102843230B (en) | Mathematical model of four-dimensional autonomous hyper-chaos system and achieving circuit of mathematical model | |
| CN108448631B (en) | Dynamic phasor modeling method based on model order reduction and including distributed power supply access power distribution network | |
| CN104993491B (en) | Linear power flow calculation method with voltage and reactive power being taken into consideration | |
| CN105426594A (en) | Method for fast simulating wind fields on basis of stability and homogeneity of time-space field and condition interpolation | |
| CN103095191A (en) | Switch reluctance motor memory sensor model modeling method | |
| CN103970717A (en) | Unconditional stability FDTD algorithm based on Associated Hermite orthogonal function | |
| CN105631105A (en) | Cable equivalent circuit modeling method containing applied field coupling | |
| CN105635009A (en) | Self-adaptive MIMO pre-distortion method for hybrid compensation of multi-branch crosstalk and IQ imbalance | |
| CN107918292A (en) | Towards the Power Electronic Circuit transient emulation GPU accelerated methods of exponential integral | |
| CN111861789A (en) | Method and system for identifying user variable relationship based on current data integer programming solution | |
| CN103729502B (en) | A kind of method improving electromagnetic transient in power system simulation velocity | |
| CN103226637B (en) | Optimizing method for state space of switch current wavelet filter circuit | |
| CN103901255A (en) | Hybrid power filter harmonic current detection method based on generalized weighted morphological filtering theory | |
| CN110569558B (en) | Hybrid electromagnetic transient simulation method suitable for micro-grid real-time simulation | |
| CN101778291A (en) | Lifting structure-based DCT conversion structure and method thereof | |
| CN106533138B (en) | The unified approach that multi-level converter particular harmonic is eliminated | |
| CN103262067B (en) | A kind of data processing method, data processing equipment and communication system | |
| CN104536807B (en) | DC/DC real-time simulators and method based on FPGA | |
| CN107480109B (en) | Incomplete parallel state estimation realization method based on network structure dynamic partitioning | |
| CN102799765B (en) | The real-time unfixed point interpolation computing method of switch in a kind of step-length | |
| Chandrasekar et al. | Integer lifting wavelet transform based hybrid active filter for power quality improvement | |
| CN103902769A (en) | Model order reduction method of asymptotic waveform estimation algorithm | |
| CN102946239B (en) | Design method of multi-ring feedback wavelet filter for switching current-fused negative phase differentiator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |