Electromagnetic Interference Reduction Calculation Method, Device, and Program, and Its Electronic Circuit

Information

  • Patent Application
  • 20080077337
  • Publication Number
    20080077337
  • Date Filed
    August 22, 2005
    19 years ago
  • Date Published
    March 27, 2008
    16 years ago
Abstract
A calculation method of electromagnetic interference reduction through spectrum diffusion that is a means for reducing the electromagnetic interference, in which, when a center frequency of a measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in the measurement bandwidth, of a plurality of spectra generated from the frequency spectrum of an electromagnetic interference signal through the spectrum diffusion, the amplitudes of all spectra included in the measurement bandwidth are added up and the maximum total sum of the amplitudes is divided by the amplitude of the electromagnetic interference signal, hence to estimate the electromagnetic interference reduction.
Description
FIELD OF THE INVENTION

The invention relates to a method of estimating electromagnetic interference reduction, and particularly to an electromagnetic interference calculation method, device, and program, and its electronic circuit which can obtain electromagnetic interference reduction only through calculation, without calibration through the measurement of electromagnetic interference reduction, which calibration is required owing to the influence of the measurement bandwidth of a measurement device, when measuring the electromagnetic interference reduction.


DESCRIPTION OF THE RELATED ART

According to miniaturization of electronic devices, influence of inductive coupling between transmission lines is not negligible in a highly-packed electronic circuit, and there may occur electromagnetic interference (EMI) causing an unexpected malfunction of an electronic circuit. It is known that electromagnetic interference often occurs especially in the clock waveform used for a digital circuit.


In the regulation about the electromagnetic interference, a test method about the discharge amount of the electromagnetic interference and its allowance are defined and the products on the market should satisfy the corresponding EMI regulations.


When each measurement device measures electromagnetic interference, its measurement bandwidth is defined in the standard regulations (in CISPR22 of Comittie International Special des Perturbations Radioelectrique and in VCCI standard based on this international standards): it is defined as 9 kHz in the measurement of a target in a frequency range of 150 kHz to 30 MHz; as 120 kHz, in the frequency range of 30 MHz to 1 GHz; and, as 1 MHz, in the frequency exceeding 1 GHz.


In a spectrum diffusion clock technique, by modulating a clock waveform of defining an operational timing inside a computer according to a frequency modulation technique, the frequency spectrum is expanded and electric power is dispersed, hence to reduce the electromagnetic interference in a specified frequency. Reduction of the electromagnetic interference is controlled while changing a parameter for modulation.


In this method, reduction of the electromagnetic interference is controlled while changing a parameter for frequency modulation. In the conventional design method, however, the electromagnetic interference reduction effect has been evaluated, according to the properties of the electromagnetic interference after a device is completed. Therefore, when a target effect is not obtained, it is necessary to change the design repeatedly.


In the measurement of the electromagnetic interference reduction, when the bandwidth of a measurement device is wider than the modulation frequency, the measured electromagnetic interference reduction depends on the bandwidth of the measurement device disadvantageously and a level difference has to be calibrated according to the bandwidth of the measurement device.


An example of the method about solving this situation is disclosed, for example, in Patent Publication Laid-Open No. 2003-150660 (Document 1).


According to the method disclosed in the document 1, a spectrum distribution of a signal when the clock signal is frequency-modulated is calculated and a theoretical value of the electromagnetic interference reduction is derived, under a condition of given parameters for frequency modulation.


In the method of the document 1, when the bandwidth of a measurement device is wider than the modulation frequency, at first, a calculation method in the case where the bandwidth of a measurement device is narrower than the modulation frequency, in which method the measured electromagnetic interference reduction does not depend on the bandwidth of the measurement device, is used to calculate the frequency, then the frequency is measured by the measurement device, and the calculated value is calibrated by using this measured value, thereby calculating the theoretical value of the electromagnetic interference reduction.


In the method of the document 1, as for a level difference according to the bandwidth of a measurement device, a difference between the measured value and the calculated value is required and the calculated value is calibrated using the difference, hence to estimate the electromagnetic interference suppression effect through calculation.


The above conventional technique has the following problems.


According to the method of document 1, when the bandwidth of a measurement device does not have to be taken into consideration, in other words, when the bandwidth of a measurement device is narrower than the modulation frequency, the electromagnetic interference reduction can be calculated using a calculation formula.


On the other hand, when it has to be taken into consideration, in other words, when the measurement bandwidth of a measurement device is wider than the modulation frequency, the level calibration value needed at that time is obtained through measurement. Therefore, when the measurement bandwidth is wider than the modulation frequency, it requires another measurement for level calibration disadvantageously.


In the case of the conventional method including the method of the document 1, there is a problem that a calculation error occurs because no consideration is taken to a difference between the definition of the measurement bandwidth of a measurement device supposed to be ideal and the definition of an actually measured measurement bandwidth.


In the definition of the ideal measurement bandwidth, an input signal level does not change within a specified band and an input signal is eliminated outside the specified band. Specifically, when the ideal filter property is obtained within the measurement bandwidth of a measurement device, the calculated value can be used as it is.


In the actual measurement device, however, although an input signal does not change in the vicinity of the center frequency within the band, sensibility to the input signal decreases in the frequency in the vicinity of the both ends within the measurement band. Therefore, the influence of the sensibility reduction in the end frequencies within the measurement band has to be taken into account.


Further, even outside a specified band, an input signal is not completely eliminated, but it is detected by a measurement device at a constant rate. This causes an error, which is not negligible.


In the calculation of electromagnetic interference reduction through spectrum diffusion measured by a measurement device having a given measurement bandwidth. In order to solve the above problems of the conventional technique, a first object of the invention is to provide an electromagnetic interference reduction calculation method, device, calculation program, and its electronic circuit which can estimate the electromagnetic interference reduction only through the calculation, without calibration of the electromagnetic interference reduction through measurement which has been required when the measurement bandwidth of the measurement device is wider than the modulation frequency.


A second object of the invention is to provide an electromagnetic interference reduction calculation method, device, program and its electronic circuit which can calibrate the influence of sensibility reduction in a measurement device in the vicinity of the frequency at the both ends within the measurement bandwidth.


SUMMARY OF THE INVENTION

