[已解决]新人求助嵌套问题

mrdeer

问题的公式是
z =2000000*x+50000000+1000000/0.8
z= 2000000*y+50000000+1000000/0.7

z的值在什么范围时2000000/(2000000+50000000/x+1000000/0.8x) < 2000000/(2000000+50000000/y+1000000/0.7y)


我现在想用while函数嵌套z=z+10000之类的，但失败了...

我目前的垃圾想法是这样的...求大佬解答

from sympy import symbols, Eq, solve


while equation_3 < equation_4:
    z = 1000000
    z = str(z)
   
    x, y = symbols("x y")
    equation_1 = Eq(2000000*x+50000000+1000000/0.8, z)
    equation_2 = Eq(2000000*y+50000000+1000000/0.7, z)
    solution = solve((equation_1, equation_2), (x, y))

    equation_3 = Eq(2000000/(2000000+50000000/x+1000000/0.8*x))
    equation_4 = Eq(2000000/(2000000+50000000/y+1000000/0.7*y))
   
z = z + 100000

return z

最佳答案

月排行榜 / 总排行榜

ouyunfu

2023-2-23 09:21:09

1. from sympy import symbols, Eq, solve
   
2. 
   
3. z = 1000000
   
4. x, y = symbols("x y")
   
5. equation_1 = Eq(2000000*x+50000000+1000000/0.8, z)
   
6. equation_2 = Eq(2000000*y+50000000+1000000/0.7, z)
   
7. 
   
8. while True:
   
9.     # 解决方程组
   
10.     solution = solve((equation_1, equation_2), (x, y))
    
11.     x_val = solution[x]
    
12.     y_val = solution[y]
    
13.    
    
14.     # 计算方程式的值
    
15.     equation_3_value = 2000000/(2000000+50000000/x_val+1000000/0.8*x_val)
    
16.     equation_4_value = 2000000/(2000000+50000000/y_val+1000000/0.7*y_val)
    
17.    
    
18.     if equation_3_value.evalf() < equation_4_value.evalf():
    
19.         z += 100000
    
20.         equation_1 = Eq(2000000*x+50000000+1000000/0.8, z)
    
21.         equation_2 = Eq(2000000*y+50000000+1000000/0.7, z)
    
22.     else:
    
23.         break
    
24. 
    
25. print(z)

复制代码

跳转到最佳答案楼层

wp231957

xy.都未知，这东西能有答案吗

mrdeer

z的值是自己输入的，这个循环的目的就是测试出z的值在什么范围时2000000/(2000000+50000000/x+1000000/0.8x) < 2000000/(2000000+50000000/y+1000000/0.7y)。我理解就是先把z值输入x和y的公式，算出x和y后，看看符不符合2000000/(2000000+50000000/x+1000000/0.8x) < 2000000/(2000000+50000000/y+1000000/0.7y)。符合就对了。

mrdeer

wp231957 发表于 2023-2-23 08:34
xy.都未知，这东西能有答案吗

z的值是自己输入的，这个循环的目的就是测试出z的值在什么范围时2000000/(2000000+50000000/x+1000000/0.8x) < 2000000/(2000000+50000000/y+1000000/0.7y)。我理解就是先把z值输入x和y的公式，算出x和y后，看看符不符合2000000/(2000000+50000000/x+1000000/0.8x) < 2000000/(2000000+50000000/y+1000000/0.7y)。符合就对了。

mrdeer

wp231957 发表于 2023-2-23 08:34
xy.都未知，这东西能有答案吗

求大佬指点！！！

ouyunfu

1. from sympy import symbols, Eq, solve
   
2. 
   
3. z = 1000000
   
4. x, y = symbols("x y")
   
5. equation_1 = Eq(2000000*x+50000000+1000000/0.8, z)
   
6. equation_2 = Eq(2000000*y+50000000+1000000/0.7, z)
   
7. 
   
8. while True:
   
9.     # 解决方程组
   
10.     solution = solve((equation_1, equation_2), (x, y))
    
11.     x_val = solution[x]
    
12.     y_val = solution[y]
    
13.    
    
14.     # 计算方程式的值
    
15.     equation_3_value = 2000000/(2000000+50000000/x_val+1000000/0.8*x_val)
    
16.     equation_4_value = 2000000/(2000000+50000000/y_val+1000000/0.7*y_val)
    
17.    
    
18.     if equation_3_value.evalf() < equation_4_value.evalf():
    
19.         z += 100000
    
20.         equation_1 = Eq(2000000*x+50000000+1000000/0.8, z)
    
21.         equation_2 = Eq(2000000*y+50000000+1000000/0.7, z)
    
22.     else:
    
23.         break
    
24. 
    
25. print(z)

复制代码

wp231957

mrdeer 发表于 2023-2-23 08:40
求大佬指点！！！

这东西可以先数学化处理（手工计算一下）
我没太仔细算，不知道是否中间有错
大致最后会简化成以下样子

mrdeer

ouyunfu 发表于 2023-2-23 09:21

这就是大佬么！！！！太厉害了！！！

mrdeer

wp231957 发表于 2023-2-23 09:21
这东西可以先数学化处理（手工计算一下）
我没太仔细算，不知道是否中间有错
大致最后会简化成以下样 ...

谢谢大佬指点！！！