一元二次拟合

Halcom

1. Rows:=Rows*1.0/1000.0
   
2. Columns:=Columns*1.0/1000.0
   
3. 
   
4. x1:=0.00000000000
   
5. x2:=0.00000000000
   
6. x3:=0.00000000000
   
7. x4:=0.00000000000
   
8. y1:=0.00000000000
   
9. x1y1:=0
   
10. x2y1:=0
    
11. for i := 0 to |Rows|-1 by 1
    
12.     x1:=x1+Rows[i]
    
13.     x2:=x2+Rows[i]*Rows[i]
    
14.     x3:=x3+Rows[i]*Rows[i]*Rows[i]
    
15.     x4:=x4+Rows[i]*Rows[i]*Rows[i]*Rows[i]
    
16.     y1:=y1+Columns[i]
    
17.     x1y1:=x1y1+Rows[i]*Columns[i]
    
18.     x2y1:=x2y1+Rows[i]*Rows[i]*Columns[i]
    
19. endfor
    
20. create_matrix(3, 3, 0, MatrixID)
    
21. set_value_matrix(MatrixID,0,0,|Rows|)
    
22. set_value_matrix(MatrixID,0,1,x1)
    
23. set_value_matrix(MatrixID,0,2,x2)
    
24. set_value_matrix(MatrixID,1,0,x1)
    
25. set_value_matrix(MatrixID,1,1,x2)
    
26. set_value_matrix(MatrixID,1,2,x3)
    
27. set_value_matrix(MatrixID,2,0,x2)
    
28. set_value_matrix(MatrixID,2,1,x3)
    
29. set_value_matrix(MatrixID,2,2,x4)
    
30. 
    
31. create_matrix(3,1,0, MatrixID1)
    
32. set_value_matrix(MatrixID1,0,0,y1)
    
33. set_value_matrix(MatrixID1,1,0,x1y1)
    
34. set_value_matrix(MatrixID1,2,0,x2y1)
    
35. 
    
36. invert_matrix(MatrixID,'general', 0, MatrixInvID)
    
37. mult_matrix(MatrixInvID,MatrixID1,'AB', MatrixMultID)
    
38. 
    
39. get_value_matrix(MatrixMultID, 0, 0, a0)
    
40. get_value_matrix(MatrixMultID, 1, 0, a1)
    
41. get_value_matrix(MatrixMultID, 2, 0, a2)
    
42. return ()

复制代码

15s1f2ds12fd1f

666666666666