说明

采用root 作对数坐标图，通常会出现横坐标不是均等宽度分布，这是由于在定义historgram，在设置bin的个数和上下限时默认时均有分割，而不是按对数宽度均匀分割。为此，需要自定义bin的分割位置。

具体代码如下

  //--------对数坐标变化--------double low=1.e-8;double up=1e4;double nbin=192;Double_t xAxis[193];for ( int i=0; i <= nbin; i++) {double val_bin=low * std::pow(10., i * std::log10(up/low)/nbin);double exp_10=4.-int(std::log10(val_bin));double factor =std::pow(10., exp_10);val_bin=int(factor*val_bin)/factor;xAxis[i] = val_bin;}TH1D* h1 = new TH1D("h1","h1",192,xAxis);

蒙卡模拟

对这样的histogram 又该如何做随机蒙卡模拟呢？采用求逆的方法

  //求和和归一化double sum=0;// vals.clear();for(int ii=0;ii<192;ii++){// value=vals[ii]*(bins[ii]-bins[ii-1]);// sum = sum+value;// value=vals[ii];//*(bins[ii]-bins[ii-1]);sum = sum+vals[ii];// cout<<"  Sum"<<sum<<endl;// cout<<""<<vals[ii]<<endl;}nvals.clear();// double norm=0.;for(int iii=0;iii<192;iii++){norm = vals[iii]/sum+norm;cout<<"  norm "<<norm<<endl;// cout<<" vals"<<vals[iii]<<endl;nvals.push_back(norm);}//--------------------------// 求逆for(int j=0;j<1000000;j++){int numberOfBin=nvals.size();int nmin=0;int nmid=numberOfBin/2;int nmax=numberOfBin-1;double a=rndm.Uniform();while(nmin!=nmax-1){if(a>nvals[nmid])nmin=nmid;elsenmax=nmid;nmid=nmin+(nmax-nmin+1)/2;}double weight=(bins[nmid]-bins[nmid-1])/(nvals[nmid]-nvals[nmid-1]);double rand=bins[nmid-1]+weight*(a-nvals[nmid-1]);h1->Fill(rand);}

Root作对数坐标图及蒙卡随机模拟相关推荐

15 单因子利率模型蒙卡模拟
15 单因子利率模型蒙卡模拟 15.1 简介 15.1.1 单因子模型蒙卡抽样 15.1.2 由抽样结果计算债券期权价格 15.2 蒙卡模拟步骤 15.2.1 蒙卡抽样利率变化路径 15.2.2 计算 ...
第9节蒙卡模拟计算美式期权价格(b)
第9节蒙卡模拟计算美式期权价格(b) 9.1 简介 9.2 参数化执行边界算法 9.3 算法Python代码实现 9.4 计算示例 9.5 参考资料 9.1 简介使用蒙卡模拟计算美式期权价格除了最 ...
第7节蒙卡模拟计算路径依赖型期权价格
第7节蒙卡模拟计算路径依赖型期权价格 7.1 简介 7.2 蒙卡模拟计算回望/亚式期权算法 7.3 算法 Python 代码实现 7.4 计算示例 7.5 相关说明 7.5.1 由均匀分布产生正态分 ...
第8节蒙卡模拟计算美式期权价格(a)
第8节蒙卡模拟计算美式期权价格(a) 8.1 简介 8.2 最小二乘法计算美式期权价格 8.3 算法Python代码实现 8.4 计算示例 8.5 相关说明 8.5.1 美式看涨期权不会被提前行使 ...
随机投点法计算定积分java_11 随机模拟积分 | 统计计算
11.4 高维定积分上面的两种计算一元函数定积分的方法可以很容易地推广到多元函数定积分, 或称高维定积分. 设\(d\)元函数\(h(x_1, x_2, \dots, x_d)\)定义于超矩形 \[ ...
基于启发式算法与单目优化和马尔科夫模型的进出口公司的货物装运策略——整数线性规划随机模拟
目录进出口公司的货物装运策略摘要一.问题的重述 1问题的背景 2问题的提出 3目的及意义二.符号说明三.模型假设四.建模准备与问题分析 4.1 线性整数规划 4.2 三维空间分割启发式算法 ...
使用计算机模拟抛硬币实验动画课件,概率论与数理统计实验2抛硬币实验的随机模拟实验报告...
概率论与数理统计实验2抛硬币实验的随机模拟实验报告 (4页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 14.90 积分 <概率论与数理统计> ...
用计算机随机模拟掷骰子的实验,云南师范大学概率论实验报告随机事件的模拟--模拟掷均匀硬币的随机试验...
云南师范大学概率论实验报告随机事件的模拟--模拟掷均匀硬币的随机试验 (7页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 9.9 积分教学类捻杈告班级1 ...
随机模拟【1】：随机模拟的研究范围与特征
本系列同步发布于本人的知乎专栏:确定性随机最新想开始再系统学习一次数理统计和随机过程方便的知识,但是耽搁了一段时间,主要原因是不知道从哪里开始.最近,刚好在网上"随机"地买了一本 ...

Root作对数坐标图及蒙卡随机模拟

Root作对数坐标图及蒙卡随机模拟

说明

蒙卡模拟

Root作对数坐标图及蒙卡随机模拟相关推荐

最新文章

热门文章