【转载】拟蒙特卡洛方法(Quasi-Monte Carlo)积分实例

sir_lab

本帖转载于Matlab技术论坛，原帖参见http://www.matlabsky.com/thread-673-1-1.html



%使用Matlab提供的函数求积分,exp(-1/2*x^2)在(0,1)间积分


format long;
syms x
a = sym(1/2);
f = exp(-a*x^2);
ezplot(f)
disp(int(f,-1,1));
fprintf('integral result:%1.18f.\n',double(int(f,0,1)));
%disp(double(int(f,0,1)));


%使用拟蒙特卡洛方法积分
%得到拟蒙特卡洛序列,即低偏差序列,halton法
%如果有相关的工具箱的话,可以用Matlab里面的haltonset,faureset,sobolset函数实现




x=halton(10000,2,5577);
n=length(x);
mju=0;
for i=1:n
    mju=mju + exp(-0.5*x(i)^2);
end
mju=mju/n;
fprintf('Quasi-Monte Carlo result:%1.18f.\n',mju);
%disp(mju);
%使用蒙特卡洛方法积分
%得到Uniform序列,
x=random('unif',0,1,10000,1);
n=length(x);
mju=0;
for i=1:n
    mju=mju + exp(-0.5*x(i)^2);
end
mju=mju/n;
fprintf('Monte Carlo result:%1.18f.\n',mju);

%=============生成HALTON序列========================
function result = halton( m,base,seeder )
%生成HALTON序列
% Check inputs
if nargin < 3
seeder = 0;
if nargin < 2
      error('MATLAB:Halton:NotEnoughInputs',...
             'Not enough input arguments. See Halton.');
end
end
res=0;
n=length(base);
for i=1:m
   
    for j=1:n
        element=0;
        temp=seeder+i;
        k=1;
        while temp>0
            element(k)=rem(temp,base(j));
            temp=fix(temp/base(j));
            k=k+1;
        end
        res(i,j)= 0;
        for k=1:length(element)
            res(i,j)=res(i,j)+element(k)/(base(j)^k);
        end
    end
   
end
result=res;

jkwlyn1013

专门从事大宗商品交易的商业活动。零售的对称。是商品流通中不可缺少的一个环节。通常有两种情况：①商业企业将商品批量售给其他商业企业用作转卖。②商业企业将用作再加工的生产资料供应给生产企业。
　　服装批发是随着商品经济的发展而产生的。商品生产和商品交换的发展，使商品购销量增大，流通范围扩展，生产者相互之间、生产者与零售商之间直接进行商品交换，常有困难或不方便，于是产生了专门向生产者直接购进商品，然后再转卖给其他生产者或零售商的服装批发商业，商业部门内部有了服装批发和零售之间的分工。服装批发业务一般由服装批发企业来经营，每次批售的商品数量较大，并按服装批发价格出售。商品的服装批发价格低于零售价格，即存在着批零差价，其差额由零售企业所耗费的流通费用、税金和利润构成。商业服装批发是生产与零售之间的中间环节。通过商业服装批发活动，   社会产品从生产领域进入流通领域，起到组织和调动地区之间商品流通的作用。还可通过商品储存发挥“蓄水池”作用，平衡商品供求。
　　在电子商务兴起后,服装批发业经常被称为:B2B

chengchneg5l8

写的不错，还要再加油啊!!!


















注册香港公司
香港公司注册