99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話：13121318867

登錄
首頁精彩閱讀Python金融大數(shù)據(jù)分析-蒙特卡洛仿真
Python金融大數(shù)據(jù)分析-蒙特卡洛仿真
2017-12-12
收藏

Python金融大數(shù)據(jù)分析-蒙特卡洛仿真

1.簡單的例子

了解一點金融工程的對這個公式都不會太陌生，是用現(xiàn)在股價預測T時間股價的公式，其背后是股價符合幾何布朗運動，也就是大名鼎鼎的BSM期權定價模型的基礎。
我們假設現(xiàn)在一個股票的價值是100，那么兩年后是多少呢？
[python]view plaincopy

    import numpy as np  
    import pandas as pd  
    import matplotlib.pyplot as plt  
    S0 = 100  
    r = 0.05  
    sigma = 0.25  
    T = 2.0  
    I = 10000  
    ST1 = S0*np.exp((r - 0.5*sigma**2)*T+sigma*np.sqrt(T)*np.random.standard_normal(I))  
    plt.hist(ST1,bins = 50)  
    plt.xlabel('price')  
    plt.ylabel('ferquency')  

    運行的結(jié)果如下所示：

很明顯，是一個lognormal分布，因為這樣的假設下，價格符合lognormal分布，收益率符合正態(tài)分布。
2.簡單的蒙特卡洛路徑

上面是一步到位的，那么如果我們中間分很多個小時間段來仿真呢？可以知道，物理問題是一樣的，結(jié)果也不會有差異。
[python]view plaincopy

    import numpy as np  
    import pandas as pd  
    import matplotlib.pyplot as plt  
    import scipy.stats as scs  
    S0 = 100  
    r = 0.05  
    sigma = 0.25  
    T = 2.0  
    I = 10000  
    #ST1 = S0*np.exp((r - 0.5*sigma**2)*T+sigma*np.sqrt(T)*np.random.standard_normal(I))  
    #plt.hist(ST1,bins = 50)  
    #plt.xlabel('price')  
    #plt.ylabel('ferquency')  
      
    M = 50  
    dt = T/M  
    S = np.zeros((M + 1,I))  
    S[0] = S0  
    print S[0]  
    for t in range(1,M+1):  
        S[t] = S[t-1]*np.exp((r-0.5*sigma**2)*dt+sigma*np.sqrt(dt)*np.random.standard_normal(I))  
    plt.hist(S[-1],bins = 50)  
    plt.xlabel('price')  
    plt.ylabel('frequency')  
    plt.show()  
    plt.plot(S[:,:],lw = 1.5)  
    plt.xlabel('time')  
    plt.ylabel('price')  
    plt.show()  

我們不僅可以得到最終的分布，也可以知道價格路徑，而這一價格路徑，才是真正代表了蒙特卡洛的精髓。

如果我們繪制得路徑更加多一點，就是這樣的一個效果：

從側(cè)面看，其實就是一個lognormal分布。

數(shù)據(jù)分析咨詢請掃描二維碼

若不方便掃碼，搜微信號：CDAshujufenxi

數(shù)據(jù)分析師資訊
更多

OK
客服在線
立即咨詢
客服在線
立即咨詢
') } function initGt() { var handler = function (captchaObj) { captchaObj.appendTo('#captcha'); captchaObj.onReady(function () { $("#wait").hide(); }).onSuccess(function(){ $('.getcheckcode').removeClass('dis'); $('.getcheckcode').trigger('click'); }); window.captchaObj = captchaObj; }; $('#captcha').show(); $.ajax({ url: "/login/gtstart?t=" + (new Date()).getTime(), // 加隨機數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調(diào)用 initGeetest 進行初始化 // 參數(shù)1：配置參數(shù) // 參數(shù)2：回調(diào)，回調(diào)的第一個參數(shù)驗證碼對象，之后可以使用它調(diào)用相應的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須，不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗服務器是否宕機 new_captcha: data.new_captcha, // 用于宕機時表示是新驗證碼的宕機 product: "float", // 產(chǎn)品形式，包括：float，popup width: "280px", https: true // 更多配置參數(shù)說明請參見：http://docs.geetest.com/install/client/web-front/ }, handler); } }); } function codeCutdown() { if(_wait == 0){ //倒計時完成 $(".getcheckcode").removeClass('dis').html("重新獲取"); }else{ $(".getcheckcode").addClass('dis').html("重新獲取("+_wait+"s)"); _wait--; setTimeout(function () { codeCutdown(); },1000); } } function inputValidate(ele,telInput) { var oInput = ele; var inputVal = oInput.val(); var oType = ele.attr('data-type'); var oEtag = $('#etag').val(); var oErr = oInput.closest('.form_box').next('.err_txt'); var empTxt = '請輸入'+oInput.attr('placeholder')+'！'; var errTxt = '請輸入正確的'+oInput.attr('placeholder')+'！'; var pattern; if(inputVal==""){ if(!telInput){ errFun(oErr,empTxt); } return false; }else { switch (oType){ case 'login_mobile': pattern = /^1[3456789]\d{9}$/; if(inputVal.length==11) { $.ajax({ url: '/login/checkmobile', type: "post", dataType: "json", data: { mobile: inputVal, etag: oEtag, page_ur: window.location.href, page_referer: document.referrer }, success: function (data) { } }); } break; case 'login_yzm': pattern = /^\d{6}$/; break; } if(oType=='login_mobile'){ } if(!!validateFun(pattern,inputVal)){ errFun(oErr,'') if(telInput){ $('.getcheckcode').removeClass('dis'); } }else { if(!telInput) { errFun(oErr, errTxt); }else { $('.getcheckcode').addClass('dis'); } return false; } } return true; } function errFun(obj,msg) { obj.html(msg); if(msg==''){ $('.login_submit').removeClass('dis'); }else { $('.login_submit').addClass('dis'); } } function validateFun(pat,val) { return pat.test(val); }