版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zyj8119

木虫 (著名写手)

模拟EPI: 10
应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【转帖】使用matlab对某正态函数进行蒙特卡罗模拟

在如下图的公式中

g(x)=f(m,n,L), 其中,m,n,L均服从正态分布,分布情况也在所给的图中.
请使用matlab,用蒙特卡罗模拟法对该函数进行模拟,得出g(x)大于0的概率.
请将所需程序告知我,并给出结果.

CODE:

function func_baidu_monte_carlo

randn('state', sum(100*clock));

%利用时钟设置随机种子，这样每次产生的随机数就不同了

N = 1000000;

% 设置样本个数

m = 0.43 + sqrt(0.008) * randn(N,1);

n = 0.19 + sqrt(0.008) * randn(N,1);

L = 0.5 + sqrt(0.008) * randn(N,1);

%随机生成m, n, L，根据他们各自的分布

g = 16.4445 .* (m.*m.*L./(m.*m+n.*n))-8.2;

%计算g

fprintf('the probability of g>0 is %.2f\n', sum(g>0)/length(g));

%统计 g>0的概率，结果为0。18

已经测试通过，可以运行，结果为

the probability of g>0 is 0.18

回复此楼

» 猜你喜欢

焦虑已经有10人回复
298-一志愿中国农业大学-求调剂已经有3人回复
265求调剂已经有3人回复
297求调剂已经有3人回复
302求调剂已经有8人回复
环境工程调剂已经有7人回复
268求调剂已经有6人回复
311求调剂已经有11人回复
工科材料085601 279求调剂已经有4人回复
301求调剂已经有9人回复

好好学习，天天向上。

1楼 2010-09-17 09:27:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 zyj8119 的主题更新