微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

Python - 使用 t 作为参数和变量的数值积分

如何解决Python - 使用 t 作为参数和变量的数值积分

我正在尝试实现数值积分。具体来说,

Integral

我使用梯形法则得到如下解析表达式

Analytical

当我尝试用数字实现它时。

gamma=0.1 #parameter

Tmax=100 

N=100000 #number of steps

t=np.linspace(0,Tmax,N)

B=np.random.rand(len(t))

def A(t):
    A=np.zeros(len(t))
    trapezi=0
    for i in range(1,N-1):
        trapezi+=(np.exp(-gamma*(Tmax-t[i]))*B[i])
    h=Tmax/N
    A=h*(trapezi+1/2*(B[-1]+B[0]*np.exp(-gamma*Tmax)))
    return A

在这里我创建了一个随机 B 矩阵。

问题是现在输出 A(t) 只是一个值(即对于这个 t,我得到 A(t)=0.47),但我需要的是从 t=0 开始的函数数组到 t=t。这是我的主要问题:如何将变量 t 用作积分中的参数和变量

非常欢迎任何建议。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。