数模
关注
分享
扫一扫
默默小艺
乐
展开
-
R学习之统计实验(一)--电梯问题（R语言编程）-----数模
问题有r个人在一楼进入电梯楼上共有n层设每个乘客在任何层楼出电梯的概率相同.试建立一个概率模型,求直到电梯中的人下完为止,电梯需停次数的数学期望并对r=10,n=7进行模拟验证直接R模拟（大数定律）注意括号各种英文输入，R语言编写时候看不出明显区别elevator <- function(r,n,m){ #r个人，n为1楼以上楼层数，m为数据个数 y=0 for...原创 2018-07-29 16:37:40 · 2292 阅读 · 0 评论 -
matlab学习之蒙特卡罗 --渡口模型和赶火车问题（matlab编程）-----数模
1. 渡口模型问题描述:问题分析:模型建立:模拟程序设计:matlab代码其中函数说明结果2. 赶火车DONE！！！1. 渡口模型问题描述:一个渡口的渡船营运者拥有一只甲板长32米,可以并排停放两列车辆的渡船.他在考虑怎样在甲板上安排过河车辆的位置,才能安全地运过最多数量的车辆. 分析:怎样安排过河车辆,关心一次可以运多少辆各类车....原创 2018-08-07 17:58:03 · 5213 阅读 · 2 评论 -
R学习之统计算法与R优化包（极大似然估计练习） --（R语言编程）-----数模
题目R代码#l=0data=rcauchy(n=1000,location = 0,scale = 10)f<-function(pa){ result= -log(pi*pa*(1+((data)/pa)^2)) -sum(result)}ML.op=optimize(f,interval = c(0,100))#不需要参数初值，interval定义参数...原创 2018-08-06 19:15:13 · 4411 阅读 · 2 评论 -
R学习之统计算法与R优化包（极大似然估计） --（R语言编程）-----数模
library(MASS);attach(geyser);hist(waiting,freq = F)#mnf&amp;amp;lt;-function(pa,data){ x&amp;amp;lt;-dnorm(data,pa[2],sqrt(pa[4])) y&amp;amp;lt;-dnorm(data,pa[3],sqrt(pa[5])) pdf=pa[1]*x+(1-pa[1])*y l=sum(log(...原创 2018-08-06 15:42:17 · 20113 阅读 · 3 评论 -
R学习之统计算法与R优化包（矩估计法也即Newton） --（R语言编程）-----数模
矩估计说明：R代码结果DONE！！！矩估计说明：就是把求方程解的Newton法过渡到统计学里面，达到参数估计的目的，其实也是对于方程的求解。 直接上例题 R代码#矩估计#下面是之前的Newtons代码Newtons<-function(fun,x,ep=1e-5,it_max=100){#itmax表示最大迭代次数 index<...原创 2018-08-06 10:25:11 · 4029 阅读 · 2 评论 -
R学习之统计算法与R优化包（Newton法） --（R语言编程）-----数模
Newton法说明：Newton法R程序：其中部分说明练习题1代码结果练习题2结果DONE！！！Newton法说明： Newton法R程序：例 {x21+x22−5=0(x1+1)∗x2−(3∗x1+1)=0{x12+x22−5=0(x1+1)∗x2−(3∗x1+1)=0\left\{\begin{array...原创 2018-08-05 17:59:21 · 12286 阅读 · 2 评论 -
R学习之统计算法与R优化包（二分法） --（R语言编程）-----数模
二分法其中函数说明结果作业题：结果DONE！！！二分法 R语言编程：#（二分法）dichotomy<-function(f,a,b,eps=1e-5){ if(f(a)*f(b)>0) list(fail="find root is fail!") else{ repeat{ if(abs(b-a)&...原创 2018-08-05 11:19:50 · 8218 阅读 · 2 评论 -
R学习之蒙特卡罗积分 --（R语言编程）-----数模
蒙特卡罗方法介绍问题1其中函数说明结果问题2结果问题3结果问题4结果问题5结果问题6结果DONE！！！蒙特卡罗方法介绍 问题1∫10sin(x)xdx∫01sin(x)xdx \int_{0}^{1} \frac {sin(x)}xdxR语言编程：#0到1 sin(x)/x积分（投点法...原创 2018-08-04 19:42:49 · 21353 阅读 · 6 评论 -
R学习之统计实验（十）--新乒乓问题（R语言编程）-----数模
问题盒中有12个乒乓球．其中9个是新的，第一次比赛时从中选出3个，比赛后仍放回盒中，第二次比赛时再从盒中选出3个．(1）求第二次取出的球都是新球的概率；(2）又已知第二次取出的球都是新球，求第一次取出的球都是新球的概率．分别作理论和模拟分析． 直接R模拟（大数定律）newpingpang<-function(n){ m<-0;z<-0 for(i in 1...原创 2018-07-30 20:03:10 · 1148 阅读 · 0 评论 -
R学习之统计实验（九）--和的偶数概率问题（R语言编程）-----数模
问题从0-9中随机不放回选4个数，求组成偶数的概率．分别作理论和模拟分析．直接R模拟（大数定律）eventest<-function(n){ m<-0 for(i in 1:n){ x<-sample(x = 0:9,size = 4,replace = F) if(sum(x)%%2==0){m=m+1} } rt<-c...原创 2018-07-30 18:09:31 · 1309 阅读 · 0 评论 -
R学习之统计实验（八）--统计次数概率问题（R语言编程）-----数模
问题从1~10随机有放回选7个数，求10恰好出现2次的概率．分别作理论和模拟分析直接R模拟（大数定律）doubleten<-function(n){ m<-0 for (i in 1:n) { u<-runif(7) x<-floor(u*10)+1 if(sum(x==10)==2){m=m+1} } rt<...原创 2018-07-30 17:36:01 · 6370 阅读 · 0 评论 -
R学习之统计实验(七）--大于1问题（R语言编程）-----数模
问题令U1U1U_1,U2U2U_2，⋅⋅⋅······是一个i.i.d序列，且UiUiU_i~UUU(0,1).令N =min {n:∑ni=1Ui∑i=1nUi\sum_{i=1}^n U_i>1} 用统计模拟方法求E(N）直接R模拟（大数定律）注意括号各种英文输入，R语言编写时候看不出明显区别morethanone <-function(n){ ...原创 2018-07-30 13:50:00 · 3743 阅读 · 5 评论 -
R学习之统计实验（六）--生日问题（R语言编程）-----数模
问题美国数学家伯格米尼曾经做过一个别开生面的实验：在个盛况空前、人山人海的世界杯赛场上，他随机地在某号看台上召唤了22个球迷，请他们分别写下自自己的生日，结果竟发现其中有两人同生日．怎么会这么揍巧呢？通过分析求出理论概率，并运用模拟方法求出概率的近似值．人数分别为40, 50, 64时，又如何？ 直接R模拟（大数定律）注意括号各种英文输入，R语言编写时候看不出明显区别birth...原创 2018-07-30 11:45:03 · 4952 阅读 · 1 评论 -
R学习之统计实验（五）--赶火车问题（R语言编程）-----数模
问题问题一列火车从A站开往B站，某人每天赶往B站赶火车．他己了解到火车从A站到B站的运行时间是服从均值为30分钟．标准差为2分钟的 正态分布．火车大约13点离开A站．此人大约13:30到达B站．火车离开A站 的时刻及概率和此人到达B站的时刻及率如表所示．用模拟实验的方法 求他赶上火车的概率． 直接R模拟（大数定律）train<-function(n){ m=0 ...原创 2018-07-30 11:22:26 · 3862 阅读 · 2 评论 -
R学习之统计实验（四）--蒲丰投针（R语言编程）-----数模
问题法国数学家蒲丰于1777年完成的投针实验是历史上第一个统计模拟实验，可以利用实验结果估计圆周率二．该实验可以表述为：向一簇距离为2a的平行线构成的平面中投掷一根长度为2L的针，求针与直线相交的概率．这里假设L直接R模拟（大数定律）注意括号各种英文输入，R语言编写时候看不出明显区别...原创 2018-07-30 10:05:55 · 7810 阅读 · 2 评论 -
R学习之统计实验（三）--报童问题（R语言编程）-----数模
问题设某报每日的潜在卖报数数服从参数为lambda的Poisson分布．每卖出一份报可得报酬a 元，卖不掉退回则侮每份赔偿b元．若某日该报童买 进n份报，试求其期望所得，并对a=1.5, b=0.6, lambda=120,买进报数n1= 100.n2 = 140分别作模拟分析． 直接R模拟（大数定律）newspaper&amp;lt;-function(n,lambda,m,a,b){...原创 2018-07-30 09:02:42 · 2604 阅读 · 2 评论 -
R学习之统计实验（二）--掷色子问题（R语言编程）-----数模
问题问题一颗股子投4次至少得到一个六点与两颗般子投24次至少得到 一个双六点，这两个事件哪一个概率更大？ 直接R模拟（大数定律）注意括号各种英文输入，R语言编写时候看不出明显区别sice&lt;-function(n){ p1=0;p2=0; for (i in 1:n) { u1&lt;-runif(4);x1&lt;-floor(u1*6)+1; ...原创 2018-07-29 21:30:00 · 7695 阅读 · 2 评论 -
matlab学习之优化模型 --（matlab编程）-----数模
这里主要对于matlab优化工具箱进行整理linprog求解线性规划：intlinprog求解求解整数或混合整数规划。例题1：例题2例题3例题4例题5DONE！！！这里主要对于matlab优化工具箱进行整理官方文档有例子，用到哪个根据目录检索，再官方文档查阅linprog求解线性规划：求解一个线性不等式、线性等式和边界的简单线性规划。 调用...原创 2018-08-08 20:32:39 · 20377 阅读 · 3 评论