[已解决]请问如何取得sympy solve解出的方程后的返回值？

veilong

1. import math
   
2. from sympy import *
   
3. def solve_c(θ,α):
   
4.     """
   
5.     求解曲线系数c
   
6.     θ为过渡角度,用弧度
   
7.     α为斜线段角度
   
8.     """
   
9.     c1,c2,c3,c4,c5,c6,c7=symbols('c1:8')
   
10.     eq1 = c1*(1/4 - 1/(2 * math.pi))-c4
    
11.     eq2 = c1/(4*math.radians(1/4*θ))-c3
    
12.     eq3 = c1*math.pi-16*c2*math.radians(1/4*θ)**2
    
13.     eq4 = c2*math.radians(1/2*θ)**2+c3*math.radians(1/2*θ)+c4+c6/(2 * math.pi)-c7
    
14.     eq5 = 2*c2*math.radians(1/2*θ)+c3-c5/(4*math.radians(1/4*θ))
    
15.     eq6 = 16*c2*math.radians(1/4*θ)**2-c6*math.pi
    
16.     eq7 = 5.298-c5/4-c7
    
17.     sol = solve((eq1,eq2,eq3,eq4,eq5,eq6,eq7),(c1,c2,c3,c4,c5,c6,c7))
    
18.     return sol
    
19. c = solve_c(19,85)
    
20. #c1=c[c1]

复制代码

                               
登录/注册后可看大图

运行后得{c1: 2.06084003807692, c2: 58.8751760609687, c3: 6.21460191754672, c4: 0.187217130537599, c5: 8.53515996192307, c6: 2.06084003807692, c7: 3.16421000951923}
加上c1=c[c1]这句则报错：name 'c1' is not defined
尝试将symbols中的c1换为x,y,z中之一则正确取得值。
请高手不吝解答，谢谢！

最佳答案

月排行榜 / 总排行榜

XiaoPaiShen

2019-10-16 07:43:51

1. import math
   
2. from sympy import *
   
3. def solve_c(θ,α):
   
4.     """
   
5.     求解曲线系数c
   
6.     θ为过渡角度,用弧度
   
7.     α为斜线段角度
   
8.     """
   
9.     c1,c2,c3,c4,c5,c6,c7=symbols('c1:8')
   
10.     eq1 = c1*(1/4 - 1/(2 * math.pi))-c4
    
11.     eq2 = c1/(4*math.radians(1/4*θ))-c3
    
12.     eq3 = c1*math.pi-16*c2*math.radians(1/4*θ)**2
    
13.     eq4 = c2*math.radians(1/2*θ)**2+c3*math.radians(1/2*θ)+c4+c6/(2 * math.pi)-c7
    
14.     eq5 = 2*c2*math.radians(1/2*θ)+c3-c5/(4*math.radians(1/4*θ))
    
15.     eq6 = 16*c2*math.radians(1/4*θ)**2-c6*math.pi
    
16.     eq7 = 5.298-c5/4-c7
    
17.     sol = solve((eq1,eq2,eq3,eq4,eq5,eq6,eq7),(c1,c2,c3,c4,c5,c6,c7))
    
18.     return sol
    
19. 
    
20. c = solve_c(19,85)
    
21. 
    
22. c1,c2,c3,c4,c5,c6,c7=symbols('c1:8')
    
23. print(c[c1])

复制代码

跳转到最佳答案楼层

XiaoPaiShen

1. import math
   
2. from sympy import *
   
3. def solve_c(θ,α):
   
4.     """
   
5.     求解曲线系数c
   
6.     θ为过渡角度,用弧度
   
7.     α为斜线段角度
   
8.     """
   
9.     c1,c2,c3,c4,c5,c6,c7=symbols('c1:8')
   
10.     eq1 = c1*(1/4 - 1/(2 * math.pi))-c4
    
11.     eq2 = c1/(4*math.radians(1/4*θ))-c3
    
12.     eq3 = c1*math.pi-16*c2*math.radians(1/4*θ)**2
    
13.     eq4 = c2*math.radians(1/2*θ)**2+c3*math.radians(1/2*θ)+c4+c6/(2 * math.pi)-c7
    
14.     eq5 = 2*c2*math.radians(1/2*θ)+c3-c5/(4*math.radians(1/4*θ))
    
15.     eq6 = 16*c2*math.radians(1/4*θ)**2-c6*math.pi
    
16.     eq7 = 5.298-c5/4-c7
    
17.     sol = solve((eq1,eq2,eq3,eq4,eq5,eq6,eq7),(c1,c2,c3,c4,c5,c6,c7))
    
18.     return sol
    
19. 
    
20. c = solve_c(19,85)
    
21. 
    
22. c1,c2,c3,c4,c5,c6,c7=symbols('c1:8')
    
23. print(c[c1])

复制代码

veilong

非常感谢！但还有点不明白，为什么符号变量如果为x，y，z可直接取值，而非x，y，z如c1等需再次symbols('c1:8')？

XiaoPaiShen

c1 是方法solve_c中的局部变量，出了该方法就无法再访问了，在方法外重新运行symbols方法获得它的值就可以访问结果了

veilong

XiaoPaiShen 发表于 2019-10-16 22:45
c1 是方法solve_c中的局部变量，出了该方法就无法再访问了，在方法外重新运行symbols方法获得它的值就可以 ...

谢谢！