蒙特卡洛求二重积分matlab

---

归弦15532699437问： 如何使用MATLAB进行蒙特卡洛模拟 小小知识站 - ？
兰溪市贝感回答： 贴一个蒙特卡洛方法的matlab程序,供大家使用. {3 x& K/ i1 i( D8 C0 c$ O % Example Monte Carlo Simulation in Matlab 0 O5 \; P＂ t# t7 v8 c& @ % Function: y = x2^2/x1 5 Z0 W4 e9 q, d5 B+ c % % Generate n samples from a normal distribution...

归弦15532699437问： matlab求一重和二重积分的近似值 - ？
兰溪市贝感回答： syms x y int(exp(-x^2)*cos(x^2),x,0,pi); vpa(ans,5) ans =.68850-.64202e-5*i int(int(cos(x^2-y^2),y,0,sqrt(pi)),x,0,sqrt(pi)); vpa(ans,5) ans =1.2401

归弦15532699437问： 如何用matlab实现蒙特卡洛法求定积分 - ？
兰溪市贝感回答： 利用蒙特卡洛法计算圆周率的代码,即用unifrnd()随机函数和mean()平均函数来计算 fun=inline('4./(1+x.*x)'); 定义积分函数 N=100000;x=unifrnd(0,1,N,1);mean(fun(x))

归弦15532699437问： 求助,matlab编写二重积分程序 - ？
兰溪市贝感回答： 应该是曲面积分∫∫f(x,y)dxdy,积分区域为以(5,0)为圆心,以4为半径的圆形区域

归弦15532699437问： 求,利用matlab计算二重数值积分f=exp((x - 2)^2+(y - 1)^2))在圆域x^2+y^2=1内的积分 - ？
兰溪市贝感回答： 提供两种方法1. 直角坐标系符号积分 syms x y int_y = int(exp((x-2)^2+(y-1)^2),y,-sqrt(1-x^2),sqrt(1-x^2)); % 先对y积分 I = vpa(int(int_y,x,-1,1))2. integral2()数值积分 f=@(x,y)exp(((x-2).^2+(y-1).^2)); y_l = @(x) -sqrt(1-x.^2); y_u = @(x) sqrt(1-x.^2); I = integral2(f,-1,1,y_l,y_u) 计算结果为:I = 6.4181e+03 就是6418.1

归弦15532699437问： 求二重积分,级数的解 - ？
兰溪市贝感回答： 用matlab计算对二重级数进行二重积分可以用循环语句来实现.求解思路:当m、n为某值,用二重积分函数计算其积分值,然后用求和函数累加.下面的计算值为(m、n各为1~100) S=5313.34975572556(耗时 1643.778 s)

归弦15532699437问： 怎么用Matlab计算这个二重积分??求程序代码 - ？
兰溪市贝感回答： 程序:fun=sin(0.5*pi*x./y);%% a=int(int(fun,y,sqrt(x),x),x,1,2); b=simple(a) %化简 I=vpa(b,4) %得到4位近似解,也可以任意N位解 结果:I =0.2719

归弦15532699437问： matlab中蒙特卡洛算法代码求解释 - ？
兰溪市贝感回答： 随机在正方形区域生成10000个点,find函数的作用:统计出在给定曲线y=2-x^2和曲线y^3=x^2之间的点数. 这样求S的目的是计算比值,比值即为所求.

归弦15532699437问： Matlab的Monte carlo法求定积分 区域( - π/2,π/2),y=cosx所围面积 - ？
兰溪市贝感回答： function y = MonteCarlo(fun,a,b,N) A=unifrnd(a,b,N); B=fun(A); y=(b-a)/N^2*sum(sum(B)); 自己编了个MonteCarlo求一维定积分的函数.在command里如下调用 fun=@(x) cos(x); MonteCarlo(fun,-pi/2,pi/2,500) % N越大,结果可信度越高,速度也越慢.用quad(fun,-pi/2,pi/2) 验证结果准确度.

归弦15532699437问： matlab求2重积分 - ？
兰溪市贝感回答： 既然为2重积分,那为什么有三个变量,Matlab中dblquad的定义为q = dblquad(fun,xmin,xmax,ymin,ymax),只有两个变量.若为三重积分的话用triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax).

---