为什么我这代码算积分要好几秒才出结果？

cyrus白

1. import random
   
2. import math
   
3. 
   
4. class Function:
   
5.     global dx
   
6.     dx = 0.000001
   
7.     def __init__(self, cor_law="pow(x,2)"):
   
8.         self.cor_law = cor_law
   
9. 
   
10.     def getValue(self,x):
    
11.         try:
    
12.             y = eval(self.cor_law)
    
13.             return y
    
14.         except:
    
15.             return "函数在该点处没意义"
    
16. 
    
17.     def derivate(self,x0):
    
18.         try:
    
19.             y0 = self.getValue(x0)
    
20.             x1 = x0 + dx
    
21.             y1 = self.getValue(x1)
    
22.             return  (y1 - y0)/(x1 - x0)
    
23.         except:
    
24.             return "函数不可导"
    
25. 
    
26.     def integral(self,up,down):
    
27.         try:
    
28.             s = 0
    
29.             if up > down:
    
30.                 while down < up:
    
31.                     yi = self.getValue(down)
    
32.                     s += yi*dx
    
33.                     down += dx
    
34. 
    
35.             else:
    
36.                 while up < down:
    
37.                     yi = self.getValue(up)
    
38.                     s -= yi*dx
    
39.                     up += dx
    
40.             return s
    
41.         except:
    
42.             return "函数不可积"
    
43. 
    
44. 
    
45. if __name__ == "__main__":
    
46.     fun = Function(cor_law='math.sin(x)')
    
47.     print(fun.getValue(3.1415926))
    
48.     print(fun.derivate(0))
    
49.     print(fun.integral(3.14159, 0))

复制代码

qiuyouzhi

效率不行呗
Python本来弱势就是效率

cyrus白

qiuyouzhi 发表于 2020-4-7 20:03
效率不行呗
Python本来弱势就是效率

  刚开始还以为电脑陷入了死循环呢。

BngThea

用numpy不香吗

cyrus白

BngThea 发表于 2020-4-7 20:26
用numpy不香吗

还在摸索中 很多都不懂

BngThea

cyrus白 发表于 2020-4-7 20:28
还在摸索中 很多都不懂

numpy和scipy底层都是c，速度很快
搞科学计算必备的

NZND

为毛不编译一下呢？

NZND

Pyinstaller他不香吗？