BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    The present invention relates a method for compensating timing errors of real-time clocks. 
         [0003]    2. Description of the Related Art 
         [0004]    As shown in  FIG. 1 , a minimum unit of real-rime clocks is “second”, which principle is to divide a 32768 Hz square wave, which is an output frequency of a crystal oscillator, into 1 Hz square waves. The 1 Hz square waves generate the unit of time “second” and calendar logic circuits. Accordingly, calendar logics of minute, hour, day, month, and week are generated. The output standard frequency of the crystal oscillator is 32768 Hz; however, errors of the output frequency of the crystal oscillator always exist due to the changes of temperature and manufacturing errors of the quartz crystals. 
         [0005]    Generally, the output frequency of the crystal oscillator is compared with the relative errors of the standard frequency 32768 Hz for deriving an error measuring standard, as the expression shown, the unit is PPM (Parts Per Million). 
         [0000]    
       
         
           
             ERR 
             = 
             
               
                 
                   f 
                   - 
                   32768 
                 
                 32768 
               
               * 
               
                 10 
                 6 
               
                
               PPM 
             
           
         
       
     
         [0006]    In present technology, the errors between the output frequency of the crystal oscillator and the standard frequency 32768 Hz are usually compensated by means of making compensations on the relative errors at regular times. According to the regular times, for example, find out the clock numbers of compensation that need to be compensated in accordance with the relative errors (ERR) and change the value of the frequency dividing counter (CNT) in every 10 seconds. As shown in  FIG. 2 , while the compensation is executed at every 10 seconds, 10 seconds, which means the total numbers of the frequency dividing counter and the clock numbers of reference frequency to be 32768*10=327680. Therefore, the 10-seconds frequency dividing counter CNT adding 1 is equal to a compensation of +3 PPM to the reference frequency. The frequency dividing counter CNT which fails to add one time is equal to a compensation of −3 PPM to the reference frequency. In general, the sampling frequency of ERR is 0.1 Hz, the ERR effective value of compensatory data is 3 PPM, the maximum errors after the compensation is 1.5 PPM. 
         [0007]    The low of the sampling frequency of ERR derived by the compensatory method at regular times is, the high the effective value of the compensatory data is. However, the low sampling frequency causes aliasing errors of frequency spectrum. As a result, the sampling frequency of ERR attained by the method aforementioned cannot be too high and low. A contradictory relationship exists between the sampling frequency apparatus of ERR and the last compensatory accuracy. 
       SUMMARY OF THE INVENTION 
       [0008]    The purpose of the present invention is to provide a method for compensating timing errors of real-time clocks with characteristics of an adjustable sampling frequency of relative error ERR and a higher compensatory accuracy. 
         [0009]    To achieve the purpose abovementioned, the solution of the present invention comprises a calculating step and a compensating step: 
         [0010]    wherein the calculating steps comprising: 
         [0011]    step 1, subject a cyclic subtraction register M2 to an assignment operation, a cyclic subtraction register M1=the cyclic subtraction register 
         [0000]    
       
         
           
             
               
                 M 
                  
                 
                     
                 
                  
                 2 
               
               = 
               
                 
                   
                     10 
                     6 
                   
                   * 
                   S 
                    
                   
                       
                   
                    
                   1 
                   * 
                   S 
                    
                   
                       
                   
                    
                   2 
                 
                 
                   S 
                    
                   
                       
                   
                    
                   3 
                 
               
             
             ; 
           
         
       
     
         [0000]    a periodic number of accumulative errors register M3 is assigned to be 0, wherein the S1 is used to adjust an ERR effective value, the S2 is used to adjust an operating frequency of compensatory circuits, and the S3 is used to adjust a calculation of times, then execute step 2; 
         [0012]    step 2, make a compensatory flag register COM assigned to be 0 when a rising edge of a 
         [0000]    
       
         
           
             32768 
             
               S 
                
               
                   
               
                
               2 
             
           
         
       
     
         [0000]    Hz clock is arrived, the M2 is executed by a subtraction, M2=M1−|ERR*S1|, M1=M2, then execute step 3; 
         [0013]    step 3, judge the M2, if M2&lt;0, an assignment to M2 is executed, 
         [0000]    
       
         
           
             
               
                 M 
                  
                 
                     
                 
                  
                 2 
               
               = 
               
                 
                   M 
                    
                   
                       
                   
                    
                   1 
                 
                 + 
                 
                   
                     
                       10 
                       6 
                     
                     * 
                     S 
                      
                     
                         
                     
                      
                     1 
                     * 
                     S 
                      
                     
                         
                     
                      
                     2 
                   
                   
                     S 
                      
                     
                         
                     
                      
                     3 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    M1=M2, and the compensatory flag register COM is assigned to be 1, then execute step 4; otherwise, execute the step 2; and 
         [0014]    step 4, judge ERR*S1, if ERR*S1&gt;0, execute 
         [0000]    
       
         
           
             
               
                 M 
                  
                 
                     
                 
                  
                 3 
               
               = 
               
                 
                   M 
                    
                   
                       
                   
                    
                   3 
                 
                 + 
                 
                   
                     S 
                      
                     
                         
                     
                      
                     2 
                   
                   
                     S 
                      
                     
                         
                     
                      
                     3 
                   
                 
               
             
             ; 
           
         
       
     
         [0000]    otherwise, execute 
         [0000]    
       
         
           
             
               M 
                
               
                   
               
                
               3 
             
             = 
             
               
                 M 
                  
                 
                     
                 
                  
                 3 
               
               - 
               
                 
                   S 
                    
                   
                       
                   
                    
                   2 
                 
                 
                   S 
                    
                   
                       
                   
                    
                   3 
                 
               
             
           
         
       
     
         [0000]    and execute said step 2; and 
         [0015]    wherein the compensating step comprising: 
         [0016]    step 1, make a frequency dividing counter CNT assigned to be 0, then execute step 2; 
         [0017]    step 2, make a compensatory flag register FLAG assigned to be 1 when a rising edge of a 1 Hz clock is arrived, then execute step 3; 
         [0018]    step 3, judge the FLAG and the M3, if FLAG=1 and M3&lt;0, execute step 4 while waiting until CNT=S4; if FLAG=1, CNT=0, and M3&gt;0, execute step 5; otherwise, execute said step 2; 
         [0019]    step 4, execute an assignment operation, CNT=0, M3=M3+S4, FLAG=0, and restart the step 2, wherein the S4 is a maximum compensatory periodic number in 1 second; and 
         [0020]    step 5, execute an assignment operation, CNT=S4, M3=M3-S4, FLAG=0, and restart the step two. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0021]      FIG. 1  is a schematic view showing a principle of real-time clocks; 
           [0022]      FIG. 2  is a schematic view showing a compensatory method being executed at regular times in present technology; 
           [0023]      FIG. 3  is a schematic view showing a basic relationship between a standard frequency and an actual frequency; 
           [0024]      FIG. 4  is a schematic view showing an order of calculating accumulative errors; 
           [0025]      FIG. 5  is a schematic view showing the order of calculating accumulative errors in execution; 
           [0026]      FIG. 6  is a schematic view showing a method for compensating errors; 
           [0027]      FIG. 7  is a schematic view showing an order of the method for compensating errors in execution. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
       [0028]    The present invention is illustrated in detail with embodiments. 
         [0029]    A basic relationship of the standard frequency (32768 Hz) and the actual frequency (32768*(1−ERRppm) can be obtained according to the basic definition of the relative error ERR, 
         [0000]    
       
         
           
             ERR 
             = 
             
               
                 
                   f 
                   - 
                   32768 
                 
                 32768 
               
               * 
               
                 10 
                 6 
               
                
               
                 PPM 
                 . 
               
             
           
         
       
     
         [0030]    The relationship between the standard period and the actual period can be calculated by the basic relationship of the standard frequency (32768 Hz) and the actual frequency (32768*(1−ERRppm). As shown in  FIG. 3 , the difference between the abovementioned frequencies is 1*ERR/10 6  periods after passing 1 period, is 2*ERR/10 6  periods after passing 2 periods, and is (10 6 /ERR)*ERR/10 6 =1 periods after passing 10 6 /ERR periods. Accordingly, after passing N periods, the relationship of the standard frequency and the actual frequency is different from each other by W periodic numbers, as shown in Table 1: 
         [0000]    
       
         
               
               
               
             
           
               
                   
                 TABLE 1 
               
               
                   
                   
               
               
                   
                 N 
                 W 
               
               
                   
                   
               
             
             
               
                   
                 
                   
                     
                       
                         
                           10 
                           6 
                         
                         ERR 
                       
                     
                   
                 
                 1 
               
               
                   
                   
               
             
          
         
       
     
         [0031]    That is, the N: W=10 6 /ERR: 1, which needs to be compensated 1 period after passing 10 6 /ERR periods. In the actual compensatory operation, 10 6 /ERR periods are usually converted into a cyclic subtraction, and 1 period is compensated when the value of 10 6 /ERR periods is subtracted to 0 or a negative number. Because the calculating times of the 10 6  is much bigger and the relative error ERR is smaller, accordingly, the value of the cyclic subtraction becomes big, which causes big power dissipation of the circuits. The value of the relative error ERR is smaller, for instance, 1.001, which has a decimal point cannot be calculated by the digital circuits. Only ERR=1 can be subjected the calculation, and the value 0.001 cannot be calculated, accordingly, the compensatory accuracy is affected. Therefore, the relationship of the numbers of the N and the W is converted into a Table 2 as below: 
         [0000]    
       
         
               
               
               
             
           
               
                   
                 TABLE 2 
               
               
                   
                   
               
               
                   
                 N 
                 W 
               
               
                   
                   
               
             
             
               
                   
                 
                   
                     
                       
                         
                           
                             
                               10 
                               6 
                             
                             × 
                             S 
                              
                             
                                 
                             
                              
                             1 
                           
                           
                             ERR 
                             × 
                             S 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                         × 
                         S 
                          
                         
                             
                         
                          
                         
                           2 
                           ÷ 
                           S 
                         
                          
                         
                             
                         
                          
                         3 
                       
                     
                   
                 
                 1 × S2 ÷ S3 
               
               
                   
                   
               
             
          
         
       
     
         [0032]    wherein the meanings of the S1, S2, and s3 shown as Table 3 represents: 
         [0000]    
       
         
               
               
             
           
               
                 TABLE 3 
               
               
                   
               
               
                 parameter 
                 controlling value 
               
               
                   
               
             
             
               
                 S1 
                 ERR effective 
               
               
                   
                 value 
               
               
                 S2 
                 compensatory 
               
               
                   
                 circuit operating 
               
               
                   
                 frequency 
               
               
                 S3 
                 Calculation of 
               
               
                   
                 times 
               
               
                   
               
             
          
         
       
     
         [0033]    The function of the S1 is to raise the compensatory accuracy. The compensatory accuracy is raised as a multiple of the S1, the S1 is a positive integer, for instance, 1, 2, 3 and N. However, if ERR=1.001, of which value 0.001 cannot be calculated by the digital circuits, consequently, only ERR=1 can be calculated, and the value of 0.001 cannot be calculated. In 
         [0000]    
       
         
           
             
               
                 
                   10 
                   6 
                 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
               
                 ERR 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
             
             , 
           
         
       
     
         [0000]    while S1 is 1000, ERR×S1=1001. As a result, the value of 1001 can be calculated by the digital circuits, and the compensatory accuracy is raised 1000 times. 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                   
               
               
                 Table 4, interpretations on the examples of the parameter S1: 
               
             
          
           
               
                   
                   
                 compensating 
               
               
                 N 
                 W 
                 accuracy 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             
                               10 
                               6 
                             
                             × 
                             1000 
                           
                           
                             ERR 
                             × 
                             1000 
                           
                         
                         = 
                         
                           
                             10 
                             9 
                           
                           
                             ERR 
                             ′ 
                           
                         
                       
                     
                   
                 
                 1 
                 
                   
                     
                       
                         
                           
                             1 
                              
                             
                                 
                             
                              
                             ppm 
                           
                           1000 
                         
                         = 
                         
                           0.001 
                            
                           
                               
                           
                            
                           ppm 
                         
                       
                     
                   
                 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             
                               10 
                               6 
                             
                             × 
                             
                               S 
                                
                               1 
                             
                           
                           
                             ERR 
                             × 
                             S1 
                           
                         
                         = 
                         
                           
                             
                               10 
                               6 
                             
                             × 
                             S1 
                           
                           
                             ERR 
                             ′ 
                           
                         
                       
                     
                   
                 
                 1 
                 
                   
                     
                       
                         
                           1 
                            
                           
                               
                           
                            
                           ppm 
                         
                         S1 
                       
                     
                   
                 
               
               
                   
               
             
          
         
       
     
         [0034]    In the expression 
         [0000]    
       
         
           
             
               
                 
                   10 
                   6 
                 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
               
                 ERR 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
             
             , 
           
         
       
     
         [0000]    the S1 can be used to adjust the compensatory accuracy. As shown in Table 4, the parameters are set to be S1=1000, S2=1 and S3=1, and the compensatory accuracy is accordingly raised 1000 times. 
         [0035]    In accordance with Table 4, 1 period will be compensated after passing 
         [0000]    
       
         
           
             
               10 
               9 
             
             
               ERR 
               ′ 
             
           
         
       
     
         [0000]    reference frequency periods. The unit of ERR is PPM, and the unit of ERR′ is PPM/1000=0.001 PPM. The unit of the ERR′ is 0.001 PPM, so the compensatory data ERR′ is accurate and makes the compensatory accuracy to be 0.001 PPM. The merit of the S1 which is adjustable is: the compensatory accuracy can be raised. The taking range of the value of the S1 to be: the S1 is a positive integer, for example, 1, 2, 3 and N. 
         [0036]    The S2 is used to adjusting the operating frequency of compensatory circuits for lowering the power dissipation of the circuits, Table 5 shows the interpretations on the examples of parameter S2: 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                   
               
               
                 Table 5, interpretations on the parameter S2: 
               
             
          
           
               
                   
                   
                 operating 
               
               
                   
                   
                 frequency of 
               
               
                   
                   
                 compensatory 
               
               
                 N 
                 W 
                 circuits 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             10 
                             6 
                           
                           ERR 
                         
                         × 
                         8 
                       
                     
                   
                 
                 1 × 8 
                 
                   
                     
                       
                         
                           1 
                           8 
                         
                          
                         f 
                       
                     
                   
                 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             10 
                             6 
                           
                           ERR 
                         
                         × 
                         S 
                          
                         
                             
                         
                          
                         2 
                       
                     
                   
                 
                 1 × S2 
                 
                   
                     
                       
                         
                           1 
                           
                             S 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                          
                         f 
                       
                     
                   
                 
               
               
                   
               
             
          
         
       
     
         [0037]    The S2 can adjust the operating frequency of compensatory circuits. The expression 
         [0000]    
       
         
           
             
               10 
               6 
             
             ERR 
           
         
       
     
         [0000]    is conducted by a division. In the actual execution of the circuits, the division is converted into a subtraction to execute the calculation for diminishing the area of the chips and saving the manufacturing cost. The period of the cyclic subtraction is S2×T, the T is the reference frequency period, the minued of the cyclic subtraction is the difference derived from the last subtraction, and the subtrahend is ERR. The initial value of the minued is 10 6 , the time have passed 
         [0000]    
       
         
           
             
               
                 10 
                 6 
               
               ERR 
             
             × 
             S 
              
             
                 
             
              
             2 
             × 
             T 
           
         
       
     
         [0000]    periods when the difference is a negative number, and one cycle of compensatory calculation is finished, S2×T are needed to be compensated. 
         [0038]    For instance, regarding Table 5, S1=1, S2=8, S3=1, the compensatory circuits are executed the cyclic subtraction of which period is 8 T. In other words, the operating frequency of compensatory circuits is 
         [0000]    
       
         
           
             
               
                 1 
                 8 
               
                
               f 
             
             , 
           
         
       
     
         [0000]    the f is the reference frequency. Therefore, the operating frequency of compensatory circuits is lowered, and the power dissipation of the circuits is accordingly decreased. The taking range of the value of the S2 is: the S2 is a positive integer, for example, 1, 2, 3, . . . and N. 
         [0039]    The function of the S3 is to adjust the calculating times of one compensatory cycle for minimize the minued and the calculating numbers calculated by the subtractor. Therefore, the area of the compensatory calculating circuits is decreased, and the manufacturing cost is saved. Table 6 shows the interpretations on the example of the parameter S3: 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                   
               
               
                 Table 6, the interpretations on the parameter S3: 
               
             
          
           
               
                   
                   
                 The initial 
               
               
                   
                   
                 value of the 
               
               
                 N 
                 W 
                 minued 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 10 
                                 6 
                               
                               × 
                               1 
                             
                             
                               ERR 
                               × 
                               1 
                             
                           
                           × 
                           
                             16 
                             ÷ 
                             16 
                           
                         
                         = 
                         
                           
                             62500 
                             ERR 
                           
                           × 
                           16 
                         
                       
                     
                   
                 
                 1 × 16 ÷ 16 = 1 
                 10 6  × 1 ÷ 16 = 62500 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 10 
                                 6 
                               
                               × 
                               S 
                                
                               
                                   
                               
                                
                               1 
                             
                             
                               ERR 
                               × 
                               S 
                                
                               
                                   
                               
                                
                               1 
                             
                           
                           × 
                           S 
                            
                           
                               
                           
                            
                           
                             2 
                             ÷ 
                             S 
                           
                            
                           
                               
                           
                            
                           3 
                         
                         = 
                         
                           
                             
                               
                                 10 
                                 6 
                               
                               × 
                               S 
                                
                               
                                   
                               
                                
                               
                                 1 
                                 ÷ 
                                 S 
                               
                                
                               
                                   
                               
                                
                               3 
                             
                             
                               ERR 
                               ′ 
                             
                           
                           × 
                           S 
                            
                           
                               
                           
                            
                           2 
                         
                       
                     
                   
                 
                 1 × S2 ÷ S3 = 1 
                 10 6  × S1 ÷ S3 
               
               
                   
               
             
          
         
       
     
         [0040]    The S3 is adjusted after taking the value of the S1 and S2, and the calculating times of 1 compensatory cycle is accordingly adjusted. 
         [0041]    For instance, as Table 6 shown, S1=1, S2=16, S3=16. According to the parameter S2, the compensatory circuits are subjected 1 subtraction every 16 T. The function of the S3 is: set the initial value of the minued of the cyclic subtraction to be 10 6 ×S1÷S3. S3=16, the initial value of the minued of the cycle subtraction is set to be 10 6 ×1+16=62500 of which calculating times of one compensatory calculating cycle are 16 times less than the calculating times of the initial value of the minued set to be 10 6 . That is, the minued, the initial value of the minued of the cyclic subtraction is minimized by 16 times, the initial value, the subtrahend ERR of the cyclic subtraction is kept the same, accordingly, the quotient, the calculating times of one compensatory cycle is minimized by 16 times. While one compensatory calculating cycle is closed, 
         [0000]    
       
         
           
             
               62500 
               ERR 
             
             × 
             16 
              
             T 
           
         
       
     
         [0000]    have been passed, and 1 T is compensated. 
         [0042]    The benefit of the S3 which can be adjusted is: to reduce the minued and the calculating numbers of the subtractor for diminishing the area of the circuits of compensatory calculation and saving the manufacturing cost. The taking range of the value of the S3 to be: the factors which is divisible by 10 6 ×S1 and 1×S2. For instance, S1=1, S2=16, the value of the S3 can be 16. Because 10 6 *1/16=15625 and 1*16/16=1 both conform to the divisibility rule, the S3 which is the factor of 10 6 ×1 and 1×16, the value of the S3 can also be 2, 4 and 8. 
         [0043]    To sum up, because 1×S2÷S3 periods are needed to be compensated after passing 
         [0000]    
       
         
           
             
               
                 
                   10 
                   6 
                 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
               
                 ERR 
                 × 
                 S 
                  
                 
                     
                 
                  
                 1 
               
             
             × 
             S 
              
             
                 
             
              
             
               2 
               ÷ 
               S 
             
              
             
                 
             
              
             3 
           
         
       
     
         [0000]    periods, and 1×S2÷S3 periods cannot be completely compensated within 1 second, the 1 second will become too long or short to be accepted. As a result, S4 periodic numbers can be compensated at most in 1 second. Consequently, finishing compensating W (1×S2÷S3) periodic numbers takes W/S4 seconds, as shown in Table 7. 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                   
               
               
                 Table 7, the interpretations on the parameter S4: 
               
             
          
           
               
                   
                   
                 the maximum compensatory 
               
               
                   
                 S4 
                 errors 
               
               
                   
                   
               
               
                   
                 8 
                 
                   
                     
                       
                         
                           8 
                           32768 
                         
                         = 
                         
                           244.140625 
                            
                           
                               
                           
                            
                           ppm 
                         
                       
                     
                   
                 
               
               
                   
                   
               
               
                   
                 taking range of the  value to be:  1, 2, 3, . . . 32768 
                 
                   
                     
                       
                         
                           S 
                            
                           
                               
                           
                            
                           4 
                         
                         32768 
                       
                     
                   
                 
               
               
                   
                   
               
             
          
         
       
     
         [0044]    The S4 is defined as the periodic numbers of reference frequency of the maximum compensation in 1 second, and the maximum compensatory errors is 
         [0000]    
       
         
           
             
               
                 S 
                  
                 
                     
                 
                  
                 4 
               
               32786 
             
             . 
           
         
       
     
         [0045]    For example, as shown in Table 7, S4=8, 8 reference frequency periods can be compensated at most in 1 second, accordingly, the maximum compensatory errors are 8 reference frequency periods, which are compensated in every second. The relative compensatory value is 
         [0000]    
       
         
           
             
               
                 8 
                 32786 
               
               = 
               
                 244.140625 
                  
                 
                     
                 
                  
                 PPM 
               
             
             , 
           
         
       
     
         [0000]    , wherein the taking range of the value of the S4 to be: 1, 2, 3, . . . and 32768. 
       Applied Examples 
       [0046]    Set the parameters to be S1=1000, S2=32768, S3=512, S4=16. 
         [0000]    
       
         
               
             
               
               
             
           
               
                   
               
               
                 Table 8 shows the relationship of N and W: 
               
             
          
           
               
                 N 
                 W 
               
               
                   
               
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 10 
                                 6 
                               
                               × 
                               1000 
                             
                             
                               ERR 
                               × 
                               1000 
                             
                           
                           × 
                           
                             32768 
                             ÷ 
                             512 
                           
                         
                         = 
                         
                           
                             
                               
                                 
                                   10 
                                   9 
                                 
                                 ÷ 
                                 512 
                               
                               
                                 ERR 
                                 ′ 
                               
                             
                             × 
                             32768 
                           
                           = 
                           
                             
                               1953125 
                               
                                 ERR 
                                 ′ 
                               
                             
                             × 
                             32768 
                           
                         
                       
                     
                   
                 
                 1 × 32768 ÷ 512 = 64 
               
               
                   
               
             
          
         
       
     
         [0047]    In accordance with Table 8, the frequency dividing counter CNT is compensated 64 periodic numbers after passing 
         [0000]    
       
         
           
             
               1953125 
               
                 ERR 
                 ′ 
               
             
             × 
             32768 
           
         
       
     
         [0000]    periods. The operating frequency of compensatory circuits is 32768 times of the reference frequency periods (32768 T), which means 1 second. The initial value of the minued of the cyclic subtraction is set to be 10 9 ÷512=1953125. S4=16, which means 16 periodic numbers can be compensated at most in 1 second. 
         [0048]    Consulting Table 8, ERR′=ERR*1000, the unit of ERR is 1 PPM, the unit of ERR′ is 0.001 PPM which means Parts Per Billion, and the compensatory accuracy is 0.001 PPM. 
         [0049]    As shown in  FIG. 4 , the schematic view showing the calculating order of the accumulative errors, after passing 
         [0000]    
       
         
           
             
               1953125 
               
                 ERR 
                 ′ 
               
             
             × 
             32768 
           
         
       
     
         [0000]    reference frequency periods, the accumulative errors are accumulated to 64 reference frequency periodic numbers. 1 second represents 32768 reference frequency periods, so 
         [0000]    
       
         
           
             1953125 
             
               ERR 
               ′ 
             
           
         
       
     
         [0000]    seconds represent the accumulative errors are accumulated to 64 reference frequency periodic numbers, and the errors are recorded. 
         [0050]    Wherein, the M1 and the M2 are two registers with 24 digits. The initial value is assigned as 1953125 at 0 second, and one subtraction operation is executed at 1 second, 2 seconds, and N seconds. The difference of the 1 second subtraction is seen as the minued, and the minued is used to subtract the absolute value of ERR′. Until the second to be 
         [0000]    
       
         
           
             
               
                 
                   1953125 
                   
                      
                     
                       ERR 
                       ′ 
                     
                      
                   
                 
                  
                 integer 
               
               + 
               1 
             
             , 
           
         
       
     
         [0000]    and M2&lt;0, the accumulative errors are accumulated to 64 reference frequency periodic numbers after passing 
         [0000]    
       
         
           
             
               1953125 
               
                 ERR 
                 ′ 
               
             
             × 
             32768 
           
         
       
     
         [0000]    reference frequency periods, and the register M3 start to add or subtract 64 to record the periodic numbers of the accumulative errors. The positive and negative signs of ERR′ represent the positive and negative errors of the crystal oscillating frequency. When the crystal oscillating frequency is smaller, ERR′ is defined as the positive sign, while the crystal oscillating frequency is bigger, ERR′ is defined as the negative sign. 
         [0000]    
       
         
           
             1953125 
             
               ERR 
               ′ 
             
           
         
       
     
         [0000]    is probably indivisible, and M2 is the negative number. The errors generated from 
         [0000]    
       
         
           
             1953125 
             
               ERR 
               ′ 
             
           
         
       
     
         [0000]    and M2 will be eliminated in the following cycles. In the next start of every cycle, the initial value of the M2 is set to be M1+1593125. 
         [0051]    As shown in  FIG. 5 , the schematic view showing the calculating order of the accumulative errors in operation, wherein the COM is the register with 1 digit, which is a compensating flag register. The COM is a high speed TTL (HTTL) and adds 64 to or subtracts 64 from the register M3 when the accumulative errors accumulate to 64 periodic numbers. After 64 is added to or taken from the register M3, the assignment of the register COM is reset to be 0. The frequency of the operated synchronic clocks of the registers M1 and M2 is 1 Hz, which means 32768 reference frequency periods. 
         [0052]      FIG. 5  shows the conversion of how many periodic numbers of reference frequency equal to the reference frequency errors ERR′, whose unit is 0.001 ppm, and the accumulative errors. The main execution of the above conversion is to set the initial value of the M1 and M2 to be 1953126, subjects one subtraction operation to the data stored in the registers M1 and M2 every second, subtracts the absolute value of ERR′ and sets the assignment of COM to be 0. Judges the difference of the last 1 second subtraction every 1 second, if the difference is larger than or equal to 0, the current 1 second continues subtracting the absolute value of ERR′ one time; if the difference are smaller than 0, the accumulative errors accumulates to 64 periods, sets the initial value of M2 to be 02&lt;=1953125+M1, M1&lt;=M2, COM=1, and simultaneously adds 64 to or subtracts 64 from the register M3. If ERR′&gt;0, M3&lt;=M3+64; if ERR′&lt;0, M3&lt;=M3-64. 
         [0053]    As  FIG. 6  shown, the schematic view of the method for compensating errors, the register M3 is a register storing periodic numbers of the accumulative errors. If M3&gt;0, the crystal oscillating frequency becomes small, and the clock slows down. While CTN=0, the compensation is to set the assignment of CNT to be 16 for fasting the clock by 16 reference frequency periods, and M3=M3-16 records the 16 periods which have been compensated. If M3&lt;0, the crystal oscillating frequency becomes big, and the clock goes fast. While CNT=16, the compensation is to reset the assignment of CNT to be 0 for slowing the clock by 16 reference frequency periods, and M3=M3+16 records 16 periods which has been compensated. 
         [0054]      FIG. 7  shows a schematic view revealing the executing order of the method for compensating errors, wherein the M3 is the register to store the periodic numbers of accumulative errors, the reference frequency of the CNT, which is 32.768 KHz, is divided as the lowest five digits of the counter of which the frequency is 1 Hz. The counter has 15 digits, the error compensation is conducted to decode the lowest five digits of the counter. Accordingly, the area needed to be decoded is diminished, and the power dissipation is simultaneously decreased. The FLAG is the flag register for keeping the maximum compensatory times to be 1 in every 1 second. 
         [0055]    The time unit of “second” is compensated in accordance with the reference frequency periodic numbers of accumulative errors M3. The compensation is conducted to interfere the counter CNT. To be more specific, while the rising edge of the 1 Hz clock is arrived, sets the flag register FLAG to be 1. Judges M3 and CNT, if M3&lt;0, CNT=16 and FLAG=1, resets the most highest digit of the counter CNT to be 0, that is to say, counter subtracts 16, and reset the FLAG to be 0; if M3&gt;0, CNT=1, and FLAG=1, the most highest digit of counter CNT is set to be 1, namely, the counter adds 16, and the FLAG is reset to be 0. The flag register FLAG is set to be 1 in the rising edge of the 1 Hz clock, and is reset to be 0 after the compensation for compensating one time at most in 1 second.
The present invention relates to a method for compensating timing errors of real-time clocks, which comprises a compensating step, wherein in step 1, assign CNT to be 0 and execute step two; in step 2, assign FLAG to be 1 when a rising edge of 1 Hz clock is arrived and execute step 3; in step 3, judge FLAG and M3, if FLAG=1 and M3<0, execute step 4 while waiting until CNT=S4; if FLAG=1, CNT=0 and M3>0, execute step 5; otherwise execute step 2; in step 4, execute an assignment operation, CNT=0, M3=M3+S4, FLAG=0 and restart step 2; in step 5, execute an assignment operation, CNT=S4, M3=M3-S4, FLAG=0, and restart step 2. A sampling frequency of relative errors ERR of the present invention is adjustable, and a compensatory accuracy is much higher.