In order to achieve the above objects, the invention is a method of calculating electromagnetic interference reduction through spectrum diffusion that is a means of reducing electromagnetic interference, comprising: adding up the amplitudes of all the spectra included in a measurement bandwidth when a center frequency of the measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in the measurement bandwidth, of a plurality of spectra generated from frequency spectrum of an electromagnetic interference signal through the spectrum diffusion; and estimating electromagnetic interference reduction by dividing the maximum total sum of the amplitudes by the amplitude of the electromagnetic interference signal.


According to the invention of Claim 2, when the measurement bandwidth of the electromagnetic interference reduction is wider than a modulation frequency that is a frequency fluctuation width set through the spectrum diffusion, the amplitudes of the spectra are added up.


According to the conventional calculation method of the electromagnetic interference reduction, when the measurement bandwidth of the electromagnetic interference reduction is wider than the modulation frequency that is a frequency fluctuation width set through the spectrum diffusion, the amplitude of the spectrum having the maximum amplitude included in the measurement bandwidth, of a plurality of spectra generated from the frequency spectrum of the electromagnetic interference signal through the spectrum diffusion, is divided by the amplitude electromagnetic interference signal, hence to estimate the electromagnetic interference reduction.


In this calculation method, since influence of the spectra other than the spectrum having the maximum amplitude is not considered, it is impossible to estimate a reliable electromagnetic interference reduction. Therefore, in the method of the invention, a difference between the calculated value according to the conventional method and the measured value measured under the same conditions is examined previously, and using the result of the difference, the calculated value is calibrated when estimating the electromagnetic interference reduction.


The difference between the calculated value according to the conventional method and the measured value is caused by the influence of sensibility reduction of a measurement device int he frequencies at the both ends of the measurement bandwidth, which is described according to Claim 2 of the invention, in addition to the influence of the spectra other than the spectrum having the maximum amplitude, which has not been a target of the calculation in the conventional method. Thus, the causes of the difference between the calculated value according to the conventional method and the measured value are two. The cause of the difference includes the influence of the spectra other than the influence of the sensibility reduction of a measurement device in the frequencies at the both ends of the measurement bandwidth.


The calculation method according to the invention can be associated with the measured value according to the conventional method as follows.


At first, the amplitudes of all the spectra included in the measurement bandwidth are added up, to estimate the electromagnetic interference reduction. The resultant value is obtained by calibrating the electromagnetic interference reduction obtained in the conventional method, with the influence taken into account, caused by the amplitudes of the spectra other than the spectrum having the maximum amplitude.


Next, the calculated electromagnetic interference reduction is calibrated according to the second calibration with the influence taken into account, caused by the sensibility reduction of a measurement device in the frequencies at the both ends of the measurement bandwidth, hence to estimate the measured value obtained in the conventional method, only through calculation.




BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a flow chart showing the flow of calculating the electromagnetic interference reduction according to a first embodiment of the invention;



FIG. 2 is a view showing one example of a frequency spectrum waveform in an envelope display for use in describing the first embodiment of the invention;



FIG. 3 is a view showing one example of a frequency spectrum waveform in a line spectrum display for use in describing the first embodiment of the invention;



FIG. 4 is a view showing an example of filter property for used in describing the first embodiment of the invention;



FIG. 5 is a view showing a comparison between the calculated value and the measured value of the electromagnetic interference reduction according to the first embodiment of the invention;



FIG. 6 is a view showing the amplitude obtained by the Bessel function, for use in describing the first embodiment of the invention;



FIG. 7 is a view showing an example of filter property of a measurement device for use in describing the first embodiment of the invention;



FIG. 8 is a flow chart for use in describing the processing performed by an electromagnetic interference reduction calculation device according to the first embodiment of the invention;



FIG. 9 is a view showing the structure of the electromagnetic interference reduction calculation device according to the first embodiment of the invention; and



FIG. 10 is a view showing an interference of a spectrum diffusion signal from an electronic circuit to a cellular phone according to a second embodiment of the invention.




DESCRIPTION OF EMBODIED CONFIGURATION OF THE INVENTION

Hereinafter, preferred embodiment of the invention will be described in detail with reference to the drawings.



FIG. 1 is a flow chart showing the flow of calculating the electromagnetic interference reduction according to the embodiment.


Referring to FIG. 1, calculation of electromagnetic interference reduction according to the invention includes the following five Steps.


At first, data input is performed. Here, clock frequency, modulation frequency, maximum frequency shift, and measurement bandwidth of a measurement device are supplied (Step 101).


This time, electromagnetic interference reduction is roughly estimated (Step 102).


The Step is performed according to the conventional technique.


A first calibration (Step 102) of the electromagnetic interference reduction and a second calibration (Step 104) of the electromagnetic interference reduction are performed according to a calculation of the invention.


By using the above results, calculation of the electromagnetic interference reduction is completed (Step 105).


Hereafter, each Step of FIG. 1 will be described more specifically.


A rough estimation of the electromagnetic interference reduction according to the conventional technique will be described at first.


It is well known that waveform when a sine wave of carrier frequency f(Hz) is frequency-modulated with a sine wave of modulation frequency fm(Hz) and the maximum frequency shift Δf(Hz)=f×modulation degree p(%) is given by the formula 1.


The formula 1 is described in, for example, “Frequency Analysis/Modulation and Noise Theory” written by Toshio Hosono et al. p. 168-178 Kourin Shoinn (1064) and “Radio Communication Technology” written by Isamu Tanimura p. 46-49 Coronasha (1960).
a=A{J0(Δffm)sin(2πft+θ0)+J1(Δffm)sin{2π(f+fm)t+θ0}-J1(Δffm)sin{2π(f-fm)t+θ0}+J2(Δffm)sin{2π(f+2fm)t+θ0}-J2(Δffm)sin{2π(f-2fm)t+θ0}+J3(Δffm)sin{2π(f+3fm)t+θ0}-J3(Δffm)sin{2π(f-3fm)t+θ0}+}Formula(1)


Wherein, Jn(x) indicates the Bessel function, θ0 indicates the initial value of each frequency. Å indicates the amplitude of a clock waveform, t indicates a time, and n indicates the order of the Bessel function.


