编程思路求解答

lzyyy

有大神可以教教我这么用python求解这个方程吗？A是已知量，怎么求出lamda的解，看了很多教程还是不太会.谢谢大神。

东水兽

MATLAB程序：
A = 1:10; % 填写A的值
for i = 1:length(A)
    syms lamda % 注意这里的lamda不是全局变量，是在循环里面的
    a = A(i);
    lamda = solve(1+0.6/lamda+(1+1.2/lamda)^a==0,lamda);
    answer{i} = double(lamda); % 这里应该要使用cell，因为有些A的值的情况下有多个lamda的解
end

python程序：
from sympy import *
import numpy as np
A = np.linspace(1,10,10) # 输入A的值
for x in A:
    lamda = symbols('lamda')
    eq = 1 + 0.6 / lamda - (1 + 1.2 / lamda) ** x
    c = solve(eq, lamda)
    print(c)
两段程序的思路是一致的

大理段戬

二分法就可以吧

lzyyy

因为A是可变的，而且比较多，所以我想用一种简单点的方法，二分法我怕跑不出来程序。

MyOpenGFM

网上教程若不行，估计一对一教学会有所帮助

lzyyy

MyOpenGFM 发表于 2021-3-4 09:50
网上教程若不行，估计一对一教学会有所帮助

您可以提供一些方向吗

lwx

个人感觉可以先解方程得出lamuda的A的关系（方程可解，用求导的方式还是比较容易的），然后编程求解即可

谷梁陌尘

import math

def func(cur_root):
    func = 1+0.6/cur_root -(1+1.2/cur_root)**2
    return func

def binary(convergence, left, right):
    print('current acceptable error: ' + str(convergence) + '\n')
    error = convergence + 1  # 循环开始条件
    cur_root = left
    count = 1
    while error > convergence:
        if abs(func(left)) < convergence:
            print('root = ' + str(left))
        elif abs(func(right)) < convergence:
            print('root = ' + str(left))
        else:
            print(str(count) + ' root = ' +str(cur_root))
            middle = (left + right) / 2
            if (func(left) * func(middle)) < 0:
                right = middle
            else:
                left = middle
            cur_root = left
        error = abs(func(cur_root))
        count += 1
convergence = float(input("your acceptable error:"))
left = float(input('left: '))
right = float(input('right: '))##输入左右区间
binary(convergence, left, right)

怎么说呢。我用python比较少，自己把A取为2，用二分法跑了一下，还是能正常运行出结果的，你可以把A加个循环什么的，或者自己手动换吧。

lzyyy

谷梁陌尘 发表于 2021-3-4 11:29
import math

def func(cur_root):

不好意思，评语太长后面没写完