如何处理不确定的功能限制?

计算科学 Python 数值分析
2021-12-10 05:58:44

接近时,如何确保我的下面的函数条件良好我得到的问题是,对于大函数返回一个不确定的形式 接近单调地从 0 增加到 1 ss00s

f(s)=xs2+1xs21xs2(xs(1xs21))2+(xs2+1xs21)2

xs=es

import numpy as np

def func(s):
    xs = np.exp(-1*s)
    num = xs**2+np.sqrt(1-xs**2)-1
    den = np.sqrt(xs**2*(xs*(np.sqrt(1-xs**2)-1))**2+(xs**2+np.sqrt(1-xs**2)-1)**2) 
    return num/den

s = np.linspace(0,20, 1000)
plot(s, func(s))

在此处输入图像描述

1个回答

该函数非常平滑,如这使您可以执行以下策略:s

  • 对于,进行通常的评估,即返回ss0f(s)

  • 对于,定义函数(即围绕 )进行泰勒展开并返回其中包含展开的前几个泰勒项。在评估这些泰勒项时,您需要计算第个导数 ,这将涉及您必须应用 l'Hopital 规则的因子s>s0g(h)=f(1/h)h=0s=g¯(h)|h=1/sg¯g(h)kg(k)(0)00