锄禾日当午wasd 2022-04-18 16:24
浏览 57
已结题

关于scipy解非线性方程的问题

问题遇到的现象和发生背景

img

我用python中的scipy模块时, 想要解这个方程,其中Vl 是未知量,其他均已知,但是解不出来,提示下图,请问有人可以帮忙解答一下原因吗?

img

(代码中的方程是经过我简化的,x是要求的未知量)

问题相关代码
from numpy import sin, cos, tan, array
from scipy.optimize import root, fsolve

def f(i, params):
    x = i[0]
    a_1, V_s1, a_2, V_s2 = params
    return array(tan((a_1 / x) * (x ** 2 - V_s1 ** 2) ** 0.5) - a_2 * ((V_s2 ** 2 - x ** 2) / (x ** 2 - V_s1 ** 2)) ** 0.5)

E_1 = 1.2503 * 10 ** 10
v_1 = 0.09
E_2 = 5.3105 * 10 ** 10
v_2 = 0.12

u_1 = ((4 * v_1 - 1) / (10 * v_1 - 2) * E_1)

V_s1 = 1150
D = 8
u_2 = ((4 * v_2 - 1) / (10 * v_2 - 2) * E_2)
V_s2 = 1720
w = 10
a_1 = D * w / V_s1
a_2 = (u_2 * V_s1) / (u_1 * V_s2)

# 传参
params = [a_1, V_s1, a_2, V_s2]

init_guess = array([0])
result_1 = fsolve(f, init_guess, args = params)
result_2 = root(f, init_guess, args = params)

  • 写回答

0条回答 默认 最新

    报告相同问题?

    问题事件

    • 系统已结题 4月26日
    • 修改了问题 4月19日
    • 创建了问题 4月18日

    悬赏问题

    • ¥100 如何寻找到黑客帮助,愿意付丰厚的酬劳
    • ¥15 java代码写在记事本上后在cmd上运行时无报错但又没生成文件
    • ¥15 关于#python#的问题:在跑ldsc数据整理的时候一直抱这种错误,要么--out识别不了参数,要么--merge-alleles识别不了参数(操作系统-linux)
    • ¥15 PPOCRLabel
    • ¥15 混合键合键合机对准标识
    • ¥100 现在不懂的是如何将当前的相机中的照片,作为纹理贴图,映射到扫描出的模型上
    • ¥15 魔霸ROG7 pro,win11.息屏后会显示黑屏,如图,如何解决?(关键词-重新启动)
    • ¥15 有没有人知道这是哪里出了问题啊?要怎么改呀?
    • ¥200 C++表格文件处理-悬赏
    • ¥15 Windows Server2016本地登录失败