The term of Jn(Δf/fm) indicates the amplitude of each spectrum and the term of sin(2p(f+nfm)t) indicates a time change of each spectrum.


In other words, the formula 1 is an aggregate of signals different in amplitude and frequency. When a signal component corresponding to the term having the largest amplitude is selected, in this aggregation, it means that the electromagnetic interference reduction is the smallest. The electromagnetic interference reduction is evaluated according to the measurement of the electromagnetic interference reduction (worst value) corresponding to the term having the largest amplitude.



FIG. 2 is a view showing one example of the frequency spectrum waveform in an envelope display for use in describing the embodiment.


With reference to FIG. 2, it is found that the spectrum 41 before modulation is narrow and that the frequency band is narrow. The spectrum 41 before modulation corresponds to the spectrum before the spectrum diffusion clock technique is adopted there. The reduction amount in FIG. 2 is shown with the amplitude before modulation set at 0 dB (decibel).


The spectrum 42 after modulation, however, expands in frequency band through frequency modulation and the electromagnetic interference is reduced.


In FIG. 2, the electromagnetic interference reduction caused by modulation is about 12 dB at the minimum.


The spectrum diffusion clock technique is a technique for reducing a level of the electromagnetic interference by using the frequency modulation, and this modulation may use a waveform other than a sine wave. The frequency spectrum component in the case of using a waveform other than a sine wave is different from that in FIG. 2, but it becomes the substantially same spectrum, and therefore, even in various kinds of modulations using a waveform other than a sine wave, the formula 1 may be used to calculate the reduction roughly.


As mentioned above, the processing of roughly calculating the electromagnetic interference reduction has been described.


A calculation method the electromagnetic interference reduction according to the embodiment will be described.


As mentioned in FIG. 1, the calculation method of the electromagnetic interference reduction according to the invention includes the first calibration of the electromagnetic interference reduction and the second calibration of the electromagnetic interference reduction.


At first, the first calibration of the electromagnetic interference reduction will be described.



FIG. 3 is a view showing one example of the frequency spectrum waveform in a line spectrum display for use in describing the embodiment.


In FIG. 3, since each time change and amplitude of a spectrum group 20 can be calculated according to sin(2p(f+nfm)t), Jn(Δf/fm) forming the formula 1 as mentioned above, the time change and amplitude of each line spectrum can be calculated easily.


That is, when the maximum frequency shift Δf and the modulation frequency fm are known and the order of the Bessel function is specified, amplitude can be required using the Bessel function. The time change can be calculated according to 2p(f+nfm)t of sin {2p(f+nfm)t}.


When the bandwidth of a band filter 22 is narrower than that of the modulation frequency 21, this calculation can be executed without problem because the spectrum detected by the measurement device is only one spectrum.


On the other hand, when it is wider than that of the modulation frequency 21, a plurality of spectra may pass the band filter 22 of the measurement device and enter the measurement device. In this case, the value of the vector sum of amplitudes of all spectra is a real measured value.


When a plurality of spectra pass the band filter 22 of the measurement device and enter the measurement device, calibration is performed while taking the influence on the measured value of the spectra into account, which is the first calibration of the electromagnetic interference reduction according to the invention.


Hereinafter, the first calibration of the electromagnetic interference reduction will be described.


At first, a spectrum having the minimum electromagnetic interference reduction, in other words, having the maximum amplitude Jn(Δf/fm) is selected.


When the measurement bandwidth of a measurement device is narrower than the bandwidth of the modulation frequency, a target spectrum for calculation is one and the absolute value (<1) of this amplitude becomes a reduction rate as it is and there is not need of calibration. Here, the amplitude before modulation is set at 1.


When the measurement bandwidth (RBW) of a measurement device is wider than the bandwidth of the modulation frequency, a target for calculation is a plurality of spectra. Therefore, calibration is required in the conventional calculation method of selecting only one spectrum having the maximum amplitude Jn(Δf/fm).


Specifically, all the amplitudes Jn(Δf/fm) of a plurality of adjacent spectra within the measurement bandwidth (RBW) are added together, and a value obtained by the absolute value of the largest one of the results is defined as the electromagnetic interference reduction.


According to the first calibration of the electromagnetic interference reduction, a level reduction rate is decreased more than in the general case of one spectrum and this decrease is caused by the calibration.


The frequency at this time is determined from the selected n easily and the frequency can be calculated according to (f+nfm) and (f−nfm).


The number of spectra Sn to be calculated is given by the formula 2. For example, when RBW=100 kHz and fm=100 kHz, spectra exist at the both ends of the measurement bandwidth and two spectra are targeted.


When fm=50 kHz, spectra exist at the both ends and in the center of the measurement bandwidth and three spectra are targeted.
Sn=Int(RBWfm)+1Formula(2)


Where, Int(x) indicates the integer after rounding x, and the RBW indicates the measurement bandwidth of the measurement device.


This time, a second calibration of the electromagnetic interference reduction shown in FIG. 1 will be described. This calibration is performed against sensibility reduction of spectrum caused by deviation of filter property from the ideal property in a measurement device.



FIG. 4 is a view showing an example of the filter property for use in describing the embodiment.


Referring to FIG. 4, the actual filter property 32 is deviated from an ideal filter property 31 represented by rectangular amplitude and it is curved, and as the result, it has a frequency band restriction. The actual filter property 32 can be seen by a spectrum analyzer that is a measurement device of the electromagnetic interference reduction.


As a concrete example, in a typical spectrum analyzer, a filter bandwidth 34 in which sensibility reduction 33 of −3 dB occurs, with respect to the sensibility in the center frequency is defined as a measurement bandwidth.


Referring to FIG. 4, in the calculation using the ideal filter property 31, four spectra may pass the filter without any influence on the amplitude, while in the calculation with the actual filter shape 32 taken into account, two spectra close to the both ends of the filter, of the four spectra, have influence on their amplitudes.


The second calibration is performed with this influence taken into account.


Specifically, after requiring the amplitude of each spectrum according to the formula 1 that is the conventional method, the spectra positioned at the both ends of the filter are calculated with calibration due to the corresponding sensibility reduction, which is the second calibration of the electromagnetic interference reduction.


An example of adopting the invention to a higher harmonic wave will be described.


