Abstract:
In radio communication system using angular modulation, the presence and quantity of cochannel interference is detected. According to the invention, when desired signal and undesired signal are received simultaneously, an envelope of the received IF output signal involves low frequency spectrum component similar to that of fading and high frequency spectrum component with several KHz. The low frequency spectrum component and high frequency spectrum component are then calculated, and the interference ratio D/U is obtained by processing said components.

Description:
BACKGROUND OF THE INVENTION 
     The present invention relates to a cochannel interference detection system in radio communicatins. The present invention is applicable to a mobile communication, and/or cordless telephone and other radio communication systems. The present invention is applicable to the radio communication system using an angle modulation system (FM or PM). 
     One of the prior systems for detecting cochannel interference uses a pilot signal for detecting interference. In this system, a transmitter transmits power modulated by a pilot frequency which is set outside of the baseband, and the frequency of the pilot is different for each zone. When a receiver receives a pilot signal which is assigned to another zone, it is judged that there is an interference. However, that system has the disadvantages that the baseband spectrum is spread, the structure of the system is complicated (a transmitter must have a pilot signal generator, and a receiver must have a pilot signal receiver), and it can only detect an interference when D/U (ratio of desired signal to undesired signal) is smaller than several dB. 
     Additional prior interference detection equipment is shown in U.S. Pat. No. 4,085,370 which uses noise in a baseband generated by interference. This system can operate correctly only when the reception level does not widely vary and the modulation index is almost constant. However, it can not detect the average of the D/U ratio. 
     Further, prior interference detection systems have the disadvantage that they can merely detect the presence of interference, they can not measure the quantity of interference. Therefore, when communication quality becomes poor becouse of interference, the automatic switching has been almost impossible. 
     SUMMARY OF THE INVENTION 
     It is an object, therefore, of the present invention to overcome the disadvantages and limitations of a prior interference detection system by providing a new and improved interference detectin system. 
     It is also an object of the present invention to provide an interference detection system which can measure interference during actual communication without taking a particular test period. 
     It is also an object of the present invention to provide an interference detection system which can measure the quantity of interference. 
     The above and other objects are attained by a cochannel interference measurement system in radio communication using angular modulation comprising; envelope detection means for envelope detection of both first signal and second signal; low frequency spectrum component detection means for detecting a low frequency spectrum component having the similar spectrum as that of fading spectrum frequency, and said low frequency spectrum component relating to sum of first signal and second signal; high frequency spectrum component detection means for detecting a high frequency spectrum component higher than a fading spectrum frequency, said high frequency spectrum component relating to product of first signal and second signal; means for providing ratio of said low frequency spectrum component to said high frequency spectrum component; and indication means for indicating interference by ratio of first signal to second signal according to said ratio of said low frequency spectrum component to said high frequency spectrum component. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The foregoing and other objects, features, and attendant advantages of the present invention will be appreciated as the same become better understood by means of the following description and accompanying drawings wherein; 
     FIG. 1 shows waveforms of output of an envelope detector, 
     FIG. 2 shows curves of the relation between f m  Δt and Z 2  /Y, 
     FIG. 3 shows curves of the error of the measured interference when the delay time Δt is large, 
     FIG. 4 shows curves between the delay time Δt and the average of cos(ψ)cos(ψ.sub.Δ), 
     FIG. 5 shows curves of the error of the measured interference when the delay time Δt is small, 
     FIG. 6 shows the block diagram of an embodiment of the interference detection system according to the present invention, 
     FIG. 7A is a flow chart showing the operation of the apparatus of FIG. 6, 
     FIG. 7B is a detailed block diagram of the calculator 15 of FIG. 6, 
     FIG. 8 is a block diagram of another embodiment of the interference detection system according to the present invention, 
     FIG. 9 is a block diagram of still another embodiment of the interference detection system according to the present invention, 
     FIG. 10 shows waveforms of the apparatus of FIG. 9, and 
     FIG. 11 shows the experimental result of the apparatus of FIG. 9. 
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     The theoretical principle of the present invention is first described. 
     It is supposed that a first desired signal e 1  which is angular modulated, and a second undesired signal e 2  which is also angular modulated are expressed by the equations (1) and (2), respectively. 
     Desired signal; 
     
         e.sub.1 =E.sub.1 (t) sin(ω.sub.1 t+(Δω.sub.1)/(P.sub.1) sin P.sub.1 t)                                            (1) 
    
     Undesired signal; 
     
         e.sub.2 =E.sub.2 (t) sin(ω.sub.2 t+φ+((Δω.sub.2)/(P.sub.2))sin(P.sub.2 t+θ)) (2) 
    
     where E 1  (t) and E 2  (t) are the amplitudes of the desired signal and the undesired signal, respectively, and it is assumed that they are in Rayleigh distribution. The ω 1  and ω 2  are the angular velocity of the desired wave and the undesired wave, respectively, and Δω 1  and Δω 2  are the frequency shifts of the desired wave and the undesired wave by a modulation signal, respectively. P 1  and the P 2  are the frequencies of the modulation signals of the desired signal and the undesired signal, respectively. φ is the phase difference between the desired wave and the undesired wave, and θ is the phase difference between the desired modulation signal and the undesired modulation signal. 
     When both the desired wave and the undesired wave are received by a receiver, the combined signal e is expressed by the equation (3). ##EQU1## 
     The envelope output of the square-detector of e is expressed by the equation (6). 
     
         R(t)=E.sub.1.sup.2 (t)+E.sub.2.sup.2 (t)+2E.sub.1 (t)E.sub.2 (t)cos(ψ(t))                                          (6) 
    
     FIG. 1 shows an example of the envelope output R(t), in which the curve (1) shows the envelope R(t), and the curve (2) shows the value E 1   2  (t)+E 2   2  (t) which is the low frequency spectrum component of the envelope output R(t). It nearly equals to the fading frequency which depends upon the velocity of a mobile, and the wavelength of electromagnetic wave. When a mobile moves at 40 Km/hour, and the carrier frequency is 900 MHz, said low frequency component is about 30 Hz. The 2E 1  (t)E 2  (t)cos(ψ(t)) in the equation (6) has the high frequency spectrum higher than Rayleigh fading frequency. That high frequency component is generated by the angular modulation, or the small difference of the carrier frequencies between the desired signal and the undesired signal. 
     According to the preferred embodiment of the present invention, the average of the low frequency component &lt;E 1   2  (t)+E 2   2  (t)&gt; and the average of the high frequency spectrum component &lt;E 1  (t)E 2  (t)&gt; are measured by using an analog-digital converter, and then, the D/U is calculated by using said two average values. 
     The analog-digital converter samples the input signal at the time (t) and the time (t+Δt) as shown in FIG. 1. The sampled values R(t) and R(t+Δt) at the time (t) and (t+Δt) are expressed by the equation (7). 
     
         R(t)=E.sub.1.sup.2 +E.sub.2.sup.2 +2E.sub.1 E.sub.2 cos(ψ) R(t+Δt)=E.sub.1.sup.2.sub.Δ +E.sub.2.sup.2.sub.Δ +2E.sub.2Δ E.sub.2Δ cos(ψ.sub.Δ)    (7) 
    
     First, the average of low frequency component value is calculated as follows. It is assumed that the average &lt;cos(ψ)=0&gt; when the reception signal is angular modulated and/or the sampling period is at random. Therefore, the average of the low frequency spectrum component value is obtained by taking a plurality of samples, as follows. ##EQU2## where &lt; &gt; shows the average, and N is sample number. 
     Next, the average of high frequency spectrum component value Y is explained. The equation (9) is expressed by taking the average of the square of the difference between R(t) and R(t+Δt). ##EQU3## 
     The sampling delay time Δt is selected so that the equation (10) is satisfied for E 1  and E 2  which vary according to Rayleigh fading, and the equation (11) is satisfied for cos(ψ), and cos(ψ.sub.Δ) which varies quickly. 
     
         E.sub.1 ≈E.sub.1Δ, E.sub.2 ≈E.sub.2Δ(10) 
    
     
         &lt;cos(ψ)cos(ψ.sub.Δ)&gt;≈0               (11) 
    
     The equation (11) means that cos(ψ) has no correlation with cos(ψ.sub.Δ). Then, the equation (9) is changed to the equation (12). 
     
         &lt;[R(t)-R(t+Δt)].sup.2 &gt;=4&lt;E.sub.1.sup.2 &gt;&lt;E.sub.2.sup.2 &gt;=Y (12) 
    
     Setting 
     
         (Γ)=(D/U)=&lt;E.sub.1.sup.2 &gt;/&lt;E.sub.2.sup.2 &gt;          (13) 
    
     the equation (14) is derived from the equations (8), (12) and (13). 
     
         (Γ).sup.2 -2K(Γ)+1=0                           (14) 
    
     Then, ##EQU4## where 
     
         K=(2X.sup.2 /Y)-1=(&lt;E.sub.1.sup.2 &gt;+&lt;E.sub.2.sup.2 &gt;)/(2&lt;E.sub.1.sup.2 &gt;&lt;E.sub.2.sup.2 &gt;)                                        (16) 
    
     Since K is equal to or larger than 1, the value of the equation (15) is a real number, from which may be calculated the interference value (Γ). 
     Next, the delay time Δt is considered. 
     It should be noted that the equation (14) is met under the condition that the equations (10) and (11) are satisfied. Therefore, when those equations (10) and (11) are not satisfied, the interference value (Γ) has an error, which is analyzed below. 
     It is assumed that the following conditions are satisfied, and it should be noted that those conditions are naturally satisfied in a mobile communication. ##EQU5## 
     It is quite clear from the equation (8) that X=&lt;R(t)&gt; does not depend upon Δt. 
     The value Y depends upon Δt, and the value Y&#39; with the influence of Δt is expressed by the equation (18). ##EQU6## In comparing Y expressed by the equation (12) with no influence of Δt with Y&#39; expressed by the equation (18) with the influence of Δt, it should be noted that Y&#39; has the additional terms &lt;Z 2  &gt;, and -8&lt;E 1   2  &gt;&lt;E 2   2  &gt;&lt;cos(ψ)cos(ψ.sub.Δ)&gt;. 
     (a) When Δt is large; 
     In this case, &lt;cos(ψ)cos(ψ.sub.Δ)&gt;≈0 is satisfied, and the equation (18) becomes; 
     
         Y&#39;=&lt;Z.sup.2 &gt;+4&lt;E.sub.1.sup.2 &gt;&lt;E.sub.2.sup.2 &gt;            (19) 
    
     where &lt;Z 2  &gt; is expressed by the equation (20). 
     
         &lt;Z.sup.2 &gt;=8b.sub.2.sup.2 [Γ.sup.2 (1-ρ.sub.1.sup.2)+(1-ρ.sub.2.sup.2)]              (20) 
    
     where b 2   2  =&lt;E 2   4  &gt;/8, and (ρ) is an auto correlation as follows. 
     
         (ρ)=J.sub.0 (2πf.sub.m Δt)                    (21) 
    
     where f m  is the fading frequency (f m  =v(velocity)/λ(wavelength), and J 0  is the Bessel function of zero&#39;th order. The equation (21) shows that Y is influenced by the product of Δt and f m . 
     FIG. 2 shows the curves of &lt;Z 2  &gt;/Y. The influence of f m  Δt depends upon the value (Γ), and the larger the value (Γ) is, the larger the influence is. 
     FIG. 3 shows the curves of the influence of f m  Δt to (Γ), in which the horizontal axis shows the true (set) value of D/U=(Γ t ), and the vertical axis shows the measured value D/U=(Γ M ). The solid lines in FIG. 3 are theoretically calculated value, and the dots show the measured values. It should be noted that the ideal line in FIG. 3 is the 45° line which has no error. The value (Γ M ) becomes theoretically and experimentally smaller than the value (Γ t ) when f m  Δt is large, and the error is large when (Γ t ) is large. It should be appreciated in FIG. 3 that the value f m  Δt must be less than 0.03 because the error becomes less than 1 dB when (Γ t ) is smaller than 15 dB. In FIG. 3, ρ 1  =ρ 2  is assumed. 
     (b) When Δt is small; 
     In this case, the equation (10) is satisfied, and Y&#39; becomes 
     
         Y&#39;=4&lt;E.sub.2.sup.2 &gt;&lt;E.sub.2.sup.2 &gt;(1-2&lt;cos(ψ)cos(ψ.sub.66)&gt;) (22) 
    
     where the value &lt;cos(ψ)cos(ψ.sub.Δ)&gt; is expressed by the equation (23). ##EQU7## where J shows the Bessel function, ##EQU8## 
     The second term and the third term become effective only when nP 2  =mP 1 , and ((2n+1)P 2  =(2m+1)P 1 , respectively are satisfied by proper integers m and n (m=n≈0). 
     FIG. 4 shows the relation between 2&lt;cos(ψ)cos(ψ 66 )&gt; (vertical axis), and the delay time Δt (horizontal axis), in which Δω 1  =1000 Hz, Δω 2  =1010 Hz, the modulation signal is a sine wave signal, P 1  =200 Hz, and P 2  is in the range between 195 Hz and 205 Hz. It should be noted in FIG. 4 that the value 2&lt;cos(ψ)cos(ψ.sub.Δ)&gt; decreases with the increase of Δt, but it does not depend upon modulation frequency P 2 . The value becomes nearly zero when Δt=0.4 msec, and after that, it decreases while vibrating. 
     It should be noted that 2&lt;cos(ψ)cos(ψ.sub.Δ)&gt; is substantially zero when Δt is larger than 0.4 msec, therefore, it is preferably that Δt is larger than 0.4 msec. 
     The value 2&lt;cos(ψ)cos(ψ.sub.Δ)&gt; shows the decrease from the value Y when Δt is small, as shown in the equation (18). 
     FIG. 5 shows the influence of Δt to the value (Γ M ) under the conditions of FIG. 4. When the value Δt is small, the value (Γ M ) becomes large, and thus the influence is large when (Γ t ) is large. The reason why (Γ M ) becomes large, is that Y becomes small when Δt is small. 
     FIG. 6 shows a block diagram of an embodiment of the present invention, in which the upper portion shows the experimental unit for testing the present apparatus, and the lower portion 90 enclosed by a dotted line is the present interference detector. In the figure, reference numeral 3 is an oscillator for the generating a desired signal (D), which is angular-modulated by the modulation signal on the terminal 3a. Oscillator 6 generates undesired signal (U), which is also angular-modulated by the modulation signal on the terminal 6a. Fading simulators 4 and 7 provide artificial fading to both the desired signal (D) and the undesired signal (U), and 5 and 8 are variable attenuators. Combiner 9 combines both the desired signal (D) and the undesired signal (U), and 10 is a receiver which provides the intermediate frequency output (IF OUT). The IF OUT signal is applied to the present interference detector 90, which has an amplifier 11, an envelope detector 12, a pair of analog-digital converters 13 and 14, a processor 15 and a delay circuit 22. Preferably, the amplifier 11 is a logarithmic amplifier, and the corresponding logarithmic expander 15a is inserted at the input portion of the processor 15 so that signal processed by the envelope detector 12 and converter 13 and/or 14 has a small dynamic range. 
     FIG. 7A shows a flow-chart of the processor 15 when the processor 15 is implemented by a programmed computer, and FIG. 7B shows a block diagram of the processor 15 when the processor is implemented by hardware. First, the sampling number N is set in the counter 15b (FIG. 7B, box 100 in FIG. 7A). The number N is usually 500-1000 in case of Rayleigh fading. The delay time Δt is also set and the value Δt is applied to the delay circuit 22 which delays the sampling trigger signal by Δt. Then, the first A/D converter 13 provides the signal value R i  (t) sampled at the time (t), and the second converter 14 provides the signal value R i  (t+Δt) sampled at the time (t+Δt). The accumulators (15c, 15d in FIG. 7B, box 108 in FIG. 7A) provide the sum of the sampled values, and the square of the difference of the two samples (R i  (t) and R i  (t+Δt)), respectively. It should be noted that the accumulators 15c and 15d, or the box 108 in FIG. 7A are reset to zero in each initial condition (i=0). 
     When the above accumulation is repeated by N times, the division is carried out by the divider 15e in FIG. 7B, or the box 112 in FIG. 7A, which divides the input value by N to provide the value X and Y in the previous equations (8) and (12). Then, the calculator 15f in FIG. 7B or the box 114 in FIG. 7A provides the value k (=(2X 2  /Y-1)) according to the equation (16). 
     The calculator 15g in FIG. 7B or the box 116 in FIG. 7A calculates the value (k 2  -1) in order to assure that said value (k 2  -1) is positive. If that value (k 2  -1) is negative, the present system is reinitialized, and the above operation is carried out again. When said value (k 2  -1) is zero or positive, the calculator 15h in FIG. 7B or the box 120 in FIG. 7A calculates the value ##EQU9## which is the resultant interference. 
     FIG. 8 shows another embodiment of the present invention, in which the calculation is accomplished by hardware without using a programmed means. In the figure, the reference numeral 20 is an antenna, 10 is a receiver, 12 is an envelope detector, 23 is a smoother circuit, 24 is a squaring circuit, 25 is a divider circuit, 26 is a level meter, 27 is a delay circuit, 28 is a differential amplifier, 29 is a square circuit, and 30 is a smoother circuit. 
     In FIG. 8, both the desired signal (D) and the undesired signal (U) enter into the receiver 10 through the antenna 20. At the output of the receiver 10, the intermediate frequency output (IF out ) is applied to the envelope detector 12. The output of the detector 12 is applied to the smoother 23, which, then, flattens the high frequency component of the envelope detected signal, therefore, that smoother 23 provides the low frequency component X expressed by the equation (8). Then, the squaring circuit 24 provides the square X 2 , as follows. ##EQU10## The output of the detector 12 is also applied to the delay circuit 27, which delays the signal by the time Δt. The delayed signal and the non-delayed signal are applied to the differential amplifier 28, the output of which is applied to the squaring circuit 29. The output of the squaring circuit 29 is smoothed by the smoother 30. Accordingly, the output of the smoother 30 is the high frequency component Y of the equation (12). The output of the divider 25 is the ratio of the low frequency component X 2  and the high frequency component Y as shown in the equation (25). ##EQU11## Therefore, the level meter 26 indicates the value Z, which then shows the ratio D/U=&lt;E 1   2  &gt;/&lt;E 2   2  &gt; by plotting the scale of the meter 26 properly. 
     FIG. 9 shows still another embodiment of the present invention, in which the reference numeral 20 is an antenna, 10 is a receiver, 43 is an AGC (Automatic Gain Control) circuit, 44 is an envelope detector, 45 is an amplifier, 46 is a high-pass filter, 47 is a detector, 48 is a low-pass filter, 49 is a smoother circuit, 50 is another smoother circuit, and 51 is a level meter. FIG. 10 shows the waveforms of the apparatus of FIG. 9, and 62 is the waveform of the output of the detector 44, 63 is the waveform of the high-pass filter 46, 64 is the waveform of the detector 47, and 65 is the waveform of the smoother circuit 50. 
     In FIG. 9, both the desired signal D and the undesired signal U are applied to the receiver 10 through the antenna 20. The intermediate frequency output IF OUT  of the receiver 10 is applied to the envelope detector 44 through the AGC amplifier 43. The time constant of the loop circuits of 43, 44 and 49 is very large so that that loop circuit removes only the variation of the median value level which is usually 0.1-1.0 Hz. The output of the envelope detector 44 is shown in the equation (6) and the waveform 62 in FIG. 10. The output of the detector 44 is applied to the high-pass filter 46 through the amplifier 45. The high-pass filter 46 removes the low frequency component (E 1   2  +E 2   2 ), and then the output of the high-pass filter 46 is the waveform 63 of FIG. 10. The waveform 63 is processed by the envelope detector 47 and the low-pass filter 48, then, the product of the desired signal E 1  and the undesired signal E 2  is obtained as shown in the following formula. 
     
         Y&#34;=2E.sub.1 E.sub.2 
    
     The value Y&#34; is an instantaneous value and varies by fading. Therefore, the waveform is processed as shown by the waveform 65 in FIG. 10 by the smoother 50, and the square root of the square average of Y&#34; is indicated by the level meter 51. The indicated value W by the meter 51 is; 
     
         W=2(&lt;E.sub.1.sup.2 &gt;&lt;E.sub.2.sup.2 &gt;).sup.1/2              (26) 
    
     On the other hand, the value X(=&lt;E 1   2  &gt;+&lt;E 2   2  &gt;) is constant because of the presence of the AGC circuit 43. 
     
         X=&lt;E.sub.1.sup.2 &gt;+&lt;E.sub.2.sup.2 &gt;=C (constant)           (27) 
    
     The following equation is derived from equations (26) and (27). ##EQU12## Therefore, the D/U ratio (=Γ) is measured by properly plotting Γ on the meter 51. 
     FIG. 11 shows the experimental curve of the apparatus of FIG. 9. In FIG. 11, the horizontal axis shows the set value D/U ratio, and the vertical axis shows the measured value W, and the solid line shows the theoretical curve, and the dots show the measured value. The measured value coincides well with the theoretical solid curve. 
     It should be appreciated of course that some modifications of the present invention are possible to those skilled in the art within the spirit of the present invention. For instance, an envelope detection may be not only a square detection, but also a 4&#39;th order detection, or any order detection. 
     It may be advantageous to provide a pilot signal outside of the signal band for the easy detection of the beat between the desired signal and the undesired signal. That pilot signal is preferably modulated in angular modulation. Alternatively, the carrier frequency of the first transmitter may be a little different from that of the second transmitter so that the beat frequency between the two carriers is easily detected. 
     As described above in detail, the present invention detects a cochannel interference very easily, and can be used for instance in a mobile communication, or cordless telephone. 
     From the foregoing, it will now be apparent that a new and improved cochannel interference detection system has been found. It should be understood of course that the embodiments disclosed are merely illustrative and are not intended to limit the scope of the invention. Reference should be made to the appended claims, therefore, rather than the specification as indicating the scope of the invention.