【LinRegArray2】

1. inputs:
   
2.         IndepArray[XSize]( numericarray ),                       
   
3.         DepArray[YSize]( numericarray ),                       
   
4.         Size( numericsimple ),                                                           
   
5.                                                              
   
6.         TgtPos( numericsimple ),                                                      
   
7.         oLRSlope( numericref ),
   
8.         oLRAngle( numericref ),
   
9.         oLRIntercept( numericref ),                                                   
   
10.         oLRValueRaw( numericref ) ;                                             
    
11.                                                                                      
    
12.                                                           
    
13. 
    
14. variables:
    
15.         var0( 0 ),
    
16.         var1( 0 ),
    
17.         var2( 0 ),
    
18.         var3( 0 ),
    
19.         var4( 0 ),
    
20.         var5( 0 ) ;
    
21. 
    
22. if Size > 1 then
    
23.         begin
    
24.         var0 = 0 ;
    
25.         var1 = 0 ;
    
26.         var2 = 0 ;
    
27.         var3 = 0 ;
    
28.         for var4 = 1 to Size
    
29.                 begin
    
30.                 var0 = var0 + IndepArray[var4] ;
    
31.                 var3 = var3 + Square( IndepArray[var4] ) ;
    
32.                 var1 = var1 + DepArray[var4] ;
    
33.                 var2 = var2 + ( IndepArray[var4] * DepArray[var4] ) ;
    
34.                 end ;
    
35.         var5 = Square( var0 ) - Size * var3 ;
    
36.         oLRSlope = ( Size * var2 - var0 * var1) / var5 ;
    
37.         oLRAngle = ArcTangent( oLRSlope ) ;
    
38.         oLRIntercept = ( var1 + oLRSlope * var0 ) / Size ;
    
39.         oLRValueRaw = oLRIntercept + oLRSlope * -TgtPos ;
    
40.         LinRegArray2 = 1 ;
    
41.         end
    
42. else
    
43.         LinRegArray2 = -1 ;

复制代码