The above description has been made in the case where an electromagnetic interference signal with the time change shown in sine (sinusoidal) wave is modulated by sine wave. While, in the modulation according to the spectrum diffusion clock technique, an electromagnetic interference signal is modulated by clock wave, that is, rectangle wave.


The rectangle wave is a composite wave of a basic sine wave and the higher harmonic wave, the electromagnetic interference reduction is important also in the higher harmonic wave, and the electromagnetic interference reduction calculation is required in the higher harmonic wave included in the rectangle wave.


In the case of the higher harmonic wave, it is basically the same as in the case of sine wave, and the same calculation is enabled by handling the variables concerned about the modulation of the N-order higher harmonic wave as follows.


When calculating the N-order higher harmonic wave of the sine wave that is a basic wave of clock, carrier frequency fN used for calculation as a carrier is represented as follows, using the clock frequency (=carrier frequency)f.

fN=Nf(Hz)  formula (3)


When calculating the N-order higher harmonic wave, the modulation frequency fmN used for calculation as a modulation frequency is represented as follows, using the modulation frequency fm for clock frequency as it is.

fmN=fm(Hz)  formula (4)


When calculating the N-order higher harmonic wave, the maximum frequency shift ΔfN used for calculation as a maximum frequency shift is represented as follows, using the maximum frequency shift Δf of clock frequency.

ΔfNNΔf (Hz)  formula (5)


In the case of the N-order higher harmonic wave, the number of spectra SnN handled here is represented as follows, using the modulation frequency fm for clock frequency and the measurement bandwidth of a measurement device.


That is, the same formula as the formula 2 is used.
SnN=Int(RBWfm)+1Formula(6)


Referring to the flow chart of FIG. 1, an example of adopting it to the higher harmonic wave will be specifically described. The corresponding step number in FIG. 1 is attached to the description corresponding to the flow chart of FIG. 1 at the end of the sentence in parentheses. A comparison between the calculation result in FIG. 1 and the measurement result is described.


At first, clock frequency (carrier frequency), modulation frequency, maximum frequency shift (=clock frequency×modulation degree (%)/100), and measurement bandwidth of a measurement device are supplied (Step 101).


Next, the electromagnetic interference reduction is roughly calculated. The input numeric values are applied to the formula 1, to calculate the amplitude given by the Bessel function. The absolute value of the maximum amplitude is searched for and the order n of the Bessel function Jn(Δf/fm) which gives the maximum amplitude is defined as nmax (Step 102).


Next, the bandwidth of a measurement device is compared with the modulation frequency using the formula 2 and the number of spectra within the measurement bandwidth of a measurement device is calculated.


According to the number of spectra calculated in the above, a combination having the maximum total sum of the coefficients of the amplitudes is selected as a combination of a group of adjacent spectra, from the vicinity of the order nmax of the Bessel function, for example, nmax−1, nmax−2, or nmax+1, nmax+2, . . . .


Concretely, when the number of spectra obtained by the formula 2 is two, one combination of the larger total sum is selected from the two combinations: the total sum of the absolute value of the Bessel function given by the order nmax and the absolute value of nmax−1, and the total sum of the absolute value of the Bessel function given by the order nmax and the absolute value of nmax+1.


The total sum obtained in this way corresponds to the first calibration of the electromagnetic interference reduction (Step 103).


A comparison between the result obtained thus and the measured value will be described.



FIG. 5 is a view showing a comparison between the calculated value and the measured value of the electromagnetic interference reduction according to the embodiment.


In FIG. 5, the measured values of the electromagnetic interference reduction effect are plotted when the carrier frequency is set at 50 MHz, the modulation degree is set 1.25%, and the modulation frequency is set at the horizontal axis as a variable. Measurement is performed in the case of RBW=3 kHz (indicated by ●) and in the case of RBW=100 kHz (indicated by ¦).


The calculation results “no calibration” (indicated by □) are obtained through a rough calculation of electromagnetic interference reduction according to the conventional technique described in Step 102 of FIG. 1, and they agree with the measurement values of RBW=3 kHz. This shows that calibration is not necessary when the measurement bandwidth is narrower than the modulation frequency fm.


According to the measurement results of RBW32 100 kHz, it is found that a large deviation occurs in the range where the modulation frequency fm is narrower than the measurement bandwidth and that this portion needs calibration.


When the first calibration of the electromagnetic interference reduction is performed (indication by □) according to the embodiment, it is found that they substantially agree with the calculation results until the modulation frequency becomes 60 kHz. However, a deviation from the measured value is found at 60 kHz and more.


Then, the second calibration of the electromagnetic interference reduction (indication by o) is performed according to the embodiment. The calibrated values agree with the measured values and it is found that the deviation found in the first calibration of the electromagnetic interference reduction is caused by the sensibility reduction in the ends of the filter bandwidth 34.


As mentioned above, by combining the first calibration of the electromagnetic interference reduction with the second calibration of the electromagnetic interference reduction according to the embodiment, it is possible to estimate the electromagnetic interference reduction accurately through calculation when the measurement bandwidth is wider than the modulation frequency.


A method of simplifying the calculation processing according to the embodiment will be described.


Though that the Bessel function Jn(Δf/fm) represents amplitude, the amplitude concerns the electromagnetic interference reduction level, and that the order n giving the maximum amplitude is searched for, has been described, a method of searching for the maximum value easier will be described.



FIG. 6 is a view showing the amplitude obtained by the Bessel function for use in describing the embodiment. The horizontal axis of FIG. 6 is the n-value of the Bessel function.


Referring to FIG. 6, convergence to zero occurs in the vicinity of n=Δf/fm in the Bessel function Jn(ΔF/Fm). Based on the property of the Bessel function, at first, the vicinity of n=Δf/fm is searched for.


The Bessel function becomes the maximum in the vicinity of n=Δf/fm.


As mentioned above, as a method requiring the maximum amplitude easier, n=Δf/fm is first calculated, and then, taking the Sn obtained by the formula 2 into account, amplitude may be calculated in the range of n=Δf/fm−Sn to n=Δf/fm+Sn.


In the above method of searching for the maximum value easier, the order giving the maximum amplitude can be found in the order of the Bessel function from n=Δf/fm−Sn to n=Δf/fm+Sn. In the range of the order, a combination of the maximum spectrum amplitude corresponding to the number of the spectra included in the measurement bandwidth can be found easier.


