python程序设计sgnx_Python中符号函数的数值积分
我是Python新手,所以我的一些問題或想法可能是愚蠢的,但是。。。在
我想畫一個分布D(x)。m和s2是一些給定的實數(shù)。有人告訴我,繪制D(x)圖的最好方法是寫一個函數(shù),它可以解出每個x的積分(在函數(shù)D(x)內(nèi))
其中chi2是這樣定義的:
所以,根據(jù)我對數(shù)學(xué)的了解,我應(yīng)該先積分,然后才能解出每個x,如果我錯了,請糾正我。在
我也被告知要用數(shù)值計算積分,但我不知道怎么做,因為函數(shù)包含符號。
我已經(jīng)嘗試過使用符號積分(盡管我被告知了什么),但是內(nèi)核從來沒有結(jié)束積分的過程,當(dāng)我試圖用數(shù)值計算它的時候。當(dāng)我嘗試數(shù)字積分時,我當(dāng)然使用了lamdify。在
下面是我的代碼:
1試解符號積分from sympy import symbols, integrate, sqrt, exp, oo
s2= 0.0628777415586
m= 5.02422436191
x, n, z=symbols ('x, n, z')
integrate(exp(-n/(z+1) * (x-m)**2/2*s2) * 1/2 / sqrt(z+1), (z, 0, oo))
內(nèi)核,永不停止工作。[我在配方中加入了1/2而不是chi2,打算稍后更改]
另一種方法是嘗試求解數(shù)值積分(從函數(shù)D(x)調(diào)用它)
^{pr2}$
錯誤來自帶有return的行,它說:TypeError: unsupported operand type(s) for *: 'float' and 'rv_frozen'
我想問題是因為chi2被凍住了,但是如何使它工作呢?我的代碼有什么錯誤嗎?我不知道我寫的是不是對的,以及如何解決這個問題。。。
我已經(jīng)為此工作了很長時間,有點絕望,所以任何幫助都是歡迎的。在
總結(jié)
以上是生活随笔為你收集整理的python程序设计sgnx_Python中符号函数的数值积分的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: idefo功能模型图_利用好预后预测模型
- 下一篇: websocket python爬虫_p