第七色在线视频,2021少妇久久久久久久久久,亚洲欧洲精品成人久久av18,亚洲国产精品特色大片观看完整版,孙宇晨将参加特朗普的晚宴

為了賬號安全,請及時(shí)綁定郵箱和手機(jī)立即綁定
已解決430363個(gè)問題,去搜搜看,總會有你想問的

用python模擬核密度的問題

用python模擬核密度的問題

手掌心 2021-08-14 17:47:03
我想使用內(nèi)核密度進(jìn)行模擬觀察,但我有以下錯(cuò)誤類型:TypeError:不支持的操作數(shù)類型為 +:'float' 和 'NoneType'如何解決?這是我使用的代碼:    from matplotlib.pyplot import *    from math import *    from array import *    import numpy as np     from numpy.random import *    from scipy.misc import *    from scipy.stats import *    from scipy import *    from random import *    N=1000    n=30    lamb=0.5    X=lamb*tan(pi*(np.reshape(rand(n,1),n)-0.5)) #loi de Cauchy      x=1    alpha=0.45    def k_gaussien(x):        sigma=1        if(sigma<=0):            return((1/(sigma*sqrt(2*pi)))*exp(-(x**2/(2*sigma**2))))    def h(n,alpha):        for i in range(1,n):               return (i**(1-alpha))    def f_PR(x,X,alpha):        global F;        F = ones((n,1))        h_n = h(n,alpha)        for k in range(2,n):            for i in  range(1,k):                F[k] = F[k-1] + k_gaussien((x-X[i])*i**alpha)            F[k] = F[k-1] * h_n        return F    # Almost sure convergence f_n(x)--> f(x) ps    figure(figsize=(20,10))    fPR=f_PR(x,X,alpha)    T=linspace(1,n,n);    plot(cumsum(fPR)/T)    plot(T,(1/pi)*(lamb/(lamb**2 + y**2))*linspace(1,1,N),"r--",lw=3)#with Cauchy density    grid(True)    title("convergence presque sure",fontsize=20,color="blue")    #Convergence in mean N(0,e2f(x))    hist(fPR,bins=linspace(-10,10,50),normed=True)    y=linspace(-10,10,100);    v=(1/pi)*(lamb/(lamb**2 + y**2))    plot(y,(1/sqrt(2*pi)*v)*exp ((-(y*y)/(2*v**2)))*linspace(1,1,n),'r--') #with cauchy density    title("convergence asymptotique", fontsize=20,color="blue")TypeError Traceback(最近一次調(diào)用最后一次)<ipython-input-76-13bc86608417> in <module>()         38 # Almost sure convergence f_n(x)--> f(x) ps         39 figure(figsize=(20,10))    ---> 40 fPR=f_PR(x,X,alpha)         41 T=linspace(1,n,n);         42 plot(cumsum(fPR)/T)類型錯(cuò)誤:不支持 + 的操作數(shù)類型:'float' 和 'NoneType'
查看完整描述

1 回答

?
慕的地8271018

TA貢獻(xiàn)1796條經(jīng)驗(yàn) 獲得超4個(gè)贊

它在你的定義中:

def k_gaussien(x):
    sigma=1
    if(sigma<=0):
           return((1/(sigma*sqrt(2*pi)))*exp(-(x**2/(2*sigma**2))))

您對 進(jìn)行了硬編碼sigma=1,但您的函數(shù)只會返回一些sigma<=0永遠(yuǎn)不會發(fā)生的情況。所以k_gaussien((x-X[i])*i**alpha)會回來None。因此F[k] = F[k-1] + k_gaussien((x-X[i])*i**alpha)嘗試求和 float 和 None 類型,這是行不通的。


查看完整回答
反對 回復(fù) 2021-08-14
  • 1 回答
  • 0 關(guān)注
  • 179 瀏覽
慕課專欄
更多

添加回答

舉報(bào)

0/150
提交
取消
微信客服

購課補(bǔ)貼
聯(lián)系客服咨詢優(yōu)惠詳情

幫助反饋 APP下載

慕課網(wǎng)APP
您的移動學(xué)習(xí)伙伴

公眾號

掃描二維碼
關(guān)注慕課網(wǎng)微信公眾號