In the N-order higher harmonic wave, it may be obtained from n=Δf nf/m−NSn to n=Δf n/fm+NSn.


In this way, search for the maximum value can be simplified extremely.


In the above embodiment, the example of calculating electromagnetic interference reduction according to the given parameters has been described. The reduction amount of electromagnetic interference according to a combination of the parameters including modulation frequency, modulation degree, and the like can be obtained not through measurement but through only calculation and the optimum combination of parameters can be derived theoretically. This is also used as a designing means for determining a parameter of spectrum diffusion at a time of designing an electronic circuit.


A concrete operation example according to the embodiment will be described using the drawings.


A calculation example in the case of the following condition is shown.


Example of Condition



  • clock frequency: 50 MHz

  • modulation frequency: 100 kHz

  • modulation degree: 1%

  • measurement bandwidth (RBW) of a spectrum analyzer that is a measurement device: 100 kHz



Under the above conditions, each electromagnetic interference reduction in a basic wave and a second order harmonic wave is required. Referring to FIG. 1 the corresponding step number int he flow chart of FIG. 1 is attached to the corresponding description a the end of sentence in parentheses.


At first, the electromagnetic interference reduction of FIG. 1 is roughly estimated (Step 102).


When Jn(Δf/fm) showing the maximum electromagnetic interference reduction is calculated, the absolute value of the coefficient becomes the maximum value, 0.391, when n=4 (at the frequency of 49.6 MHz and 50.4). That is, in the vicinity of 50 MHz, the level reduction at the frequency of the least electromagnetic interference reduction is 0.391 times as much as that before the modulation.


This time, the first calibration of the electromagnetic interference reduction is performed (Step 103). According to the formula 2, the number of spectra Sn captured at once within the RBW of a measurement device is two.


Here, the method of simplifying the calculation processing according to the embodiment is used.


The respective coefficients indicated by J4−1(Δf/fm) and J4+1(Δf/fm) in the vicinity of n=4 are 0.365 and 0.261.


According to the result, it s found that the maximum amplitude taking the influence by the RBW of a measurement device into account is obtained by a combination of the spectrum of the coefficient 0.391 first obtained and the adjacent spectrum of the coefficient 0.365 (49.6 MHz and 49.7 MHz).


The second calibration of the electromagnetic interference reduction is performed (Step 104). In the second calibration, influence by the detection sensibility reduction of a measurement device is calibrated.



FIG. 7 is a view showing an example of the filter property of a measurement device for use in describing the embodiment.


Referring to FIG. 7, the detection sensibility of a measurement device is highest at the outer frequency and it gets lower when it is deviated from the center frequency.


The second calibration of the electromagnetic interference reduction described later is performed by using FIG. 7.


Since the spectra used for calculating the maximum amplitude are positioned at the both ends of the measurement bandwidth, each sensibility reduction in the measurement device is required from FIG. 7 and calibrated by −3 dB.


As the result, 0.391 and 0.365 become 0.277 and 0.258, 0.535 in total, and the electromagnetic interference reduction is 0.535 times as much as the level before the modulation.


With respect for the electromagnetic interference reduction in the second order harmonic wave at the frequency 100 MHz, the setting parameters are replaced with the parameters corresponding to the second order harmonic wave of a clock frequency according to the formula 3, the formula 4, the formula 5, and the formula 6 and calculated in the same way. As the result, the respective coefficients handled in the stage of the first calibration become 0.318 and 0.291, the total sum after the second calibration becomes 0.431, and in the second order harmonic wave, the electromagnetic interference can be reduced 0.431 times.


The above mentioned method of the embodiment is a calculation method of the electromagnetic interference reduction through the spectrum diffusion that is a means of reducing the electromagnetic interference. Of a plurality of spectra generated from the frequency spectrum of an electromagnetic interference signal according to the spectrum diffusion, the center frequency of the measurement bandwidth is changed so that the spectrum having the maximum amplitude may be included in the measurement bandwidth, the amplitudes of all the spectra included in the measurement bandwidth are added together, and the maximum total sum of the amplitudes is divided by the amplitude of the electromagnetic interference signal, hence to calculate the electromagnetic interference reduction.


A calculation device of electromagnetic interference reduction which provided the method of the invention will be described.



FIG. 8 is a flow chart for use in describing the processing performed by the electromagnetic interference reduction calculation device according to the embodiment. In FIG. 8, the right frames in the flow chart are to describe the corresponding steps roughly.


Referring to FIG. 8, numerical data is supplied at first by the input means (Step 801).


As the parameters necessary for calculation, clock frequency f, modulation frequency fm, maximum frequency shift Δf, and measurement bandwidth RBW of a measurement device are entered from an input unit such as a keyboard. When the maximum frequency shift has to be required from the modulation degree (%) given to a clock frequency, the modulation degree is entered, hence to require the maximum frequency shift from Δf=f×(modulation degree/100).


A temporary order of the Bessel function having the maximum amplitude is required by using the relationship n=Δf/fm (Step 802).


Here, the temporary order of the Bessel function indicates n. As mentioned above, the Bessel function gets the maximum amplitude in the order of the vicinity of the temporary order n.


Calculation is performed to find the Bessel function having the maximum amplitude from the order of n=Δf/fm−Sn to the order of n=Δf/fm+Sn, hence to determine the order giving the maximum amplitude (Step 803).


The amplitude corresponding to the order n giving the maximum amplitude and the amplitudes corresponding to the orders before and after n are added together, and a combination of the maximum total sum of the amplitudes is determined, hence to estimate the first calibration value for calibrating the amplitude of the order n (Step 804).


This time, the second calibration value is estimated at the both ends of the measurement bandwidth (Step 805).


The first calibration value and the second calibration value are added, hence to estimate the amplitude (Step 806).


By using thus calculated amplitude, the electromagnetic interference reduction is estimated and displayed (Step 807).


According to the above processing, it is possible to estimate the electromagnetic interference reduction when a user enters data.



FIG. 9 is a view showing the structure of an electromagnetic interference reduction calculation device according to the embodiment.


The electromagnetic interference reduction calculation device according to the embodiment comprises an input unit 110, a calculating unit 120, a storing unit 130 and an output unit 140.


