蒙特卡洛方法matlab程序

hattrick110

贴一个蒙特卡洛方法的matlab程序，供大家使用。
祝大家比赛都能取得好成绩
* L' ]& P+ z4 i! p) A$ c

% Example Monte Carlo Simulation in Matlab
8 J6 `6 G& f. }/ d  o% Function: y = x2^2/x1
%
% Generate n samples from a normal distribution
0 z+ f2 }, B7 w/ f% r = ( randn(n,1) * sd ) + mu
: d# S# M' @$ E) H" E2 P, }1 l+ r% mu : mean
% sd : standard deviation
) T" B) [8 E" j* E6 M0 P1 z%
) {; g& d- x% m/ G% Generate n samples from a uniform distribution
% r = a + rand(n,1) * (b-a)
% a : minimum
* d3 \1 h, s1 G- ]; ], e% b : maximum
) B8 t' m* U" ^/ u' n# Gn = 100000; % The number of function evaluations
' A2 X# O5 }& f1 x, ]' ?% --- Generate vectors of random inputs
% x1 ~ Normal distribution N(mean=100,sd=5)
% x2 ~ Uniform distribution U(a=5,b=15)
x1 = ( randn(n,1) * 5 ) + 100;
8 r; v1 B1 X5 V' m# |- o# tx2 = 5 + rand(n,1) * ( 15 - 5 );
% --- Run the simulation
% Note the use of element-wise multiplication
y = x2.^2 ./ x1;
8 {$ ]6 ]5 L' d% O! M3 j% O  r6 Z% --- Create a histogram of the results (50 bins)
hist(y,50);
3 K, X# h: n8 t$ C2 B8 u% E8 B% --- Calculate summary statistics
4 N2 y$ E! q5 {; V" t, m' cy_mean = mean(y)
y_std = std(y)
6 Q: V1 r7 S5 C$ M* g/ @. Q% s+ ?y_median = median(y)

mma

多谢拉！！！！！！！

sjtucat5

thankyouverymuch

tzhggx

顶啊！不错啊！

BigTou

谢谢啦...................

rever99

个盲目反卷积程序是用模拟退火法做的，想请教

jafy2009

哈哈哈。。。谢谢！！！

eyehappy

thank you!!!!

cherry_s1218

看看。。。。。。

shuxinjian24134

xinga!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!