The input unit 110 enters the setting parameters for spectrum diffusion. When calculating the electromagnetic interference reduction of clock waveform, it enters clock frequency, modulation frequency, maximum frequency shift, and measurement bandwidth.


The calculating unit 120 calculates electromagnetic interference reduction through spectrum diffusion.


At this stage, the temporary order of the Bessel function giving the maximum amplitude can be calculated and the order of the Bessel function giving the maximum amplitude can be calculated.


Calibration value in the case of a plurality of spectra can be calculated and the calibration values at the both ends of the measurement bandwidth can be calculated.


The calculating unit 120 calculates amplitude, hence to estimate the electromagnetic interference reduction.


The storing unit 130 stores the data necessary for calculating by the calculating unit 120.


The output unit 140 outputs the electromagnetic interference reduction calculated by the calculating unit 120.


According to the above mentioned embodiment, when a measurement device having a given measurement bandwidth calculates the electromagnetic interference reduction through spectrum diffusion, when the measurement bandwidth of the electromagnetic interference reduction which has to be calibrated by measurement is wider than the modulation frequency, the electromagnetic interference reduction is estimated from the maximum total sum of the amplitudes of all the spectra included in the measurement bandwidth, thereby making it unnecessary to calibrate the electromagnetic interference level, which has been required due to the influence to estimate the electromagnetic interference reduction only through calculation.


With respect for the electromagnetic interference reduction calculation device of the invention, its operation may be realized not only in a way of hardware but also in a way of software by running a calculation program (application) 200 for executing the above calculations, on a computer, that is, an electromagnetic interference reduction calculation device. The calculation program 200 is stored in a magnetic disk, a semiconductor memory, and the other recording medium, and loaded from the recording medium into the electromagnetic interference reduction calculation device to execute the operation, thereby estimating the above numerical information.


Second Embodiment

This time, a second embodiment of the invention will be described using the drawings.


An example of an electronic circuit with the maximum electromagnetic interference reduction will be described in this embodiment, based on the calculation method according to the first embodiment.


According to the calculation method described in the first embodiment, it is also possible to set the spectrum diffusion of an electronic circuit, using the respective values of the clock frequency, the maximum frequency shift, the measurement bandwidth, and the modulation frequency, after estimating the modulation frequency giving the maximum electromagnetic interference reduction, when the clock frequency, the maximum frequency shift, and the measurement bandwidth are constant.


Referring to FIG. 5, it is found that when the influence by the measurement bandwidth of a measurement device is not taken into account, the calculated value (graph of RBW=3 kHz, the value indicated by ●) of the reduction effect decreases according as the modulation frequency increases.


While, when the calibration is necessary, the calculated value (graph of RBW=100 kHz, the value indicated by ¦) of the reduction effect increases in the effect according as the modulation frequency increases.


These two graphs approach at the frequency where the modulation frequency becomes equal to the measurement bandwidth of a measurement device. An intersection of these two graphs depends on the filter property as mentioned below.


The intersection of the two graphs appears at fm=RBW when an ideal rectangular filter property is obtained such that the sensibility reduction rapidly falls down in the band exceeding the measurement bandwidth. While, in the case of the filter property deviated from the actual rectangle as shown in FIG. 7, the intersection is deviated to the right side, that is, toward the larger modulation frequency fm.


According to the result, it is found that the modulation frequency fm having the maximum reduction effect of the electromagnetic interference appears when it becomes substantially equal to the frequency of the measurement bandwidth of a measurement device. By setting the modulation frequency fm at the same value of RBW, the substantially maximum effect can be obtained.


When a measurement bandwidth is determined in the regular standards of a target measurement (in CISPR22 of Comittie International Special des Perturbations Radioelectrique and in VCCI standard based on this international standards), as mentioned in Background Art, it is preferable that, by using this effect, the modulation frequency is set at 9 kHz according to the spectrum diffusion technique when the maximum reduction effect is expected in the frequency range from 150 kHz to 30 MHz, the modulation frequency is set at 120 kHz when the maximum reduction effect is expected in the frequency range from 30 MHz to 1 GHz, and the modulation frequency is set at 1 MHz when the maximum reduction effect is expected in the frequency exceeding 1 GHz.


Due to the influence of the above filter property, the maximum reduction effect of the electromagnetic interference can be obtained when the modulation frequency fm exceeds the measurement bandwidth. Taking the filter property into account, the modulation frequency fm having the maximum reduction effect of the electromagnetic interference is in the range from the measurement bandwidth to the frequency exceeding the measurement bandwidth.



FIG. 10 is a view showing an interference of spectrum diffusion signal 320 from an electronic circuit 310 to a cellular phone 330, according to the embodiment.


In order to minimize the interference from the electronic circuit 310 to the cellular phone 330, the modulation frequency fm of the spectrum diffusion signal 320 has to be equal to the exclusive bandwidth of the cellular phone, 12.5 kHz or 25 kHz. When measuring the electromagnetic interference, it is defined that the measurement bandwidth has to be equal to the exclusive bandwidth of a device to be interfered.


As mentioned above, although the invention has been described using the preferred embodiments, the invention is not always restricted to the above embodiments but various modifications can be made within the scope of the technical spirit.


According to the electromagnetic interference reduction calculation method, device, calculation program, and its electronic circuit of the invention, the following effects can be achieved.


First, in the calculation of the electromagnetic interference reduction through the spectrum diffusion, which reduction is measured by a measurement device having a given measurement bandwidth, when the measurement bandwidth of the electromagnetic interference reduction requiring the calibration through measurement is wider than the modulation frequency, the electromagnetic interference reduction is estimated from the maximum total sum of the amplitudes of all the spectra included in the measurement bandwidth, thereby making the calibration through measurement unnecessary and making it possible to estimate the electromagnetic interference reduction only through the calculation.


Second, when adding the amplitudes of all the spectra included in the measurement bandwidth, since the influence of the sensibility reduction of a measurement device is calibrated on the spectra in the vicinity of the frequency at the ends of the measurement bandwidth, the electromagnetic interference reduction can be estimated through calculation at higher precision.


Third, since the spectrum diffusion is performed on an electronic circuit under the condition that the estimated electromagnetic interference reduction gets the maximum value, an electronic circuit device having the maximum electromagnetic interference can be manufactured.

Claims
  • 1. A method of calculating electromagnetic interference reduction through spectrum diffusion that is a means of reducing electromagnetic interference, comprising adding up amplitudes of all spectra included in a measurement bandwidth when a center frequency of said measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in said measurement bandwidth, of a plurality of spectra generated from frequency spectrum of an electromagnetic interference signal through said spectrum diffusion, and estimating electromagnetic interference reduction by dividing the maximum total sum of said amplitudes by the amplitude of said electromagnetic interference signal.
  • 2. The calculation method of electromagnetic interference reduction according to claim 1, in which when said measurement bandwidth of the electromagnetic interference reduction is wider than a modulation frequency that is a frequency fluctuation width set through said spectrum diffusion, said amplitudes of the spectra are added up.
  • 3. The calculation method of electromagnetic interference reduction according to claim 1 or claims 2, in which when adding up the amplitudes of all the spectra included in said measurement bandwidth, the spectra in the vicinity of frequencies at both ends within the said measurement bandwidth is calibrated considering influence of sensibility reduction of a measurement device in said end frequencies.
  • 4. The calculation method of electromagnetic interference reduction according to one of claim 1 to claim 3, in which when calculating a plurality of spectra generated from the frequency spectrum of said electromagnetic interference signal by Bessel function, calculation of the Bessel function is performed on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting a value obtained by adding 1 to an integral part of a value obtained by dividing said bandwidth by said modulation frequency from a value obtained by dividing a maximum frequency shift that is an expansion of frequencies set through said spectrum diffusion by said modulation frequency, to the order of the Bessel function corresponding to the value obtained by adding 1 to said integral part of said value obtained by dividing said maximum frequency shift by said modulation frequency, and the amplitude of the spectrum having the maximum amplitude is calculated from the amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 5. The calculation method of electromagnetic interference reduction according to claim 4, in which the maximum total sum of said amplitudes is calculated by using the amplitudes of said different orders of the Bessel functions.
  • 6. The calculation method of electromagnetic interference reduction according to one of claim 1 to claim 3, in which when calculating a spectrum of harmonic wave using the Bessel function, which harmonic wave is obtained by performing the integral multiple of the frequency of a basic wave included in said electromagnetic interference signal generated from the frequency spectrum of said electromagnetic interference signal. calculation of the Bessel function is performed on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting said integral multiple value of a value obtained by adding 1 to the integral part of the value obtained by dividing said measurement bandwidth by said modulation frequency from a value obtained by dividing by said modulation frequency, said integral multiple value of the maximum frequency shift that is the expansion of the frequencies set through said spectrum diffusion of said basic wave, to the order of the Bessel function corresponding to the value obtained by adding said integral multiple value to the value obtained by dividing said integral multiple of said maximum frequency shift by said modulation frequency, and the amplitude of the spectrum having the maximum amplitude is calculated from said amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 7. The calculation method of electromagnetic interference reduction according to claim 1 or claim 2, in which when generating a plurality of spectra from the frequency spectrum of said electromagnetic interference signal through said spectrum diffusion, clock frequency, modulation frequency, maximum frequency shift, and measurement bandwidth are used as input data to set said spectrum diffusion.
  • 8. A device for calculating electromagnetic interference reduction through spectrum diffusion that is a means of reducing electromagnetic interference, comprising an input unit which enters a setting condition of the spectrum diffusion, a calculating unit which calculates electromagnetic interference reduction through the spectrum diffusion, and an output unit which supplies calculation result of said electromagnetic interference reduction, in which when a center frequency of a measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in said measurement bandwidth, of a plurality of spectra generated from a frequency spectrum of an electromagnetic interference signal through said spectrum diffusion, said calculating unit adds up the amplitudes of all spectra included in said measurement bandwidth and calculates said electromagnetic interference reduction by dividing the maximum total sum of said amplitudes by the amplitude of said electromagnetic interference signal.
  • 9. The calculation device of electromagnetic interference reduction according to claim 8, in which when said measurement bandwidth of the electromagnetic interference reduction is wider than a modulation frequency that is a frequency fluctuation width set through said spectrum diffusion, said calculating unit adds up the amplitudes of said spectra.
  • 10. The calculation device of electromagnetic interference reduction according to claim 8 or claim 9, in which when adding up the amplitudes of all the spectra included in said measurement bandwidth, said calculating unit calibrates the spectra in the vicinity of frequencies at both ends within the said measurement bandwidth, considering influence of sensibility reduction of a measurement device in said end frequencies.
  • 11. The calculation device of an electromagnetic interference reduction according to one of claim 8 to claim 10, in which when calculating a plurality of spectra generated from the frequency spectrum of said electromagnetic interference signal by Bessel function, said calculating unit calculates the Bessel function on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting a value obtained by adding 1 to an integral part of a value obtained by dividing said bandwidth by said modulation frequency from a value obtained by dividing a maximum frequency shift that is an expansion of frequencies set through said spectrum diffusion by said modulation frequency, to the order of the Bessel function corresponding to the value obtained by adding 1 to said integral part of said value obtained by dividing said maximum frequency shift by said modulation frequency, and calculating the amplitude of the spectrum having the maximum amplitude from said amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 12. The calculation device of electromagnetic interference reduction according to claim 11, in which said calculating unit calculates the maximum total sum of said amplitudes by using the amplitudes of said plurality of different orders of the Bessel functions.
  • 13. The calculation device of electromagnetic interference reduction according to one of claim 8 to claim 10, in which when calculating a spectrum of harmonic wave using the Bessel function, which harmonic wave is obtained by performing the integral multiple of the frequency of a basic wave included in said electromagnetic interference signal generated from the frequency spectrum of said electromagnetic interference signal. said calculating unit calculates the Bessel function on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting said integral multiple value of a value obtained by adding 1 to the integral part of the value obtained by dividing said measurement bandwidth by said modulation frequency from a value obtained by dividing by said modulation frequency, said integral multiple value of the maximum frequency shift that is the expansion of the frequencies set through said spectrum diffusion of said basic wave, to the order of the Bessel function corresponding to the value obtained by adding said integral multiple value to the value obtained by dividing said integral multiple of said maximum frequency shift by said modulation frequency, and calculates the amplitude of the spectrum having the maximum amplitude, from said amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 14. The calculation device of electromagnetic interference reduction according to one of claim 8 or claim 13, in which said input unit enters basic frequency, modulation frequency, maximum frequency shift, and measurement bandwidth of an electromagnetic interference signal.
  • 15. A program run on an electromagnetic interference reduction calculation device which calculates electromagnetic interference reduction through spectrum diffusion that is a means of reducing the electromagnetic interference, comprising a function of adding up amplitudes of all spectra included in a measurement bandwidth when a center frequency of said measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in said measurement bandwidth, of a plurality of spectra generated from frequency spectrum of an electromagnetic interference signal through said spectrum diffusion, and a function of estimating electromagnetic interference reduction by dividing the maximum total sum of said amplitudes by the amplitude of said electromagnetic interference signal.
  • 16. The calculation program according to claim 15, comprising a function of adding up the amplitudes of said spectra when said measurement bandwidth of the electromagnetic interference reduction is wider than a modulation frequency that is a frequency fluctuation width set through said spectrum diffusion.
  • 17. The calculation program according to claim 15 or claim 16, comprising a function of calibrating the spectra in the vicinity of frequencies at both ends within the said measurement bandwidth, considering influence of sensibility reduction of a measurement device in said end frequencies when adding up the amplitudes of all the spectra included in said measurement bandwidth.
  • 18. The calculation program according to one of claim 15 or claim 17, comprising a function of, when calculating a plurality of spectra generated from the frequency spectrum of said electromagnetic interference signal by Bessel function, calculating the Bessel function on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting a value obtained by adding 1to an integral part of a value obtained by dividing said bandwidth by said modulation frequency from a value obtained by dividing a maximum frequency shift that is an expansion of frequencies set through said spectrum diffusion by said modulation frequency, to the order of the Bessel function corresponding to the value obtained by adding 1 to said integral part of said value obtained by dividing said maximum frequency shift by said modulation frequency, and calculating the amplitude of the spectrum having the maximum amplitude from the amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 19. The calculation program according to claim 18, comprising a function of calculating the maximum total sum of said amplitudes by using the amplitudes of said different orders of the Bessel functions.
  • 20. The calculation program according to one of claim 15 to claim 17, comprising when calculating a spectrum of harmonic wave using the Bessel function, which harmonic wave is obtained by performing the integral multiple of the frequency of a basic wave included in said electromagnetic interference signal generated from the frequency spectrum of said electromagnetic interference signal, a function of calculating the Bessel function on a plurality of different orders, from the order of the Bessel function corresponding to a value which is obtained by subtracting said integral multiple value of a value obtained by adding 1 to the integral part of the value obtained by dividing said measurement bandwidth by said modulation frequency from a value obtained by dividing by said modulation frequency, said integral multiple value of the maximum frequency shift that is the expansion of the frequencies set through said spectrum diffusion of said basic wave, to the order of the Bessel function corresponding to the value obtained by adding said integral multiple value to the value obtained by dividing said integral multiple of said maximum frequency shift by said modulation frequency, and calculating the amplitude of the spectrum having the maximum amplitude from said amplitude of the order of the Bessel function having the maximum amplitude, of the plurality of different orders of the Bessel functions.
  • 21. An electronic circuit using a spectrum diffusion clock with a clock signal spectrum diffused, in which according to a calculation method of calculating electromagnetic interference reduction, comprising: adding up the amplitudes of all spectra included in a measurement bandwidth, when a center frequency of said measurement bandwidth is changed so that a spectrum having the maximum amplitude is included in said measurement bandwidth, of a plurality of spectra generated from frequency spectrum of an electromagnetic interference signal through said spectrum diffusion; and dividing the maximum total sum of said amplitudes by the amplitude of said electromagnetic interference signal, a modulation frequency having the maximum electromagnetic interference reduction is calculated when clock frequency, maximum frequency shift, and measurement bandwidth are respectively constant, and said spectrum diffusion is set based on the respective values of said clock frequency, said maximum frequency shift, said measurement bandwidth, and said modulation frequency.
  • 22. the electronic circuit according to claim 21, in which when the measurement bandwidth of said electromagnetic interference reduction is wider than the modulation frequency that is a frequency fluctuation width set through said spectrum diffusion, the amplitudes of said spectra are added up.
  • 23. The electronic circuit, according to claim 21 or claim 22, in which according to the result, derived from said calculation method, that said electromagnetic interference reduction becomes the maximum when said modulation frequency is equal to said measurement bandwidth and that said modulation frequency having the maximum electromagnetic interference reduction becomes the frequency exceeding said measurement bandwidth owing to influence of sensibility reduction of a measurement device in the frequencies at both ends of said measurement bandwidth, the minimum of said modulation frequency is defined as said measurement bandwidth, the maximum of said modulation frequency is defined as a constant frequency exceeding said measurement bandwidth, and said modulation frequency is set in a range of said minimum modulation frequency to said maximum modulation frequency.
  • 24. The electronic circuit according to claim 23, in which based on a regulation of measurement bandwidth in each clock frequency range, in the standard “CISPR22” of the Comittie International Special des Perturbations Radioelectrique, when the clock frequency is in the range from not less than 150 kHz to less than 30 MHz, the minimum of said modulation frequency is 9 kHz and the maximum of said modulation frequency is 11.7 kHz, when said clock frequency is in the range from not less than 30 MHz to less than 1 GHz, the minimum of said modulation frequency is 120 kHz and the maximum of said modulation frequency is 156 kHz, and when said clock frequency is 1 GHz and more, the minimum of said modulation frequency is 1 MHz and the maximum of said modulation frequency is 1.3 MHz.
  • 25. The electronic circuit according to claim 23, in which based on setting the measurement bandwidth equal to an exclusive bandwidth of a device to be interrupted, when measuring an electromagnetic interference, said modulation frequency is set in a vicinity of 12.5 kHz or 25 kHz that is the exclusive bandwidth of a cellular phone.
Priority Claims (1)
Number Date Country Kind
2004240651 Aug 2004 JP national
PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/JP05/15593 8/22/2005 WO 3/16/2007